Random: unsigned int gsl_ran_hypergeometric(const gsl_rng * r, unsigned int n1, unsigned int n2, unsigned int t)
This function returns a random integer from the hypergeometric
distribution. The probability distribution for hypergeometric
random variates is,
where @math{C(a,b) = a!/(b!(a-b)!)}. The domain of @math{k} is
@math{max(0,t-n_2), ..., max(t,n_1)}.
Function: double gsl_ran_hypergeometric_pdf(unsigned int k, unsigned int n1, unsigned int n2, unsigned int t)
This function computes the probability @math{p(k)} of obtaining k
from a hypergeometric distribution with parameters n1, n2,
n3, using the formula given above.