Conservative methods for the Toda lattice equations

We are concerned with the numerical integration of the Todalattice equations by using different conservative methods. Numericalexperiments suggest that the global error for isospectral schemesdecreases exponentially with time but it is almost constantfor either symplectic or more general integrators. We providea theoretical explanation for these experimental findings.