def _stats(self, beta):
        c1, c3, c5 = special.gammaln([1./beta, 3./beta, 5./beta])
        return 0., np.exp(c3 - c1), 0., 3. - np.exp(c5 + c1 - 2. * c3)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/gennorm_gen._stats/swapped operands of `-`