function y = f2(x)

alfa = 0.5;
beta = -0.1;
gama = 20;
Da   = 0.1;
CAe  = 1.0;
n    = 1.0;

y(1,1) = 1 - x(1) - Da*CAe^(n-1)*x(1)^n.*exp(gama*x(2)./(1+x(2)));
y(2,1) = x(2) - beta - alfa*(1 - x(1));
