public static double Ierfc(double y) { double s = Normal.Inverse(-0.5 * y + 1); double r = s * Math.Sqrt(2) / 2.0; return r; }