function c = fbatelada (t,x)
global mixmax ks alfa yxs yps kp w
% equações (modelos) cinéticas

%mix=mixmax*(x(2)/(ks+x(2)))*(kp/(kp+x(3)));
mix=mixmax*(x(2)/(ks+x(2)+(x(2)^2/w)));
mip=alfa*mix;

% equaçao para x (dx/dt)
if(x(2)<=0);
c(1,1)=0;
else
c(1,1)=mix*x(1);
end

% equaçao para s (ds/dt)
if(x(2)<=0);
c(2,1)=0;

else
c(2,1)=(-(1/yxs)*mix*x(1))-((1/yps)*mip*x(1));
end;
% equaçao para p (dp/dt)
c(3,1)=mip*x(1);