pub fn poisson(r: &mut Rng, mu: f64) -> u32
This function returns a random integer from the Poisson distribution with mean mu. The probability distribution for Poisson variates is,
p(k) = {\mu^k \over k!} \exp(-\mu)
for k >= 0.