我是一个蟒蛇初学者,我想计算pi.我尝试使用Chudnovsky算法,因为我听说它比其他算法更快.
这是我的代码:
from math import factorial
from decimal import Decimal, getcontext
getcontext().prec=100
def calc(n):
t= Decimal(0)
pi = Decimal(0)
deno= Decimal(0)
k = 0
for k in range(n):
t = ((-1)**k)*(factorial(6*k))*(13591409+545140134*k)
deno = factorial(3*k)*(factorial(k)**3)*(640320**(3*k))
pi += Decimal(t)/Decimal(deno)
pi = pi * Decimal(12)/Decimal(640320**(1.5))
pi = 1/pi
return pi
print calc(25)
出于某种原因,与可接受的值相比,此代码产生pi的vakue,最多只有15位小数.我试图通过提高精度值来解决这个问题; 这会增加位数,但只有前15位仍然准确.我尝试改变计算算法的方式,但它也没有用.所以我的问题是,是否可以对此代码进行一些操作以使其更准确或者我是否必须使用其他算法？我很感激这方面的帮助,因为我不知道如何在python中使用这么多位数.我希望能够控制程序确定和显示的(正确)数字的数量 - 无论是10,100,1000等.