Random: unsigned int gsl_ran_logarithmic(const gsl_rng * r, double p)
This function returns a random integer from the logarithmic
distribution. The probability distribution for logarithmic random variates
is,
for @c{$k \ge 1$}
@math{k >= 1}.
Function: double gsl_ran_logarithmic_pdf(unsigned int k, double p)
This function computes the probability @math{p(k)} of obtaining k
from a logarithmic distribution with probability parameter p,
using the formula given above.