> restart;

> g:=x->c*x+((d-c)/2)*(abs(x+1)-abs(x-1)):

> a:=15: b:=25.58: c:=-5/7: d:=-8/7:

> chua:=diff(x(t),t)=a*(y(t)-x(t)-g(x(t))),
diff(y(t),t)=x(t)-y(t)+z(t),
diff(z(t),t)=-b*y(t):

> warpocz:=x(0)=.2,y(0)=.3,z(0)=.1:

> with(DEtools):

> with(plots):

Warning, the name changecoords has been redefined

> F:=dsolve({chua,warpocz},numeric,output=listprocedure);

> odeplot(F,[x(t),y(t)],0..150);

> odeplot(F,[x(t),z(t)],0..150);

> odeplot(F,[z(t),y(t)],0..150);

> odeplot(F,[x(t),y(t),z(t)],0..150);

>