Construction of classical superintegrable systems with higher order integrals of motion from ladder operators

We construct integrals of motion for multidimensional classical systems from ladder operators of one-dimensional systems. This method can be used to obtain new systems with higher order integrals. We show how these integrals generate a polynomial Poisson algebra. We consider a one-dimensional system with third order ladders operators and found a family of superintegrable systems with higher order integrals of motion. We obtain also the polynomial algebra generated by these integrals. We calculate numerically the trajectories and show that all bounded trajectories are closed.