Multi-value Numerical Methods for Hamiltonian Systems