Function erfinv 

pub fn erfinv<F>(p: F) -> F
where
    F: Float,

The x with erf(x) = p. ±1 gives ±∞, and a magnitude beyond 1 gives NaN.