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

Question

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

Answer 1

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

然而，在特殊情况下，上不完全伽马函数与SciPy 中可用a=0的指数积分一致：exp1

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

（与Wolfram Alpha 相比）。

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

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