软件下载 | 资讯教程 | 最近更新 | 下载排行 | 一键转帖 | 发布投稿
您的位置:最火下载站 > 电脑教程 > 编程开发 > 工具与技巧 > 计算物理:蒙特卡洛法求圆周率

计算物理:蒙特卡洛法求圆周率

计算物理:蒙特卡洛法求圆周率,在linux下用 g++ pi.cpp -o pi -lm 编译。

//蒙特卡洛法求圆周率
#include <iostream>
#include <cmath>
#include <cstdlib>
#include <ctime>

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<<"运行了"<<i<<"次运算,完成了总运算的"<<(double)i/(double)n*100.0<<"%"<<endl;
}
}
return 4.0*((double)sum/(double)n);
}

int main()
{
long n = 0;
cout<<"输入模拟次数:";
cin>>n;
if (n <= 0)
{
cout<<"输入错误"<<endl;
exit(1);
}
double pi = FindPi(n);
cout.precision(50);
cout<<"模拟"<<n<<"次,pi = "<<pi<<endl;
return 0;
}

    相关阅读
    网友评论
    栏目导航
    推荐软件