def EvalPoissonPmf(k, lam):
<EOL>return lam ** k * math.exp(-lam) / math.factorial(k)<EOL>
Computes the Poisson PMF.

    k: number of events
    lam: parameter lambda in events per unit time

    returns: float probability