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.

This content is only available via PDF.
You do not currently have access to this content.