Example RE9-3 Linear Stability  |POLVER05_0
d(ca)/d(t) = vo*(cao-ca)/V-k*ca #
d(T)/d(t) = (UA*(Ta-T)-rho*vo*Cp*(T-To)-dhrxn*k*ca*V)/(rho*V*Cp) #
d(x)/d(t) = (-L*x+N/J*y)/tau #
d(y)/d(t) = (J*(1-L)*x-(M*N)*y)/tau #
vo = 400 #
cao = .5 #
V = 48 #
UA = 150*250 #
Ta = 530 #
rho = 50 #
Cp = .75 #
To = 530 #
dhrxn = -30000 #
L = 7.1 #
N = 14.21 #
M = 3.5 #
E = 30000 #
y1 = y/628 #
x1 = x/.0684 #
k = 2*7.08*10^11*exp(-E/1.987/T) #
tau = V/vo #
J = -dhrxn/rho/Cp #
t(0)=0
ca(0)=.0684
T(0)=628.645
x(0)=-.02
y(0)=-2
t(f)=.2