plot3(0,0,0);
xlabel('x');
ylabel('y');
zlabel('z');
axis([-3 3 -3 3 -3 3]);
hold on
Delta_t = 360;
for t=0:(Delta_t):12*3600    
    p=SatellitePos(getGPSSatellite( 1 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3);
    %plot(p1,p2);
    hold on;
    pause(0.01);
end
for t=0:(Delta_t):12*3600    
    p=SatellitePos(getGPSSatellite( 5 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3,'k');  
    hold on;
    pause(0.01);
end
for t=0:(Delta_t):12*3600     
    p=SatellitePos(getGPSSatellite( 9 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3,'r');
    hold on;
    pause(0.01);
end
%17 %19