Several numerical methods for treating stochastic differential equations are considered. Both the convergence in the mean square limit and the convergence of the moments is discussed and the generation of appropriate random numbers is treated. The necessity of simulations at various time steps with an extrapolation to time step zero is emphasized and demonstrated by a simple example.