function dx = mysys(t,x,params)


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

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

% ToDo
dx=[-alfa*x(1)+alfa*x(2)-alfa*(-Rg1*x(1)+Rg3*x(1)^3);...
    x(1)-x(2)+x(3)-gcs*(x(2)-x(5));...
    -beta*x(2);...
    -alfa*x(4)+alfa*x(5)-alfa*(-Rg1*x(4)+Rg3*x(4)^3);...
    x(4)-x(5)+x(6)+gcs*(x(2)-x(5));...
    -beta*x(5)...
    ];

     
     
     
     
     
     
     
     
     
     
     
     
 
