A simple, very accurate algorithm for numerical simulation of stochastic differential equations is described. Its relationship to colored noise is elucidated and exhibited by explicit results. The especially delicate problem of mean first passage times is highlighted and highly accurate agreement between the numerical simulations and analytic results are shown.