% Problem 4.1 - Transient response to initial conditions

A1=[0 1 0; 0 0 1; 0 0 0];
A2=[0 1 0; 0 0 0; 0 0 0];
A3=[0 0 0; 0 0 0; 0 0 0];
x0=[1; 1; 1];

% Add B, C, and D matrices
B=[1; 1; 1];
C=[1 1 1];
D=0;

% 4.1(a)
[y1,x1,t]=initial(A1,B,C,D,x0);
plot(t,x1)
grid
xlabel('time [secs]')
ylabel('x1,x2,x3')

% 4.1(b)
[y2,x2,t]=initial(A2,B,C,D,x0)
figure(2)
plot(t,x2)
grid
xlabel('time [secs]')
ylabel('x1,x2,x3')

% 4.1(c)
[y3,x3,t]=initial(A3,B,C,D,x0)
figure(3)
plot(t,x3)
grid
xlabel('time [secs]')
ylabel('x1,x2,x3')