Abstract

This paper is concerned with the numerical approximation of stochastic ordinary differential equations, which satisfy a global monotonicity condition. This condition includes several equations with super-linearly growing drift and diffusion coefficient functions such as the stochastic Ginzburg-Landau equation and the 3/2-volatility model from mathematical finance. Our analysis of the mean-square error of convergence is based on a suitable generalization of the notions of C-stability and B-consistency known from deterministic numerical analysis for stiff ordinary differential equations. An important feature of our stability concept is that it does not rely on the availability of higher moment bounds of the numerical one-step scheme. While the convergence theorem is derived in a somewhat more abstract framework, this paper also contains two more concrete examples of stochastically C-stable numerical one-step schemes: the split-step backward Euler method from Higham et al. (2002) and a newly proposed explicit variant of the Euler-Maruyama scheme, the so called projected Euler-Maruyama method. For both methods the optimal rate of strong convergence is proven theoretically and verified in a series of numerical experiments.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call