%     Reactor descontínuo isotérmico com reacção em série reversível A <--> B <-->C

%    Limpar todo o workspace
clc; clear; clf;

%     as concentrações iniciais são:
Ao=1;Bo=0;Co=0;
%     as constantes de velocidade são:

global k
k=[0.01 0.01;0.01 0];
% A->B A<-B
% B->C B<-C

%     Construção e resolução dos balanços de massa feitos na subrotina BMser.m
%     Serão resolvidos pelo metodo ode45 do Matlab entre t=0 e 1000 com condições iniciais C(0)=Co
[t,C]=ode45(@BMser,[0 1000],[Ao Bo Co]);
%     C é uma matriz com 3 colunas (uma por espécie) e N linhas, N determinado pelo dt determinado pelo método ode45.

plot(t,C);
%     Melhoramentos ao gráfico

title('Reaccao A->B->C')
legend('A', 'B', 'C')
ylabel('C (M)')
grid on
xlabel('t (min.)')

%    subrotina BMser.m
function dydt = BMser(t,y)
global k
dAdt = -k(1,1)*y(1)+k(1,2)*y(2);                                     %     dA/dt=-k₁A + k_-1B
dBdt=k(1,1)*y(1)-k(1,2)*y(2)-k(2,1)*y(2)+k(2,2)*y(3);    %     dB/dt=k₁A - k_-1B - k₂B + k_-2C
dCdt=k(2,1)*y(2)-k(2,2)*y(3);                                          %     dC/dt=k₂B + k_-2C
dydt=[dAdt; dBdt; dCdt];
end