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