function f=seh(t,x)
% sistema com 1 gdl sob excitação harmônica, c.i. não nulas

m=100; c=100; k=2000; % sistema
f0=150; wf=6; % força

f=[x(2); -(c/m)*x(2)-(k/m)*x(1)+(1/m)*f0*sin(wf*t)];