T=0:0.01:5;U=0.25*ones(size(T));[Y,X]=lsim(A-B*K,B,C,D,U,T);