def _logpdf(self, x, c):
        return (log(x+1) - (1-x)**2 / (2.0*x*c**2) - log(2*c) -
                0.5*(log(2*pi) + 3*log(x)))
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/fatiguelife_gen._logpdf/swapped operands of `-`