The Negative Binomial Distribution

Random: unsigned int gsl_ran_negative_binomial (const gsl_rng * r, double p, double n)
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^n (1-p)^k

Note that n is not required to be an integer.

Function: double gsl_ran_negative_binomial_pdf (unsigned int k, double p, double n)
This function computes the probability p(k) of obtaining k from a negative binomial distribution with parameters p and n, using the formula given above.