clear all

if false
  my_time = 0:0.01:10;
else
  my_time = linspace(0, 10, 1024);
end

i = 1;
for t = my_time
  x_1(i) = (exp(-5*t)*(4*exp(5*t)+45*exp(4*t)-9))/2.0E+1;
  x_2(i) = -(exp(-5*t)*(9*exp(4*t)-9))/4.0E+0;
  i = i + 1;
end

clf
plot(my_time, x_1)
hold
plot(my_time, x_2)
grid on
legend("x_1", "x_2")
xlabel("time [s]")
ylabel("state vector")
title("time response")

if false
  print("time_response.png", "-dpng");
end
