def _ppf(self, q, a):
        fac = special.gammainccinv(a, 1-abs(2*q-1))
        return where(q > 0.5, fac, -fac)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/dgamma_gen._ppf/original