function dx = rigid(t,x)
%Definição das Variáveis
%m=0.1; l=0.2; k=0.015/(0.2^2); g=9.81;
% % g=9.81; m=g^-2; l=g; k=m; 
% % 
% dy = zeros(2,1);    % a column vector
% dy(1) = y(2);
% dy(2) = -(g/l)*sin(y(1))-(k/m)*y(2);
dx = zeros(2,1);    % a column vector
dx(1) = x(2) - x(1)^7*(x(1)^4 + 2*x(2)^2 -10);
dx(2) = -x(1)^3 -3*x(2)^5*(x(1)^4 + 2*x(2)^2 -10);