function f=oded5(~,x)
global zz
r=norm(x([1,2]));
f=[x(3);x(4);-x(1)/r^3;-x(2)/r^3];
zz=zz+1;
