Process Performance Determination in C# Part 4

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace ExcelNORMSDIST
{
    public class Excel
    {
        static void Main(string[] args)
        {
            Console.WriteLine("Value for NORMSDIST(0) is: {0}.", NORMSDIST(0));
            Console.WriteLine("Value for NORMSDIST(0.2) is: {0}.", NORMSDIST(0.2));
            Console.WriteLine("Value for NORMSDIST(-0.2) is: {0}.", NORMSDIST(-0.2));
            Console.WriteLine("Value for NORMSDIST(-1) is: {0}.", NORMSDIST(-1));
            Console.WriteLine("Value for NORMSDIST(-2) is: {0}.", NORMSDIST(-2));
            Console.WriteLine("Value for NORMSDIST(-3) is: {0}.", NORMSDIST(-3));
            Console.WriteLine("Value for NORMSDIST(-4) is: {0}.", NORMSDIST(-4));
            Console.WriteLine("Value for NORMSDIST(-5) is: {0}.", NORMSDIST(-5));
            Console.WriteLine("Value for NORMSDIST(-7) is: {0}.", NORMSDIST(-7));
            Console.ReadKey();
        }

        public static double NORMSDIST(double z)
        {
            double sign = 1;
            if (z < 0) sign = -1;
            return 0.5 * (1.0 + sign * erf(Math.Abs(z) / Math.Sqrt(2)));
        }

        private static double erf(double x)
        {
            //A&S formula 7.1.26
            double a1 = 0.254829592;
            double a2 = -0.284496736;
            double a3 = 1.421413741;
            double a4 = -1.453152027;
            double a5 = 1.061405429;
            double p = 0.3275911;
            x = Math.Abs(x);
            double t = 1 / (1 + p * x);
            //Direct calculation using formula 7.1.26 is absolutely correct
            //But calculation of nth order polynomial takes O(n^2) operations
            //return 1 - (a1 * t + a2 * t * t + a3 * t * t * t + a4 * t * t * t * t + a5 * t * t * t * t * t) * Math.Exp(-1 * x * x);

            //Horner's method, takes O(n) operations for nth order polynomial
            return 1 - ((((((a5 * t + a4) * t) + a3) * t + a2) * t) + a1) * t * Math.Exp(-1 * x * x);
        }
    }//end class program
}