function [hist_x,y] = simula_DLTI_StateSpace_discreto_mor(A,B,C,D,u,x0)
%

hist_x = [];
N = length(u);
y = zeros(N,1);
x = x0;
for i=1:N
  y(i) = C * x  +  D * u(i);
  x = A * x  +  B * u(i);
  hist_x = [hist_x x];
end


