function up=ler(t,u)
% sistema linear, excitação pulso retangular, c.i. não nulas

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

% força
F=1500; 
t0=1.5; 
t1=5;

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