% Simulação de Equações de Balanço para Reator Contínuo sem reciclo de
% células

clear all
global mixmax ks ki alfa beta yxs yps D sal kp
%parâmetros cinéticos empregados nos modelos
mixmax=0.042; % 0.042
ks=0.03; % 0.03
ki=800; % 800
kp=0.01;% 0.01
alfa=0.05; % 0.05
beta=0.005; % 0.005

% vazão de alimentação
D=0.03; % 0.04
% Concentração da Solução de Alimentação
sal=4.5; % 33
% fatores de conversão empregado nos modelos
yxs=0.23; % 0.23
yps=0.40; % 0.4
% Condições iniciais das variáveis dependentes
x0=4.6;
s0=0.575;
p0=0.8065;
% intervalo da variável independente, tempo
tspan=[0,5000];
y0=[x0;s0;p0];
% Chamando a rotina Runge-Kutta 4/5 ordem para integrar y0 num tempo span
[t,xcn]=ode45('fcontinuo',tspan,y0);
% Como a minha saída é a matriz xcn, preciso identificar as colunas
x=xcn(:,1);
s=xcn(:,2);
p=xcn(:,3);
% montando os gráficos
subplot(221), plot(t,x)
%axis([0 500 0 6])
title('X(g/L)');
subplot(222), plot(t,s)
%axis([0 500 0 6])
title('S(g/L)');
subplot(223), plot(t,p)
%axis([0 500 0 6])
title('P(g/L)');