%matlab time series of flight through model(neuts)
     clf
     load neutfoxTA36  -ascii;    
        tim=neutfoxTA36  (:,1);
        den1=neutfoxTA36 (:,2);
        den2=neutfoxTA36 (:,3);
        den3=neutfoxTA36 (:,4);
semilogy(tim,den1,'b-x')
     plot(tim,den1,'b-','linewidth',2)
     axis([-10,10,-25,5])
     xlabel ('Time from CA along TA (minutes)','fontsize',16)
     ylabel ('Denstiy cm^-^3','fontsize',16)
     title( 'Yung Model of Neutrals in Titans Ionosphere','fontsize',16)
hold on
     plot(tim,den2,'r-','linewidth',2)
hold on     
     plot(tim,den3,'k-','linewidth',2)
     
            legend('METAM','C_4N_2','C_2H_3CN')
