function J=J_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;

J(1,1)=-teta-ep;
J(1,2)=-ep*x(1)*gama/(x(2)*x(2));
J(2,1)=Bo*ep;
J(2,2)=-teta-beta+Bo*ep*x(1)*gama/(x(2)*x(2));
