pub fn qfunc(x: f64) -> f64
Gaussian tail probability Q(x) = P(N(0,1) > x) = ½ erfc(x/√2).
Q(x) = P(N(0,1) > x) = ½ erfc(x/√2)