#蒙特卡洛方法计算圆周率
from random import random
s=1000*1000#撒点的个数
hist=0.0#圆内的点的初值
for i in range(0,s+1):
x=random()
y=random()
if pow(x**2+y**2,0.5)<=1:
hist+=1
print("pi is :",4*hist/s)
运行结果:
== RESTART: C:/Users/lenovo/AppData/Local/Programs/Python/Python36/calpi.py ==
pi is : 3.14204