查看: 5611|回复: 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=OTY4MzN8NjMzYmY4OWY1MTYzZDJiYzNhZTA5MTAyYmY1MzhiZDB8MTc4NzQ2MzEzNw%3D%3D&request=yes&_f=.nb

zan

圆周率

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

相关帖子

使用道具举报

变动的永恒

0 主题	12 听众	534 积分

升级 78%

TA的每日心情

	奋斗 2015-4-18 18:22

签到天数: 163 天

[LV.7]常住居民III

自我介绍: 数学水平停留在高中水平了。

7^#

发表于 2014-8-22 20:00 |只看该作者 |招呼Ta 关注Ta

楼主我要了谢谢你了撒

我也说一句

发表

使用道具举报

蔡骏鹏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 9 E C5 Y4 Y: D J& r
不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。 t' w& K) Z' [- ~: h
Module[
. p* T+ x3 F3 b# m{i,x, ...

我也说一句

发表

使用道具举报

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
) S) }- c+ m2 {( k# y1 J+ J这种方法是比较粗略的，计算不精确

不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。
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]
]