function up=lse(t,u)
% sistema linear, 1 gdl, sem excitação (vibração livre), c.i. não nulas

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

up=[u(2); -(c/m)*u(2)-(k/m)*u(1)];