Stability of numerical solution of systems of ordinary differential equations