%Ec mov
function yp=ecdifmov(t,y)
global GEO mu wt
t;
[GU U]=Gpot(y(1),y(2),y(3),GEO,mu);
yp=[y(4),y(5),y(6),2*wt*y(5)+wt^2*y(1)+GU(1),-2*wt*y(4)+wt^2*y(2)+GU(2),GU(3)]';
end