Contractivity-preserving explicit Hermite–Obrechkoff ODE solver of order 13

Author: Nguyen-Ba Truong   Desjardins Steven   Sharp Philip   Vaillancourt Rémi  

Publisher: Springer Publishing Company

ISSN: 0923-2958

Source: Celestial Mechanics and Dynamical Astronomy, Vol.117, Iss.4, 2013-12, pp. : 423-434

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

A new optimal, explicit, Hermite–Obrechkoff method of order 13, denoted by HO(13), that is contractivity-preserving (CP) and has nonnegative coefficients is constructed for solving nonstiff first-order initial value problems. Based on the CP conditions, the new 9-derivative HO(13) has maximum order 13. The new method usually requires significantly fewer function evaluations and significantly less CPU time than the Taylor method of order 13 and the Runge–Kutta method DP(8,7)13M to achieve the same global error when solving standard $$N$$ N -body problems.