Numerical integration of stochastic differential equations

A procedure for numerical integration of a stochastic differential equation, by extension of the Runge-Kutta method, is presented. The technique produces results which are statistically correct to a given order in the time step. Second- and third-order approximations are explicitly displayed.