// 
// pi -- critical
//
#include <omp.h>
#include <stdio.h>
#include <math.h>

const int n=10000;
 
inline double f(double x)
{
	return 4/(x*x+1);
}

int main()
{
	double a=0.0, b=1.0;
	double h=(b-a)/n;
	double mypi=0.0;
	int i;
	
	#pragma omp parallel for num_threads(4) private(i) shared(a,h,mypi)
		for(i=1; i<n; i++)
		{
			#pragma omp critical
			mypi = mypi + f(a+i*h);
		}
	// printf("mypi=%.10f\n", mypi);
		
	mypi = mypi + (f(a) + f(b))/2;
	mypi = h*mypi;
	
	printf("mypi=%.10f, err=%.2e\n", mypi, fabs(mypi-4*atan(1)));

	return 0;
}