pub fn erfcinv(y: f64) -> f64
Inverse of erfc on (0, 2), by bisection (±∞ at the ends).
erfc
(0, 2)
±∞