1. 程式人生 > >蒙特卡洛方法計算圓周率

蒙特卡洛方法計算圓周率

#蒙特卡洛方法計算圓周率
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