pi_mc.c
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include <stdlib.h>
	#include <stdio.h>
	#include <math.h>
	#include <string.h>
	#define SEED 35791246

	int main(int argc, char* argv)
	{
	int niter=0;
	double x,y;
	int i,count=0; /* # of points in the 1st quadrant of unit circle */
	double z;
	double pi;

	niter=99999999;

	/* initialize random numbers */
	srand(SEED);
	count=0;
	for ( i=0; i<niter; i++) {
	x = (double)rand()/RAND_MAX;
	y = (double)rand()/RAND_MAX;
	z = xx+yy;
	if (z<=1) count++;
	}

	pi=(double)count/niter*4;
	printf("# of trials= %d , estimate of pi is %g \n",niter,pi);

	exit (0);
	}