function dy = mysys(t,y,params)


% 2 Chua csatolas
if nargin < 3; params= [9.5 15 8/7 4/63 0.1]; end
if nargin < 2; y = 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
dy=[-alfa*y(1)+alfa*y(2)-alfa*(-Rg1*y(1)+Rg3*y(1)^3);...
    y(1)-y(2)+y(3)-gcs*(y(2)-y(5));...
    -beta*y(2);...
    -alfa*y(4)+alfa*y(5)-alfa*(-Rg1*y(4)+Rg3*y(4)^3);...
    y(4)-y(5)+y(6)+gcs*(y(2)-y(5));...
    -beta*y(5)...
    ];

     
     
     
     
     
     
     
     
     
     
     
     
 
