查看: 5347|回复: 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=OTY4MzN8NTUyZmM0NGY4YTM4NTc2NTRhYmRjMTIyOTlmZDk1NDV8MTc2NTQ4Mjg1Ng%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
( C& D' M& J) L( j; \" {找了好久呢

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

我也说一句

发表

使用道具举报

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
( d, G6 W& `7 u9 [* y/ R8 F这种方法是比较粗略的，计算不精确

不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。
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 " u- w: ^3 u' l6 Y5 Q8 x; M& _
不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。) d" q' h& U, A7 _. r9 E
Module[
# @/ m2 g0 u6 ^" j0 ]( l7 k" W{i,x, ...