function t = plots_TC11_lon_ekf(t, parOEM, Z, U,syekf, sxekf, sxekfstd, method);


XL = t(end) + 5;

% Plot time histories of measured and estimated observation variables and control inputs
figure(1)
subplot(311),plot(t,Z(:,1), t,syekf(:,1),'r-.');grid;ylabel('{\alpha} (°)');
title('EKF: Time histories of output variables (measured and estimated); input variables')
subplot(312),plot(t,Z(:,2), t,syekf(:,2),'r-.');grid;ylabel('q (°/s)');
subplot(313),plot(t,U(:,1));grid;ylabel('{\delta_e} (°)'); xlabel('Time in sec'); 

% Plot time histories of estimated parameters
figure(2)
subplot(811),plot(t,sxekf(:,3),'k', t,sxekf(:,3)+sxekfstd(:,3),'r',...
                                    t,sxekf(:,3)-sxekfstd(:,3),'r',...
                                    [0 XL],[parOEM(1,1) parOEM(1,1)]);grid;ylabel('Z_0');
title('EKF: Convergence of parameter estimates with error bounds')
subplot(812),plot(t,sxekf(:,4),'k', t,sxekf(:,4)+sxekfstd(:,4),'r',...
                                    t,sxekf(:,4)-sxekfstd(:,4),'r',...
                                    [0 XL],[parOEM(1,2) parOEM(1,2)]);grid;ylabel('Z_{\alpha}');
subplot(813),plot(t,sxekf(:,5),'k', t,sxekf(:,5)+sxekfstd(:,5),'r',...
                                    t,sxekf(:,5)-sxekfstd(:,5),'r',...
                                    [0 XL],[parOEM(1,3) parOEM(1,3)]);grid;ylabel('Z_q');
subplot(814),plot(t,sxekf(:,6),'k', t,sxekf(:,6)+sxekfstd(:,6),'r',...
                                    t,sxekf(:,6)-sxekfstd(:,6),'r',...
                                    [0 XL],[parOEM(1,4) parOEM(1,4)]);grid;ylabel('Z_{\deltae}');
subplot(815),plot(t,sxekf(:,7),'k', t,sxekf(:,7)+sxekfstd(:,7),'r',...
                                    t,sxekf(:,7)-sxekfstd(:,7),'r',...
                                    [0 XL],[parOEM(1,5) parOEM(1,5)]);grid;ylabel('M_0');
subplot(816),plot(t,sxekf(:,8),'k', t,sxekf(:,8)+sxekfstd(:,8),'r',...
                                    t,sxekf(:,8)-sxekfstd(:,8),'r',...
                                    [0 XL],[parOEM(1,6) parOEM(1,6)]);grid;ylabel('M_{\alpha}');
subplot(817),plot(t,sxekf(:,9),'k', t,sxekf(:,9)+sxekfstd(:,9),'r',...
                                    t,sxekf(:,9)-sxekfstd(:,9),'r',...
                                    [0 XL],[parOEM(1,7) parOEM(1,7)]);grid;ylabel('M_q');
subplot(818),plot(t,sxekf(:,10),'k',t,sxekf(:,10)+sxekfstd(:,10),'r',...
                                    t,sxekf(:,10)-sxekfstd(:,10),'r',...
                                    [0 XL],[parOEM(1,8) parOEM(1,8)]);grid;ylabel('M_{\deltae}');
                                    xlabel('Time in sec');  

% Plot time histories of estimated parameters
figure(3)
subplot(321),plot(t,sxekf(:,4),'m',[0 XL],[parOEM(1,2) parOEM(1,2)]);
                  axis([0 XL,-2 2]);grid;legend('Z_{\alpha} EKF','    MLE',1);
title('EKF: Convergence of parameter estimates')
subplot(323),plot(t,sxekf(:,5),'m',[0 XL],[parOEM(1,3) parOEM(1,3)]);
                  axis([0 XL,-2 2]);grid;legend('Z_q EKF','    MLE',1);
subplot(325),plot(t,sxekf(:,6),'m',[0 XL],[parOEM(1,4) parOEM(1,4)]);
                  axis([0 XL,-2 2]);grid;legend('Z_{\deltae} EKF','     MLE',1);
                  xlabel('Time in sec')
subplot(322),plot(t,sxekf(:,8),'m',[0 XL],[parOEM(1,6) parOEM(1,6)]);
                  axis([0 XL, -10 0]);grid;legend('M_{\alpha} EKF','     MLE',1);
subplot(324),plot(t,sxekf(:,9),'m',[0 XL],[parOEM(1,7) parOEM(1,7)]);
                  axis([0 XL,-5 0]);grid;legend('M_q EKF','     MLE',1);
subplot(326),plot(t,sxekf(:,10),'m',[0 XL],[parOEM(1,8) parOEM(1,8)]);
                  axis([0 XL,-10 0]);grid;legend('M_{\deltae} EKF','     MLE',1);
                  xlabel('Time in sec')