function f=reator(x)

gama=20;
teta=0.1;
Caf=1;
Da=0.134;
Tf=1;
Tc=0.8;
Bo=0.36;
beta=0.8;

ep=exp(-gama*(1/x(2)-1))*x(1)*Da;

f(1,1)=teta*(Caf-x(1))-ep;
f(2,1)=teta*Tf+beta*Tc-(teta+beta)*x(2)+Bo*ep;
