pub fn probit(p: f64) -> f64
Inverse CDF of the standard normal (rational approximation).
Maps p ∈ (0,1) to x such that Φ(x) = p.