#include <stdio.h>
#include <omp.h>
#include <stdio.h>
#include <stdlib.h>

static long num = 100000000;


int main() {
  double x,v, width, sum=0.0;
  int i;
  double area;

  int tid; 
  int numthreads;

  width=1./(double)num;
  
#pragma omp parallel for private(x,v) reduction(+:sum) 
  for (i=0; i<num; i++) {
    x = (i+0.5)*width;
    v = 4.0/(1+x*x);
    sum+=v;
  }
  area = width*sum;
  
#pragma omp parallel private(numthreads, tid)
  {
  tid = omp_get_thread_num();
  if (tid ==0) numthreads = omp_get_num_threads();
  else numthreads = 44;
  
if (tid ==0 ) {
    printf("I am %d and I computed the value of PI: %12.10g with %d threads\n",
	 tid,area,numthreads);
  }
  }
}