

function F=chua(t,vec,flag,c3)

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

 x=vec(1);
 y=vec(2);
 z=vec(3);

 g = m1*x+(m0-m1)*(abs(x+1)-abs(x-1))/2;

 F(1)= c1*(y-x-g);
 F(2)= c2*(x-y+z);
 F(3)= -c3*y ;

 F=transpose(F);

