H. Scott Fogler Example 6.6 -polymath-

Fogler ex-6.6

Polymath

Don't forget to select Stiff method!!!!
polymath file....

https://hotfile.com/dl/202153300/34a8e79/66.pol.html

++++++++++++++++++++ODE Report (STIFF)+++++++++++++++++++

Differential equations as entered by the user
[1] d(CH)/d(to) = r1H+r2H
[2] d(Cm)/d(to) = r1m
[3] d(Cx)/d(to) = r1x+r2x
[4] d(Ct)/d(to) = r2t
[5] d(Cme)/d(to) = r1me+r2me

Explicit equations as entered by the user
[1] k1 = 55.2
[2] k2 = 30.2
[3] r1m = -k1*Cm*(CH^(0.5))
[4] r2t = k2*Cx*(CH^0.5)
[5] r1H = r1m
[6] r2H = -r2t
[7] r1x = -r1m-r2t
[8] r2x = -r2t
[9] r1me = -r1m
[10] r2me = r2t

Comments
[1] d(CH)/d(to) = r1H+r2H
mol/ft3
[2] d(Cm)/d(to) = r1m
mol/ft3
[3] d(Cx)/d(to) = r1x+r2x
mol/ft3
[8] k1 = 55.2
(ft3/lbmol )^0.5/h

Independent variable
variable name : to
initial value : 0
final value : 0.5

Precision
Independent variable accuracy. eps = 0,00001
First stepsize guess. h1 = 0,0001
Minimum allowed stepsize. hmin = 0,00000001
Good steps = 1
Bad steps = 0

General
number of differential equations: 5
number of explicit equations: 10