pub fn poisson_pmf(k: u64, lambda: f64) -> f64
Poisson probability mass function: P(k;λ) = λᵏ · e⁻λ / k!