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