function val=energyval(z)
% Compute Energy at z
x=z(1); y=z(2); v=z(3); w=z(4);
val=(v^2+w^2)/2-1/sqrt(x^2+y^2);