Верхняя Неполная Гамма-функция порядка 0 в scipy

Question

Я пытаюсь реализовать верхнюю неполную гамма-функцию нулевого порядка в Python. Обычно мы используем функцию gammaincc, но в соответствии с документами она определена только для положительных значений a . Есть ли способ реализовать это в python для случая = 0? Спасибо.

Answer 1

SciPy реализует регуляризованную неполную гамма-функцию *1001* с функцией деления по гамме (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')