pub fn beta_Qinv(Q: f64, a: f64, b: f64) -> f64
This function computes the cumulative distribution functions P(x), Q(x) and their inverses for the beta distribution with parameters a and b.