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

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

% força
f0=150; wf=6; 

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