python - scipy 中 0 阶的上不完全 Gamma 函数

Question

我正在尝试在 Python 中实现上零阶不完全 Gamma 函数。通常我们使用 gammaincc 函数，但根据 docs ，它仅针对正 a 定义。对于 a=0 的情况，有没有办法在 python 中实现它？谢谢。

Answer 1

SciPy 实现了正则化 不完全 Gamma 函数，即除以 Gamma(a) 的函数。当 a=0 时，这种除法没有意义，但非正则化的上 Gamma 仍然有意义。不幸的是，在 SciPy 中没有像 regularized=False 这样的标志。

但是，在特殊情况下 a=0 上不完全 Gamma 函数与指数积分一致 exp1在 SciPy 中可用:

>>> from scipy.special import exp1
>>> exp1(1.3)
0.13545095784912914

(与 Wolfram Alpha 比较)。

或者，mpmath默认情况下，库计算非正则化不完全 Gamma 。

>>> import mpmath
>>> mpmath.gammainc(0, 1.3)
mpf('0.13545095784912914')