蒙特卡洛方法计算圆周率

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/qq_38669138/article/details/79823031
#蒙特卡洛方法计算圆周率
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

猜你喜欢

转载自blog.csdn.net/qq_38669138/article/details/79823031