function z=moto(t,x)
% Es 4.12 di Ascher-Petzold 
% Serve ad apprezzare controllo dinamico passo
% prova [t,x]==RK(@moto, [0, 17.1], x0, 1.e-3, GRK4, bRK4);
% e [t1,x1]=ode45(@moto, [0, 17.1], x0)
% dove x0=[0.994; 0; 0; -2.00158510]
% fai plot(x(:,1),x(:,2)) e plot(x1(:,1), x1(:,2))
% Osserva diverso numero punti di integrazione!
z=zeros(4,1);
mu=0.012277471;
mu1=1-mu;
D1=( (x(1)+mu)^2 + x(2)^2 )^(3/2);
D2= ( (x(1)-mu1)^2 + x(2)^2)^(3/2);
z(1)=x(3);
z(2)=x(4);
z(3)=x(1)+2*x(4)-mu1*(x(1)+mu)/D1-mu*(x(1)-mu1)/D2;
z(4)=x(2)-2*x(3)-mu1*x(2)/D1-mu*x(2)/D2;
