pub fn inv_beta_reg(a: f64, b: f64, x: f64) -> f64
Computes the inverse of the regularized incomplete beta function