python计算的实际pi值是

Question

我的代码输出了不正确的pi值。。。你知道吗

$ python pi_Chudnovsky.py 100 1000
3.1415926535896751768742638018492028796422794574832770797807531304291682072583579638804596101871423663813700249451933537852

我认为我的算法是正确的。为了确认一下，请告诉我是不是错了。你知道吗

#!/usr/bin/python

import sys
import gmpy2
from gmpy2 import mpfr

n = int(sys.argv[1])
l = int(sys.argv[2])

gmpy2.get_context().precision = n * 4

p0 = mpfr(1)
q0 = mpfr(1)
t0 = mpfr(0)

p = mpfr(0)
q = mpfr(0)
t = mpfr(0)
a = mpfr(0)

for loop in range(1, l+1, 1):
    p = (2*loop - 1)*(6*loop - 1)*(6*loop - 5)
    q = (loop**3)*(640320**3)/24
    a = (13591409 + 545140134*loop)
    if loop % 2 == 0:
        a = -1 * a
    else:
        a = 1 * a

    t = t0*q + a*p

    p = p*p0
    q = q*q0

    p0 = p
    q0 = q
    t0 = t

pi = 640320 * gmpy2.sqrt(640320)*q
pi = pi/(12*(t + 13591409*q))

print(pi)

预期产量：

3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066486