python - 欧拉-马歇罗尼常数

Question

在编程中，我只使用整数。不过这次要进行一些计算。我需要计算Euler-Mascheroni Constant γ .最多 n 位小数。{虽然 n ∈ [30, 150]对我来说已经足够了。

[x] = gif(x) = math.floor(x)

但是，我怀疑精度
Numerical Algorithm
我需要使用 Python 获得更高的准确度。

Answer 1

来自 French Wikipedia discussion page , 近似到 6 位小数:

import math as m
EulerMascheroniApp = round( (1.-m.gamma(1+1.e-8))*1.e14 )*1.e-6
print(EulerMascheroniApp)
# 0.577216 

这个常量也可以在 sympy 模块中使用，名称为 EulerGamma :

>>> import sympy
>>> sympy.EulerGamma
EulerGamma
>>> sympy.EulerGamma.evalf()
0.577215664901533
>>> - sympy.polygamma(0,1)
EulerGamma
>>> sympy.stieltjes(0)
EulerGamma
>>> sympy.stieltjes(0, 1)
EulerGamma

文档:

math.gamma ;

sympy.EulerGamma ;

sympy.functions.special ;

sympy: numerical evaluation .

在最后一个文档链接上，您可以找到有关如何更精确地评​​估常量的更多信息，如果默认值为 .evalf()是不足够的。
如果您仍然想自己计算常数作为练习，我建议将您的结果与 sympy 的常数进行比较，以检查准确性和正确性。