A weak trapezoidal method for a class of stochastic differential equations

We present a numerical method for the approximation of solutions for the class of stochastic differential equations driven by Brownian motions which induce stochastic variation in fixed directions. This class of equations arises naturally in the study of population processes and chemical reaction kinetics. We show that the method constructs paths that are second order accurate in the weak sense. The method is simpler than many second order methods in that it neither requires the construction of iterated Ito integrals nor the evaluation of any derivatives. The method consists of two steps. In the first an explicit Euler step is used to take a fractional step. This fractional point is then combined with the initial point to obtain a higher order, trapezoidal like, approximation. The higher order of accuracy stems from the fact that both the drift and the quadratic variation of the underlying SDE are approximated to second order.

[1]  W. Wasow,et al.  On the Approximation of Linear Elliptic Differential Equations by Difference Equations with Positive Coefficients , 1952 .

[2]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[3]  E. Helfand.  Numerical integration of stochastic differential equations , 1979, The Bell System Technical Journal.

[4]  T. Kurtz.  Representations of Markov Processes as Multiparameter Time Changes , 1980 .

[5]  D. Talay,et al.  Discretization and simulation of stochastic differential equations , 1985 .

[6]  J. Norris.  Simplified Malliavin calculus , 1986 .

[7]  Denis R. Bell.  The Malliavin Calculus , 1987 .

[8]  D. Talay,et al.  Expansion of the global error for numerical schemes solving stochastic differential equations , 1990 .

[9]  Eckhard Platen,et al.  On weak implicit and predictor-corrector methods , 1995 .

[10]  Jessica G. Gaines,et al.  Variable Step Size Control in the Numerical Solution of Stochastic Differential Equations , 1997, SIAM J. Appl. Math..

[11]  Desmond J. Higham,et al.  Mean-Square and Asymptotic Stability of the Stochastic Theta Method , 2000, SIAM J. Numer. Anal..

[12]  Jonathan C. Mattingly,et al.  Ergodicity for SDEs and approximations: locally Lipschitz vector fields and degenerate noise , 2002 .

[13]  K. Burrage,et al.  Predictor-Corrector Methods of Runge-Kutta Type for Stochastic Differential Equations , 2002, SIAM J. Numer. Anal..

[14]  G. N. Milstein,et al.  Numerical Integration of Stochastic Differential Equations with Nonglobally Lipschitz Coefficients , 2005, SIAM J. Numer. Anal..

[15]  S. Ethier,et al.  Markov Processes: Characterization and Convergence , 2005 .

[16]  Jonathan C. Mattingly,et al.  An adaptive Euler–Maruyama scheme for SDEs: convergence and stability , 2006, math/0601029.

[17]  Nicola Bruti-Liberati,et al.  Strong Predictor-Corrector Euler Methods for Stochastic Differential Equations , 2008 .

[18]  David F. Anderson.  Incorporating postleap checks in tau-leaping. , 2007, The Journal of chemical physics.

[19]  Weak convergence in the Prokhorov metric of methods for stochastic differential equations , 2007, 0707.4466.