function dx = mysys(t,x,params)
    
% Chua
if nargin < 3; params= [9.5 15 8/7 4/63]; end
if nargin < 2; x = zeros(3,1); end
if nargin < 1; t = 0; end

alfa=params(1);
beta=params(2);
Rg1=params(3);
Rg3=params(4);

dx = [ -alfa*x(1)+alfa*x(2)-alfa*(-Rg1*x(1)+Rg3*x(1)^3);...
       x(1)-x(2)+x(3);...
       -beta*x(2)...
     ];
          

 
 
 
 
 
 
 
 
 