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

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

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

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

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