def _stats(self, s):
        p = exp(s*s)
        mu = sqrt(p)
        mu2 = p*(p-1)
        g1 = sqrt((p-1))*(2+p)
        g2 = np.polyval([1, 2, 3, 0, -6.0], p)
        return mu, mu2, g1, g2
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/lognorm_gen._stats/original