def _pmf(self, k, lambda_):
        fact = (exp(-lambda_)-1)
        return fact*exp(-lambda_*k)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/planck_gen._pmf/swapped operands of `-`