clear
clf
r=load('kepler.dat');
hold on;
axis equal;
mkr=6378.;
grid on;

z0=r(:,1)*0;
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,z0/mkr,'color','black');
comet3(r(:,2)./mkr,r(:,3)./mkr,z0/mkr,0.01);
