Among all statistical distributions, standard normal distribution has been the most important and practical distribution in which calculation of area under probability density function and cumulative distribution function are required. Unfortunately, the cumulative distribution function of this is, in general, expressed as a definite integral with no closed form or analytical solution. Consequently, it has to be approximated. In this paper, attempts have been made for Winitzki's approximation to be proved by a new approach. Then, the approximation is improved with some modifications and shown that the maximum error resulted from this is less than 0.0000584. Finally, an inverse function for computation of normal distribution quantiles has been derived.