function up=lep(t,u)
% sistema linear, 1 gdl, excitação periódica, c.i. nulas

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

% força periódica
F1=150;
w1=3;
F2=100;
w2=6;

up=[u(2); -(c/m)*u(2)-(k/m)*u(1)+...
           (F1/m)*sin(w1*t)+(F2/m)*sin(w2*t)];