Example 9-6 Falling Off the Upper Steady State   |POLVER05_0
d(Ca)/d(t) =  1/tau*(Ca0-Ca)+ra  #
d(Cb)/d(t) =  1/tau*(Cb0-Cb)+rb  #
d(Cc)/d(t) =  1/tau*(0-Cc)+rc  #
d(Cm)/d(t) =  1/tau*(Cm0-Cm)  #
d(T)/d(t) =  (Q-Fa0*ThetaCp*(T-T0)+(-36000)*ra*V)/NCp  #
Fa0 =  80  #
T0 =  70  #
V =  (1/7.484)*500  #
UA =  16000  #
Ta1 =  60  #
k =  16.96e12*exp(-32400/1.987/(T+460))  #
Fb0 =  1000  #
Fm0 =  100  #
mc =  1000  #
ra =  -k*Ca  #
rb =  -k*Ca  #
rc =  k*Ca  #
Nm =  Cm*V  #
Na =  Ca*V  #
Nb =  Cb*V  #
Nc =  Cc*V  #
ThetaCp =  35+Fb0/Fa0*18+Fm0/Fa0*19.5  #
v0 =  Fa0/0.923+Fb0/3.45+Fm0/1.54  #
Ta2 =  T-(T-Ta1)*exp(-UA/(18*mc))  #
Ca0 =  Fa0/v0  #
Cb0 =  Fb0/v0  #
Cm0 =  Fm0/v0  #
Q =  mc*18*(Ta1-Ta2)  #
tau =  V/v0  #
NCp =  Na*35+Nb*18+Nc*46+Nm*19.5  #
X = 1-Ca/(80/v0) #
t(0)=0
Ca(0)=0.03789
Cb(0)=2.12
Cc(0)=0.143
Cm(0)=0.2265
T(0)=138.53
t(f)=4