function up=lei2(t,u)
% sistema linear, excitação tipo pulso, c.i. não nulas, 2 gdl

% sistema
M=[2 0; 0 1]; 
C=[0.3 -0.05; -0.05 0.05]; 
K=[3 -1; -1 1];
A=[zeros(2) eye(2); -inv(M)*K -inv(M)*C];

% excitação
F=[0; 1];
t1=1; 
t2=1.1;
f=inv(M)*F;

up=A*u+[0; 0; f]*(stepfun(t,t1)-stepfun(t,t2));