function xdot=reac2(t,x,a0,b0,k,r,d)
xdot(1)=-k*x(1)^2+r*x(1)*x(2);
xdot(2)=k*x(1)^2-r*x(1)*x(2)-d*x(2);

