// MODIFIED FOR TEST CODE
double MonteCarlo_integrate(int Num_samples, int print)
{
Random R = new_Random_seed(SEED);
int under_curve = 0;
int count;
for (count=0; count<Num_samples; count++)
{
double x= Random_nextDouble(R);
double y= Random_nextDouble(R);
if ( x*x + y*y <= 1.0)
under_curve ++;
}
Random_delete(R);
// ADDED FOR TEST CODE print results
if( print )
{
printf("\n%d\n", under_curve);
}
return ((double) under_curve / Num_samples) * 4.0;
}
double *RandomVector(int N, Random R)
{
int i;
double *x = (double *) malloc(sizeof(double)*N);
for (i=0; i<N; i++)
x[i] = Random_nextDouble(R);
return x;
}
