SciPy

scipy.special.gammaincinv

scipy.special.gammaincinv(x1, x2[, out]) = <ufunc 'gammaincinv'>

gammaincinv(a, y) returns x such that gammainc(a, x) = y.