我写了这个程序来计算π。它给了我50个小数位。我怎样才能更精确?
# pi.py - arctan(1) * 4 = pi
from math import *
from decimal import *
err = 0.0000000000000001
def arctan(n, err):
"""Uses Gregory's formula for calculating atan."""
temp = n
atan = 0
i = 3
while (abs(atan - n) > err):
atan = n
n = n - (pow(temp, i)/i) + ((pow(temp, i + 2)) / (i + 2))
i += 4
return n
def euler(a, b, err):
"""Uses Euler's formula and fibonacci numbers."""
euler = 0
temp = 5
while (abs(temp - euler) > err):
temp = euler
euler += arctan(1/b, err)
a = b + a
b = b + a
return euler
pi = euler(1, 2, err) * 4
print(Decimal.from_float(pi))
必须将Decimal prec设置为更高的值。见this example thread。official python site有很多例子。
此外,您应该使用十进制进行所有计算,而不仅仅是最后一步。否则你就得不到所需的精度。
或者看看支持任意精度浮点运算的mpmath。
而且,有些算法会无限期地一次生成一位数的pi(而不需要非常高的精度要求)。
相关问题 更多 >
编程相关推荐