pub fn binomial_pmf(n: u32, k: u32, p: f32) -> f32
Binomial probability P(X = k) for n trials with success rate p