Numerical solution of Hamiltonian systems by multi-value methods