python - 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')