xjwjwsn 2014-12-14 23:23 采纳率: 100%
浏览 2214
已采纳

用蒙特卡洛方法求解圆周率

//蒙特卡洛法求圆周率
#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<<"运行了"<<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;
}

这个代码为什么算出来的总是在3左右,而不是3.14?

展开全部

  • 写回答

1条回答 默认 最新

  • devmiao 2014-12-14 23:31
    关注

    选取的样本是不是太小了?改大点看看。

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
编辑
预览

报告相同问题?

悬赏问题

  • ¥15 VAE代码如何画混淆矩阵
  • ¥15 求遗传算法GAMS代码
  • ¥15 雄安新区高光谱数据集的下载网址打不开
  • ¥66 android运行时native和graphics内存详细信息获取
  • ¥100 求一个c#通过CH341读取数据的Demo,能够读取指定地址值的功能
  • ¥15 rk3566 Android11 USB摄像头 微信
  • ¥15 torch框架下的强化学习DQN训练奖励值浮动过低,希望指导如何调整
  • ¥35 西门子博图v16安装密钥提示CryptAcquireContext MS_DEF_PROV Error of containger opening
  • ¥15 mes系统扫码追溯功能
  • ¥40 selenium访问信用中国