clear
r=load('kepler.dat');
mkr=6378.;

nt=length(r(:,1));
t=1:10000:nt;
%f1=figure(1);
%set(gcf,'Position',[100 100 1000 1000]);
hold on;
axis equal;
im= imread('earth1.jpg');
[x,y,z]=sphere(50)
surf(x,y,-z,im,'FaceColor','texturemap','Edgecolor','none');
plot3(r(:,2)./mkr,r(:,3)./mkr,r(:,4)./mkr);
quiver (r(t,2)./mkr,r(t,3)./mkr,r(t,6),r(t,7),'Color','black');
quiver (r(t,2)./mkr,r(t,3)./mkr,r(t,4),r(t,5),'Color','red');
view(-45,2);
legend('Maa','traj','a','v');
grid on;
