scipy.special.gammainc

scipy.special.gammainc(x1, x2[, out])

y=gammainc(a,x) returns the incomplete gamma integral defined as 1 / gamma(a) * integral(exp(-t) * t**(a-1), t=0..x). Both arguments must be positive.