% 本程序用來實現CSTR模型的仿真
% 
%   完成日期：2007年4月22日

clear
N = 100;

Q = 100.0;
Qc = 15.0;

Ca(1) = 0.037;
T(1) = 402.35;
Tc(1) = 345.44;
H(1) = 6.0; % 單位是分米

for i = 1:N
    [Ca(i+1), T(i+1), Tc(i+1), H(i+1)] = CSTR_RK4(Q, Qc, Ca(i), T(i), Tc(i), H(i));
end

subplot(4,1,1)
plot(Ca)
grid on

subplot(4,1,2)
plot(T)
grid on

subplot(4,1,3)
plot(Tc)
grid on

subplot(4,1,4)
plot(H)
grid on