A sharp Pólya-based approximation to the normal cumulative distribution function

We study an expansion of the cumulative distribution function of the standard normal random variable that results in a family of closed form approximations that converge at 0. One member of the family that has only five explicit constants offers the absolute error of 5.79·10−6 across the entire range of real numbers. With its simple form and applicability for all real numbers, our approximation surpasses either in computational efficiency or in relative error, and most often in both, other approximation formulas based on numerical algorithms or ad-hoc approximations. An extensive overview and classification of the existing approximations from the literature is included.