A new higher-order weak approximation scheme for stochastic differential equations and the Runge–Kutta method

Abstract

The authors report on the construction of a new algorithm for the weak approximation of stochastic differential equations. In this algorithm, an ODE-valued random variable whose average approximates the solution of the given stochastic differential equation is constructed by using the notion of free Lie algebras. It is proved that the classical Runge–Kutta method for ODEs is directly applicable to the ODE drawn from the random variable. In a numerical experiment, this is applied to the problem of pricing Asian options under the Heston stochastic volatility model. Compared with some other methods, this algorithm is significantly faster.