|
Accurate methods used to evaluate the inverse of the standard normal cumulative distribution function at probability p commonly used today are too cumbersome and/or slow to obtain a large number of evaluations reasonably quickly, e.g., as required in certain Monte Carlo applications. Previously reported simple approximations all have a maximum absolute error em>10 to the -4 for a p-range of practical concern such as Min[p,1-p]>=10 to the -6. An 11-term polynomial-based approximation is presented for which em<10 to the -6 in this range.
|
|