hamilton_eqs receives a Hamiltonian and returns a sequence with Hamilton's equations and a list with the p's and q's involved.

Some useful conventions were adopted to represent the p's and q's. All p's and q's must appear as pn or qn where n is a positive integer, as in p1, p2, and the time dependence need not be explicit, as in pn or qn instead of pn(t) or qn(t). The Hamilton equations will be automatically returned using pn(t) or qn(t).

This function is part of the DEtools package, and so it can be used in the form hamilton_eqs(..) only after executing the command with(DEtools). However, it can always be accessed through the long form of the command by using DEtools[hamilton_eqs](..).