pub fn weibull(r: &mut Rng, a: f64, b: f64) -> f64
This function returns a random variate from the Weibull distribution. The distribution function is,
p(x) dx = {b \over a^b} x^{b-1} \exp(-(x/a)^b) dx
for x >= 0.