pub fn binomial_cdf(k: u32, n: u32, p: f64) -> f64
Binomial CDF via direct sum (small n only; for large use normal approx in caller).