我安切洛蒂就是饿死:谁会用蒙特卡洛法求Pi啊。帮帮小弟我啊~谢谢。

来源:百度文库 编辑:高校问答 时间:2024/05/05 08:34:42
可以用C/C++编写,谢谢了。

double PI( DOWRD dwCount/*测试次数*/ )
{
double R = 32767.0;
DWORD count = 0;
srand( (unsigned)time( NULL ) ); // 这段代码当时没写,因为没有msdn,只记得一个rand(),连32767也是试出来的。
for( DWORD i=0;i<dwCount;i++ )
{
int x = rand();
int y = rand();
double r = sqrt( x*x + y*y );
if( r<R )
{
count++;
}
}
return 4.0*count/dwCount;
}