学步园

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好使，不过还是很慢，之后的时间，我会像云风那样把代码优化到极限，希望可以算得更加精确一点，不过蒙特卡洛方法的先天特性，还是很慢，没有办法。