使用python和蒙特卡洛方法算π

来源:互联网 发布:js中的事件冒泡 编辑:程序博客网 时间:2024/06/01 14:36
Python语言: Codee#24529
01 # -*- coding: cp936 -*-
02 #monta cola蒙特卡洛方法算PI
03 #younger.x.shen@gmail.com
04 #2011/12/08
05
06 import random
07
08 def getThePi(randsum):
09
10     inner = 0.0
11     outer = 0.0
12
13     for i in range(randsum) :
14         temp1 = random.random()
15         temp2 = random.random()
16
17         if temp1 ** 2 + temp2 ** 2 <= 1 :
18             inner += 1
19
20     perc = inner / randsum
21     print(inner)
22     circle = 4 * perc
23     # s = PI * R * R
24     # PI = s / (R*R)
25
26     return circle
不知道为何高亮这么不好用,插件又装不上去,很郁闷,写多了java ee ,这点代码我也是敲了10多分钟的,算法功力大退,虽然是差不多第一次使用python吧,感觉用python描述算法确实是比c好使,不过还是很慢,之后的时间,我会像云风那样把代码优化到极限,希望可以算得更加精确一点,不过蒙特卡洛方法的先天特性,还是很慢,没有办法。
原创粉丝点击