function xcn = fbatalim (t,x)
global mixmax ks alfa yxs yps va sal kp
% equações (modelos) cinéticas
mix=mixmax*x(2)/(ks+x(2))*exp(-kp*x(3));
mip=alfa*mix;
% equações (modelos) de balanço no reator 
% equaçao para x (dx/dt)
xcn(1,1)=mix*x(1)-va/x(4)*x(1);
if(x(2)<=0);
xcn(2,1)=0;
else
% equaçao para s (ds/dt)
xcn(2,1)=va/x(4)*sal-(1/yxs*mix*x(1))-(1/yps*mip*x(1))-va/x(4)*x(2);
end;
% equaçao para p (dp/dt)
xcn(3,1)=mip*x(1)-va/x(4)*x(3);
% equaçao para v (dv/dt)
xcn(4,1)=va;