计算圆周率近似值
方法:
1.拉马努金法计算圆周率近似值:
“数学家拉马努金(Srinivasa Ramanujan)找到了一个无限序列,可以用来生成π的数值近似值:
编写一个函数 ,使用这个公式计算并返回π的近似估计。它应当使用一个循环来计算求和的每一项,直到最后一项的值小于1e-15 (这是Python对的标记法)。可以通过和math.pi 比较来检查计算的结果。”
代码实现
#拉马努金计算圆周率
import math
#求阶乘的函数
def factorial(n):
if n==0:
return 1
else:
return n*factorial(n-1)
#计算π值的函数
def pi():
sum =0
k=0
f=2*(math.sqrt(2))/9801
while True:
fz = (26390*k + 1103)*factorial(4*k) #求和项分子
fm = (396**(4*k))*((factorial(k))**4) #求和项分母
t = f*fz/fm
sum += t
if t<1e-15: #最后一项小于10^(-15)时跳出循环
break
k += 1 #更新k值
return 1/sum
pri