function up=lec(t,u)
% sistema linear, excitação constante, c.i. nulas

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

% força
mc=10
g=9.81;
F=mc*g; 

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