How do you Program to calculate the value of pi using Monte Carlo method?

Answer 1

//Program to calculate the value of pi using monte carlo method:

/*

NADEEM AHMAD

MCA 2012

NITISH HOODA

M TECH 2012

SHARDA UNIVERSITY

GREATER NOIDA*/

#include<stdlib.h>

#include<stdio.h>

#include<math.h>

#include<string.h>

#include<conio.h>

#define SEED 35791246

void 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;

printf("Enter the number of iterations used to estimate pi: ");

scanf("%d",&niter);

/* initialize random numbers */

srand(SEED);

count=0;

for ( i=0; i<niter; i++) {

x = (double)rand()/RAND_MAX;

y = (double)rand()/RAND_MAX;

z = x*x+y*y;

if (z<=1) count++;

}

pi=(double)count/niter*4;

printf("# of trials= %d , estimate of pi is %g \n",niter,pi);

getch();

}