查看: 5443|回复: 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=OTY4MzN8M2E1ZWQyZDg2NTBkZTA1MjhkMDkwMjEwYTBkNDJmNWF8MTc3NTcxOTA2NA%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 - [5 ~4 G5 `& r( U5 @* ?
找了好久呢

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

我也说一句

发表

使用道具举报

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 4 p5 U1 G/ U4 P. }2 ~0 n
这种方法是比较粗略的，计算不精确

不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。
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 9 ~; z* @2 `2 C; O2 j: i
不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。
5 Q; N# V0 g: Z- ^* D+ ^: X$ B# uModule[2 s/ M2 C+ H6 }. p. i! R2 Z' |2 ~
{i,x, ...