Runge-Kutta-Nyström triples (Q1096341)

It is shown how the embedded Runge-Kutta-Nyström (RKN) process for the numerical solution of a special second-order initial value problem can be extended by the addition of a so-called dense output formulae for the dependent variable and its derivative. It results in RKN triples yielding non-mesh point solutions and thereby enhancing the practicability of the integration process. Continuity and order of interpolants are discussed. The equivalence of two possible modes of implementation of dense RKN formulae analogous to the RK methods is analyzed. Some new formulae of orders 4 and 6 and numerical tests are presented.