clear
a=load('lorentz_prooton.dat');
r_maa=6371;
hold on
axis equal
im= imread('earth1.jpg');
[x,y,z]=sphere(50);
surf(x*r_maa,y*r_maa,-z*r_maa,im,'FaceColor','texturemap','Edgecolor','none');
plot3(a(:,1),a(:,2),a(:,3));
grid on
