Multi-step methods for random ODEs driven by Itô diffusions

Linear multi-step methods are derived for random ordinary differential equations (RODEs) driven by the solutions of Itô stochastic differential equations (SODEs) via strong Itô–Taylor schemes for SODEs. Due to the special structure of the RODE–SODE pair it is not necessary to restrict the intensity of the noise. Pathwise convergence is established as well as the B-stability of implicit multi-step methods. Numerical comparisons are provided for explicit schemes applied to a low dimensional RODE and implicit schemes applied to a high dimensional RODE obtained with the method of lines by spatially discretizing a random partial differential equation with finite difference quotients.