Error growth in the numerical integration of periodic orbits by multistep methods, with application to reversible systems

We study the growth with time of (the coefficients of the asymptoticexpansion of) the error in the numerical integration with linearmultistep methods of periodic solutions of systems of ordinarydifferential equations. Particular attention is devoted to reversiblesystems. It turns out that symmetric linear multistep methodscannot be recommended in spite of the fact that they mimic thereversibility of the true flow. For reversible second-ordersystems, linear multistep methods without parasitic double rootsare useful.