function dy=dfcstr(x)
alfa=0.5;
beta=-0.1;
gama=20;
Da=0.1;
aux=Da*exp(gama*x./(1+x))./(1+Da*exp(gama*x./(1+x)));
dy=1-alfa*gama*aux.*(1-aux)./(1+x).^2;
