function up=leh(t,u)
% sistema linear, excitação harmônica (c.i. nulas)

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

% força harmônica
F=150; 
wf=3;

up=[u(2); -(c/m)*u(2)-(k/m)*u(1)+(F/m)*sin(wf*t)];