%Program 9_1 
A=[2 3;0 -1];
B=[0 1; 1 0];
C=[1 1; 0 -1];
D=[1 0; 1 0];
x0=[2 -1];
dt=0.01;
t=0:dt:2;
f(:,1)=ones(length(t),1);
f(:,2)=exp(-3*t)';
sys=ss(A,B,C,D);
y=lsim(sys,f,t,x0);
subplot(2,1,1);
plot(t,y(:,1),'r');
ylabel('y1(t)');
xlabel('t');
subplot(2,1,2);
plot(t,y(:,2));
ylabel('y2(t)');
xlabel('t');