Strong convergence of the symmetrized Milstein scheme for some CEV-like SDEs

Abstract

In this paper, we study the rate of convergence of a symmetrized version of the Milstein scheme applied to the solution of the one dimensional SDE \[X_{t}=x_{0}+\int_{0}^{t}{b(X_{s})\,ds}+\int_{0}^{t}{\sigma\vert X_{s}\vert^{\alpha}\,dW_{s}},\qquad x_{0}>0,\sigma>0,\alpha\in[\frac{1}{2},1).\] Assuming $b(0)/\sigma^{2}$ big enough, and $b$ smooth, we prove a strong rate of convergence of order one, recovering the classical result of Milstein for SDEs with smooth diffusion coefficient. In contrast with other recent results, our proof does not relies on Lamperti transformation, and it can be applied to a wide class of drift functions. On the downside, our hypothesis on the critical parameter value $b(0)/\sigma^{2}$ is more restrictive than others available in the literature. Some numerical experiments and comparison with various other schemes complement our theoretical analysis that also applies for the simple projected Milstein scheme with same convergence rate.