SciPy

scipy.special.gammainc

scipy.special.gammainc(a, x) = <ufunc 'gammainc'>

Incomplete gamma function

Defined as:

1 / gamma(a) * integral(exp(-t) * t**(a-1), t=0..x)

a must be positive and x must be >= 0.