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

Question

我正在尝试在Python中实现零阶的上部不完整伽玛函数。通常我们使用printf("%p %p %p %p %p",p,p+1,p+2,p+3,p+4); 函数，但是根据docs，它仅定义为正 a 。对于a = 0的情况，有什么办法可以在python中实现吗？谢谢。

Answer 1

SciPy实现了 regularized 不完全伽玛函数，该函数除以Gamma（a）。当a = 0时，这种除法是没有意义的，但是非正规化的上伽玛仍然有意义。不幸的是，SciPy中没有像regularized=False这样的标志。

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

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

（与Wolfram Alpha比较）。

或者，mpmath库默认情况下计算非正则不完整伽玛。

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