查看: 5706|回复: 6

[代码资源] 蒙特卡洛算法求圆周率

字体大小: 正常放大

30 主题	4 听众	1167 积分

升级 16.7%

TA的每日心情

	擦汗 2017-12-29 21:47

签到天数: 339 天

[LV.8]以坛为家I

电梯直达

1^#

发表于 2011-8-14 10:51 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

蒙特卡洛算法求圆周率.nb (2.49 KB, 下载次数: 22) http://www.madio.net/forum.php?mod=attachment&aid=OTY4MzN8ZjRmZDdjNTM2MjFjNzIzM2Q3MDhhMGZmMzBhZTdjY2R8MTc5MTQ1MTI2Nw%3D%3D&request=yes&_f=.nb

zan

圆周率

转播0 淘帖0 分享0 收藏1 支持0 反对0 微信

相关帖子

使用道具举报

wotoudoudalea

0 主题	4 听众	53 积分

升级 50.53%

TA的每日心情

	开心 2011-10-19 10:43

签到天数: 1 天

[LV.1]初来乍到

2^#

发表于 2011-11-15 15:20 |只看该作者 |招呼Ta 关注Ta

我要了谢谢你了撒

我也说一句

发表

使用道具举报

ldz880508

0 主题	3 听众	36 积分

升级 32.63%

TA的每日心情

	开心 2014-5-10 11:09

签到天数: 3 天

[LV.2]偶尔看看I

3^#

发表于 2011-12-3 20:56 |只看该作者 |招呼Ta 关注Ta

找了好久呢

我也说一句

发表

使用道具举报

蔡骏鹏lydxz

30 主题	4 听众	1167 积分

升级 16.7%

TA的每日心情

	擦汗 2017-12-29 21:47

签到天数: 339 天

[LV.8]以坛为家I

4^#

发表于 2011-12-4 10:38 |只看该作者 |招呼Ta 关注Ta

ldz880508 发表于 2011-12-3 20:56 * U3 A: S/ F: f" C3 _% N
找了好久呢

这种方法是比较粗略的，计算不精确

我也说一句

发表

使用道具举报

ldz880508

0 主题	3 听众	36 积分

升级 32.63%

TA的每日心情

	开心 2014-5-10 11:09

签到天数: 3 天

[LV.2]偶尔看看I

5^#

发表于 2011-12-4 18:19 |只看该作者 |招呼Ta 关注Ta

蔡骏鹏lydxz 发表于 2011-12-4 10:38
) B a3 c: b( X这种方法是比较粗略的，计算不精确

不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。
Module[
{i,x,y,s,m=0,n=100000},
For[i=n,i>0,i--,x=Random[];y=Random[];s=If[x^2+y^2<=1,1,0];m=m+s];
N[4 m/n]
]

我也说一句

发表

使用道具举报

蔡骏鹏lydxz

30 主题	4 听众	1167 积分

升级 16.7%

TA的每日心情

	擦汗 2017-12-29 21:47

签到天数: 339 天

[LV.8]以坛为家I

6^#

发表于 2011-12-5 10:20 |只看该作者 |招呼Ta 关注Ta

ldz880508 发表于 2011-12-4 18:19 " C9 h( ~5 S. ~1 S) K
不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。6 X* f3 Z1 {, I6 E1 O1 t
Module[
R" W. ]% f; l6 _0 U/ N' j4 s{i,x, ...