def _cdf(self, x, a, c):
        xc = x**c
        val1 = special.gammainc(a, xc)
        val2 = special.gammaincc(a, xc)
        return np.where(c > 0, val1, val2)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/gengamma_gen._cdf/original