%This script produces a plot of x_1 and its estimate in Example 3.12
% The full-order estimator is in the loop

% written by S. H. Zak on Nov. 30, 2000


tspan=[0 3];
x0=[1 .2 -.1 0 0 0]';
[t,x]=ode23('comp_fo',tspan,x0);
plot(t,x(:,1),t,x(:,4),'--')
set(gca,'Fontsize',12)
xlabel('Time (sec)','Fontsize',12)
title('A plot of x_1 and its estimate (full-order estimator) in Example 3.12.','Fontsize',12)
legend('{\it x}_1','{\it x}_1-estimate')
grid
