def _pdf(self, x, a, b, c, z):
        Cinv = gam(a)*gam(b)/gam(a+b)*special.hyp2f1(c, a, a+b, -z)
        return 1.0/Cinv * x**(a-1.0) * (x-1.0)**(b-1.0) / (1.0+z*x)**c
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/gausshyper_gen._pdf/swapped operands of `-`