  We revisit an algorithm by Skeel et al. for computing the modified, or
shadow, energy associated with the symplectic discretization of Hamiltonian
systems. By rephrasing the algorithm as a Richardson extrapolation scheme
arbitrary high order of accuracy is obtained, and provided error estimates show
that it does capture the theoretical exponentially small drift associated with
such discretizations. Several numerical examples illustrate the theory.
