function y = f1(x)

alfa = 0.5;
beta = -0.1;
gama = 20;
Da   = 0.1;

y = x - beta - alfa*Da*exp(gama*x./(1+x))./(1+Da*exp(gama*x./(1+x)));
