function V = henon_heiles(x,y)

% henon(X,Y) returns the value of the Henon-Heiles potential at (X,Y).

V = 1/2*(x.^2+y.^2) + (x.^2).*y - 1/3*(y.^3);
