python中不完整的gamma函数？

Question

在scipy.special.gammainc不能取负值的第一个参数.在python中是否还有其他实现？我可以肯定地进行手动集成,但我想知道是否存在已经存在的良好替代方案.

Correct result: 1 - Gamma[-1,1] = 0.85

Use Scipy: scipy.special.gammainc(-1, 1) = 0

谢谢.

Answer 1

每当我需要特殊功能时,我通常会达到mpmath,而且我并不太关心性能.(尽管它在许多情况下的表现都非常好.)

例如:

>>> import mpmath
>>> mpmath.gammainc(-1,1)
mpf('0.14849550677592205')
>>> 1-mpmath.gammainc(-1,1)
mpf('0.85150449322407795')
>>> mpmath.mp.dps = 50 # arbitrary precision!
>>> 1-mpmath.gammainc(-1,1)
mpf('0.85150449322407795208164000529866078158523616237514084')