#include <iostream.h>
#include <math.h>

int main( int argc, char *argv[])
{
   int n = 1000000000;
   int i;
  
   double pi;
   double h;
   double x;

   // $INIT$
 
   h = 1.0 / (double) n;
   pi= 0.0;

   // $SUM$ DECOMP(i=1,n+1) REDUCE(pi) REDUCETYPE(double)
   
   for(i=1; i<n+1;i++)
   {
	  x = h * ((double)i - 0.5);
	  pi +=(4.0 / (1.0 + x*x));
   }
   pi = h * pi;

   // $ENDSUM$
   
   // $ROOT$
   cout<<"pi = "<<pi<<endl;
   // $ENDROOT$

   // $FIN$
 
   return 0;
}