1 import random
 2 i = 0
 3 numin = 0
 4 numall = 0
 5 while 1:
 6     i += 1
 7     numall += 1
 8     x = random.uniform(-5,5)
 9     y = random.uniform(-5,5)
10     if y*y+x*x <= 25:
11         numin += 1
12     if i % 50000 == 0:
13         print(4.0*numin/numall)

 

　　蒙特卡洛法计算圆周率计算原理直观。如图所示（图来自网络），假设有一个正方形，我们做它的内切圆，设内切圆的半径是r，那么内切圆面积是πr*r，正方形面积是4*r*r，圆的面积/正方形面积=π/4。随机向正方形平面射点，当点的数量足够多，我们可以看做落在圆上的点的数量代表圆的面积、所有的点代表正方形面积，因此，落在圆上的点的数量/所有的点数量=π/4。π因此得以算出。

　　此方法效率很低，本人在单线程的情况下运行5小时才精确到3.14159。