数学建模社区-数学中国

标题: 这个程序为啥求不出pi [打印本页]

作者: 三世清华 时间: 2013-8-27 10:59
标题: 这个程序为啥求不出pi
flag=0;
N=1000000;
x=rand(1,N)*2-1;
y=rand(1,N)*2-1;
m=sqrt(x.^2+y.^2);
if m<=1
flag=flag+1
end

pai=4* (flag/N)

新建 Microsoft Word 文档.doc

23.5 KB, 下载次数: 0, 下载积分: 体力 -2 点

作者: 净心、精心 时间: 2013-8-27 12:57
m<=1 错误 m是一个矩阵怎么能比较大小呢

作者: 袁海亮 时间: 2013-8-27 18:22
楼主你的程序我运行出
pai =

0

作者: 三世清华 时间: 2013-8-28 07:55

袁海亮发表于 2013-8-27 18:22
楼主你的程序我运行出
pai =

亲，答案不是pai的答案，老师显示为0，这是为啥呢？秦？？？

作者: kong1234 时间: 2013-8-28 09:32
N=10000000;
x=rand(1,N)*2-1;
y=rand(1,N)*2-1;
m=sqrt(x.^2+y.^2);
f=find(m<=1);
flag=length(f);
pai=4*(flag/N);

亲测可用

作者: kong1234 时间: 2013-8-28 09:33
而且我发现结果也比较随机，，

作者: magic2728 时间: 2013-8-28 22:01
m应该得加上m（i），然后再循环

作者: 三世清华 时间: 2013-8-29 07:50
谢谢，非常高感谢！！！！！！！！！！！！！！

作者: 三世清华 时间: 2013-8-29 07:52

三世清华发表于 2013-8-28 07:55
亲，答案不是pai的答案，老师显示为0，这是为啥呢？秦？？？

亲，我也不知道，这不网上求助吗？？

作者: gaoyingbetty 时间: 2013-8-29 14:29

作者: gaoyingbetty 时间: 2013-8-29 14:30

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)