pub fn pareto(r: &mut Rng, a: f64, b: f64) -> f64
This function returns a random variate from the Pareto distribution of order a. The distribution function is,
p(x) dx = (a/b) / (x/b)^{a+1} dx
for x >= b.