pub fn pareto_P(x: 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.