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

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

% força
F=150; 
t0=0; 
t1=0.001;

up=[u(2); -(c/m)*u(2)-(k/m)*u(1)+...
           ((F/m)*(1/t1)*(stepfun(t,t0)-stepfun(t,t1)))];