Published by the Foundation for Open Access Statistics Editors-in-chief: Bettina Grün, Torsten Hothorn, Rebecca Killick, Edzer Pebesma, Achim Zeileis    ISSN 1548-7660; CODEN JSSOBK
Authors: George Marsaglia
Title: Evaluating the Normal Distribution
Abstract: This article provides a little table-free C function that evaluates the normal distribution with absolute error less than 8 x 10 -16 . A small extension provides relative error near the limit available in double precision: 14 to 16 digits, the limits determined mainly by the computer's ability to evaluate exp(-t) for large t. Results are compared with those provided by calls to erf or erfc functions, the best of which compare favorably, others do not, and all appear to be much more complicated than need be to get either absolute accuracy less than 10-15 or relative accuracy to the exp()-limited 14 to 16 digits. Also provided: A short history of the error function erf and its intended use, as well as, in the "browse files" attachment, various erf or erfc versions used for comparison.

Page views:: 29923. Submitted: 2004-06-05. Published: 2004-07-18.
Paper: Evaluating the Normal Distribution     Download PDF (Downloads: 32231)
Supplements: sources.c: C code Download (Downloads: 1466; 13KB)

DOI: 10.18637/jss.v011.i04

This work is licensed under the licenses
Paper: Creative Commons Attribution 3.0 Unported License
Code: GNU General Public License (at least one of version 2 or version 3) or a GPL-compatible license.