function up=rceh(t,u)
% sistema não linear, rigidez cúbica, excitação harmônica

% sistema
m=1;
c=0.168; % 0.2 (frequência); 0.168 (amplitude)
k=-0.5; % -1 (frequência); -0.5 (amplitude)
k3=0.5; % 1 (frequência); 0.5 (amplitude)

% força harmônica
  % alteração na frequência
% Fo=0.3; wf=1.4; % periódica
% Fo=0.3; wf=1.29; % caótica
  % alteração na amplitude (exercício)
% wf=1; Fo=0.210; % 0.100 0.177 0.178 0.197 0.210

up=[u(2); -(c/m)*u(2)-(k/m)*u(1)-(k3/m)*u(1)^3+(Fo/m)*cos(wf*t)];