pub fn erfcinv<F: Float>(q: F) -> F
The x with erfc(x) = q. 0 gives ∞, 2 gives -∞, and q outside 0 to 2 gives NaN.