为什么我的Chudnovsky算法程序给出了错误的结果？

Question

我正在尝试在Python中编码Chudnovsky算法。但是，当我运行代码时，它给了我很少的数字（-5.051212624421025E-55），这不是PI。我在中学时，我不认识有人可以帮助我。我在做什么错？

这是Chudnovsky公式的链接： https://levelup.gitconnected.com/generating-the-phalue-the-value-of-pi-pi-pi-pi-pi-pi-a-nonning-noming-number-forn-number-findecimals-decimals-decimals-place-place-place-plass-place-in-python-e93986bbbbbbbb474d

这是我的代码：

import math

def fact(exi):
    memory = exi
    for i in range(1, exi):
        memory *= i
    return memory

k = 10
s = 0

for i in range(0, k):
    a = -1^k
    b = fact(6*k)
    c = (545140134*k) + 13591409
    d = fact(3*k)
    e = (fact(k))^3
    f = (3 * k) + 3/2
    g = math.pow(640320, f)

    numerator = (a*b*c)
    denominator = (d*e*f)

    s += (numerator / denominator)

s *= 12
print(1 / s)

这是我的更新代码：

import math

def fact(exi):
    memory = exi
    for i in range(1, exi):
        memory *= i
    return memory

k = 17
s = 0

for i in range(1, k):
    a = (-1)**i
    b = fact(6*i)
    c = (545140134*i) + 13591409
    d = fact(3*i)
    e = (fact(i))**3
    f = (3 * i) + 3/2
    g = math.pow(640320, f)
    num = (a*b*c)
    den = (d*e*g)
    s += (num / den)
s *= 12
print(1 / s)

Answer 1

我看到两个错误：

1. 与 i （公式中的 q ），其中您当前在循环中使用 k 。在您的代码中 k 是 i 。
   

2. python中的启动操作员是 ** ，而不是^（它是位xor）。