% Simulação de Equações de Balanço para Reator Descontínuo Alimentado

clear all
global mixmax ks alfa yxs yps va sal kp
%parâmetros cinéticos empregados nos modelos
mixmax=0.25;
ks=0.03;
kp=0.05;
alfa=0.2;

% vazão de alimentação
va=0.08; % (l/h)
% Concentração da Solução de Alimentação
sal=150;
% fatores de conversão empregado nos modelos
yxs=0.52;
yps=0.51;

% Condições iniciais das variáveis dependentes
x0=2;
s0=200;
p0=0;
v0=2;
% intervalo da variável independente, tempo
tspan=[0,15];
y0=[x0;s0;p0;v0];
% Chamando a rotina Runge-Kutta 4/5 ordem para integrar y0 do modelo fbatalim num tempo span
[t,xcn]=ode45('fbatalim',tspan,y0);
% Como a minha saída é a matriz xcn, preciso identificar as colunas
x=xcn(:,1);
s=xcn(:,2);
p=xcn(:,3);
v=xcn(:,4);
% montando os gráficos
subplot(221), plot(t,x)
title('X(g/L)');
subplot(222), plot(t,s)
title('S(g/L)');
subplot(223), plot(t,p)
title('P(g/L)');
subplot(224), plot(t,v)
title('V(L)');