使用 Python 计算 π 值(5)

当前位置:

首页 > Python基础教程 >

使用 Python 计算 π 值(5)

输出结果：

哦，不，我们得到的是同样的精度。好吧，让我们试试第三个公式， Chudnovsky 算法，它看起来是这个样子：

wpid-826dc7788dba249ee86fc0135e06b035-2013-05-28-12-54.png

再一次，让我们看一下这个计算公式（假设我们有一个阶乘公式）。点击这里可下载用 python 实现的 Chudnovsky 公式。

下面是程序和输出结果：

def chudnovsky(n):

    pi=Decimal(0)

    k=0

    while k &lt; n:

        pi+=(Decimal(-1)**k)*(Decimal(factorial(6*k))/((factorial(k)**3)*(factorial(3*k)))*(13591409+545140134*k)/(640320**(3*k)))

        k+=1

    pi=pi*Decimal(10005).sqrt()/4270934400

    pi=pi**(-1)

    return pi

所以我们有了什么结论？花哨的算法不会使机器浮点世界达到更高标准。我真的很期待能有一个比我们用求和公式时所能得到的更好的精度。我猜那是过分的要求。如果你真的需要用PI，就只需使用math.pi变量了。然而，作为乐趣和测试你的计算机真的能有多快，你总是可以尝试第一个计算出Pi的百万位或者更多位是几。

栏目列表