function up=lei(t,u)
% sistema linear, 1 gdl, excitação impulsiva, c.i. nulas

% sistema
m=100; 
c=50; 
k=2000;

% força impulsiva
I=150; % impulso
t0=0; % início
t1=0.001; % final

% % versão original
% up=[u(2); -(c/m)*u(2)-(k/m)*u(1)+...
%            ((I/m)*(1/t1)*(stepfun(t,t0)-stepfun(t,t1)))];

% versão alternativa
up=[u(2); -(c/m)*u(2)-(k/m)*u(1)+...
           (I/m)*(1/(t1-t0))*((t>=t0)-(t>=t1))];