Numerical simulation of a linear stochastic oscillator with additive noise Original Research Article

The ability of numerical methods to reproduce long-time features of a linear stochastic oscillator is examined. It is shown that certain, widely-used, methods fail to capture the correct second moment growth rate, whereas a customized extension of the partitioned Euler method behaves well in this respect. It is also shown that the partitioned Euler method inherits an infinite-oscillation property. A weaker oscillation result is proved for a wide class of numerical methods.