Function rgsl::randist::pareto::pareto

pub fn pareto(r: &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.