Computer Evaluation of the Normal and Inverse Normal Distribution Functions

An algorithm involving polynomial approximations for evaluation of the normal distribution function is presented which may be implemented in fast and accurate computer programs of moderate length. Using this approximation and simple two-step Newton-Raphson iteration, the evaluation of the inverse normal distribution function is achieved with reasonable accuracy.