function dy = sync_chua_sys(t,y, params)


% 2 Chua csatolas

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) ];