global gamk;
gamma=6.67*10^-11;
mkm=5.97219*10^24;
mkr=6400.;
gamk=gamma*mkm*10^-9;
r0=[0.;6400.;11;0.];
options = odeset('RelTol',1e-5,'AbsTol',1e-5,'MaxStep',1000.,'InitialStep',1.);
[t,r]=ode45(@kplf,[0,400000], r0,options);
%comet (r(:,1),r(:,2),0.01);
hold on;
axis equal;
%axis([-2 2. -2. 2.]);
grid on;
z=r(:,1)*0;
sphere
plot3(r(:,1)./mkr,r(:,2)./mkr,z);
comet3(r(:,1)./mkr,r(:,2)./mkr,z,0.01);

