def _entropy(self, lam):
        def integ(p):
            return log(pow(p, lam-1)+pow(p-1, lam-1))
        return integrate.quad(integ, 0, 1)[0]
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/tukeylambda_gen._entropy/swapped operands of `-`