pub fn geometric_pmf(k: u32, p: f32) -> f32
Geometric distribution: probability of first success at trial k