ODEs solved by RK4 methodes

Question

I\'m trying to simulate free particle in 2D which is describes by Hamilton function

H = p_x^2/2m + p_y^2/2m.

At points x = -60, x = 60, y =