A new family of multistep numerical integration methods based on Hermite interpolation

In this paper, a new family of explicit and implicit multistep methods is presented both for the error-controlled and uncontrolled modes. The main concept is to replace the Newton interpolation with the Hermite interpolation, where the Hermite polynomial is fitted to the function values and its derivatives. This idea is very useful in the numerical solution of problems (e.g., orbit propagation problem) where higher-order derivatives can easily be computed. In addition to the theoretical concept, the stability regions of the proposed methods are determined. The new methods are more stable than the well-known multistep numerical integrators (i.e., Adams–Bashforth and Adams–Bashforth–Moulton) in the explicit, implicit, and predictor–corrector forms. Using the second-order derivatives gives smaller error constants in the proposed method. The new integrators are numerically tested for a few examples, and the solutions are compared with those of the well-known multistep methods. Moreover, the CPU time and absolute integration error are compared in the satellite orbit propagation problem using various integration methods. The CHAMP mission, i.e., a German small-satellite mission for geoscientific and atmospheric research and applications, is considered as a case study for comparing the achievable accuracy of the proposed method with the existing method for solving the two-body problem.