function dydt = chuaD(t,y,a,b,m0,m1)
% states [Uc1,Uc2,Il]
% a=9;b=14.286;m0=-1/7;m1=2/7;
dydt = [a*(y(2)-(m1*y(1)+.5*(m1-m0)*(abs(y(1)-1)-abs(y(1)+1))));y(1)-y(2)-y(3);b*y(2)];