一、蒙特卡洛求圆周率
向区域内随即撒点
当点的数目足够多时,落在圆的点数目与在正方形点数目成正比
即圆的面积和正方形的面积成正比
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: #在圆内部加1
hits=hits+1
pi=4*(hits/darts)
print("圆周率值是:{}".format(pi))
print("运行时间是:{:.5f}s".format(perf_counter()-start))
*****************************************************
圆周率值是:3.141804
运行时间是:1.11708s