查看: 5520|回复: 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=OTY4MzN8YTUwNTY2OWFhMjVjNjc4ZjQ1ZDgxMGFjZWUzNGQ5MmF8MTc4MzU2NzExMQ%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
7 K) G- h$ H4 v找了好久呢

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

我也说一句

发表

使用道具举报

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
3 U8 x+ a2 _6 o T, M9 @5 h4 Z这种方法是比较粗略的，计算不精确

不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。
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
' t. [8 P8 T# m, ^/ R不精确是应该的吧，蒙特卡罗算法求的就是近似值嘛，我把你的算法稍微改了下写法，嘿嘿。$ x; Z; g4 Z" k; L& L$ @. N
Module[' |" e7 c# y9 P& j L4 T
{i,x, ...