计算思维:蒙特卡罗方法求“圆周率”

#CalPiV2.py
from random import random
from time import perf_counter
DARTS = 1000*1000
hits=0.0
start = perf_counter()
for i in range(1,DARTS+1):
x,y = random(),random()
dist=pow(x**2+y**2,0.5)
if dist <=1.0:
hits=hits+1
pi=4*(hits/DARTS)
print("圆周率值是:{}".format(pi))
print("运行时间是:{:.5f}s".format(perf_counter()-start))

猜你喜欢

转载自www.cnblogs.com/oycc2000/p/11230552.html