fn beta_quantile(p: f64, a: f64, b: f64) -> f64
The x where the Beta(a, b) distribution reaches p, by bisection.
p