cosmo norminv hdr

function y = cosmo_norminv(p, mu, sd)