
% matlab 3d plotting of particle positions
clf
load oTAtraj -ascii;
x1=oTAtraj(:,1);
y1=oTAtraj(:,2);
z1=oTAtraj(:,3);
az=0.
zen=90.
plot3(x1,y1,z1,'or','Markersize',2);
xlabel('Titan Radii (x points to sun)','Fontsize',13);
ylabel('Titan Radii','Fontsize',13);
zlabel('Titan Radii','Fontsize',13);
title('TA trajectory','Fontsize',15)
axis([-12.,12.,-12.,12.,-12.,12.])
axis('square')
view(az,zen);
box on
grid off
hold on
sphere(25)
colormap(gray)
rotate3d on
