Random: unsigned int gsl_ran_binomial(const gsl_rng * r, double p, unsigned int n)
This function returns a random integer from the binomial distribution,
the number of successes in n independent trials with probability
p. The probability distribution for binomial variates is,
for @c{$0 \le k \le n$}
@math{0 <= k <= n}.
Function: double gsl_ran_binomial_pdf(unsigned int k, double p, unsigned int n)
This function computes the probability @math{p(k)} of obtaining k
from a binomial distribution with parameters p and n, using
the formula given above.