计算趋于无穷大

Question

我试图计算一个数的伽马函数，它将趋于无穷大。例如，当我使用maple时，它会返回正确的答案（2.57e1133（是的，它很大））。我尝试使用十进制，但0成功。有解决办法吗？提前谢谢

代码

import scipy as sp                             
from scipy.special import gamma
from decimal import Decimal
from scipy import special


def teste(k):
    Bk2 = Decimal(gamma((1/(2*k))+(3/4)))
    return Bk2

print(teste(0.001))

结果

无穷大

Answer 1

我使用gammaln来防止溢出，然后将结果转换为十进制。最后，我使用指数函数来恢复预期结果

import scipy as sp                             
from scipy.special import gammaln
from decimal import Decimal
from scipy import special
from numpy import exp

def teste(k):
    Bk2 = exp(Decimal(gammaln((1/(2*k))+(3/4))))
    return Bk2

print(teste(0.001))