function up=rcse(t,u)
% sistema não linear, rigidez cúbica, sem excitação

% sistema
m=1;
c=0.016;
k=-0.5;
k3=0.5;

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