function up=mqs(t,u)
% máquina de serraria, modelo com 1 gdl, excitação harmônica, c.i. nulas

% sistema
m=3410; % (kg)
c=1.97e4; % (kg/s)
k=1.14e7; % (N/m)

% força harmônica
F=1000; 
wf=30;

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