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

Question

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

Answer 1

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

但是，在特殊情况下 a=0 上不完全伽玛函数与指数积分 exp1 一致，后者在 SciPy:

中可用

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

（与 Wolfram Alpha 相比）。

或者，mpmath 库默认计算非正则化不完全伽马。

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

Upper Incomplete Gamma Function of order 0 in scipy