Second-Order Methods for Solving Stochastic Differential Equations

In this paper we discuss the numerical methods with second-order accuracy for solving stochastic differential equations. An unbiased sample approximation method for $I_n=\int ^{t_{n+1}}_{t_n}(B_u-B_{t_n})^2du$ is proposed, where {$B_u$} is a Brownian motion. Then second-order schemes are derived both for scalar cases and for system cases. The errors are measured in the mean square sense. Several numerical examples are included, and numerical results indicate that second-order schemes compare favorably with Euler's schemes and 1.5th-order schemes.