
function F=chua_synchro(t,vec,flag,c3,c)

 c1=15.6;
 c2=1;
 m0=-8/7;
 m1=-5/7;

 x1 = vec(1);
 y1 = vec(2);
 z1 = vec(3);
 x2 = vec(4);
 y2 = vec(5);
 z2 = vec(6);

 g1 = m1*x1+(m0-m1)*(abs(x1+1)-abs(x1-1))/2;
 g2 = m1*x2+(m0-m1)*(abs(x2+1)-abs(x2-1))/2;

 F(1)= c1*(y1-x1-g1) + c*(x2-x1);
 F(2)= c2*(x1-y1+z1) + c*(y2-y1);
 F(3)= -c3*y1        + c*(z2-z1);
 F(4)= c1*(y2-x2-g2) + c*(x1-x2);
 F(5)= c2*(x2-y2+z2) + c*(y1-y2);
 F(6)= -c3*y2        + c*(z1-z2);

 F=transpose(F);


