1

我们最近深入研究了微积分中的无限级数，话虽如此，我玩得很开心。我在python中导出了我自己的inverse tan infinte系列，并设置为1得到pi/4*4得到pi。我知道这不是最快的算法，所以请不要讨论我的算法。我想讨论的是如何在 python 中表示非常非常小的数字。我注意到的是，当我的程序迭代系列时，它会在小数点后 20 位（给或取）处停止。我尝试使用十进制模块并且只推到大约 509。我想要一个无限（几乎）表示。

有没有办法做这样的事情？我认为没有数据类型能够处理如此庞大的数据，但如果你能告诉我一个解决方法，我将非常感激。

4

1 回答 1

1

Python 的十进制模块要求您指定“上下文”，这会影响表示的精确程度。

对于这类事情，我可能会推荐gmpy2 - 您可以对有理数（任意精度）进行计算并在最后一步转换为十进制。

这是一个示例 - 根据需要替换您自己的算法：

import gmpy2
# See https://gmpy2.readthedocs.org/en/latest/mpfr.html
gmpy2.get_context().precision = 10000
pi = 0
for n in range(1000000):
    # Formula from http://en.wikipedia.org/wiki/Calculating_pi#Arctangent
    numer = pow(2, n + 1)
    denom = gmpy2.bincoef(n + n, n) * (n + n + 1)
    frac = gmpy2.mpq(numer, denom)
    pi += frac
    # Print every 1000 iterations
    if n % 1000 == 0:
        print(gmpy2.mpfr(pi))
于 2013-05-18T20:12:44.987 回答