The Implementation of Milstein Scheme in Two-Dimensional SDEs Using the Fourier Method

Multiple stochastic integrals of higher multiplicity cannot always be expressed in terms of simpler stochastic integrals, especially when the Wiener process is multidimensional. In this paper we describe how the Fourier series expansion of Wiener process can be used to simulate a two-dimensional stochastic differential equation (SDE) using Matlab program. Our numerical experiments use Matlab to show how our truncation of Itô’-Taylor expansion at an appropriate point produces Milstein method for the SDE