pub fn quantile(p: f64, mu: f64, sigma: f64) -> f64
Inverse normal cdf with mean mu, std-dev sigma.
mu
sigma