On two-dimensional Hamiltonian systems with sixth-order integrals of motion

We obtain 21 two-dimensional natural Hamiltonian systems with sextic invariants, which are polynomial of the sixth order in momenta. Following to Bertrand, Darboux, and Drach these results of the standard brute force experiments can be applied to construct a new mathematical theory.