Fixed step, symmetric Runge-Kutta-Nyström formulae have proved to be very efficient for the numerical integration of a large class of reversible second order systems of ordinary differential equations of the special form d 2y dt2 = f(t, y). However for some important classes of problems it is necessary, for the sake of efficiency, to allow a variable steplength of integration to be used and in ...