— Function File: beta_inv (x, a, b)

For each component of x, compute the quantile (the inverse of the CDF) at x of the Beta distribution with parameters a and b.