% CSTR模型－－微分方程組
% 系統輸入是: Q，Qc; 系統輸出是: Ca，T, Tc, H
%   完成日期：2007年4月22日

function [K, L, M, N] = CSTR_MD(Q, Qc, Ca, T, Tc, H)

   % 模型參數
   
   Qf = 100.0;
   Tf = 320.0;
   Tcf = 300.0;
   Caf = 1.0;
   A = 16.66; % 單位是平方分米
   V = 100.0;
   k0 = 7.2*10^10;
   dH = -5.0*10^4;
   pCp = 239.0;
   pcCpc = 4175.0;
   EdR = 8750.0;
   UAc = 5.0*10^4;
   Vc = 10.0;

   % CSTR模型
   
   K = -k0*Ca*exp(-EdR/T)+(Qf*Caf-Q*Ca)/(A*H);
   L = k0*(-dH)*Ca*exp(-EdR/T)/pCp+(Qf*Tf-Q*T)/(A*H)+UAc*(Tc-T)/(pCp*A*H);
   M = Qc*(Tcf-Tc)/Vc+UAc*(T-Tc)/(pcCpc*Vc);
   N = (Qf-Q)/A;
end