High order methods for the integration of the Bateman equations and other problems of the form of y′ = F(y,t)y

The statistical properties of the algorithms were also analyzed during the Monte Carlo problem. The new methods are shown to yield smaller standard deviations on final quantities as compared to the reference predictor-corrector method, by up to a factor of 1.4.