function x=ltisol(A,B,u,x0,t)
    n=length(x0);//orden del sistema
    r=length(t);//numero de puntos en el tiempo
    x=zeros(n,r);
    x(:,1)=x0;
    iA=inv(A)
    for k=2:r
        eA=expm(A*t(k));
        x(:,k)=eA*x0+iA*(eA-eye(n,n))*B*u;
    end    
endfunction