For Calogero and Toda dynamical equations two numerical methods of arbitrary high order, conserving the Hamiltonian are developed. The methods consist of modifications of conventional polynomial extrapolation with the Gragg method used as a basic method. The theoretical study is confirmed by a number of numeiical examples.