function xdot=reac6(t,x,a,b) 
xdot(1)=a-x(1)+x(1)^2*x(2);
xdot(2)=b-x(1)^2*x(2);



