%Matlab Program2-3
%Input the system parameter
F=[0 1;-0.632121 1];
G=[0 1]';
C=[0.264241 0.367879];
x0=[0;0];
N=80;
Tsam=0.01;
u=ones(N,1);
x=zeros(2,N);
y=zeros(1,N);
t=zeros(1,N);
x(:,1)=x0;
y(1)=C*x0;
t(1)=0;

%Main loop
for i=2:N
    x(:,i)=F*x(:,i-1)+G*u(i);
    y(i)=C*x(:,i);
    t(i)=Tsam*(i-1);
end
plot(t,y,'.',t,y)
grid
xlabel('T s')
ylabel('y')