在linux下用 g++ pi.cpp -o pi -lm 編譯。

//蒙特卡洛法求圓周率

#include

#include

#include

#include

using namespace std;

double FindPi(const long n)

{

srand(time(NULL));

long sum = 0;

double x = 0.0, y = 0.0;

for (long i = 0; i < n; i++)

{

x = (double)(rand())/(double)RAND_MAX;

y = (double)(rand())/(double)RAND_MAX;

if (x*x + y*y <= 1)

{

sum++;

}

if (i % 10000 == 0)

{

system("clear");

cout<

}

}

return 4.0*((double)sum/(double)n);

}

int main()

{

long n = 0;

cout<

cin>>n;

if (n <= 0)

{

cout<

exit(1);

}

double pi = FindPi(n);

cout.precision(50);

cout<

return 0;

}

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

以上是生活随笔為你收集整理的蒙特卡洛法求圆周率 c语言,c++蒙特卡洛法求圆周率的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。