def _stats(self, a):
        ea = exp(a)
        mu2 = 2.*ea/(ea-1.)**2
        mu4 = 2.*ea*(ea**2+10.*ea+1.) / (ea-1.)**4
        return 0., mu2, 0., 3. - mu4/mu2**2
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/dlaplace_gen._stats/swapped operands of `-`