Differential algebraic description for third and fifth order

The programs were written in COSY INFINITY [11]. The trajectory equations, Eqs. (1) and (2), were tracked by using an eighth order Runge-Kutta integrator [12] to produce the transfer map in fixed and rotating coordinates, respectively. Then, all the third and fifth order aberration coefficients were calculated in terms of Eqs. ................
................