    







      dde = (de+5)/57.3;

       won = ones(size(dde));

      data = lsim(a,b,c,d,[dde won],time);


     Nz = data(:,3);

    hold on

   figure(1)


   hold on
   subplot(5,1,1),plot(time,alpha,time,data(:,1))
   ylabel('alpha')
   grid

   hold on
   subplot(5,1,3),plot(time,az,time,Nz)
   ylabel('az')
   grid
 
   hold on
   subplot(5,1,5),plot(time,qq,time,data(:,2))
   ylabel('q')
   grid 


 %  data(:,1)+4 is alpha
 %  Nz is normal acceleration
 %  data(:,2) is pitch rate 



