pub fn pareto_Qinv(Q: f64, a: f64, b: f64) -> f64
This function computes the cumulative distribution functions P(x), Q(x) and their inverses for the Pareto distribution with exponent a and scale b.