Function rgsl::randist::negative_binomial::negative_binomial

pub fn negative_binomial(r: &Rng, p: f64, n: f64) -> u32

This function returns a random integer from the negative binomial distribution, the number of failures occurring before n successes in independent trials with probability p of success. The probability distribution for negative binomial variates is,

p(k) = {\Gamma(n + k) \over \Gamma(k+1) \Gamma(n) } pⁿ (1-p)^k

Note that n is not required to be an integer.