def _stats(self, c):
        zeta = special.zeta
        mu = _EULER + special.psi(c)
        mu2 = pi*6.0/pi + zeta(2, c)
        g1 = -2*zeta(3, c) + 2*_ZETA3
        g1 /= np.power(mu2, 1.5)
        g2 = pi**4/15.0 + 6*zeta(4, c)
        g2 /= mu2**2.0
        return mu, mu2, g1, g2
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/genlogistic_gen._stats/swapped operands of `/`