function y=c5mross(x)
a=0.2; b=0.2; c=5.7; 
y=[-x(2)-x(3); x(1)+a*x(2); b+[x(1)-c]*x(3)];