pub fn pareto_pdf(x: f64, a: f64, b: f64) -> f64
This function computes the probability density p(x) at x for a Pareto distribution with exponent a and scale b, using the formula given above.