function xcn = fcontinuo (t,x)
global mixmax ks ki alfa beta yxs yps D sal kp



if (t > 500)
  D = 0.02;
end
  
  

% equações (modelos) cinéticas
mix=mixmax*x(2)/(ks+x(2)+(x(2)^2/ki))*exp(-kp*x(3));
%mix=mixmax*x(2)/(ks+x(2));

mip=alfa*mix+beta;
% equações (modelos) de balanço no reator 
% equaçao para x (dx/dt)
xcn(1,1)=mix*x(1)-D*x(1);
% equaçao para s (ds/dt)
xcn(2,1)=D*sal-(1/yxs*mix*x(2))-(1/yps*mip*x(1)*x(2))-D*x(2);
% equaçao para p (dp/dt)
xcn(3,1)=mip*x(1)-D*x(3);