function z=E(t,y)
M=1;L=4;
z=[y(2),-M/y(1)^2+L^2/y(1)^3-3*M*L^2/y(1)^4,L/y(1)^2]';
