def _pdf(self, x, a, b):
        x2 = x*x
        trm = _norm_pdf(a + b * log(x + sqrt(x2+1)))
        return b*sqrt(x2+1.0)*trm/1.0
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/johnsonsu_gen._pdf/swapped operands of `/`