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