EN
Initial value problems for systems of ordinary differential equations (ODEs) are solved numerically by using a combination of (a) the θ-method, (b) the sequential splitting procedure and (c) Richardson Extrapolation. Stability results for the combined numerical method are proved. It is shown, by using numerical experiments, that if the combined numerical method is stable, then it behaves as a second-order method.