function xdot=reac3(t,x,a0,b0,k1,k2)
xdot(1)=-k1*x(1);
xdot(2)=k1*x(1)-k2*x(2);


