Numerical Analysis Seminar

Wednesday, October 15, 2008

Speaker: Jean-Marie Linhart, Texas A&M University
Title: Computing the Logarithm of the Normal Distribution
Time: 3:00-4:00 pm
Place: Blocker 627

Abstract

We present and compare three C functions to compute the logarithm of the cumulative standard normal distribution. The first is a new algorithm derived from Algorithm 304's calculation of the standard normal distribution via a series or continued fraction approximation, and it is good to the accuracy of the machine. The second is based on Algorithm 715's calculation of the standard normal distribution via rational Chebyshev approximation. This is related to and an improvement on the algorithm for the logarithm of the normal distribution available in the software package R. The third is a new and simple algorithm which uses the compiler's implementation of the error function and complement of the error function to compute the log of the normal distribution.

Last revised: 09/11/08 By: abnersg@math