![8c95a6620a82f0f67c25cf3d4898b380.png](https://i-blog.csdnimg.cn/blog_migrate/04aabb1bbf19f868647cd2c363752d43.jpeg)
在 上图 中, 圆的 半径 为 1, 则 右上角 正方形的 面积 为 1, 右上角 1/4 圆的 面积 为 1/4 Π, 假设 x, y 都 在 [0, 1) 内 均匀随机 取值, 如果 x*x + y*y <=1, 则点 落入 1/4 的圆的范围内, 设 p= 点 x,y 落入 1/4 圆内的 概率, 则有 p=1/4圆的面积 /正方形面积 = 1/4 Π/1, 则有 1/4 Π=p, 也就是 Π=4P, 就是 说, 随机 在 [0, 1) 取值 x,y, 计算 x*x + y*y <=1 的 概率, 再乘以4, 就可以 算出 圆周率的 值。
下边 是 python 程序。
import random
def compute_pi( ):
count=0
incirclecount=0
divisor=10
while(1):
count+=1
x=random.random( )
y=random.random( )
if(x*x+y*y<=1):
incirclecount+=1
if(count%divisor==0):
pi=4*incirclecount/count
print ("n now the loop count %d pi is %f"%(count,pi))
divisor*=10
compute_pi( )
下边是程序运行结果
now the loop count 10 pi is 3.200000
now the loop count 100 pi is 3.280000
now the loop count 1000 pi is 3.084000
now the loop count 10000 pi is 3.144000
now the loop count 100000 pi is 3.142560
now the loop count 1000000 pi is 3.144048
now the loop count 10000000 pi is 3.142086
now the loop count 100000000 pi is 3.141474
可以看出,随着 循环次数的 增加 ,逐渐接近 圆周率,经过 1亿次循环以后, 误差已经在小数点后 4位,也就是小于 1/30000 。 这个方法,是不是很有意思