Monte-Carlo simulation techniques are widely used to approximate quantities of interest of the solution of random ordinary differential equations.
Over the last decades, methods based on spectral expansions of the solution process have drawn great interest.