数学方法,套用数学公式:
pi=0
N=100
for k in range(N):
pi += 1/pow(16,k)*(\
4/(8*k+1)-2/(8*k+4)-\
1/(8*k+5)-1/(8*k+6))
print("圆周率值是:{}".format(pi))
注:在一行中代码写不开,可以分成多行来写,只需用“\”来连接即可
蒙特卡罗方法计算圆周率:
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)
#用点到圆心的距离和1相比来确定,点是否在圆内
if dist<=1.0:
hits=hits+1
pi=4*(hits/DARTS)
print("圆周率的值是:{}".format(pi))
print("运行时间是:{:.5f}s".format(perf_counter()-start))