Monte Carlo Generation of Pi | Compiled Thoughts by John Ciacia

Monte Carlo Generation of Pi

May 21, 2009 · C++ · Snippets · No Comments

[code lang="cpp"]#include
#include
#include
#include
int main ()
{
int in = 0;
srand(time(NULL));

for(int i = 0; i < 1000000; i++) {
float x = (rand()%1000000)/1000000.0;
float y = (rand()%1000000)/1000000.0;
double r = sqrt((pow(x,2)+pow(y,2)));
if(r <= 1) in++;
}
printf("percent error=%fn", ((((4*in/1000000.0)-3.141592)/3.141592)*100));
printf("pi=%fn", (4*in/1000000.0));
return 0;
}
[/code]

C++ cpp monte carlo pi

Related posts:

Euler’s Constant Calculate Euler’s constant. [code lang="cpp"]#include /** * This code will calculate e * as the sum of the infinite series...

Leave a Reply