pub fn negative_binomial_pmf(n: u32, r: u32, p: f32) -> f32
Negative binomial: P(r-th success on n-th trial)