Estimation of global errors in the numerical solution of ordinary differential equations Free

The aim of this work is to present a new method to estimate errors propagated in the numerical integration of ordinary differential equations. The method has proved to be very efficient in many numerical applications, and its main advantages are shown in a set of examples concerning some cases of the gravitational N‐body problem.