在 Python 2 中计算 PI，如何停止截断答案 [重复]

Question

请放轻松，我已经学习 Python 一个星期了！ 我想我会尝试使用 Rumanujan 公式计算 Pi。我相信我能够正确编码。 我的答案是截断，我希望它用 200 dp 表示。在 C 语言中，我可能会使用 malloc 来执行此操作，但我知道 Python 不能那样工作。 我想从中学到的一点是：截断是由表示浮点数的限制引起的，如果是的话，是否可以修复？

谢谢。

import math
from decimal import *

getcontext().prec = 200
def iterate(n):
    sum = 0
    Decimal(sum)
    sum = (math.factorial(4*n))
    sum = (sum/math.pow(math.factorial(n), 4))
    sum = sum*((26390*n +1103)/math.pow(396, (4*n)))
    return sum
ans=0
Decimal(ans)    
print "Choose the number of iterations:\n"
itnum = int(raw_input())

for n in range (0, itnum+1):
    this_iteration = 0
    Decimal(this_iteration)
    this_iteration = iterate(n)
    ans = ans + this_iteration

ans =  ans*(math.pow(8, 0.5)/9801)
ans = 1/ans
print "%.200f" % ans

Answer 1

你的sn-p

sum = 0
Decimal(sum)

将sum 设置为int 0，并计算并丢弃Decimal 等效项。改为使用 assignment 语句：

sum = Decimal(0)

接下来，您需要确保每个中间结果也正确转换为 Decimal（默认情况下 floats 不是）。

就个人而言，我建议改用gmpy2，但是，我有偏见:-)。