[code lang="cpp"]#include
#include
#include
#include
{
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]
21. May 2009
No Comments »