Julia has a great ODE package, providing a vast number of algorithms. The DiffEqPhys.jl package allows to use them for the numerical solution of Hamiltonian systems.
DiffEqPhys.jl