Excel Solutions To The Chemical Engineering Problem Set

Download as xls, pdf, or txt
Download as xls, pdf, or txt
You are on page 1of 59

EXCEL SOLUTIONS TO THE CHEMICAL ENGINEERING PROBLEM SET

Mathematical Software - Sessions 16 and 116*

Mordechai Shacham, Department of Chemical Engineering, Ben-Gurion University


of the Negev, Beer Sheva, Israel 84105 ([email protected])

Michael B. Cutlip, Department of Chemical Engineering, Box U-3222, University


of Connecticut, Storrs, CT 06269-3222 ([email protected])

No. File name - Description

Linear algebraic equations

1 Demo 1 - Steady-State Material Balances on a Separation Train

Nonlinear algebraic equations

2 Demo 2 - Molar Volume and Compressibility factor from Van Der Waals Equation (Pr=0.5)
3 Demo 3 - Three Pase Equilibrium - Bubble Point
4 Demo 4- Terminal velocity of falling particles
5 Demo 5 - Expediting the Solution of Sytems of Nonlinear Algebraic Equations

Polynomial regression

6 Demo 6a - Fitting Equations to Vapor Pressure Data - Polynomials

Multiple linear regression

7 Demo 6b - Fitting Equations to Vapor Pressure Data - Clapeyron and Riedel equations

Multiple nonlinear regression

8 Demo 6c - Fitting Equations to Vapor Pressure Data - Antoine equation

Ordinary differential equations

9 Demo 7- Heat Exchange in a Series of Tanks


10 Demo 8 - Shooting Method for Solving Two-Point Boundary Value Problems.
Solved by Explicit Euler's method and goal seek
11 Demo 9 - Reversible , Exothermic Gas Phase Reaction in a Catalytic Reactor
12 Demo 10 - Dynamics of a heated tank

Differential algebraic equations

13 Demo 11 - Binary Batch Distillation


Solved by Implicit Euler's Method

Partial diffential equation

14 Demo 12- Method of lines for Partial Differential Equations


Demo 1 - Steady-State Material Balances on a Separation Train

Original A matrix b (base case) -1%Xylene-2%xylene -1%Styrene -2%Styrene


0.07 0.18 0.15 0.24 10.5 10.395 10.29 10.5 10.5
0.04 0.24 0.1 0.65 17.5 17.5 17.5 17.325 17.15
0.54 0.42 0.54 0.1 28 28 28 28 28
0.35 0.16 0.21 0.01 14 14 14 14 14
Inverse A matrix solution (base case)
-2.810484 1.221774 -1.8024194 6.06048387 D1 26.25 26.55 26.84 26.04 25.82
94.3629 -29.83065 -36.685484 41.1370968 B1 17.5 7.59 -2.32 22.72 27.94
-66.04435 20.24597 30.528226 -36.205645 D2 8.75 15.68 22.62 5.21 1.66
-24.50806 9.362903 8.9596774 -9.9919355 B2 17.5 20.07 22.65 15.86 14.22

Molal Flow Mole %


Xylene 4.9875 11.4
D Styrene 5.25 12
Toluene 21.525 49.2
Benzene 11.9875 27.4
Total 43.75 100

Molal Flow Mole %


Xylene 5.5125 21
B Styrene 12.25 46.67
Toluene 6.475 24.67
Benzene 2.0125 7.67
Total 26.25 100
Demo 2 - Molar Volume and Compressibility factor from Van Der Waals Equation
Part a - Solution for P=56 (atm) and T = 450 K

Equations Initial value Solution1

Constants P = 56 56 56 56
R = 0.08206 0.08206 0.08206 0.08206
T = 450 450 450 450
Tc = 405.5 405.5 405.5 405.5
Pc = 111.3 111.3 111.3 111.3
Functions Pr = P/Pc 0.50314 0.50314 0.50314
a = 27*(R^2*Tc^2/Pc)/64 4.19695 4.19695 4.19695
b = R*Tc/(8*Pc) 0.03737 0.03737 0.03737
Unknown V 0.7 0.7 0.57489
Equations Z = P*V/(R*T) 1.06155 1.06155 0.87183
f(V) = (P+a/(V^2))*(V-b)-R*T = 0 5.85576 5.85576 8.4999E-07

1
Solution is obtained by Goal Seek(see under Tools dropdown menu) by setting the value
of the cell b14 at zero while changing cell b4

Demo 2 - Molar Volume and Compressibility factor from Van Der Waals Equation
Part b - Solution for Pr = 1, 2, 4, 10, 20

Parameter Pr 1 2 4 10
Function P = Pr*Pc 111.3 222.6 445.2 1113
Unknown V1 0.23351 0.07727 0.06065 0.05088
Equations Z = P*V/(R*T) 0.70381 0.46578 0.73126 1.53341
f(V) = (P+a/(V^2))*(V-b)-R*T = 0 3.9402E-06 7.6045E-07 2.2077E-06 6.1840E-08

1
Solution is obtained by Goal Seek(see under Tools dropdown menu) by setting the value
of the cell which contains the value of f(V) at zero while changing the cell which contains the value of V
20
2226
0.04618
2.78348
6.9615E-09
Demo 3 - Three Phase Equilibrium - Bubble Point

Equations
Constants A = 1.7; 1.7
B = 0.7; 0.7
z1 = 0.2; 0.2
z2 = 0.8; 0.8
Unknowns x11 0
x12 1
x21 1
x22 0
t 100
beta 0.8
Functions of the unknowns p1 = 10^(7.62231-1417.9/(191.15+t)); 565.3426070855
p2 = 10^(8.10765-1750.29/(235+t)); 763.6664853284
gamma11 = 10^(A*x21*x21/((A*x11/B+x21)^2)); 50.1187233627
gamma21 = 10^(B*x11*x11/((x11+B*x21/A)^2)); 1
gamma12 = 10^(A*x22*x22/((A*x12/B+x22)^2)); 1
gamma22 = 10^(B*x12*x12/((x12+B*x22/A)^2)); 5.0118723363
k11 = gamma11*p1/760; 37.2819075391
k21 = gamma21*p2/760; 1.0048243228
k12 = gamma12*p1/760; 0.7438718514
k22 = gamma22*p2/760; 5.0360512263
y1 = k11*x11; 0
y2 = k21*x21; 1.0048243228
f(1) = x11*(beta+(1-beta)*k11/k12)-z1 ; -0.2
f(2) = x12*k12-x11*k11 ; 0.7438718514
f(3) = x21*(beta+(1-beta)*k21/k22)-z2 ; 0.0399052463
f(4) = x22*k22-x21*k21 ; -1.0048243228
f(5) = x11-y1+x21-y2 ; -0.0048243228
f(6) = (x11-x12)+(x21-x22) ; 0
Sum of squares of errors sum=f1^2+f2^2+f3^2+f4^2+f5^2+f6^2 1.6046329538

1
Solution is obtained by Solver (see under Tools dropdown menu) by searching minimum
of the cell containing the variable sum(b32) while changing x11,x12,x21,x22,t and beta
Initial values Solution1
1.7 1.7
0.7 0.7
0.2 0.2
0.8 0.8
0 0.0226982226
1 0.6867461056
1 0.9772980597
0 0.3132508697
100 88.5377723763
0.8 0.7330014796
565.342607085 357.049430646
763.666485328 498.657698154
50.1187233627 33.3664325023
1 1.0046055616
1 1.102820177
5.0118723363 3.1342255915
37.2819075391 15.6756127995
1.0048243228 0.6591503907
0.7438718514 0.5181069951
5.0360512263 2.0564548934
0 0.3558085493
1.0048243228 0.6441863979
-0.2 -1.6720806E-06
0.7438718514 -5.8807909E-07
0.0399052463 -1.5426696E-06
-1.0048243228 -0.000000114
-0.0048243228 1.33522879E-06
0 -6.9290854E-07
1.6046329538 7.79748123E-12
Solution2
0.7054310618
0.1779556981
0.3743915788
1.06
0.552347277
2.63
5
1.5
1.5
0.1530837848
0.4201773594
0.2422216613
0.0001074354
0.0007723965
0.0005268022
8.8565921E-07
Demo 4- Terminal velocity of falling particles

Constants rho = 994.6;


g = 9.80665;
rhop = 1800;
Dp = 0.208e-3;
vis = 8.931e-4;
Parameter Acceleration (*g)
Unknown vt
Functions of the unknown Re = Dp*vt*rho/vis;
CD= if (Re<0.1) then (24/Re)else (24*(1+0.14*Re^0.7)/Re)
f = vt^2*(3*CD*rho)-4*g*(rhop-rho)*Dp;
1
Solution is obtained by Goal Seek(see under Tools dropdown menu) by setting the value
of the cell b12 at zero while changing cell b4
Part a Part b
Equations Initial values Solution1 Equations Initial values Solution1
994.6 994.6 994.6
9.80665 9.80665 9.80665
1800 1800 1800
0.000208 0.000208 0.000208
0.0008931 0.0008931 0.0008931
1 1 1 30 30 30
0.025 0.025 0.0157824487 0.025 0.025 0.206021
5.7909752547 5.7909752547 3.6558307941 5.7909752547 5.7909752547 47.72261
6.1282521833 6.1282521833 8.842276939 6.1282521833 6.1282521833 1.556618
4.8570587332 4.8570587332 0.000408065 -185.71254242 -185.71254242 1.83E-07
Demo 5 - Expediting the Solution of Sytems of Nonlinear Algebraic Equations

Equations Initial values Solution2

Constants CA0 = 1.5 1.5 1.5 1.5


CB0 = 1.5 1.5 1.5 1.5
KC1 = 1.06 1.06 1.06 1.06
KC2 = 2.63 2.63 2.63 2.63
KC3 = 5 5 5 5
Unknowns CD 0 0 0.705334413
CX 0 0 0.177792441
CZ 0 0 0.373976447
CY = CX+CZ 0 0 0.551768887
CC = CD-CY 0 0 0.153565525
CA = CA0-CD-CZ 1.5 1.5 0.420689141
CB = CB0-CD-CY 1.5 1.5 0.2428967
f(CD) = CC*CD-KC1*CA*CB = 0 -2.385 -2.385 5.176E-09
f(CX) = CX*CY-KC2*CB*CC = 0 0 0 -1.3358E-07
f(CZ) = CZ-KC3*CA*CX = 0 0 0 -2.9923E-07
Sum of squares of errors sum=f(CD)^2+f(CX)^2+f(CZ)^2 5.688225 5.688225 1.07409E-13

1
No convergence from the initial estimates given in the book
Initial estimates closer to the solution had to be used

2
Solution is obtained by Solver (see under Tools dropdown menu) by searching minimum
of the cell containing the variable sum(c20) while changing CD,CX and CZ
Initial values1 Solution2 Initial values Solution2

0 0.055582 10 0.7054310618
0.5 0.597286 10 0.1779556981
1 1.082034 10 0.3743915788
1.5 1.67932 20 0.552347277
-1.5 -1.623737 -10 0.1530837848
0.5 0.362383 -18.5 0.4201773594
0 -0.234902 -28.5 0.2422216613
0 -1.92E-05 -658.885 0.0001074354
0.75 -0.0001 -549.55 0.0007723965
-0.25 -0.000197 935 0.0005268022
0.625 4.93E-08 1610359.645725 8.8565921E-07
Demo 6. Fitting Polynomials and Correlation Equations to Vapor Pressure Data
Part a. Polynomials

T (deg. C) P (mm. Hg) TK TK2 TK3 TK4 TK5


-36.7 1 236.45 55908.6025 13219589.061 3125771834 7.391E+11
-19.6 5 253.55 64287.6025 16300121.614 4132895835 1.048E+12
-11.5 10 261.65 68460.7225 17912748.042 4686870525 1.226E+12
-2.6 20 270.55 73197.3025 19803530.191 5357845093 1.45E+12
7.6 40 280.75 78820.5625 22128872.922 6212681073 1.744E+12
15.4 60 288.55 83261.1025 24024991.126 6932411190 2E+12
26.1 100 299.25 89550.5625 26798005.828 8019303244 2.4E+12
42.2 200 315.35 99445.6225 31360177.055 9889431834 3.119E+12
60.6 400 333.75 111389.0625 37176099.609 1.2408E+10 4.141E+12
80.1 760 353.25 124785.5625 44080499.953 1.5571E+10 5.501E+12

a2 a1 a0
coeff.s 0.0861525706 -44.998001397 5862.6992453
std.dev.s 0.0090580486 5.3598069732 784.266243103
R2, SE (y) 0.9841859905 34.8824456904 #N/A
95% conf. i 0.02088786 12.3597148803 1808.5179566
Variance 1216.7850173

a3 a2 a1 a0
coeff.s 0.0007449112 -0.5709694343 146.383159187 -12540.815474
std.dev.s 4.772551E-05 0.0421282801 12.2943522468 1186.3813337
R2, SE (y) 0.9996198817 5.8414234444 #N/A #N/A
95% conf. i 0.0001128517 0.0996165312 29.0712253228 2805.3173017
Variance 34.122227857

a4 a3 a2 a1 a0
coeff.s 3.963148E-06 -0.0039170137 1.4716764942 -248.69004375 15917.8482
std.dev.s 1.755758E-07 0.0002065973 0.0906077788 17.552853021 1267.27108
R2, SE (y) 0.999996306 0.630809207 #N/A #N/A #N/A
95% conf. i 4.296163E-07 0.000505523 0.2217081739 42.950076058 3100.88561
Variance 0.3979202556

a5 a4 a3 a2 a1 a0
coeff.s -2.5051E-09 7.6474065E-06 -0.0060735 2.0995932784 -339.63942 21160.111
std.dev.s 8.142829E-09 1.1977257E-05 0.0070133773 2.0434983789 296.266634 17097.426
R2, SE (y) 0.9999963914 0.6970675689 #N/A #N/A #N/A #N/A
95% conf. i 2.093196E-08 3.0788738E-05 0.0180285876 5.2530169329 761.58301 43950.644
Variance 0.4859031956
Cpcalc (4) residual(4)
1.0477125 -0.047712495
4.51835895 0.4816410459 Vapor Pressure of Benzene - Residual Plot
10.4153728 -0.4153727642
20.7392276 -0.7392275768 1
39.1623364 0.8376635717 0.8
59.694174 0.3058259535 0.6
100.338419 -0.3384191186
0.4
200.264723 -0.2647226584
0.2

Residual
399.767921 0.2320792791
760.051755 -0.0517552375 0
-0.2 0 100 200 300 400 500 600 700 800
Degrees of t-distibution
-0.4
Freedom for 95% conf. Int
1 12.7062 -0.6
2 4.3027 -0.8
3 3.1824 -1
4 2.7764
Heat Capacity
5 2.5706
6 2.4469
7 2.3646
8 2.306
9 2.2622
10 2.2281
Demo 6- Fitting Equations to Vapor Pressure Data
Part b. Clapeyron and Riedel Equations
Clapeyron Equation
T (deg. C) TK P (mm. Hg) Trec logP logP(calc) residual
-36.7 236.45 1 0.004229224 0 0.1441392725 -0.1441392725
-19.6 253.55 5 0.003943995 0.69897 0.7246740368 -0.0257040368
-11.5 261.65 10 0.003821899 1 0.9731793734 0.0268206266
-2.6 270.55 20 0.003696174 1.30103 1.2290714177 0.0719585823
7.6 280.75 40 0.003561888 1.60206 1.502389222 0.099670778
15.4 288.55 60 0.003465604 1.7781513 1.6983588706 0.0797924294
26.1 299.25 100 0.003341688 2 1.9505696212 0.0494303788
42.2 315.35 200 0.00317108 2.30103 2.2978129703 0.0032170297
60.6 333.75 400 0.002996255 2.60206 2.6536398628 -0.0515798628
80.1 353.25 760 0.002830856 2.8808136 2.9902802527 -0.1094666527

a1 a0
coeff.s -2035.331034 8.752010003
std.dev.s 66.621127947 0.235184304 VP of Benzene- Clapeyron - residual plot
R2, SE (y) 0.9915015982 0.087129536
0.15
95% conf. int. 153.62832105 0.542335004
Variance 0.0075915561 0.1

0.05
Degrees of t-distibution
Residual

0
Freedom for 95% conf. Int
0 0.5 1 1.5 2 2.5
1 12.7062 -0.05
2 4.3027
3 3.1824 -0.1
4 2.7764
5 2.5706 -0.15
6 2.4469
-0.2
7 2.3646
8 2.306 Log(P) - measured
9 2.2622
10 2.2281
Riedels' Equation
Trec logT T2 logP logP(calc)
0.004229224 2.373739318 55908.6025 0 -0.0107226
0.003943995 2.404063615 64287.6025 0.69897 0.722590957
0.003821899 2.417720739 68460.7225 1 1.011328031
0.003696174 2.432247538 73197.3025 1.30103 1.293057333
0.003561888 2.448319765 78820.5625 1.60206 1.576922101
0.003465604 2.460221079 83261.1025 1.7781513 1.770027635
0.003341688 2.476034159 89550.5625 2 2.006505993
0.00317108 2.498792835 99445.6225 2.30103 2.312241806
0.002996255 2.523421274 111389.0625 2.60206 2.606704193
0.002830856 2.548082171 124785.5625 2.8808136 2.875459448

a3 a2 a1 a0
coeff.s 4.44453E-05 -75.73554721 -9317.202915 216.6853854
on - residual plot std.dev.s 2.04966E-05 23.93880773 1989.887826 64.08471065
R2, SE (y) 0.999750783 0.017228791 #N/A #N/A
95% conf. int. 5.02E-05 58.57586864 4869.056521 156.8088785
Variance 0.000296831

VP of Benzene- Riedel - residual plot


1.5 2 2.5 3
0.03

0.02

0.01
Residual

0
measured
0 0.5 1 1.5 2 2.5 3
-0.01

-0.02

-0.03
Log(P) - measured
residual
0.010723
-0.023621
-0.011328
0.007973
0.025138
0.008124
-0.006506
-0.011212
-0.004644
0.005354
Demo 6 - Fitting Equations to Vapor Pressure Data
Part c. - Antoine equation
T (deg. C) P (mm. Hg)
No. of data points (N) 10 10^(A+B/(TC+C))
No. of parameters (m) 3 -36.7 1
-19.6 5
-11.5 10
-2.6 20
Solution1 Initial estimate 7.6 40
A= 6.618507372 8.7520099666 15.4 60
B= -1054.9777433 -2035.3310261 26.1 100
C= 202.1405011427 273 42.2 200
Sum(Err^2)= 15.9731581096 47701.117369 60.6 400
Variance 2.2818797299 2805.9480805 80.1 760

1
Solution is obtained by using the solver to minimize Sum(Err^2)
while changing A, B and C. The initial estimate from Clapeyron's
equation were used

Nonlinear Regression Results for Antoine Equation

800
700
600
Pressure (mmHg)

500
400
300
200
100
0
-60 -40 -20 0 20 40 60 80 100
Temperature (C)
(Pv)calc Residual Squares of
10^(A+B/(TC+C)) (Pv)-(Pv)calc the Residuals
1.744724994 -0.7447249939 0.5546153166
6.903828597 -1.9038285969 3.6245633265
12.15200385 -2.1520038479 4.6311205612
21.45219427 -1.4521942706 2.1088681996
38.77825989 1.2217401057 1.4926488859
58.74055394 1.2594460612 1.5862043812
99.14884529 0.8511547128 0.7244643452
199.9186761 0.0813239002 0.0066135767
401.0729043 -1.0729042534 1.151123537
759.695146 0.3048540302 0.0929359797

Residual Plot of the Antoine Equation for Benzene


1.5
1
0.5
0
Residual

0 100 200 300 400 500 600 700 800


-0.5 (Pv)-(Pv)calc
-1
-1.5
-2
-2.5
Vapor Pressure (mmHg)
enzene

800
(Pv)-(Pv)calc
Degrees of t-distibution
Freedom for 95% conf. Int
1 12.7062
2 4.3027
3 3.1824
4 2.7764
5 2.5706
6 2.4469
7 2.3646
8 2.306
9 2.2622
10 2.2281
11 2.201
12 2.1788
13 2.1604
14 2.1448
15 2.1315
16 2.1199
17 2.1098
18 2.1009
19 2.093
20 2.086
21 2.0796
22 2.0739
23 2.0687
24 2.0639
25 2.0595
Demo 7 - Heat Exchange in a Series of Tanks Step no. (i) t T1(t) T2(t)
0 0 20 20
Definitions Equations/values 1 1 21.091682655 21.14531614
2 2 22.074551435 22.27475879
3 3 22.959452403 23.38011342
Final value (Ind. Var.) tf=200 200 4 4 23.756150536 24.45496079
Integration Stepsize h 1 5 5 24.473437487 25.49440063
Constants W=100 100 6 6 25.119228592 26.49481267
Cp=2.0 2 7 7 25.700650229 27.4536501
T0=20 20 8 8 26.224118447 28.36926158
UA=10. 10 9 9 26.695409776 29.24073801
Tsteam=250 250 10 10 27.119724965 30.06778098
M=1000 1000 11 11 27.50174638 30.85059006
Independent Variable t 0 1 12 12 27.845689667 31.58976649
Variables T1(0)=20 20 21.0916826555 13 13 28.15535028 32.28623113
T2(0)=20 20 21.1453161398 14 14 28.434145355 32.94115474
T3(0)=20 20 21.147081869 15 15 28.685151427 33.55589906
Differential equations k11=h*d(T1)/d(t)=h*(W*Cp*(T0-T1)+UA*(Tsteam-T1))/(M*Cp) 1.15 1.0353733212 16 16 28.9111383755 34.13196708
k12=h*d(T2)/d(t)=h*(W*Cp*(T1-T2)+UA*(Tsteam-T2))/(M*Cp) 1.15 1.1389100709 17 17 29.114599991 34.67096135
k13=h*d(T3)/d(t)=h*(W*Cp*(T2-T3)+UA*(Tsteam-T3))/(M*Cp) 1.15 1.1440880177 18 18 29.297781493 35.17454921
T1+k11/2 20.575 21.6093693161 19 19 29.462704311 35.64443396
T2+k12/2 20.575 21.7147711753 20 20 29.611188386 36.08233109
T3+k13/2 20.575 21.7191258779 21 21 29.744872255 36.48994897
k21 1.089625 0.9810162218 22 22 29.865231135 36.86897321
k22 1.147125 1.1308859582 23 23 29.973593198 37.22105422
k23 1.147125 1.1409689004 24 24 30.071154233 37.54779746
T1+k21/2 20.5448125 21.5821907664 25 25 30.158990835 37.85075606
T2+k22/2 20.5735625 21.7107591189 26 26 30.23807229 38.13142525
T3+k23/2 20.5735625 21.7175663192 27 27 30.309271273 38.3912385
k31 1.0927946875 0.9838699695 28 28 30.37337347 38.63156494
k32 1.1442571875 1.1285893691 29 29 30.431086256 38.85370801
k33 1.1471321875 1.1407314484 30 30 30.4830465 39.05890487
T1+k31 21.0927946875 22.075552625 31 31 30.529827586 39.24832675
T2+k32 21.1442571875 22.273905509 32 32 30.571945751 39.42307972
T3+k33 21.1471321875 22.2878133174 33 33 30.609865771 39.5842061
k41 1.0352565578 0.9320669744 34 34 30.6440061 39.73268616
k42 1.1391324641 1.1187951841 35 35 30.674743478 39.86944014
k43 1.1439768391 1.1371701526 36 36 30.702417097 39.99533047
37 37 30.727332337 40.11116413
38 38 30.749764142 40.2176952
39 39 30.769960048 40.31562733
Heat Exchange in a Series of Tanks

55
50
eg. C)

45
T1(t)
Heat Exchange in a Series of Tanks 40 40 30.78814292 40.40561631
41 41 30.804513407 40.48827265
55 42 42 30.81925216 40.56416405
43 43 30.832521822 40.63381791
50 44 44 30.844468825 40.69772368
Temperature (deg. C)

45 45
46
45
46
30.855225007
30.864909062
40.75633524
40.81007307
T1(t)
40 T2(t) 47 47 30.873627855 40.85932643
48 48 30.881477599 40.90445542
35 T3(t) 49 49 30.888544918 40.94579289
30 50
51
50
51
30.894907798
30.900636456
40.98364631
41.01829954
25 52 52 30.905794108 41.05001446
53 53 30.910437669 41.07903254
20 54 54 30.914618382 41.1055763
0 20 40 60 80 55 55 30.91838238 41.12985069
Time (min) 56 56 30.9217712 41.15204441
57 57 30.924822239 41.17233106
58 58 30.927569164 41.19087033
59 59 30.930042288 41.20780899
60 60 30.932268903 41.22328195
61 61 30.934273579 41.23741311
62 62 30.936078438 41.25031622
63 63 30.937703397 41.2620957
64 64 30.939166388 41.27284735
65 65 30.940483555 41.28265905
66 66 30.941669432 41.29161135
67 67 30.942737107 41.29977809
68 68 30.943698361 41.30722694
69 69 30.944563801 41.31401986
70 70 30.945342979 41.32021361
71 71 30.946044491 41.32586015
72 72 30.946676081 41.33100703
73 73 30.947244716 41.33569775
74 74 30.947756672 41.33997212
75 75 30.948217599 41.34386653
76 76 30.948632583 41.34741426
77 77 30.949006203 41.35064573
78 78 30.949342582 41.35358875
79 79 30.949645433 41.3562687
80 80 30.949918097 41.35870879
T3(t) k11 k12 k13 T1(t)+hk11/2 T2(t)+hk12/2 T3(t)+hk13/2 k21 k22 k23 T1(t)+hk21/2 T2(t)+hk22/2 T3(t)+hk23/2 k31
20 1.15 1.15 1.15 20.575 20.575 20.575 1.089625 1.147125 1.147125 20.5448125 20.5735625 20.5735625 1.092795
21.14708187 1.035373 1.13891 1.144088 21.60936932 21.714771175 21.71912588 0.981016 1.130886 1.140969 21.58219077 21.71075912 21.717566319 0.98387
22.28785835 0.932172 1.118605 1.137251 22.54063748 22.834061527 22.85648372 0.883233 1.106487 1.133475 22.51616797 22.82800244 22.854596028 0.885802
23.42103886 0.839257 1.091033 1.128802 23.37908115 23.92563009 23.98543999 0.795196 1.075717 1.124092 23.35705064 23.9179719 23.983084768 0.79751
24.54477075 0.755604 1.057844 1.118295 24.13395263 24.983882873 25.10391833 0.715935 1.040088 1.112477 24.11411802 24.97500457 25.101009183 0.718018
25.65685204 0.680289 1.020432 1.105471 24.81358202 26.004616472 26.20958734 0.644574 1.000873 1.098455 24.79572443 25.99483737 26.206079527 0.646449
26.75489864 0.612481 0.979968 1.090217 25.42546909 26.984796431 27.30000709 0.580326 0.959143 1.081979 25.40939147 26.97438431 27.295888089 0.582014
27.83647381 0.551432 0.937432 1.072535 25.97636609 27.922365978 28.37274144 0.522482 0.915788 1.063099 25.96189101 27.91154419 28.368023179 0.524001
28.89918646 0.496468 0.893639 1.052512 26.47235223 28.816081266 29.42544225 0.470403 0.871547 1.041937 26.45931996 28.80503492 29.420154804 0.471771
29.94076407 0.446982 0.849263 1.030294 26.91890076 29.665369751 30.45591086 0.423515 0.827026 1.018666 26.90716749 29.65425113 30.450097241 0.424747
30.9591048 0.402429 0.804855 1.006072 27.3209394 30.470208727 31.46214084 0.381301 0.782722 0.993496 27.31037565 30.45914199 31.45585284 0.382411
31.95231283 0.362317 0.760863 0.980066 27.68290469 31.231021404 32.44234591 0.343295 0.739033 0.966656 27.67339388 31.22010667 32.435640743 0.344294
32.9187204 0.326203 0.717643 0.952511 28.00879096 31.948588237 33.3949759 0.309077 0.696277 0.938386 28.00022814 31.93790516 33.387913572 0.309976
33.85689905 0.293688 0.675481 0.923649 28.30219439 32.623971506 34.31872341 0.27827 0.654702 0.908931 28.29448507 32.61358234 34.311364646 0.279079
34.76566274 0.264415 0.634593 0.893721 28.56635272 33.258451381 35.21252318 0.250533 0.614498 0.87853 28.55941184 33.24840368 35.204927838 0.251262
35.64406441 0.238059 0.595146 0.862963 28.80418098 33.853471929 36.07554598 0.225561 0.575804 0.847415 28.79793193 33.84380083 36.067771842 0.226217
36.49138787 0.21433 0.557257 0.831601 29.01830361 34.410595725 36.90718836 0.203078 0.538718 0.815805 29.01267744 34.40132598 36.89929027 0.203669
37.30713611 0.192967 0.521009 0.799847 29.211083491 34.931465878 37.70705953 0.182836 0.503304 0.783905 29.20601811 34.92261357 37.699088774 0.183368
38.09101714 0.173733 0.48645 0.767898 29.38464796 35.417774455 38.4749662 0.164612 0.469598 0.751906 29.38008748 35.40934845 38.466970133 0.165091
38.84292832 0.156416 0.453605 0.735936 29.54091233 35.871236392 39.21089628 0.148204 0.437611 0.71998 29.53680641 35.86323967 39.202918084 0.148635
39.56293969 0.140825 0.422474 0.704124 29.681601 36.293568128 39.91500191 0.133432 0.407335 0.688282 29.67790433 36.28599881 39.907080499 0.13382
40.25127696 0.126788 0.393043 0.672611 29.80826646 36.686470266 40.58758237 0.120132 0.378747 0.656951 29.80493827 36.67932261 40.579752398 0.120481
40.90830451 0.114151 0.365281 0.641525 29.9223065 37.051613677 41.22906718 0.108158 0.351811 0.626109 29.91931004 37.04487882 41.221359163 0.108472
41.53450881 0.102773 0.339149 0.610982 30.02497956 37.39062853 41.8399998 0.097377 0.326482 0.595863 30.02228177 37.3842952 41.832440242 0.09766
42.13048241 0.092529 0.314597 0.581079 30.11741864 37.705095807 42.42102195 0.087671 0.302707 0.566302 30.11498975 37.69915086 42.413633544 0.087926
42.69690872 0.083306 0.29157 0.5519 30.20064382 37.996540912 42.97285882 0.078932 0.280428 0.537504 30.19845703 37.99096986 42.965660681 0.079162
43.23454773 0.075002 0.270008 0.523515 30.2755735 38.266429044 43.49630524 0.071065 0.259582 0.509531 30.27360468 38.2612164 43.489313159 0.071272
43.74422264 0.067527 0.249847 0.49598 30.34303453 38.516162041 43.99221287 0.063981 0.240106 0.482434 30.34126196 38.51129172 43.985439564 0.064167
44.22680758 0.060796 0.231023 0.469342 30.40377136 38.747076457 44.46147843 0.057604 0.221934 0.456252 30.40217547 38.742532 44.454933785 0.057772
44.68321637 0.054736 0.213469 0.443633 30.45845423 38.960442648 44.90503291 0.051862 0.204999 0.431016 30.45701741 38.95620748 44.898724276 0.052013
45.11439228 0.04928 0.19712 0.418879 30.50768656 39.157464691 45.32383193 0.046693 0.189235 0.406744 30.50639296 39.1535223 45.317764335 0.046829
45.52129879 0.044368 0.181908 0.395096 30.552011638 39.33928097 45.71884694 0.042039 0.174577 0.383449 30.55084698 39.33561508 45.713023376 0.042161
45.90491141 0.039946 0.167771 0.372292 30.5919186 39.506965318 46.09105754 0.037849 0.160961 0.361135 30.59087002 39.50355997 46.08547915 0.037959
46.26621028 0.035964 0.154645 0.350469 30.62784782 39.661528568 46.44144455 0.034076 0.148324 0.339801 30.62690376 39.65836824 46.43611087 0.034175
46.60617382 0.032379 0.142469 0.32962 30.66019578 39.803920446 46.770984 0.030679 0.136608 0.319439 30.65934582 39.80099013 46.765893178 0.030769
46.92577303 0.029152 0.131183 0.309738 30.68931945 39.935031711 47.08064195 0.027621 0.125754 0.300036 30.68855421 39.93231695 47.075790912 0.027702
47.2259667 0.026246 0.120732 0.290807 30.7155402 40.055696472 47.37136997 0.024868 0.115706 0.281576 30.71485124 40.05318341 47.366754596 0.024941
47.5076972 0.02363 0.111061 0.272808 30.73914739 40.166694634 47.64410131 0.02239 0.106412 0.264039 30.7385271 40.16437003 47.639716615 0.022455
47.77188703 0.021275 0.102118 0.255721 30.76040152 40.268754414 47.89974772 0.020158 0.097821 0.247402 30.75984306 40.26660567 47.895587994 0.020216
48.01943584 0.019154 0.093855 0.239522 30.77953715 40.362554899 48.13919682 0.018149 0.089885 0.23164 30.77903435 40.36057006 48.135255749 0.018201
48.25121808 0.017245 0.086225 0.224184 30.79676542 40.448728603 48.36330995 0.01634 0.08256 0.216725 30.79631274 40.44689633 48.359580738 0.016387
48.4680811 0.015526 0.079183 0.209679 30.81227645 40.527864006 48.57292047 0.014711 0.075802 0.20263 30.81186889 40.52617361 48.569395971 0.014754
48.67084363 0.013979 0.072688 0.195978 30.82624142 40.600508047 48.76883254 0.013245 0.069571 0.189323 30.82587449 40.59894945 48.765505323 0.013283
48.86029475 0.012585 0.066701 0.183051 30.83881443 40.667168559 48.95182017 0.011924 0.063829 0.176776 30.83848406 40.66573228 48.948682615 0.011959
49.03719307 0.011331 0.061186 0.170867 30.85013421 40.728316632 49.12262662 0.010736 0.05854 0.164956 30.84983678 40.72699377 49.119671006 0.010767
49.20226635 0.010201 0.056107 0.159396 30.86032569 40.784388891 49.28196412 0.009666 0.053672 0.153833 30.86005791 40.78317111 49.279182674 0.009694
49.35621122 0.009185 0.051433 0.148605 30.86950134 40.835789686 49.43051379 0.008702 0.049192 0.143375 30.86926024 40.83466918 49.427898733 0.008728
49.49969331 0.008269 0.047134 0.138465 30.87776239 40.882893188 49.56892573 0.007835 0.045072 0.133552 30.87754533 40.88186266 49.56646937 0.007858
49.6333474 0.007445 0.04318 0.128944 30.88520003 40.92604539 49.69781944 0.007054 0.041285 0.124333 30.8850046 40.92509804 49.695514153 0.007075
49.75777787 0.006703 0.039546 0.120013 30.89189631 40.965566005 49.81778418 0.006351 0.037805 0.115689 30.89172036 40.96469549 49.8156225 0.006369
49.8735592 0.006035 0.036208 0.111641 30.89792514 41.001750267 49.92937965 0.005718 0.034609 0.10759 30.89776673 41.00095068 49.927354277 0.005734
49.98123664 0.005433 0.033142 0.1038 30.90335304 41.034870637 50.03313669 0.005148 0.031674 0.100008 30.90321042 41.03413648 50.031240497 0.005163
50.08132699 0.004892 0.030328 0.096462 30.90823992 41.065178407 50.12955805 0.004635 0.02898 0.092914 30.90811151 41.06450459 50.127784113 0.004648
50.17431938 0.004404 0.027745 0.0896 30.91263969 41.092905213 50.21911923 0.004173 0.026509 0.086283 30.91252409 41.092287 50.217460876 0.004185
50.26067617 0.003965 0.025376 0.083187 30.91660092 41.11826446 50.30226949 0.003757 0.024242 0.080088 30.91649683 41.11769746 50.300720246 0.003768
50.34083393 0.00357 0.023204 0.077198 30.92016731 41.141452651 50.37943269 0.003382 0.022164 0.074305 30.9200736 41.14093279 50.377986349 0.003392
50.41520436 0.003214 0.021212 0.071608 30.92337821 41.162650639 50.45100835 0.003045 0.02026 0.068909 30.92329384 41.16217416 50.449658951 0.003054
50.48417529 0.002894 0.019387 0.066395 30.92626907 41.182024795 50.51737264 0.002742 0.018514 0.063878 30.92619311 41.18158822 50.516114461 0.00275
50.54811169 0.002605 0.017716 0.061535 30.92887178 41.199728093 50.57887935 0.002468 0.016916 0.05919 30.9288034 41.19932819 50.577706933 0.002476
50.60735671 0.002346 0.016184 0.057008 30.93121507 41.215901137 50.63586093 0.002222 0.015452 0.054825 30.9311535 41.21553494 50.63476907 0.002229
50.66223262 0.002112 0.014782 0.052794 30.93332479 41.230673096 50.6886295 0.002001 0.014112 0.050761 30.93326935 41.23033786 50.687613224 0.002007
50.71304183 0.001901 0.013499 0.048872 30.93522422 41.244162597 50.73747779 0.001801 0.012885 0.046981 30.93517431 41.24385578 50.736532379 0.001807
50.7600679 0.001712 0.012325 0.045224 30.93693432 41.256478537 50.78268015 0.001622 0.011763 0.043466 30.93688939 41.25619781 50.781801122 0.001627
50.80357645 0.001541 0.01125 0.041834 30.93847397 41.267720844 50.82449347 0.00146 0.010737 0.0402 30.93843351 41.26746405 50.823676585 0.001464
50.84381611 0.001388 0.010268 0.038684 30.93986015 41.277981188 50.86315813 0.001315 0.009798 0.037167 30.93982373 41.27774635 50.86239937 0.001319
50.88101945 0.001249 0.009369 0.035759 30.94110817 41.287343629 50.89889888 0.001184 0.00894 0.03435 30.94107538 41.28712892 50.898194442 0.001187
50.91540383 0.001125 0.008548 0.033044 30.94223179 41.295885227 50.9319257 0.001066 0.008155 0.031736 30.94220226 41.29568897 50.931271994 0.001069
50.94717227 0.001013 0.007797 0.030525 30.94324341 41.303676599 50.96243463 0.000959 0.007438 0.029312 30.94321683 41.30349724 50.96182828 0.000962
50.97651425 0.000912 0.007111 0.028189 30.9441542 41.310782442 50.9906086 0.000864 0.006783 0.027064 30.94413027 41.31061857 50.990046419 0.000866
51.00360652 0.000821 0.006484 0.026023 30.9449742 41.317262007 51.01661817 0.000778 0.006185 0.024981 30.94495266 41.31711231 51.016097165 0.00078
51.02861383 0.000739 0.005912 0.024017 30.94571247 41.323169546 51.04062229 0.0007 0.005638 0.023052 30.94569307 41.32303283 51.04013964 0.000702
51.05168968 0.000665 0.005389 0.022159 30.94637716 41.328554717 51.06276898 0.00063 0.005139 0.021265 30.94635969 41.32842989 51.062322042 0.000632
51.07297696 0.000599 0.004912 0.020438 30.94697559 41.333462961 51.08319602 0.000568 0.004684 0.019611 30.94695986 41.333349 51.082782314 0.000569
51.09260865 0.000539 0.004476 0.018846 30.94751437 41.337935852 51.10203158 0.000511 0.004268 0.01808 30.94750021 41.33783183 51.10164878 0.000512
51.11070841 0.000486 0.004079 0.017373 30.94799945 41.342011414 51.11939483 0.00046 0.003889 0.016665 30.9479867 41.34191649 51.119040755 0.000461
51.12739121 0.000437 0.003716 0.016011 30.94843618 41.345724415 51.1353965 0.000414 0.003543 0.015356 30.9484247 41.3456378 51.135069119 0.000415
51.14276385 0.000394 0.003385 0.014751 30.94882937 41.349106641 51.15013946 0.000373 0.003227 0.014146 30.94881904 41.34902763 51.149836859 0.000374
51.15692549 0.000354 0.003083 0.013587 30.94918338 41.352187142 51.16371919 0.000336 0.002939 0.013028 30.94917408 41.35211508 51.163439587 0.000337
51.16996818 0.000319 0.002807 0.012512 30.9495021 41.354992466 51.17622429 0.000302 0.002676 0.011996 30.94949372 41.35492675 51.175966026 0.000303
51.18197731 0.000287 0.002556 0.011519 30.94978905 41.357546862 51.18773694 0.000272 0.002436 0.011042 30.94978151 41.35748694 51.187498467 0.000273
51.19303208 0.000259 0.002327 0.010603 30.9500474 41.359872483 51.19833333 0.000245 0.002218 0.010162 30.95004061 41.35981785 51.198113202 0.000246
k32 k33 T1(t)+hk31 T2(t)+hk32 T3(t)+hk33 k41 k42 k43
1.144257 1.147132 21.09279469 21.1442572 21.14713219 1.035257 1.139132 1.143977
1.128589 1.140731 22.07555262 22.2739055 22.28781332 0.932067 1.118795 1.13717
1.104677 1.133068 22.9603538 23.3794353 23.42092601 0.839163 1.091195 1.128746
1.074318 1.123573 23.75696209 24.4544314 24.54461215 0.755519 1.057981 1.118259
1.039036 1.111894 24.47416814 25.4939971 25.65666524 0.680212 1.020547 1.10545
1.000115 1.097845 25.11988642 26.4945152 26.75469742 0.612412 0.980065 1.090208
0.958629 1.08137 25.70124249 27.4534415 27.83626882 0.55137 0.937513 1.072536
0.915477 1.062512 26.22465167 28.3691271 28.89898579 0.496412 0.893707 1.052519
0.871403 1.041387 26.69588985 29.2406649 29.9405737 0.446932 0.849319 1.030306
0.82702 1.018165 27.12015719 30.0677584 30.95892898 0.402383 0.804901 1.006088
0.782828 0.99305 27.50213552 30.8506086 31.95215445 0.362276 0.7609 0.980085
0.739228 0.966268 27.84604002 31.5898183 32.91858122 0.326166 0.717673 0.952531
0.696543 0.93806 28.15566571 32.2863093 33.85677999 0.293655 0.675504 0.923669
0.655022 0.908665 28.43442935 32.9412535 34.765564 0.264385 0.634611 0.893741
0.614859 0.878323 28.68540711 33.5560135 35.64398568 0.238032 0.595159 0.862983
0.576194 0.847264 28.91136857 34.1320932 36.49132845 0.214306 0.557267 0.83162
0.539129 0.815707 29.11480724 34.6710956 37.30709499 0.192945 0.521016 0.799865
0.503727 0.783857 29.29796809 35.1746887 38.09099314 0.173713 0.486454 0.767915
0.470027 0.751903 29.46287231 35.6445764 38.84292012 0.156398 0.453607 0.735951
0.43804 0.720018 29.61133964 36.0824744 39.56294589 0.140809 0.422474 0.704138
0.407761 0.688356 29.74500843 36.4900916 40.25129612 0.126774 0.393041 0.672623
0.379165 0.657058 29.86535374 36.8691139 40.90833522 0.114138 0.365278 0.641536
0.352219 0.626245 29.97370358 37.2211919 41.53454968 0.102761 0.339145 0.610991
0.326877 0.596023 30.07125361 37.5479314 42.1305321 0.092518 0.314593 0.581087
0.303088 0.566484 30.15908031 37.8508856 42.69696597 0.083297 0.291565 0.551907
0.280794 0.537703 30.23815285 38.1315499 43.23461134 0.074994 0.270003 0.523521
0.259933 0.509744 30.3093438 38.391358 43.74429149 0.067519 0.249842 0.495985
0.240441 0.482658 30.37343877 38.6316791 44.22688066 0.060789 0.231018 0.469345
0.222252 0.456485 30.43114505 38.8538166 44.68329273 0.05473 0.213464 0.443636
0.2053 0.431255 30.48309943 39.059008 45.11447107 0.049275 0.197114 0.418881
0.189519 0.406987 30.52987524 39.2484243 45.52137925 0.044363 0.181903 0.395098
0.174845 0.383694 30.57198865 39.4231719 45.90499285 0.039941 0.167766 0.372293
0.161213 0.361381 30.6099044 39.5842929 46.26629209 0.03596 0.15464 0.350469
0.148562 0.340045 30.64404088 39.7327678 46.60625546 0.032376 0.142463 0.32962
0.136831 0.31968 30.67477479 39.8695168 46.92585404 0.029149 0.131178 0.309737
0.125962 0.300274 30.70244529 39.9954023 47.22604668 0.026243 0.120727 0.290805
0.115901 0.281809 30.72735772 40.1112313 47.5077758 0.023627 0.111056 0.272807
0.106594 0.264267 30.74978699 40.217758 47.77196396 0.021272 0.102114 0.25572
0.097991 0.247624 30.76998062 40.3156859 48.01951086 0.019152 0.093851 0.23952
0.090044 0.231855 30.78816144 40.4056709 48.25129099 0.017243 0.086221 0.224182
0.082707 0.216934 30.80453008 40.4883235 48.46815174 0.015524 0.079179 0.209676
0.075939 0.202831 30.81926717 40.5642113 48.67091188 0.013977 0.072685 0.195975
0.069698 0.189517 30.83253534 40.6338618 48.86036052 0.012584 0.066698 0.183048
0.063947 0.176962 30.84448099 40.6977644 49.0372563 0.011329 0.061183 0.170865
0.058649 0.165134 30.85523596 40.756373 49.20232699 0.0102 0.056104 0.159393
0.053773 0.154003 30.86491893 40.8101081 49.35626928 0.009184 0.051431 0.148603
0.049286 0.143538 30.87363674 40.8593588 49.49974877 0.008268 0.047131 0.138462
0.045159 0.133707 30.8814856 40.9044854 49.63340029 0.007444 0.043178 0.128942
0.041365 0.124481 30.88855212 40.9458206 49.75782822 0.006702 0.039544 0.12001
0.037879 0.115829 30.89491428 40.9836719 49.87360706 0.006034 0.036206 0.111638
0.034677 0.107723 30.90064229 41.0183232 49.98128206 0.005433 0.03314 0.103798
0.031737 0.100133 30.90579936 41.0500363 50.08137004 0.004891 0.030326 0.09646
0.029038 0.093033 30.9104424 41.0790526 50.17436012 0.004404 0.027744 0.089597
0.026562 0.086395 30.91462264 41.1055948 50.26071468 0.003965 0.025375 0.083184
0.024291 0.080194 30.91838621 41.1298677 50.34087029 0.003569 0.023203 0.077195
0.022209 0.074405 30.92177465 41.1520601 50.41523865 0.003214 0.021211 0.071606
0.020301 0.069003 30.92482535 41.1723455 50.48420758 0.002893 0.019386 0.066393
0.018553 0.063967 30.92757196 41.1908836 50.54814209 0.002605 0.017714 0.061533
0.016951 0.059274 30.93004481 41.2078212 50.60738529 0.002345 0.016183 0.057007
0.015484 0.054903 30.93227117 41.2232932 50.66225945 0.002112 0.014781 0.052792
0.014141 0.050834 30.93427562 41.2374234 50.71306701 0.001901 0.013498 0.04887
0.012913 0.04705 30.93608028 41.2503257 50.76009151 0.001712 0.012324 0.045223
0.011788 0.043531 30.93770505 41.2621044 50.80359857 0.001541 0.01125 0.041833
0.01076 0.04026 30.93916788 41.2728553 50.84383681 0.001387 0.010267 0.038683
0.009819 0.037223 30.9404849 41.2826664 50.88103881 0.001249 0.009369 0.035758
0.008959 0.034402 30.94167064 41.2916181 50.91542193 0.001125 0.008547 0.033043
0.008173 0.031785 30.94273819 41.2997842 50.94718917 0.001012 0.007796 0.030524
0.007454 0.029358 30.94369934 41.3072326 50.97653002 0.000912 0.007111 0.028188
0.006798 0.027107 30.94456468 41.314025 51.00362123 0.000821 0.006484 0.026022
0.006198 0.025021 30.94534377 41.3202183 51.02862755 0.000739 0.005911 0.024016
0.005651 0.023089 30.94604521 41.3258645 51.05170245 0.000665 0.005389 0.022158
0.005151 0.021299 30.94667672 41.331011 51.07298885 0.000599 0.004912 0.020437
0.004694 0.019643 30.9472453 41.3357014 51.09261971 0.000539 0.004476 0.018845
0.004278 0.01811 30.94775719 41.3399754 51.11071871 0.000485 0.004078 0.017372
0.003897 0.016692 30.94821807 41.3438696 51.12740078 0.000437 0.003716 0.01601
0.003551 0.015382 30.94863301 41.347417 51.14277274 0.000394 0.003385 0.014751
0.003234 0.01417 30.94900658 41.3506483 51.15693374 0.000354 0.003083 0.013587
0.002945 0.01305 30.94934292 41.3535911 51.16997584 0.000319 0.002807 0.012512
0.002682 0.012016 30.94964574 41.3562708 51.18198442 0.000287 0.002556 0.011519
0.002442 0.011061 30.94991837 41.3587107 51.19303867 0.000259 0.002327 0.010602
0.002223 0.01018 30.95016383 41.360932 51.20321198 0.000233 0.002119 0.009756
Demo 8 - Shooting Method for Solving Two-Point Boundary Value Problems
Solved by Explicit Euler's method and goal seek

Definitions Equations/values Final values Step No.

Final value ind.var. zf=0.001 0.001 1


Integration stepsize h 0.00001 2
Constants k = 0.001 0.001 3
DAB = 1.2E-9 1.20E-09 4
Independent variable z 0 0.00001 0.001 5
Dependent variables CA 0.2 0.19868 0.13770 6
y -131.913 -130.247 -7.53E-14 7
Differential equations f1=d(CA)/d(z) = y -131.913 -130.247 -7.53E-14 8
f2=d(y)/d(z) = k*CA/DAB 1.6667E+05 165567.388 1.1475E+05 9
10
Note Solution is obtained by using goal seek to set y(0.001)=0 11
while changing the value of y(0) 12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
z CA y f1=dCA/dz f2=dy/dz

0 0.2 -131.9134 -131.9134057 1.66667E+05


0.00001 0.198681 -130.2467 -130.246739 1.65567E+05
0.00002 0.197378 -128.5911 -128.5910651 1.64482E+05
0.00003 0.196092 -126.9462 -126.9462451 1.63410E+05
0.00004 0.194823 -125.3121 -125.3121411 1.62353E+05
0.00005 0.19357 -123.6886 -123.6886158 1.61308E+05
0.00006 0.192333 -122.0755 -122.0755333 1.60278E+05
0.00007 0.191112 -120.4728 -120.4727582 1.59260E+05
0.00008 0.189908 -118.8802 -118.880156 1.58256E+05
0.00009 0.188719 -117.2976 -117.2975932 1.57266E+05
0.0001 0.187546 -115.7249 -115.7249371 1.56288E+05
0.00011 0.186389 -114.1621 -114.1620558 1.55324E+05
0.00012 0.185247 -112.6088 -112.6088182 1.54372E+05
0.00013 0.184121 -111.0651 -111.0650941 1.53434E+05
0.00014 0.18301 -109.5308 -109.5307541 1.52508E+05
0.00015 0.181915 -108.0057 -108.0056696 1.51596E+05
0.00016 0.180835 -106.4897 -106.4897126 1.50696E+05
0.00017 0.17977 -104.9828 -104.982756 1.49808E+05
0.00018 0.17872 -103.4847 -103.4846736 1.48933E+05
0.00019 0.177685 -101.9953 -101.9953398 1.48071E+05
0.0002 0.176665 -100.5146 -100.5146297 1.47221E+05
0.00021 0.17566 -99.04242 -99.0424192 1.46383E+05
0.00022 0.17467 -97.57858 -97.57858492 1.45558E+05
0.00023 0.173694 -96.123 -96.12300418 1.44745E+05
0.00024 0.172733 -94.67555 -94.67555499 1.43944E+05
0.00025 0.171786 -93.23612 -93.23611605 1.43155E+05
0.00026 0.170854 -91.80457 -91.80456673 1.42378E+05
0.00027 0.169936 -90.38079 -90.3807871 1.41613E+05
0.00028 0.169032 -88.96466 -88.96465784 1.40860E+05
0.00029 0.168142 -87.55606 -87.55606032 1.40118E+05
0.0003 0.167266 -86.15488 -86.15487652 1.39389E+05
0.00031 0.166405 -84.76099 -84.76098906 1.38671E+05
0.00032 0.165557 -83.37428 -83.37428117 1.37964E+05
0.00033 0.164724 -81.99464 -81.99463669 1.37270E+05
0.00034 0.163904 -80.62194 -80.62194008 1.36586E+05
0.00035 0.163097 -79.25608 -79.25607634 1.35915E+05
0.00036 0.162305 -77.89693 -77.89693111 1.35254E+05
0.00037 0.161526 -76.54439 -76.54439055 1.34605E+05
0.00038 0.16076 -75.19834 -75.19834139 1.33967E+05
0.00039 0.160008 -73.85867 -73.85867094 1.33340E+05
0.0004 0.15927 -72.52527 -72.52526702 1.32725E+05
0.00041 0.158545 -71.19802 -71.19801798 1.32121E+05
0.00042 0.157833 -69.87681 -69.87681272 1.31527E+05
0.00043 0.157134 -68.56154 -68.56154062 1.30945E+05
0.00044 0.156448 -67.25209 -67.2520916 1.30374E+05
0.00045 0.155776 -65.94836 -65.94835603 1.29813E+05
0.00046 0.155116 -64.65022 -64.65022481 1.29264E+05
0.00047 0.15447 -63.35759 -63.35758928 1.28725E+05
0.00048 0.153836 -62.07034 -62.07034127 1.28197E+05
0.00049 0.153215 -60.78837 -60.78837306 1.27680E+05
0.0005 0.152608 -59.51158 -59.51157738 1.27173E+05
0.00051 0.152012 -58.23985 -58.23984739 1.26677E+05
0.00052 0.15143 -56.97308 -56.97307671 1.26192E+05
0.00053 0.15086 -55.71116 -55.71115934 1.25717E+05
0.00054 0.150303 -54.45399 -54.45398973 1.25253E+05
0.00055 0.149759 -53.20146 -53.20146272 1.24799E+05
0.00056 0.149227 -51.95347 -51.95347354 1.24356E+05
0.00057 0.148707 -50.70992 -50.70991782 1.23923E+05
0.00058 0.1482 -49.47069 -49.47069155 1.23500E+05
0.00059 0.147705 -48.23569 -48.23569111 1.23088E+05
0.0006 0.147223 -47.00481 -47.00481322 1.22686E+05
0.00061 0.146753 -45.77795 -45.77795498 1.22294E+05
0.00062 0.146295 -44.55501 -44.5550138 1.21913E+05
0.00063 0.14585 -43.33589 -43.33588746 1.21541E+05
0.00064 0.145416 -42.12047 -42.12047403 1.21180E+05
0.00065 0.144995 -40.90867 -40.90867193 1.20829E+05
0.00066 0.144586 -39.70038 -39.70037986 1.20488E+05
0.00067 0.144189 -38.4955 -38.49549685 1.20157E+05
0.00068 0.143804 -37.29392 -37.29392221 1.19837E+05
0.00069 0.143431 -36.09556 -36.09555552 1.19526E+05
0.0007 0.14307 -34.9003 -34.90029666 1.19225E+05
0.00071 0.142721 -33.70805 -33.70804577 1.18934E+05
0.00072 0.142384 -32.5187 -32.51870323 1.18653E+05
0.00073 0.142059 -31.33217 -31.3321697 1.18382E+05
0.00074 0.141746 -30.14835 -30.14834605 1.18121E+05
0.00075 0.141444 -28.96713 -28.96713343 1.17870E+05
0.00076 0.141154 -27.78843 -27.78843316 1.17629E+05
0.00077 0.140876 -26.61215 -26.61214682 1.17397E+05
0.00078 0.14061 -25.43818 -25.43817619 1.17175E+05
0.00079 0.140356 -24.26642 -24.26642323 1.16963E+05
0.0008 0.140113 -23.09679 -23.09679012 1.16761E+05
0.00081 0.139882 -21.92918 -21.92917922 1.16569E+05
0.00082 0.139663 -20.76349 -20.76349304 1.16386E+05
0.00083 0.139455 -19.59963 -19.5996343 1.16213E+05
0.00084 0.139259 -18.43751 -18.43750585 1.16050E+05
0.00085 0.139075 -17.27701 -17.2770107 1.15896E+05
0.00086 0.138902 -16.11805 -16.11805201 1.15752E+05
0.00087 0.138741 -14.96053 -14.96053307 1.15618E+05
0.00088 0.138591 -13.80436 -13.8043573 1.15493E+05
0.00089 0.138453 -12.64943 -12.64942825 1.15378E+05
0.0009 0.138327 -11.49565 -11.49564955 1.15272E+05
0.00091 0.138212 -10.34292 -10.34292498 1.15177E+05
0.00092 0.138109 -9.191158 -9.191158377 1.15090E+05
0.00093 0.138017 -8.040254 -8.040253684 1.15014E+05
0.00094 0.137936 -6.890115 -6.890114921 1.14947E+05
0.00095 0.137867 -5.740646 -5.740646178 1.14889E+05
0.00096 0.13781 -4.591752 -4.591751613 1.14842E+05
0.00097 0.137764 -3.443335 -3.443335434 1.14803E+05
0.00098 0.13773 -2.295302 -2.295301901 1.14775E+05
0.00099 0.137707 -1.147555 -1.147555313 1.14756E+05
0.001 0.137695 -7.53E-14 -7.52731E-14 1.14746E+05
Demo 9 - Reversible , Exothermic Gas Phase Reaction in a Catalytic Reactor

Equations Specified values\equations


Final value (Ind. Var.) WF 20
Integration Stepsize h 0.2
Constants Ta = 500 500
delH = -40000 -40000
CPA = 40 40
FA0 = 5 5
Independent Variable W0 0 0.2
Dependent Variables x 0 0.00146882
T 450 451.50882
y 1 0.9985
Functions k = .5*exp(5032*(1/450-1/T)) 0.5 0.5190374569
CA = .271*(1-x)*(450/T)/(1-.5*x)*y 0.271 0.2694910419
CC = .271*.5*x*(450/T)/(1-.5*x)*y 0 0.000198208
Kc = 25000*exp(delH/8.314*(1/450-1/T)) 25000 24122.567892
rA = -k*(CA^2-CC/Kc) -0.0367205 -0.0376953099
Differential equations f1=d(x)/d(W) = -rA/FA0 0.0073441 0.007539062
f2=d(T)/d(W) = (.8*(Ta-T)+rA*delH)/(CPA*FA0) 7.5441 7.7330266987
f3=d(y)/d(W) = -0.015*(1-.5*x)*(T/450)/(2*y) -0.0075 -0.0075309168

Conversion, Reduced Pressure and Temperature Profiles

1.2

0.8
x
0.6 T*10-3
y(W)
0.4

0.2

-1.66533453693773E-16

-0.2
Catalysts Weight (W)

Concentration Profiles in Catalytic Reactor


0.3

0.25

0.2
Concentration

0.15 CA
CC
0.1

0.05

0
0 2 4 6 8 10 12 14 16 18 20
0.25

0.2

Concentration
0.15 CA
CC
0.1

0.05

0
0 2 4 6 8 10 12 14 16 18 20
Catalyst's Weight (W)
Step no. (i) W x T(W) y(W) k
0 0 0 450 1 0.5
1 0.2 0.00146882 451.50882 0.9985 0.519037457
2 0.4 0.0029766324 453.05542534 0.996993817 0.539164914
3 0.6 0.0045252056 454.64155424 0.995481333 0.560466976
4 0.8 0.0061164229 456.26905827 0.993962426 0.583036566
5 1 0.0077522918 457.93991189 0.992436971 0.606975922
6 1.2 0.0094349548 459.65622301 0.990904833 0.632397725
7 1.4 0.0111667015 461.42024466 0.989365877 0.6594264
8 1.6 0.0129499809 463.23438788 0.987819958 0.688199599
9 1.8 0.0147874165 465.10123603 0.986266925 0.718869905
10 2 0.0166818222 467.02356072 0.984706623 0.751606799
11 2.2 0.0186362198 469.0043395 0.983138886 0.786598935
12 2.4 0.0206538594 471.04677556 0.981563544 0.824056765
13 2.6 0.0227382411 473.15431983 0.979980416 0.864215597
14 2.8 0.0248931403 475.33069566 0.978389312 0.90733916
15 3 0.0271226357 477.57992646 0.976790036 0.953723763
16 3.2 0.0294311401 479.90636689 0.975182378 1.003703176
17 3.4 0.0318234362 482.31473789 0.973566119 1.057654377
18 3.6 0.0343047163 484.81016628 0.971941027 1.116004329
19 3.8 0.0368806278 487.39822957 0.970306859 1.179238015
20 4 0.0395573236 490.08500679 0.968663357 1.247907994
21 4.2 0.0423415211 492.87713629 0.967010248 1.322645806
22 4.4 0.0452405681 495.78188159 0.965347244 1.40417565
23 4.6 0.0482625187 498.80720673 0.963674039 1.493330839
24 4.8 0.0514162203 501.9618625 0.961990308 1.591073704
25 5 0.0547114129 505.2554856 0.960295705 1.698519765
26 5.2 0.0581588445 508.69871284 0.958589863 1.816967216
27 5.4 0.0617704036 512.30331303 0.956872389 1.947933091
28 5.6 0.0655592731 516.08233983 0.955142863 2.093197828
29 5.8 0.0695401085 520.05030932 0.953400836 2.254860489
30 6 0.0737292462 524.22340682 0.951645827 2.435407547
31 6.2 0.0781449466 528.61972853 0.949877318 2.63779907
32 6.4 0.0828076786 533.25956467 0.948094754 2.865577348
33 6.6 0.0877404535 538.16573193 0.946297533 3.123004633
34 6.8 0.0929692192 543.36396503 0.944485009 3.415238909
35 7 0.0985233241 548.88337876 0.942656478 3.748559654
36 7.2 0.1044360663 554.75701433 0.940811181 4.130659753
37 7.4 0.1107453439 561.02248631 0.93894829 4.571025536
38 7.6 0.1174944249 567.72274927 0.937066904 5.081434949
39 7.8 0.12473286 574.90700622 0.935166041 5.676615152
40 8 0.1325175628 582.63178341 0.933244628 6.375116545
41 8.2 0.1409140828 590.96219795 0.931301489 7.200482227
42 8.4 0.1499980985 599.97344389 0.929335334 8.182822689
43 8.6 0.1598571508 609.75251744 0.92734475 9.360948191
44 8.8 0.1705926235 620.40018811 0.925328183 10.7852699
45 9 0.1823219491 632.03319362 0.923283927 12.52175938
46 9.2 0.1951809631 644.78658098 0.921210111 14.65735717
47 9.4 0.209326227 658.81601567 0.919104686 17.30734174
48 9.6 0.2249369721 674.29970799 0.916965422 20.62529119
49 9.8 0.2422160214 691.43931744 0.914789905 24.81633728
50 10 0.2613885912 710.45873581 0.912575551 30.15429325
51 10.2 0.2826971746 731.59895224 0.910319644 37.00265771
52 10.4 0.3063897256 755.10622407 0.908019395 45.83793209
53 10.6 0.3326971613 781.20957478 0.905672059 57.27027577
54 10.8 0.3617951033 810.08254915 0.903275105 72.0501507
55 11 0.3937447174 841.78409716 0.900826458 91.03965751
56 11.2 0.4284103571 876.17630963 0.898324829 115.1158147
57 11.4 0.4653603552 912.82536668 0.895770093 144.9690229
58 11.6 0.503774271 950.90902222 0.893163661 180.7828264
59 11.8 0.5424031767 989.17720068 0.890508727 221.8538548
60 12 0.579645632 1026.0283142 0.887810229 266.3275534
61 12.2 0.6137835042 1059.7453637 0.885074427 311.3023348
62 12.4 0.6433460751 1088.8601383 0.882308113 353.4452282
63 12.6 0.6674644841 1112.5074592 0.879517692 389.9269062
64 12.8 0.6860363429 1130.5893121 0.876708471 419.1793385
65 13 0.6996189505 1143.6674483 0.87388436 441.0660479
66 13.2 0.7091406031 1152.6741669 0.871047975 456.4933989
67 13.4 0.7156021861 1158.6136105 0.868200944 466.8244455
68 13.6 0.7198842304 1162.368764 0.865344231 473.4205828
69 13.8 0.7226731504 1164.627789 0.862478384 477.4126921
70 14 0.7244648449 1165.8877812 0.859603702 479.6471537
71 14.2 0.7256007916 1166.4910177 0.856720337 480.7189117
72 14.4 0.7263092461 1166.6662794 0.853828348 481.0305361
73 14.6 0.726740084 1166.5637843 0.850927742 480.8482813
74 14.8 0.7269907227 1166.2811719 0.848018488 480.3459374
75 15 0.7271241536 1165.8815779 0.845100535 479.6361391
76 15.2 0.7271810248 1165.4057439 0.842173817 478.7916504
77 15.4 0.727187565 1164.8799594 0.839238258 477.859441
78 15.6 0.7271607185 1164.321209 0.836293775 476.8698525
79 15.8 0.727111453 1163.7404865 0.83334028 475.8425179
80 16 0.7270468841 1163.1449252 0.830377681 474.790169
81 16.2 0.7269716391 1162.5391643 0.827405883 473.7210828
82 16.4 0.7268887347 1161.9262286 0.824424788 472.6406534
83 16.6 0.726800141 1161.3080938 0.821434296 471.5524041
84 16.8 0.7267071462 1160.6860526 0.818434305 470.4586407
85 17 0.7266105931 1160.0609506 0.81542471 469.3608736
86 17.2 0.7265110318 1159.4333406 0.812405406 468.2600923
87 17.4 0.7264088199 1158.803582 0.809376283 467.1569432
88 17.6 0.7263041875 1158.1719068 0.806337231 466.0518467
89 17.8 0.7261972806 1157.5384623 0.803288137 464.9450731
90 18 0.7260881884 1156.9033393 0.800228887 463.8367932
91 18.2 0.7259769626 1156.2665908 0.797159363 462.7271108
92 18.4 0.7258636293 1155.6282442 0.794079447 461.6160852
93 18.6 0.7257481974 1154.9883098 0.790989018 460.5037452
94 18.8 0.7256306641 1154.3467859 0.787887951 459.3900993
95 19 0.7255110183 1153.7036627 0.78477612 458.2751418
96 19.2 0.7253892432 1153.0589246 0.781653398 457.1588571
97 19.4 0.725265318 1152.4125522 0.778519654 456.0412229
98 19.6 0.7251392187 1151.7645229 0.775374753 454.922212
99 19.8 0.7250109194 1151.114812 0.77221856 453.8017935
100 20 0.7248803921 1150.4633929 0.769050936 452.6799339
CA CC Kc rA f1 f2 f3 T*10-3
0.271 0 25000 -0.036721 0.007344 7.5441 -0.0075 0.45
0.269491 0.000198 24122.57 -0.037695 0.007539 7.733027 -0.007531 0.451509
0.267963 0.0004 23260.86 -0.038714 0.007743 7.930645 -0.007562 0.453055
0.266416 0.000606 22414.85 -0.03978 0.007956 8.13752 -0.007595 0.454642
0.264848 0.000815 21584.53 -0.040897 0.008179 8.354268 -0.007627 0.456269
0.263259 0.001028 20769.87 -0.042067 0.008413 8.581556 -0.007661 0.45794
0.261648 0.001246 19970.87 -0.043294 0.008659 8.820108 -0.007695 0.459656
0.260014 0.001468 19187.51 -0.044582 0.008916 9.070716 -0.00773 0.46142
0.258356 0.001695 18419.79 -0.045936 0.009187 9.334241 -0.007765 0.463234
0.256674 0.001926 17667.69 -0.04736 0.009472 9.611623 -0.007802 0.465101
0.254966 0.002163 16931.22 -0.04886 0.009772 9.903894 -0.007839 0.467024
0.25323 0.002404 16210.37 -0.050441 0.010088 10.21218 -0.007877 0.469004
0.251467 0.002652 15505.14 -0.05211 0.010422 10.53772 -0.007916 0.471047
0.249674 0.002905 14815.55 -0.053872 0.010774 10.88188 -0.007956 0.473154
0.24785 0.003164 14141.59 -0.055737 0.011147 11.24615 -0.007996 0.475331
0.245994 0.003429 13483.28 -0.057713 0.011543 11.6322 -0.008038 0.47758
0.244105 0.003701 12840.63 -0.059807 0.011961 12.04185 -0.008081 0.479906
0.242179 0.00398 12213.66 -0.062032 0.012406 12.47714 -0.008125 0.482315
0.240217 0.004267 11602.38 -0.064398 0.01288 12.94032 -0.008171 0.48481
0.238216 0.004561 11006.83 -0.066917 0.013383 13.43389 -0.008218 0.487398
0.236173 0.004864 10427.01 -0.069605 0.013921 13.96065 -0.008266 0.490085
0.234087 0.005175 9862.964 -0.072476 0.014495 14.52373 -0.008315 0.492877
0.231956 0.005496 9314.717 -0.075549 0.01511 15.12663 -0.008366 0.495782
0.229776 0.005826 8782.302 -0.078843 0.015769 15.77328 -0.008419 0.498807
0.227546 0.006167 8265.754 -0.08238 0.016476 16.46812 -0.008473 0.501962
0.225261 0.006519 7765.111 -0.086186 0.017237 17.21614 -0.008529 0.505255
0.222919 0.006883 7280.414 -0.090289 0.018058 18.023 -0.008587 0.508699
0.220517 0.007259 6811.702 -0.094722 0.018944 18.89513 -0.008648 0.512303
0.218051 0.007649 6359.021 -0.099521 0.019904 19.83985 -0.00871 0.516082
0.215516 0.008054 5922.415 -0.104728 0.020946 20.86549 -0.008775 0.52005
0.212908 0.008474 5501.929 -0.110393 0.022079 21.98161 -0.008843 0.524223
0.210222 0.00891 5097.608 -0.116568 0.023314 23.19918 -0.008913 0.52862
0.207453 0.009365 4709.498 -0.123319 0.024664 24.53084 -0.008986 0.53326
0.204595 0.009839 4337.643 -0.130719 0.026144 25.99117 -0.009063 0.538166
0.201642 0.010334 3982.083 -0.138853 0.027771 27.59707 -0.009143 0.543364
0.198586 0.010852 3642.855 -0.147819 0.029564 29.36818 -0.009226 0.548883
0.19542 0.011394 3319.991 -0.157732 0.031546 31.32736 -0.009314 0.554757
0.192136 0.011964 3013.516 -0.168727 0.033745 33.50131 -0.009407 0.561022
0.188724 0.012563 2723.443 -0.180961 0.036192 35.92128 -0.009504 0.567723
0.185174 0.013194 2449.776 -0.194618 0.038924 38.62389 -0.009607 0.574907
0.181475 0.013861 2192.499 -0.209913 0.041983 41.65207 -0.009716 0.582632
0.177615 0.014567 1951.58 -0.2271 0.04542 45.05623 -0.009831 0.590962
0.17358 0.015316 1726.961 -0.246476 0.049295 48.89537 -0.009953 0.599973
0.169356 0.016112 1518.551 -0.268387 0.053677 53.23835 -0.010083 0.609753
0.164928 0.016961 1326.227 -0.293233 0.058647 58.16503 -0.010221 0.6204
0.160278 0.017869 1149.817 -0.321475 0.064295 63.76694 -0.010369 0.632033
0.155389 0.018842 989.0994 -0.353632 0.070726 70.14717 -0.010527 0.644787
0.150243 0.019888 843.7869 -0.390269 0.078054 77.41846 -0.010696 0.658816
0.144822 0.021015 713.5196 -0.431976 0.086395 85.69805 -0.010878 0.6743
0.13911 0.022232 597.8525 -0.479314 0.095863 95.09709 -0.011072 0.691439
0.133093 0.02355 496.2447 -0.532715 0.106543 105.7011 -0.01128 0.710459
0.126762 0.024979 408.0496 -0.592314 0.118463 117.5364 -0.011501 0.731599
0.120116 0.02653 332.5079 -0.657686 0.131537 130.5168 -0.011737 0.755106
0.113168 0.028211 268.7458 -0.727449 0.14549 144.3649 -0.011985 0.78121
0.105948 0.030031 215.7804 -0.79874 0.159748 158.5077 -0.012243 0.810083
0.098513 0.031991 172.5341 -0.866641 0.173328 171.9611 -0.012508 0.841784
0.090949 0.034084 137.8614 -0.92375 0.18475 183.2453 -0.012774 0.876176
0.083383 0.036289 110.5853 -0.960348 0.19207 190.4183 -0.013032 0.912825
0.075978 0.038567 89.54129 -0.965723 0.193145 191.3409 -0.013275 0.950909
0.068932 0.040854 73.62334 -0.931061 0.186212 184.2556 -0.013492 0.989177
0.062458 0.043063 61.82281 -0.853447 0.170689 168.5852 -0.013679 1.026028
0.056753 0.045097 53.25453 -0.739064 0.147813 145.5739 -0.013832 1.059745
0.051956 0.04686 47.16683 -0.60296 0.120592 118.2366 -0.013952 1.08886
0.048119 0.048292 42.9386 -0.464296 0.092859 90.40926 -0.014046 1.112507
0.045192 0.049374 40.06914 -0.339565 0.067913 65.39068 -0.014121 1.130589
0.043049 0.050133 38.16597 -0.238041 0.047608 45.03359 -0.014182 1.143667
0.041529 0.050626 36.93182 -0.16154 0.032308 29.69722 -0.014235 1.152674
0.040469 0.050914 36.14998 -0.107051 0.02141 18.77577 -0.014284 1.158614
0.039732 0.051055 35.66827 -0.069723 0.013945 11.29513 -0.014329 1.162369
0.039216 0.051096 35.38305 -0.044792 0.008958 6.299961 -0.014373 1.164628
0.038845 0.051068 35.22543 -0.028399 0.00568 3.016183 -0.014417 1.165888
0.03857 0.050995 35.15034 -0.017711 0.003542 0.876308 -0.01446 1.166491
0.038356 0.050894 35.12857 -0.010771 0.002154 -0.512475 -0.014503 1.166666
0.038182 0.050772 35.1413 -0.006266 0.001253 -1.413062 -0.014546 1.166564
0.038033 0.050639 35.17643 -0.003336 0.000667 -1.99797 -0.01459 1.166281
0.0379 0.050496 35.2262 -0.001422 0.000284 -2.37917 -0.014634 1.165882
0.037778 0.050348 35.28561 -0.000164 3.27E-05 -2.628922 -0.014678 1.165406
0.037663 0.050196 35.35142 0.000671 -0.000134 -2.793752 -0.014722 1.16488
0.037552 0.050041 35.42156 0.001232 -0.000246 -2.903612 -0.014767 1.164321
0.037443 0.049884 35.49467 0.001614 -0.000323 -2.977807 -0.014813 1.16374
0.037336 0.049725 35.56989 0.001881 -0.000376 -3.028804 -0.014859 1.163145
0.03723 0.049565 35.64663 0.002073 -0.000415 -3.064679 -0.014905 1.162539
0.037124 0.049403 35.72454 0.002215 -0.000443 -3.090674 -0.014952 1.161926
0.037019 0.049241 35.80336 0.002325 -0.000465 -3.110206 -0.015 1.161308
0.036913 0.049077 35.88294 0.002414 -0.000483 -3.12551 -0.015048 1.160686
0.036807 0.048913 35.96318 0.002489 -0.000498 -3.13805 -0.015097 1.160061
0.036701 0.048748 36.04401 0.002555 -0.000511 -3.148793 -0.015146 1.159433
0.036595 0.048582 36.12538 0.002616 -0.000523 -3.158376 -0.015195 1.158804
0.036489 0.048415 36.20728 0.002673 -0.000535 -3.167222 -0.015245 1.158172
0.036382 0.048247 36.28968 0.002727 -0.000545 -3.175615 -0.015296 1.157538
0.036274 0.048078 36.37258 0.002781 -0.000556 -3.183743 -0.015348 1.156903
0.036167 0.047909 36.45598 0.002833 -0.000567 -3.191733 -0.0154 1.156267
0.036059 0.047738 36.53986 0.002886 -0.000577 -3.199672 -0.015452 1.155628
0.03595 0.047567 36.62425 0.002938 -0.000588 -3.20762 -0.015505 1.154988
0.035841 0.047395 36.70913 0.002991 -0.000598 -3.215616 -0.015559 1.154347
0.035732 0.047222 36.79452 0.003044 -0.000609 -3.22369 -0.015614 1.153704
0.035622 0.047048 36.88041 0.003098 -0.00062 -3.231862 -0.015669 1.153059
0.035511 0.046873 36.96683 0.003152 -0.00063 -3.240147 -0.015725 1.152413
0.035401 0.046697 37.05376 0.003207 -0.000641 -3.248555 -0.015781 1.151765
0.035289 0.04652 37.14123 0.003263 -0.000653 -3.257096 -0.015838 1.151115
0.035177 0.046343 37.22923 0.00332 -0.000664 -3.265775 -0.015896 1.150463
Demo 10 - Dynamics of a heated tank

Equations Specified values\equations


Final value (Ind. Var.) tf 60
Integration Stepsize h 0.5
Constants WC = 500 500
rhoVCp = 4000 4000
Tr = 80 80
Kc = 0 0
Kr = 250 250
tauI = 2 2
taum = 5 5
taud = 1 1
Independent Variable t 0 0.5
Dependent Variables T 80 80
Tm 80 80
T0 80 80
errsum 0 0
Functions step = if (t<10) then (0) else 0 0
Ti = 60-20*step 60 60
q = 10000+Kc*(Tr-Tm)+Kc/tauI*errsum 10000 10000
Differential equations f1=d(T)/d(t) = (WC*(Ti-T)+q)/rhoVCp 0 0
f2= d(Tm)/d(t) = (T0-Tm)/taum 0 0
f3=d(T0)/d(t) = (T-T0-(taud/2)*dTdt)*2/taud 0 0
f4 = d(errsum)/d(t) = Tr-Tm 0 0

Open Loop Response to Step Down in Inlet

80
Temperature (deg C)

75

70 T(t)
Tm(t)
T0(t)
65

60
0 10 20 30 40 50 60
time (t)
Step no. (i) t T(t) Tm(t) T0(t) errsum(t)
0 0 80 80 80 0
1 0.5 80 80 80 0
2 1 80 80 80 0
3 1.5 80 80 80 0
4 2 80 80 80 0
5 2.5 80 80 80 0
6 3 80 80 80 0
7 3.5 80 80 80 0
8 4 80 80 80 0
9 4.5 80 80 80 0
10 5 80 80 80 0
11 5.5 80 80 80 0
12 6 80 80 80 0
13 6.5 80 80 80 0
14 7 80 80 80 0
15 7.5 80 80 80 0
16 8 80 80 80 0
17 8.5 80 80 80 0
18 9 80 80 80 0
19 9.5 80 80 80 0
20 10 80 80 80 0
21 10.5 80 80 80 0
22 11 78.75 80 81.25 0
23 11.5 77.578125 80.125 79.921875 0
24 12 76.4794922 80.1046875 78.676757813 -0.0625
25 12.5 75.4495239 79.961894531 77.509460449 -0.11484375
26 13 74.4839287 79.716651123 76.415119171 -0.095791016
27 13.5 73.5786831 79.386497928 75.389174223 0.045883423
28 14 72.7300154 78.986765557 74.427350834 0.352634459
29 14.5 71.9343895 78.530824085 73.525641407 0.85925168
30 15 71.1884901 78.030305817 72.680288819 1.593839638
31 15.5 70.4892095 77.495304117 71.887770768 2.578686729
32 16 69.8336339 76.934550782 71.144785095 3.83103467
33 16.5 69.2190318 76.355574214 70.448236026 5.363759279
34 17 68.6428423 75.764840395 69.795221275 7.185972172
35 17.5 68.1026647 75.167878483 69.183019945 9.303551975
36 18 67.5962481 74.569392629 68.609081198 11.71961273
37 18.5 67.1214826 73.973361486 68.071013624 14.43491642
38 19 66.6763899 73.3831267 67.566575272 17.44823568
39 19.5 66.2591156 72.801471557 67.093664318 20.75667233
40 20 65.8679209 72.230690833 66.650310298 24.35593655
41 20.5 65.5011758 71.67265278 66.234665904 28.24059113
42 21 65.1573523 71.128854092 65.844999285 32.40426474
43 21.5 64.8350178 70.600468611 65.47968683 36.8398377
44 22 64.5328292 70.088390433 65.137206403 41.53960339
45 22.5 64.2495274 69.59327203 64.816131003 46.49540817
46 23 63.9839319 69.115557927 64.515122815 51.69877216
47 23.5 63.7349362 68.655514416 64.232927639 57.14099319
48 24 63.5015026 68.213255738 63.968369662 62.81323599
49 24.5 63.2826587 67.788767131 63.720346558 68.70660812
50 25 63.0774926 67.381925074 63.487824898 74.81222455
51 25.5 62.8851493 66.992515056 63.269835842 81.12126201
52 26 62.7048274 66.620247135 63.065471102 87.62500449
53 26.5 62.5357757 66.264769531 62.873879158 94.31488092
54 27 62.3772897 65.925680494 62.694261711 101.1824962
55 27.5 62.2287091 65.602538616 62.525870354 108.2196559
56 28 62.0894148 65.294871789 62.368003457 115.4183866
57 28.5 61.9588264 65.002184956 62.22000324 122.7709507
58 29 61.8363997 64.723966785 62.081253038 130.2698582
59 29.5 61.7216248 64.45969541 61.951174723 137.9078748
60 30 61.6140232 64.208843341 61.829226303 145.6780271
61 30.5 61.5131468 63.970881637 61.714899659 153.5736055
62 31 61.4185751 63.74528344 61.60771843 161.5881646
63 31.5 61.3299141 63.531526939 61.507236028 169.7155229
64 32 61.2467945 63.329097848 61.413033777 177.9497595
65 32.5 61.1688699 63.13749144 61.324719166 186.2852105
66 33 61.0958155 62.956214213 61.241924218 194.7164648
67 33.5 61.027327 62.784785213 61.164303954 203.2383577
68 34 60.9631191 62.622737088 61.091534957 211.8459651
69 34.5 60.9029241 62.469616874 61.023314022 220.5345965
70 35 60.8464914 62.324986589 60.959356896 229.2997881
71 35.5 60.7935857 62.18842362 60.89939709 238.1372948
72 36 60.7439866 62.059520967 60.843184772 247.043083
73 36.5 60.6974874 61.937887347 60.790485723 256.0133225
74 37 60.6538944 61.823147185 60.741080366 265.0443789
75 37.5 60.613026 61.714940503 60.694762843 274.1328053
76 38 60.5747119 61.612922737 60.651340165 283.275335
77 38.5 60.5387924 61.51676448 60.610631405 292.4688736
78 39 60.5051179 61.426151172 60.572466942 301.7104914
79 39.5 60.473548 61.340782749 60.536687758 310.9974158
80 40 60.4439513 61.26037325 60.503144773 320.3270244
81 40.5 60.4162043 61.184650403 60.471698225 329.6968378
82 41 60.3901915 61.113355185 60.442217086 339.1045126
83 41.5 60.3658046 61.046241375 60.414578518 348.547835
84 42 60.3429418 60.983075089 60.388667361 358.0247143
85 42.5 60.3215079 60.923634316 60.364375651 367.5331768
86 43 60.3014137 60.86770845 60.341602172 377.0713596
87 43.5 60.2825753 60.815097822 60.320252037 386.6375054
88 44 60.2649144 60.765613244 60.300236284 396.2299565
89 44.5 60.2483572 60.719075548 60.281471517 405.8471499
90 45 60.2328349 60.675315144 60.263879547 415.4876121
91 45.5 60.2182827 60.634171585 60.247387075 425.1499545
92 46 60.20464 60.595493134 60.231925383 434.8328687
93 46.5 60.19185 60.559136359 60.217430047 444.5351222
94 47 60.1798594 60.524965727 60.203840669 454.255554
95 47.5 60.1686182 60.492853222 60.191100627 463.9930711
96 48 60.1580796 60.462677962 60.179156838 473.7466445
97 48.5 60.1481996 60.43432585 60.167959535 483.5153055
98 49 60.1389371 60.407689218 60.157462064 493.2981426
99 49.5 60.1302535 60.382666503 60.147620685 503.094298
100 50 60.1221127 60.359161921 60.138394392 512.9029647
101 50.5 60.1144807 60.337085168 60.129744743 522.7233838
102 51 60.1073256 60.316351126 60.121635697 532.5548412
103 51.5 60.1006178 60.296879583 60.114033465 542.3966656
104 52 60.0943292 60.278594971 60.106906374 552.2482258
105 52.5 60.0884336 60.261426111 60.100224726 562.1089284
106 53 60.0829065 60.245305973 60.09396068 571.9782153
107 53.5 60.0777248 60.230171443 60.088088138 581.8555623
108 54 60.072867 60.215963113 60.082582629 591.7404766
109 54.5 60.0683128 60.202625065 60.077421215 601.632495
110 55 60.0640433 60.19010468 60.072582389 611.5311825
111 55.5 60.0600406 60.17835245 60.06804599 621.4361302
112 56 60.056288 60.167321804 60.063793115 631.3469539
113 56.5 60.05277 60.156968935 60.059806045 641.263293
114 57 60.0494719 60.147252646 60.056068168 651.1848086
115 57.5 60.0463799 60.138134199 60.052563907 661.1111823
116 58 60.0434812 60.129577169 60.049278663 671.0421152
117 58.5 60.0407636 60.121547319 60.046198747 680.9773266
118 59 60.0382159 60.114012462 60.043311325 690.9165529
119 59.5 60.0358274 60.106942348 60.040604367 700.8595467
120 60 60.0335882 60.10030855 60.038066594 710.8060755
step Ti q f1 f2 f3 f4
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
0 60 10000 0 0 0 0
1 40 10000 -2.5 0 2.5 0
1 40 10000 -2.34375 0.25 -2.65625 0
1 40 10000 -2.197266 -0.040625 -2.490234 -0.125
1 40 10000 -2.059937 -0.285586 -2.334595 -0.104687
1 40 10000 -1.93119 -0.490487 -2.188683 0.038105
1 40 10000 -1.810491 -0.660306 -2.05189 0.283349
1 40 10000 -1.697335 -0.799465 -1.923647 0.613502
1 40 10000 -1.591252 -0.911883 -1.803419 1.013234
1 40 10000 -1.491799 -1.001037 -1.690705 1.469176
1 40 10000 -1.398561 -1.070003 -1.585036 1.969694
1 40 10000 -1.311151 -1.121507 -1.485971 2.504696
1 40 10000 -1.229204 -1.157953 -1.393098 3.065449
1 40 10000 -1.152379 -1.181468 -1.30603 3.644426
1 40 10000 -1.080355 -1.193924 -1.224403 4.23516
1 40 10000 -1.012833 -1.196972 -1.147877 4.832122
1 40 10000 -0.949531 -1.192062 -1.076135 5.430607
1 40 10000 -0.890185 -1.18047 -1.008877 6.026639
1 40 10000 -0.834549 -1.16331 -0.945822 6.616873
1 40 10000 -0.782389 -1.141561 -0.886708 7.198528
1 40 10000 -0.73349 -1.116076 -0.831289 7.769309
1 40 10000 -0.687647 -1.087597 -0.779333 8.327347
1 40 10000 -0.644669 -1.056771 -0.730625 8.871146
1 40 10000 -0.604377 -1.024156 -0.684961 9.399531
1 40 10000 -0.566604 -0.990237 -0.642151 9.91161
1 40 10000 -0.531191 -0.955428 -0.602016 10.40673
1 40 10000 -0.497991 -0.920087 -0.56439 10.88444
1 40 10000 -0.466867 -0.884517 -0.529116 11.34449
1 40 10000 -0.437688 -0.848977 -0.496046 11.78674
1 40 10000 -0.410332 -0.813684 -0.465043 12.21123
1 40 10000 -0.384687 -0.77882 -0.435978 12.61807
1 40 10000 -0.360644 -0.744536 -0.408729 13.00748
1 40 10000 -0.338103 -0.710955 -0.383184 13.37975
1 40 10000 -0.316972 -0.678178 -0.359235 13.73523
1 40 10000 -0.297161 -0.646284 -0.336783 14.07432
1 40 10000 -0.278589 -0.615334 -0.315734 14.39746
1 40 10000 -0.261177 -0.585374 -0.296 14.70513
1 40 10000 -0.244853 -0.556436 -0.2775 14.99782
1 40 10000 -0.22955 -0.528543 -0.260157 15.27603
1 40 10000 -0.215203 -0.501704 -0.243897 15.5403
1 40 10000 -0.201753 -0.475923 -0.228653 15.79116
1 40 10000 -0.189143 -0.451196 -0.214362 16.02912
1 40 10000 -0.177322 -0.427513 -0.200965 16.25472
1 40 10000 -0.166239 -0.404858 -0.188405 16.46847
1 40 10000 -0.155849 -0.383213 -0.176629 16.6709
1 40 10000 -0.146109 -0.362554 -0.16559 16.86251
1 40 10000 -0.136977 -0.342858 -0.155241 17.04379
1 40 10000 -0.128416 -0.324096 -0.145538 17.21521
1 40 10000 -0.12039 -0.30624 -0.136442 17.37726
1 40 10000 -0.112866 -0.289261 -0.127914 17.53038
1 40 10000 -0.105811 -0.273126 -0.11992 17.67501
1 40 10000 -0.099198 -0.257805 -0.112425 17.81158
1 40 10000 -0.092998 -0.243267 -0.105398 17.94048
1 40 10000 -0.087186 -0.22948 -0.098811 18.06211
1 40 10000 -0.081737 -0.216413 -0.092635 18.17685
1 40 10000 -0.076628 -0.204036 -0.086845 18.28506
1 40 10000 -0.071839 -0.192317 -0.081418 18.38708
1 40 10000 -0.067349 -0.181227 -0.076329 18.48324
1 40 10000 -0.06314 -0.170737 -0.071558 18.57385
1 40 10000 -0.059194 -0.160819 -0.067086 18.65922
1 40 10000 -0.055494 -0.151446 -0.062893 18.73963
1 40 10000 -0.052026 -0.14259 -0.058962 18.81535
1 40 10000 -0.048774 -0.134228 -0.055277 18.88664
1 40 10000 -0.045726 -0.126333 -0.051822 18.95376
1 40 10000 -0.042868 -0.118882 -0.048583 19.01692
1 40 10000 -0.040188 -0.111852 -0.045547 19.07637
1 40 10000 -0.037677 -0.105221 -0.0427 19.13229
1 40 10000 -0.035322 -0.098969 -0.040032 19.1849
1 40 10000 -0.033114 -0.093075 -0.03753 19.23439
1 40 10000 -0.031045 -0.087521 -0.035184 19.28092
1 40 10000 -0.029104 -0.082287 -0.032985 19.32468
1 40 10000 -0.027285 -0.077357 -0.030923 19.36583
1 40 10000 -0.02558 -0.072714 -0.028991 19.40451
1 40 10000 -0.023981 -0.068341 -0.027179 19.44086
1 40 10000 -0.022482 -0.064225 -0.02548 19.47503
1 40 10000 -0.021077 -0.060351 -0.023888 19.50715
1 40 10000 -0.01976 -0.056704 -0.022395 19.53732
1 40 10000 -0.018525 -0.053273 -0.020995 19.56567
1 40 10000 -0.017367 -0.050045 -0.019683 19.59231
1 40 10000 -0.016282 -0.047009 -0.018453 19.61733
1 40 10000 -0.015264 -0.044154 -0.017299 19.64084
1 40 10000 -0.01431 -0.041468 -0.016218 19.66291
1 40 10000 -0.013416 -0.038943 -0.015204 19.68365
1 40 10000 -0.012577 -0.036569 -0.014254 19.70312
1 40 10000 -0.011791 -0.034338 -0.013363 19.72141
1 40 10000 -0.011054 -0.03224 -0.012528 19.73857
1 40 10000 -0.010363 -0.030269 -0.011745 19.75469
1 40 10000 -0.009716 -0.028417 -0.011011 19.76983
1 40 10000 -0.009108 -0.026676 -0.010323 19.78404
1 40 10000 -0.008539 -0.025041 -0.009678 19.79737
1 40 10000 -0.008005 -0.023504 -0.009073 19.8099
1 40 10000 -0.007505 -0.022061 -0.008506 19.82165
1 40 10000 -0.007036 -0.020706 -0.007974 19.83268
1 40 10000 -0.006596 -0.019433 -0.007476 19.84303
1 40 10000 -0.006184 -0.018237 -0.007009 19.85275
1 40 10000 -0.005797 -0.017114 -0.00657 19.86187
1 40 10000 -0.005435 -0.01606 -0.00616 19.87042
1 40 10000 -0.005095 -0.01507 -0.005775 19.87845
1 40 10000 -0.004777 -0.01414 -0.005414 19.88599
1 40 10000 -0.004478 -0.013268 -0.005076 19.89306
1 40 10000 -0.004199 -0.012448 -0.004758 19.89969

0
Demo 11- Binary Batch Distillation
Solved by Implicit Euler's Method

Definitions Equations/values

Final value ind.var. x2(f)= 0.8


Integration stepsize h 0.02
Independent variable x2 0.4
Dependent variables L= 100
T= 95.583
Explicit equations x1 = 1-x2 0.6
k1 = 10^(6.90565-1211.033/(T+220.79))/(760*1.2) 1.3116
k2 = 10^(6.95464-1344.8/(T+219.482))/(760*1.2) 0.53250
Differential equations f1=d(L)/d(x2) = L/(k2*x2-x2) -534.754
f2=f(T)=(1-k1*x1-k2*x2)=0 6.8122E-05
Sum of squares of errors [Ln-(Ln-1+h/2(f1n+f1n-1))]^2+f2n+1^2 4.6406E-09
Step no. (i) x2 L T x1
0 0.4 100 95.5851 0.6
1 0.42 89.97911 96.12759 0.58
uations/values 2 0.44 81.19123 96.70554 0.56
3 0.46 73.44217 97.25485 0.54
4 0.48 66.57683 97.81125 0.52
5 0.5 60.45657 98.44261 0.5
0.42 6 0.52 54.97335 99.03052 0.48
89.976 7 0.54 50.04422 99.63007 0.46
96.142 8 0.56 45.59569 100.2537 0.44
0.58 9 0.58 41.56225 100.8966 0.42
1.3321 10 0.6 37.89499 101.5128 0.4
0.54185 11 0.62 34.55202 102.1519 0.38
-467.601 12 0.64 31.48953 102.8349 0.36
-2.1100E-04 13 0.66 28.67771 103.4824 0.34
4.4521E-08 14 0.68 26.08691 104.1937 0.32
15 0.7 23.69168 104.8927 0.3
16 0.72 21.47178 105.605 0.28
17 0.74 19.40994 106.312 0.26
18 0.76 17.48838 107.0412 0.24
19 0.78 15.6934 107.7916 0.22
20 0.8 14.00585 108.5954 0.2

1
For x1=0.4 find the solution for f1=0 by changing T using GoalSeek
For other values of x1 use the solver sequentially, row after row to
minimize sum (sum=0) by changing L and T.
k1 k2 f1 f2 sum1
1.311644 0.532535 -3.65E-07
1.331584 0.541608 0.0007340863 0.000206 5.814704E-07
1.353083 0.551408 -0.0007997676 -0.000346 7.594085E-07
1.373765 0.560853 -7.887225E-07 0.000175 3.046521E-08
1.394961 0.570551 3.482989E-05 0.000756 5.721233E-07
1.419318 0.581717 0.000203198 -0.000517 0.000000309
1.442292 0.592269 0.000317507 -0.00028 1.793036E-07
1.466015 0.603187 -0.000827971 -8.79E-05 6.932619E-07
1.491009 0.614711 0.0001815011 -0.000282 1.125947E-07
1.517118 0.626775 -0.0001960199 -0.000719 5.555583E-07
1.542468 0.638512 -9.328065E-05 -9.4E-05 1.754085E-08
1.569106 0.650869 0.0004268424 0.000201 2.224154E-07
1.597964 0.664285 -0.0006675949 -0.000409 0.000000613
1.625696 0.677204 -0.0001374698 0.000309 1.142985E-07
1.656583 0.691623 -0.0006856512 -0.00041 6.384922E-07
1.687373 0.70603 0.0001257131 -0.000433 0.000000203
1.719198 0.720952 0.0001105735 -0.000461 2.248664E-07
1.751236 0.736009 0.0004496446 3.21E-05 2.032108E-07
1.784762 0.751799 -0.0008703936 0.00029 8.415323E-07
1.819772 0.768326 0.0005894741 0.000356 4.741664E-07
1.857852 0.786345 0.0003154381 -0.000646 5.171343E-07

ng GoalSeek
Demo 12 - Method of lines for Partial Differential Equations
Solved by Explicit Euler's method

Definitions Equations/values

Final value ind.var. tf=6000 6000


Integration Stepsize h 30
Constants T1 = 0 0
alpha = 2.e-5 2.00E-05
deltax = .10 0.1
alpha/deltax^2 2.00E-03
Independent Variable t0=0 0 30
Variables T2 100 94
T3 100 100
T4 100 100
T5 100 100
T6 100 100
T7 100 100
T8 100 100
T9 100 100
T10 100 100
T11 = (4*T10-T9)/3 100 100
Differential equations f1=d(T2)/d(t) = alpha/deltax^2*(T3-2*T2+T1) -0.2 -0.176
f2=d(T3)/d(t) = alpha/deltax^2*(T4-2*T3+T2) 0 -0.012
f3=d(T4)/d(t) = alpha/deltax^2*(T5-2*T4+T3) 0 0
f4=d(T5)/d(t) = alpha/deltax^2*(T6-2*T5+T4) 0 0
f5=d(T6)/d(t) = alpha/deltax^2*(T7-2*T6+T5) 0 0
f6=d(T7)/d(t) = alpha/deltax^2*(T8-2*T7+T6) 0 0
f7=d(T8)/d(t) = alpha/deltax^2*(T9-2*T8+T7) 0 0
f8=d(T9)/d(t) = alpha/deltax^2*(T10-2*T9+T8) 0 0
f9=d(T10)/d(t) = alpha/deltax^2*(T11-2*T10+T9) 0 0

Temperature Profiles for a One-Dimentional Slab


120

100
Temperature (C)

80
T2
60 T3
T4
40 T5
20

0
0 1000 2000 3000 4000 5000 6000
Time (s)
Step no. (i) t T2 T3 T4
0 0 100 100 100
1 30 94 100 100
2 60 88.72 99.64 100
3 90 84.052 99.0064 99.9784
4 120 79.906144 98.167456 99.921376
5 150 76.20745408 97.17701248 99.82078048
6 180 72.89318034 96.07746506 99.6725561
7 210 69.9106466 94.90211344 99.4755425
8 240 67.21549582 93.67703117 99.2305837
9 270 64.77025819 92.4225522 98.93986467
10 300 62.54318034 91.15445331 98.60642231
11 330 60.5072659 89.88489507 98.23378868
12 360 58.63948769 88.62317093 97.82573409
13 390 56.92013943 87.37630373 97.3860853
14 420 55.33230092 86.14952077 96.91859993
15 450 53.86139605 84.94663232 96.42688255
16 480 52.49482647 83.77033316 95.91433159
17 510 51.22166728 82.62244267 95.38410863
18 540 50.03241377 81.5040961 94.83912385
19 570 48.91876988 80.41589683 94.28203286
20 600 47.8734713 79.35803737 93.71524144
21 630 46.89013699 78.33039565 93.14091547
22 660 45.96314429 77.33261132 92.56099429
23 690 45.08752365 76.36414628 91.97720583
24 720 44.25886959 75.42433249 91.39108272
25 750 43.47326519 74.51240973 90.8039786
26 780 42.72721795 73.62755519 90.21708401
27 810 42.01760511 72.76890669 89.63144176
28 840 41.3416269 71.9355807 89.04796138
29 870 40.69676651 71.12668631 88.46743265
30 900 40.08075571 70.3413359 87.89053805
31 930 39.49154518 69.57865322 87.31786423
32 960 38.92727895 68.8377794 86.74991249
33 990 38.38627224 68.11787736 86.1871082
34 1020 37.86699221 67.4181349 85.62980943
35 1050 37.36804124 66.73776681 85.07831461
36 1080 36.8881423 66.07601614 84.5328695
37 1110 36.42612619 65.43215492 83.99367337
38 1140 35.98092034 64.8054843 83.46088454
39 1170 35.55153896 64.19533448 82.93462537
40 1200 35.13707435 63.6010642 82.41498661
41 1230 34.73668928 63.02206015 81.90203139
42 1260 34.34961018 62.45773618 81.39579863
43 1290 33.97512113 61.90753236 80.8963062
44 1320 33.61255853 61.37091412 80.40355363
45 1350 33.26130636 60.84737115 79.91752457
46 1380 32.92079186 60.33641647 79.43818888
47 1410 32.59048183 59.83758534 78.96550464
48 1440 32.26987913 59.35043429 78.49941972
49 1470 31.95851969 58.8745401 78.03987336
50 1500 31.65596973 58.40949887 77.58679743
51 1530 31.3618233 57.95492504 77.14011762
52 1560 31.07570001 57.51045049 76.69975443
53 1590 30.79724303 57.0757237 76.2656241
54 1620 30.52611729 56.65040888 75.83763939
55 1650 30.26200775 56.23418522 75.41571025
56 1680 30.00461793 55.82674607 74.9997445
57 1710 29.75366854 55.42779829 74.58964827
58 1740 29.50889622 55.0370615 74.18532655
59 1770 29.27005236 54.65426749 73.78668356
60 1800 29.03690213 54.27915954 73.39362312
61 1830 28.80922344 53.91149191 73.00604893
62 1860 28.58680615 53.55102923 72.6238649
63 1890 28.36945116 53.19754598 72.24697531
64 1920 28.15696978 52.85082605 71.87528509
65 1950 27.94918297 52.51066222 71.50869991
66 1980 27.74592075 52.17685573 71.1471264
67 2010 27.5470216 51.84921587 70.79047222
68 2040 27.35233196 51.52755959 70.43864621
69 2070 27.1617057 51.21171113 70.09155844
70 2100 26.97500369 50.90150164 69.7491203
71 2130 26.79209334 50.59676889 69.41124459
72 2160 26.61284827 50.2973569 69.0778455
73 2190 26.43714789 50.0031157 68.74883871
74 2220 26.26487709 49.71390101 68.42414139
75 2250 26.0959259 49.429574 68.1036722
76 2280 25.93018923 49.150001 67.78735137
77 2310 25.76756658 48.87505332 67.47510063
78 2340 25.60796179 48.60460695 67.16684325
79 2370 25.45128279 48.33854242 66.86250406
80 2400 25.2974414 48.07674454 66.56200938
81 2430 25.14635311 47.81910224 66.26528707
82 2460 24.99793687 47.56550839 65.97226648
83 2490 24.85211495 47.31585958 65.68287845
84 2520 24.70881273 47.07005603 65.39705529
85 2550 24.56795856 46.82800139 65.11473077
86 2580 24.42948362 46.58960258 64.83584005
87 2610 24.29332174 46.35476969 64.56031974
88 2640 24.15940931 46.12341582 64.2881078
89 2670 24.02768514 45.89545695 64.01914355
90 2700 23.89809034 45.67081184 63.75336765
91 2730 23.77056821 45.4494019 63.49072205
92 2760 23.64506414 45.23115108 63.23114999
93 2790 23.52152551 45.0159858 62.97459597
94 2820 23.3999016 44.80383479 62.72100569
95 2850 23.28014349 44.59462906 62.47032607
96 2880 23.16220402 44.38830174 62.22250519
97 2910 23.04603764 44.18478809 61.97749229
98 2940 22.93160041 43.98402531 61.73523771
99 2970 22.81884988 43.78595256 61.4956929
100 3000 22.70774505 43.59051082 61.25881038
101 3030 22.59824629 43.39764285 61.02454371
102 3060 22.49031531 43.20729311 60.79284746
103 3090 22.38391506 43.0194077 60.56367719
104 3120 22.27900971 42.83393431 60.33698947
105 3150 22.1755646 42.65082214 60.11274176
106 3180 22.07354618 42.47002187 59.89089248
107 3210 21.97292195 42.29148556 59.67140093
108 3240 21.87366045 42.11516667 59.45422732
109 3270 21.7757312 41.94101993 59.23933268
110 3300 21.67910465 41.76900138 59.02667889
111 3330 21.58375217 41.59906822 58.81622865
112 3360 21.48964601 41.43117889 58.60794544
113 3390 21.39675922 41.26529291 58.40179353
114 3420 21.30506569 41.10137092 58.19773794
115 3450 21.21454006 40.93937463 57.9957444
116 3480 21.12515773 40.77926674 57.79577941
117 3510 21.03689481 40.62101096 57.59781011
118 3540 20.94972809 40.46457194 57.40180438
119 3570 20.86363503 40.30991526 57.20773071
120 3600 20.77859374 40.15700737 57.01555827
121 3630 20.69458294 40.00581561 56.82525687
122 3660 20.61158192 39.85630812 56.6367969
123 3690 20.52957058 39.70845388 56.45014939
124 3720 20.44852934 39.56222261 56.26528592
125 3750 20.36843918 39.41758481 56.08217867
126 3780 20.28928156 39.27451171 55.90080037
127 3810 20.21103848 39.13297522 55.72112429
128 3840 20.13369237 38.99294796 55.54312422
129 3870 20.05722617 38.8544032 55.36677448
130 3900 19.98162322 38.71731485 55.19204989
131 3930 19.90686732 38.58165746 55.01892577
132 3960 19.83294269 38.44740615 54.84737791
133 3990 19.75983394 38.31453665 54.67738257
134 4020 19.68752606 38.18302524 54.50891647
135 4050 19.61600445 38.05284876 54.34195677
136 4080 19.54525484 37.92398458 54.17648107
137 4110 19.47526334 37.79641059 54.01246741
138 4140 19.40601637 37.67010516 53.84989421
139 4170 19.33750072 37.54504718 53.68874033
140 4200 19.26970346 37.42121598 53.528985
141 4230 19.202612 37.29859137 53.37060785
142 4260 19.13621405 37.1771536 53.21358889
143 4290 19.07049758 37.05688334 53.05790848
144 4320 19.00545087 36.9377617 52.90354737
145 4350 18.94106247 36.81977019 52.75048663
146 4380 18.87732118 36.70289072 52.59870768
147 4410 18.81421608 36.58710556 52.4481923
148 4440 18.75173649 36.4723974 52.29892257
149 4470 18.68987195 36.35874925 52.15088091
150 4500 18.62861227 36.24614451 52.00405003
151 4530 18.56794747 36.13456691 51.85841298
152 4560 18.50786779 36.02400051 51.71395307
153 4590 18.44836368 35.9144297 51.57065393
154 4620 18.38942582 35.80583919 51.42849947
155 4650 18.33104508 35.69821401 51.28747388
156 4680 18.27321251 35.59153946 51.14756162
157 4710 18.21591938 35.48580118 51.00874742
158 4740 18.15915712 35.38098504 50.87101627
159 4770 18.10291737 35.27707724 50.73435343
160 4800 18.04719192 35.17406422 50.59874439
161 4830 17.99197274 35.07193269 50.46417489
162 4860 17.93725197 34.97066963 50.33063093
163 4890 17.88302191 34.87026225 50.19809873
164 4920 17.82927502 34.77069801 50.06656473
165 4950 17.7760039 34.67196464 49.93601561
166 4980 17.72320131 34.57405005 49.80643827
167 5010 17.67086015 34.47694242 49.67781983
168 5040 17.61897348 34.38063013 49.55014761
169 5070 17.56753447 34.28510178 49.42340913
170 5100 17.51653644 34.19034618 49.29759215
171 5130 17.46597284 34.09635236 49.1726846
172 5160 17.41583724 34.00310952 49.0486746
173 5190 17.36612334 33.91060709 48.92555047
174 5220 17.31682497 33.81883467 48.80330074
175 5250 17.26793605 33.72778205 48.68191408
176 5280 17.21945065 33.63743921 48.56137938
177 5310 17.17136292 33.54779631 48.44168568
178 5340 17.12366715 33.45884367 48.32282221
179 5370 17.07635771 33.37057179 48.20477836
180 5400 17.02942909 33.28297134 48.08754368
181 5430 16.98287588 33.19603314 47.97110791
182 5460 16.93669277 33.10974819 47.85546092
183 5490 16.89087453 33.02410763 47.74059275
184 5520 16.84541604 32.93910275 47.62649359
185 5550 16.80031228 32.854725 47.5131538
186 5580 16.75555831 32.77096596 47.40056386
187 5610 16.71114927 32.68781738 47.28871442
188 5640 16.6670804 32.60527111 47.17759625
189 5670 16.62334702 32.52331918 47.06720028
190 5700 16.57994453 32.44195372 46.95751757
191 5730 16.53686841 32.361167 46.84853931
192 5760 16.49411422 32.28095142 46.74025684
193 5790 16.4516776 32.20129951 46.6326616
194 5820 16.40955426 32.12220392 46.5257452
195 5850 16.36773998 32.04365742 46.41949933
196 5880 16.32623063 31.96565289 46.31391584
197 5910 16.28502213 31.88818333 46.20898667
198 5940 16.24411047 31.81124186 46.10470392
199 5970 16.20349173 31.7348217 46.00105977
200 6000 16.16316202 31.65891618 45.89804653
T5 T6 T7 T8 T9 T10 T11
100 100 100 100 100 100 100
100 100 100 100 100 100 100
100 100 100 100 100 100 100
100 100 100 100 100 100 100
99.998704 100 100 100 100 100 100
99.99414208 99.99992224 100 100 100 100 100
99.98408719 99.9995801 99.99999533 100 100 100 100
99.9663249 99.99867544 99.9999707 99.99999972 100 100 100
99.93881899 99.99681212 99.99989473 99.999998 99.99999998 100 100
99.89980446 99.99351749 99.99971597 99.99999192 99.99999986 100 100
99.84783085 99.98826662 99.99936061 99.99997584 99.9999994 99.99999999 100
99.78177249 99.98050611 99.99873189 99.99994034 99.99999802 99.99999997 100
99.70081748 99.96967564 99.99771085 99.99987129 99.99999467 99.99999989 100
99.60444396 99.95522626 99.99615836 99.99974907 99.99998758 99.99999968 100
99.49238938 99.93663525 99.99391788 99.99954794 99.999974 99.9999992 100
99.36461677 99.91341746 99.99081872 99.9992357 99.99994995 99.99999819 100
99.22128075 99.88513349 99.98667967 99.99877353 99.99990999 99.99999626 100
99.06269497 99.8513951 99.98131253 99.99811609 99.99984698 99.99999281 100
98.8893018 99.81186813 99.9745257 99.99721173 99.99975187 99.99998698 100.0001
98.7016451 99.76627361 99.9661274 99.99600298 99.99961357 99.99997757 100.0001
98.50034608 99.71438713 99.95592871 99.99442708 99.99941878 99.99996301 100.0001
98.28608226 99.65603716 99.94374612 99.99241668 99.99915193 99.99994124 100.0002
98.05956955 99.5911024 99.92940381 99.98990056 99.99879517 99.99990967 100.0003
97.821547 99.51950851 99.91273553 99.98680443 99.99832836 99.99986509 100.0004
97.57276422 99.44122445 99.89358605 99.98305173 99.99772913 99.99980362 100.0005
97.31397095 99.35625853 99.87181229 99.97856444 99.99697296 99.99972064 100.0006
97.04590866 99.2646545 99.84728419 99.97326382 99.99603331 99.99961073 100.0008
96.76930393 99.16648753 99.81988519 99.96707121 99.99488178 99.99946764 100.001
96.48486322 99.06186037 99.78951249 99.95990868 99.9934883 99.9992842 100.0012
96.19326894 98.95089967 99.75607714 99.95169969 99.99182128 99.99905237 100.0015
95.8951766 98.83375248 99.71950384 99.94236963 99.98984785 99.99876312 100.0017
95.59121284 98.710583 99.67973071 99.93184638 99.98753407 99.99840651 100.002
95.28197414 98.58156966 99.63670878 99.9200607 99.98484516 99.99797162 100.0023
94.96802617 98.44690227 99.59040155 99.90694665 99.98174568 99.99744656 100.0027
94.64990366 98.30677966 99.5407843 99.89244189 99.97819979 99.99681852 100.003
94.32811056 98.16140738 99.48784348 99.87648791 99.97417144 99.99607377 100.0034
94.00312062 98.01099574 99.43157598 99.85903025 99.96962457 99.99519768 100.0037
93.67537806 97.85575805 99.37198842 99.84001865 99.96452329 99.99417475 100.0041
93.34529857 97.69590907 99.30909641 99.81940712 99.9588321 99.9929887 100.0044
93.01327036 97.53166368 99.24292381 99.79715397 99.952516 99.99162243 100.0047
92.67965526 97.36323569 99.17350201 99.77322189 99.94554066 99.99005817 100.0049
92.34478997 97.19083684 99.10086923 99.74757782 99.93787259 99.98827747 100.0051
92.00898727 97.01467597 99.0250698 99.72019299 99.9294792 99.98626128 100.0052
91.67253727 96.83495828 98.94615356 99.69104277 99.92032895 99.98399 100.0052
91.33570867 96.65188474 98.8641752 99.66010659 99.91039144 99.98144355 100.0051
90.99874993 96.4656516 98.77919365 99.6273678 99.89963748 99.97860147 100.0049
90.66189051 96.27645002 98.69127158 99.59281353 99.88803913 99.97544291 100.0046
90.32534198 96.08446575 98.6004748 99.55643455 99.87556983 99.97194676 100.0041
89.98929917 95.88987886 98.50687184 99.51822508 99.86220432 99.96809168 100.0034
89.65394118 95.69286366 98.41053346 99.47818264 99.84791881 99.96385619 100.0025
89.31943246 95.4935885 98.31153222 99.43630786 99.83269088 99.95921869 100.0014
88.98592372 95.29221576 98.20994214 99.3926043 99.81649957 99.95415758 100
88.65355288 95.08890182 98.10583828 99.34707829 99.79932534 99.94865126 99.99843
88.32244591 94.88379707 97.9992965 99.29973871 99.78115007 99.94267822 99.99652
87.99271767 94.67704597 97.89039306 99.25059686 99.76195708 99.9362171 99.9943
87.66447267 94.4687871 97.77920447 99.19966625 99.74173106 99.9292467 99.99175
87.33780579 94.25915327 97.66580713 99.14696243 99.72045811 99.92174607 99.98884
87.01280296 94.04827165 97.55027722 99.09250285 99.69812565 99.91369455 99.98555
86.6895418 93.83626387 97.43269042 99.03630668 99.67472242 99.9050718 99.98185
86.36809221 93.62324614 97.3131218 98.97839465 99.65023843 99.89585782 99.97773
86.04851693 93.40932944 97.19164563 98.91878891 99.62466497 99.88603304 99.97316
85.73087205 93.19461966 97.06833526 98.85751287 99.59799449 99.87557832 99.96811
85.41520752 92.97921774 96.94326298 98.79459111 99.57022062 99.86447497 99.96256
85.10156757 92.76321984 96.81649995 98.7300492 99.54133811 99.85270479 99.95649
84.78999118 92.54671751 96.6881161 98.66391358 99.51134278 99.84025013 99.94989
84.48051239 92.32979785 96.55818003 98.59621148 99.48023147 99.82709383 99.94271
84.17316077 92.11254365 96.42675899 98.52697079 99.44800201 99.81321934 99.93496
83.86796168 91.8950336 96.29391878 98.45621996 99.41465318 99.79861065 99.9266
83.56493663 91.67734239 96.15972374 98.38398788 99.38018463 99.78325235 99.91761
83.26410355 91.45954093 96.02423671 98.31030384 99.34459689 99.76712964 99.90797
82.96547708 91.24169643 95.88751899 98.23519739 99.30789127 99.75022833 99.89767
82.66906884 91.02387263 95.74963034 98.15869832 99.27006986 99.73253485 99.88669
82.37488761 90.80612986 95.61062895 98.08083653 99.23113547 99.71403625 99.875
82.08293962 90.58852527 95.47057146 98.00164202 99.19109158 99.69472022 99.8626
81.7932287 90.3711129 95.32951293 97.92114476 99.14994232 99.67457507 99.84945
81.50575652 90.15394385 95.18750683 97.8393747 99.10769243 99.65358976 99.83556
81.2205227 89.93706639 95.04460513 97.75636169 99.06434721 99.63175387 99.82089
80.93752504 89.72052609 94.9008582 97.67213543 99.01991248 99.6090576 99.80544
80.65675964 89.50436596 94.7563149 97.58672542 98.97439456 99.5854918 99.78919
80.37822103 89.28862651 94.6110226 97.50016094 98.92780025 99.56104791 99.77213
80.10190235 89.07334595 94.46502713 97.41247099 98.88013675 99.535718 99.75425
79.82779538 88.85856021 94.31837289 97.32368431 98.83141168 99.50949475 99.73552
79.55589077 88.64430308 94.17110282 97.23382927 98.78163302 99.48237143 99.71595
79.28617805 88.43060632 94.02325842 97.1429339 98.7308091 99.45434189 99.69552
79.01864577 88.21749975 93.87487982 97.05102589 98.67894856 99.42540058 99.67422
78.75328158 88.00501132 93.72600578 96.95813248 98.62606032 99.3955425 99.65204
78.49007232 87.7931672 93.57667372 96.86428055 98.57215358 99.36476321 99.62897
78.22900408 87.5819919 93.42691974 96.76949652 98.51723777 99.33305883 99.605
77.97006229 87.3715083 93.27677867 96.67380639 98.46132256 99.30042598 99.58013
77.71323178 87.16173776 93.12628411 96.5772357 98.4044178 99.26686185 99.55434
77.45849684 86.95270018 92.97546843 96.47980953 98.34653351 99.23236409 99.52764
77.20584129 86.74441408 92.8243628 96.3815525 98.28767991 99.19693086 99.50001
76.9552485 86.53689663 92.67299726 96.28248876 98.22786732 99.16056082 99.47146
76.70670148 86.33016378 92.52140071 96.18264199 98.16710622 99.12325308 99.44197
76.46018289 86.12423026 92.36960097 96.08203536 98.10540718 99.08500721 99.41154
76.2156751 85.91910966 92.21762479 95.98069161 98.04278087 99.04582321 99.38017
75.97316023 85.7148145 92.06549789 95.87863296 97.97923806 99.00570151 99.34786
75.73262018 85.51135624 91.91324499 95.77588116 97.91478956 98.96464298 99.31459
75.49403667 85.30874541 91.76088984 95.67245749 97.84944626 98.92264884 99.28038
75.25739126 85.10699155 91.60845523 95.56838276 97.78321909 98.87972074 99.24522
75.02266538 84.90610335 91.45596306 95.46367729 97.71611901 98.83586067 99.20911
74.78984036 84.70608866 91.30343433 95.35836094 97.648157 98.791071 99.17204
74.55889745 84.5069545 91.15088919 95.2524531 97.57934408 98.74535444 99.13402
74.32981788 84.30870716 90.99834694 95.14597273 97.50969124 98.69871403 99.09505
74.10258279 84.11135219 90.8458261 95.03893829 97.4392095 98.65115312 99.05513
73.87717336 83.91489446 90.6933444 94.93136783 97.36790984 98.60267537 99.01426
73.65357073 83.71933819 90.54091881 94.82327895 97.29580325 98.55328475 98.97245
73.43175608 83.52468698 90.38856558 94.7146888 97.22290069 98.50298549 98.92968
73.21171062 83.33094384 90.23630026 94.60561412 97.14921306 98.4517821 98.88597
72.99341562 83.13811123 90.0841377 94.49607122 97.07475127 98.39967934 98.84132
72.77685238 82.94619108 89.93209213 94.38607601 96.99952615 98.34668222 98.79573
72.56200229 82.75518482 89.7801771 94.27564399 96.9235485 98.29279597 98.74921
72.34884683 82.56509341 89.62840557 94.16479025 96.84682908 98.23802607 98.70176
72.13736754 82.37591734 89.47678992 94.0535295 96.76937857 98.18237819 98.65338
71.92754609 82.18765671 89.32534194 93.94187607 96.6912076 98.12585821 98.60408
71.71936423 82.00031119 89.17407288 93.82984391 96.61232675 98.06847219 98.55385
71.51280386 81.81388007 89.02299344 93.71744662 96.5327465 98.01022637 98.50272
71.30784697 81.6283623 88.87211383 93.60469742 96.4524773 97.95112717 98.45068
71.10447568 81.44375647 88.72144375 93.4916092 96.3715295 97.89118118 98.39773
70.90267225 81.26006086 88.57099244 93.37819449 96.28991338 97.83039511 98.34389
70.70241907 81.07727344 88.42076867 93.2644655 96.20763915 97.76877584 98.28915
70.50369868 80.89539189 88.27078077 93.15043411 96.12471694 97.70633037 98.23353
70.30649377 80.71441363 88.12103663 93.03611188 96.04115677 97.64306584 98.17704
70.11078715 80.53433582 87.97154377 92.92151006 95.95696862 97.57898947 98.11966
69.9165618 80.35515538 87.82230927 92.80663959 95.87216236 97.51410864 98.06142
69.72380086 80.176869 87.67333985 92.69151114 95.78674777 97.44843079 98.00233
69.53248762 79.99947316 87.52464188 92.57613506 95.70073455 97.38196347 97.94237
69.34260552 79.82296415 87.37622135 92.46052144 95.61413232 97.31471431 97.88157
69.15413816 79.64733806 87.22808392 92.34468009 95.52695059 97.24669103 97.81994
68.96706932 79.47259082 87.08023494 92.22862055 95.43919878 97.17790141 97.75747
68.78138292 79.29871818 86.93267943 92.11235211 95.35088625 97.10835331 97.69418
68.59706305 79.12571574 86.78542211 91.99588379 95.26202222 97.03805463 97.63007
68.41409398 78.95357896 86.63846743 91.8792244 95.17261586 96.96701333 97.56515
68.23246011 78.78230317 86.49181954 91.76238247 95.08267622 96.89523743 97.49942
68.05214604 78.61188357 86.34548234 91.64536632 94.99221227 96.82273498 97.43291
67.87313652 78.44231524 86.19945945 91.52818404 94.90123287 96.74951407 97.36561
67.69541646 78.27359317 86.05375427 91.41084349 94.80974682 96.67558283 97.29753
67.51897094 78.10571223 85.90836996 91.29335234 94.71776278 96.60094939 97.22868
67.3437852 77.93866722 85.76330944 91.17571802 94.62528935 96.52562192 97.15907
67.16984466 77.77245283 85.61857542 91.05794779 94.53233502 96.44960862 97.0887
66.99713489 77.6070637 85.47417041 90.94004868 94.4389082 96.37291768 97.01759
66.82564163 77.44249437 85.3300967 90.82202755 94.3450172 96.2955573 96.94574
66.65535077 77.27873935 85.18635641 90.70389108 94.25067023 96.21753569 96.87316
66.48624837 77.11579306 85.04295147 90.58564575 94.15587541 96.13886107 96.79986
66.31832066 76.95364988 84.89988362 90.46729787 94.06064077 96.05954165 96.72584
66.15155401 76.79230415 84.75715445 90.34885359 93.96497425 95.97958561 96.65112
65.98593498 76.63175016 84.61476538 90.23031888 93.86888369 95.89900116 96.57571
65.82145025 76.47198216 84.47271768 90.11169956 93.77237685 95.81779646 96.4996
65.65808669 76.31299438 84.33101246 89.99300128 93.67546139 95.73597967 96.42282
65.4958313 76.154781 84.1896507 89.87422956 93.57814488 95.65355894 96.34536
65.33467126 75.9973362 84.04863325 89.75538975 93.4804348 95.57054238 96.26724
65.17459388 75.84065413 83.90796082 89.63648706 93.38233856 95.48693808 96.18847
65.01558664 75.68472892 83.76763399 89.51752658 93.28386344 95.4027541 96.10905
64.85763717 75.52955468 83.62765324 89.39851323 93.18501666 95.31799847 96.02899
64.70073322 75.37512555 83.48801893 89.27945184 93.08580537 95.2326792 95.9483
64.54486274 75.22143561 83.3487313 89.16034708 92.98623659 95.14680424 95.86699
64.39001378 75.06847898 83.20979051 89.0412035 92.88631727 95.06038154 95.78507
64.23617456 74.91624976 83.07119659 88.92202555 92.7860543 94.97341897 95.70254
64.08333344 74.76474206 82.93294952 88.80281754 92.68545446 94.88592438 95.61941
63.93147893 74.61394999 82.79504915 88.68358367 92.58452444 94.79790558 95.5357
63.78059966 74.46386767 82.65749528 88.56432805 92.48327086 94.70937034 95.4514
63.63068443 74.31448925 82.52028759 88.44505465 92.38170026 94.62032636 95.36654
63.48172215 74.16580886 82.38342571 88.32576736 92.27981909 94.53078132 95.2811
63.33370188 74.01782067 82.2469092 88.20646997 92.17763372 94.44074283 95.19511
63.18661281 73.87051885 82.11073753 88.08716615 92.07515044 94.35021846 95.10857
63.04044429 73.72389761 81.97491013 87.96785949 91.97237546 94.25921574 95.0215
62.89518577 73.57795116 81.83942634 87.84855348 91.86931492 94.16774213 94.93388
62.75082684 73.43267375 81.70428546 87.72925154 91.76597487 94.07580504 94.84575
62.60735724 73.28805964 81.56948672 87.60995698 91.66236128 93.98341183 94.7571
62.4647668 73.14410312 81.43502931 87.49067302 91.55848005 93.89056981 94.66793
62.32304552 73.00079851 81.30091236 87.37140282 91.45433702 93.79728622 94.57827
62.1821835 72.85814016 81.16713496 87.25214944 91.34993792 93.70356825 94.48811
62.04217097 72.71612245 81.03369614 87.13291588 91.24528843 93.60942304 94.39747
61.90299827 72.57473978 80.9005949 87.01370505 91.14039415 93.51485766 94.30635
61.76465589 72.4339866 80.7678302 86.89451979 91.03526062 93.41987912 94.21475
61.62713443 72.29385737 80.63540096 86.77536286 90.92989328 93.32449438 94.12269
61.49042458 72.15434661 80.50330606 86.65623697 90.82429752 93.22871033 94.03018
61.35451719 72.01544886 80.37154435 86.53714475 90.71847865 93.13253382 93.93722
61.2194032 71.87715869 80.24011464 86.41808876 90.61244193 93.03597161 93.84381
61.08507367 71.73947072 80.10901573 86.2990715 90.50619252 92.93903043 93.74998
60.95151978 71.60237959 79.97824638 86.18009542 90.39973553 92.84171691 93.65571
60.8187328 71.46588001 79.84780531 86.06116288 90.29307601 92.74403765 93.56102
60.68670414 71.3299667 79.71769125 85.94227622 90.18621892 92.64599919 93.46593
60.5554253 71.19463442 79.58790287 85.82343768 90.07916918 92.54760798 93.37042
60.42488789 71.05987798 79.45843886 85.70464948 89.97193161 92.44887043 93.27452
60.29508364 70.92569223 79.32929784 85.58591377 89.86451101 92.34979287 93.17822
60.16600436 70.79207205 79.20047846 85.46723265 89.75691209 92.2503816 93.08154
60.037642 70.65901237 79.07197933 85.34860817 89.6491395 92.15064282 92.98448
59.90998856 70.52650817 78.94379904 85.23004232 89.54119782 92.05058269 92.88704
59.7830362 70.39455444 78.81593618 85.11153705 89.43309158 91.95020729 92.78925
59.65677714 70.26314625 78.68838933 84.99309427 89.32482525 91.84952266 92.69109
59.53120371 70.13227869 78.56115704 84.87471583 89.21640323 91.74853477 92.59258
59.40630835 70.00194689 78.43423787 84.75640355 89.10782988 91.6472495 92.49372
59.28208357 69.87214604 78.30763035 84.63815919 88.99910948 91.54567272 92.39453
59.158522 69.74287135 78.18133302 84.51998447 88.89024626 91.44381019 92.295
59.03561635 69.61411809 78.05534441 84.40188109 88.78124439 91.34166763 92.19514
58.91335944 69.48588156 77.92966303 84.28385069 88.67210798 91.2392507 92.09496
58.79174415 69.35815712 77.8042874 84.16589487 88.56284111 91.13656499 91.99447
58.67076348 69.23094016 77.67921603 84.0480152 88.45344777 91.03361604 91.89367
58.5504105 69.10422611 77.55444743 83.9302132 88.34393191 90.93040931 91.79257
58.4306784 68.97801046 77.4299801 83.81249038 88.23429743 90.82695021 91.69117
f1 f2 f3 f4 f5 f6 f7 f8 f9
-0.2 0 0 0 0 0 0 0 0
-0.176 -0.012 0 0 0 0 0 0 0
-0.1556 -0.02112 -0.00072 0 0 0 0 0 0
-0.138195 -0.027965 -0.001901 -4.32E-05 0 0 0 0 0
-0.12329 -0.033015 -0.003353 -0.000152 -2.59E-06 0 0 0 0
-0.110476 -0.036652 -0.004941 -0.000335 -1.14E-05 -1.56E-07 0 0 0
-0.099418 -0.039178 -0.006567 -0.000592 -3E-05 -8.21E-07 -9.33E-09 0 0
-0.089838 -0.040836 -0.008165 -0.000917 -6.21E-05 -2.53E-06 -5.75E-08 -5.6E-10 0
-0.081508 -0.041816 -0.009691 -0.0013 -0.00011 -6E-06 -2E-07 -3.94E-09 -2.24E-11
-0.074236 -0.04227 -0.011115 -0.001732 -0.000175 -1.18E-05 -5.36E-07 -1.56E-08 -1.79E-10
-0.067864 -0.042319 -0.012421 -0.002202 -0.000259 -2.1E-05 -1.18E-06 -4.59E-08 -8E-10
-0.062259 -0.042057 -0.013602 -0.002699 -0.000361 -3.4E-05 -2.3E-06 -1.11E-07 -2.6E-09
-0.057312 -0.041562 -0.014655 -0.003212 -0.000482 -5.17E-05 -4.07E-06 -2.36E-07 -7E-09
-0.052928 -0.040893 -0.015583 -0.003735 -0.00062 -7.47E-05 -6.7E-06 -4.53E-07 -1.61E-08
-0.04903 -0.040096 -0.016391 -0.004259 -0.000774 -0.000103 -1E-05 -8E-07 -3.36E-08
-0.045552 -0.03921 -0.017085 -0.004778 -0.000943 -0.000138 -1.54E-05 -1.33E-06 -6.43E-08
-0.042439 -0.038263 -0.017674 -0.005286 -0.001125 -0.000179 -2.19E-05 -2.1E-06 -1.15E-07
-0.039642 -0.037278 -0.018166 -0.00578 -0.001318 -0.000226 -3E-05 -3.17E-06 -1.94E-07
-0.037121 -0.036273 -0.01857 -0.006255 -0.00152 -0.00028 -4E-05 -4.61E-06 -3.13E-07
-0.034843 -0.035262 -0.018893 -0.00671 -0.00173 -0.00034 -5.25E-05 -6.49E-06 -4.85E-07
-0.032778 -0.034255 -0.019144 -0.007142 -0.001945 -0.000406 -6.7E-05 -8.89E-06 -7.26E-07
-0.0309 -0.033259 -0.019331 -0.00755 -0.002164 -0.000478 -8.39E-05 -1.19E-05 -1.05E-06
-0.029187 -0.032282 -0.01946 -0.007934 -0.002386 -0.000556 -0.000103 -1.56E-05 -1.49E-06
-0.027622 -0.031327 -0.019537 -0.008293 -0.002609 -0.000638 -0.000125 -2E-05 -2E-06
-0.026187 -0.030397 -0.01957 -0.008626 -0.002832 -0.000726 -0.00015 -2.52E-05 -2.77E-06
-0.024868 -0.029495 -0.019563 -0.008935 -0.003053 -0.000818 -0.000177 -3.13E-05 -3.66E-06
-0.023654 -0.028622 -0.019521 -0.00922 -0.003272 -0.000913 -0.000206 -3.84E-05 -4.77E-06
-0.022533 -0.027778 -0.019449 -0.009481 -0.003488 -0.001012 -0.000239 -4.64E-05 -6.11E-06
-0.021495 -0.026963 -0.019351 -0.00972 -0.003699 -0.001115 -0.000274 -5.56E-05 -7.73E-06
-0.020534 -0.026178 -0.01923 -0.009936 -0.003905 -0.001219 -0.000311 -6.58E-05 -9.64E-06
-0.01964 -0.025423 -0.019089 -0.010132 -0.004106 -0.001326 -0.000351 -7.71E-05 -1.19E-05
-0.018809 -0.024696 -0.018932 -0.010308 -0.0043 -0.001434 -0.000393 -9E-05 -1.45E-05
-0.018034 -0.023997 -0.01876 -0.010465 -0.004489 -0.001544 -0.000437 -0.000103 -1.75E-05
-0.017309 -0.023325 -0.018577 -0.010604 -0.004671 -0.001654 -0.000483 -0.000118 -2.09E-05
-0.016632 -0.022679 -0.018383 -0.010726 -0.004846 -0.001765 -0.000532 -0.000134 -2.48E-05
-0.015997 -0.022058 -0.018182 -0.010833 -0.005014 -0.001876 -0.000582 -0.000152 -2.92E-05
-0.015401 -0.021462 -0.017973 -0.010925 -0.005175 -0.001986 -0.000634 -0.00017 -3.41E-05
-0.01484 -0.020889 -0.01776 -0.011003 -0.005328 -0.002096 -0.000687 -0.00019 -4E-05
-0.014313 -0.020338 -0.017542 -0.011068 -0.005475 -0.002206 -0.000742 -0.000211 -4.55E-05
-0.013815 -0.019809 -0.017321 -0.011121 -0.005614 -0.002314 -0.000798 -0.000233 -5.21E-05
-0.013346 -0.0193 -0.017099 -0.011162 -0.005747 -0.002421 -0.000855 -0.000256 -5.94E-05
-0.012903 -0.018811 -0.016874 -0.011193 -0.005872 -0.002527 -0.000913 -0.00028 -6.72E-05
-0.012483 -0.01834 -0.01665 -0.011215 -0.005991 -0.002631 -0.000972 -0.000305 -7.57E-05
-0.012085 -0.017887 -0.016425 -0.011228 -0.006102 -0.002733 -0.001031 -0.000331 -8.49E-05
-0.011708 -0.017451 -0.016201 -0.011232 -0.006208 -0.002833 -0.001091 -0.000358 -9.47E-05
-0.01135 -0.017032 -0.015978 -0.011229 -0.006307 -0.002931 -0.001152 -0.000387 -0.000105
-0.01101 -0.016628 -0.015756 -0.011218 -0.006399 -0.003027 -0.001213 -0.000416 -0.000117
-0.010687 -0.016238 -0.015536 -0.011201 -0.006486 -0.00312 -0.001274 -0.000446 -0.000129
-0.010379 -0.015863 -0.015318 -0.011179 -0.006567 -0.003211 -0.001335 -0.000476 -0.000141
-0.010085 -0.015501 -0.015103 -0.01115 -0.006643 -0.0033 -0.001396 -0.000508 -0.000155
-0.009805 -0.015152 -0.014889 -0.011117 -0.006712 -0.003386 -0.001457 -0.00054 -0.000169
-0.009537 -0.014816 -0.014679 -0.011079 -0.006777 -0.00347 -0.001518 -0.000572 -0.000184
-0.009282 -0.014491 -0.014471 -0.011037 -0.006837 -0.003551 -0.001578 -0.000606 -0.000199
-0.009038 -0.014177 -0.014266 -0.010991 -0.006892 -0.00363 -0.001638 -0.00064 -0.000215
-0.008804 -0.013874 -0.014064 -0.010941 -0.006942 -0.003706 -0.001698 -0.000674 -0.000232
-0.00858 -0.013581 -0.013866 -0.010889 -0.006988 -0.00378 -0.001757 -0.000709 -0.00025
-0.008365 -0.013298 -0.01367 -0.010833 -0.007029 -0.003851 -0.001815 -0.000744 -0.000268
-0.008159 -0.013025 -0.013477 -0.010775 -0.007067 -0.00392 -0.001873 -0.00078 -0.000287
-0.007961 -0.01276 -0.013288 -0.010715 -0.007101 -0.003986 -0.00193 -0.000816 -0.000307
-0.007772 -0.012504 -0.013102 -0.010653 -0.007131 -0.004049 -0.001987 -0.000852 -0.000327
-0.007589 -0.012256 -0.012919 -0.010588 -0.007157 -0.00411 -0.002043 -0.000889 -0.000348
-0.007414 -0.012015 -0.012739 -0.010522 -0.00718 -0.004169 -0.002097 -0.000926 -0.00037
-0.007245 -0.011783 -0.012563 -0.010455 -0.0072 -0.004225 -0.002151 -0.000963 -0.000392
-0.007083 -0.011557 -0.01239 -0.010386 -0.007217 -0.004279 -0.002205 -0.001 -0.000415
-0.006926 -0.011339 -0.01222 -0.010316 -0.007231 -0.004331 -0.002257 -0.001037 -0.000439
-0.006775 -0.011127 -0.012052 -0.010245 -0.007242 -0.004381 -0.002308 -0.001074 -0.000462
-0.00663 -0.010921 -0.011888 -0.010173 -0.00725 -0.004428 -0.002358 -0.001112 -0.000487
-0.00649 -0.010722 -0.011728 -0.010101 -0.007256 -0.004473 -0.002408 -0.001149 -0.000512
-0.006354 -0.010528 -0.01157 -0.010028 -0.00726 -0.004516 -0.002456 -0.001186 -0.000537
-0.006223 -0.01034 -0.011415 -0.009954 -0.007261 -0.004557 -0.002504 -0.001224 -0.000563
-0.006097 -0.010158 -0.011263 -0.00988 -0.007261 -0.004596 -0.00255 -0.001261 -0.00059
-0.005975 -0.00998 -0.011113 -0.009806 -0.007258 -0.004633 -0.002595 -0.001298 -0.000617
-0.005857 -0.009808 -0.010967 -0.009732 -0.007253 -0.004669 -0.00264 -0.001335 -0.000644
-0.005742 -0.00964 -0.010823 -0.009657 -0.007247 -0.004702 -0.002683 -0.001372 -0.000672
-0.005632 -0.009478 -0.010682 -0.009582 -0.007239 -0.004734 -0.002726 -0.001408 -0.0007
-0.005525 -0.009319 -0.010544 -0.009508 -0.007229 -0.004763 -0.002767 -0.001445 -0.000728
-0.005421 -0.009165 -0.010408 -0.009433 -0.007218 -0.004792 -0.002808 -0.001481 -0.000757
-0.00532 -0.009015 -0.010275 -0.009359 -0.007205 -0.004818 -0.002847 -0.001517 -0.000786
-0.005223 -0.008869 -0.010145 -0.009285 -0.007191 -0.004843 -0.002885 -0.001553 -0.000815
-0.005128 -0.008727 -0.010016 -0.009211 -0.007176 -0.004867 -0.002923 -0.001589 -0.000844
-0.005036 -0.008588 -0.009891 -0.009137 -0.00716 -0.004888 -0.00296 -0.001624 -0.000874
-0.004947 -0.008453 -0.009767 -0.009063 -0.007142 -0.004909 -0.002995 -0.001659 -0.000904
-0.004861 -0.008322 -0.009646 -0.00899 -0.007123 -0.004928 -0.00303 -0.001694 -0.000934
-0.004777 -0.008193 -0.009527 -0.008918 -0.007104 -0.004946 -0.003064 -0.001729 -0.000965
-0.004695 -0.008068 -0.009411 -0.008845 -0.007083 -0.004962 -0.003096 -0.001763 -0.000995
-0.004616 -0.007947 -0.009296 -0.008774 -0.007061 -0.004978 -0.003128 -0.001797 -0.001026
-0.004539 -0.007828 -0.009184 -0.008702 -0.007039 -0.004992 -0.003159 -0.001831 -0.001057
-0.004464 -0.007712 -0.009074 -0.008631 -0.007016 -0.005005 -0.00319 -0.001864 -0.001088
-0.004391 -0.007599 -0.008965 -0.008561 -0.006992 -0.005016 -0.003219 -0.001897 -0.001119
-0.00432 -0.007488 -0.008859 -0.008491 -0.006968 -0.005027 -0.003248 -0.001929 -0.00115
-0.004251 -0.00738 -0.008755 -0.008422 -0.006943 -0.005037 -0.003275 -0.001962 -0.001181
-0.004183 -0.007275 -0.008652 -0.008353 -0.006917 -0.005046 -0.003302 -0.001994 -0.001212
-0.004118 -0.007172 -0.008552 -0.008285 -0.006891 -0.005053 -0.003328 -0.002025 -0.001244
-0.004054 -0.007072 -0.008453 -0.008217 -0.006864 -0.00506 -0.003354 -0.002057 -0.001275
-0.003992 -0.006974 -0.008356 -0.00815 -0.006837 -0.005066 -0.003378 -0.002088 -0.001306
-0.003931 -0.006878 -0.008261 -0.008084 -0.00681 -0.005071 -0.003402 -0.002118 -0.001337
-0.003872 -0.006784 -0.008167 -0.008018 -0.006782 -0.005075 -0.003425 -0.002148 -0.001369
-0.003815 -0.006692 -0.008075 -0.007953 -0.006754 -0.005079 -0.003447 -0.002178 -0.0014
-0.003758 -0.006602 -0.007985 -0.007888 -0.006725 -0.005081 -0.003469 -0.002208 -0.001431
-0.003703 -0.006515 -0.007896 -0.007824 -0.006696 -0.005083 -0.00349 -0.002237 -0.001462
-0.00365 -0.006429 -0.007809 -0.007761 -0.006667 -0.005084 -0.003511 -0.002265 -0.001493
-0.003598 -0.006345 -0.007723 -0.007698 -0.006638 -0.005085 -0.00353 -0.002294 -0.001524
-0.003547 -0.006263 -0.007639 -0.007636 -0.006608 -0.005085 -0.003549 -0.002322 -0.001555
-0.003497 -0.006182 -0.007556 -0.007575 -0.006578 -0.005084 -0.003568 -0.002349 -0.001585
-0.003448 -0.006104 -0.007475 -0.007514 -0.006549 -0.005083 -0.003586 -0.002377 -0.001616
-0.003401 -0.006027 -0.007395 -0.007453 -0.006519 -0.005081 -0.003603 -0.002404 -0.001646
-0.003354 -0.005951 -0.007316 -0.007394 -0.006488 -0.005078 -0.00362 -0.00243 -0.001677
-0.003309 -0.005877 -0.007239 -0.007335 -0.006458 -0.005076 -0.003636 -0.002456 -0.001707
-0.003264 -0.005805 -0.007163 -0.007277 -0.006428 -0.005072 -0.003651 -0.002482 -0.001737
-0.003221 -0.005734 -0.007088 -0.007219 -0.006397 -0.005068 -0.003667 -0.002508 -0.001767
-0.003178 -0.005664 -0.007015 -0.007162 -0.006367 -0.005064 -0.003681 -0.002533 -0.001796
-0.003137 -0.005596 -0.006943 -0.007105 -0.006336 -0.005059 -0.003695 -0.002557 -0.001826
-0.003096 -0.00553 -0.006872 -0.007049 -0.006306 -0.005054 -0.003709 -0.002582 -0.001855
-0.003056 -0.005464 -0.006802 -0.006994 -0.006275 -0.005048 -0.003722 -0.002606 -0.001884
-0.003018 -0.0054 -0.006733 -0.006939 -0.006245 -0.005042 -0.003734 -0.002629 -0.001913
-0.002979 -0.005337 -0.006665 -0.006885 -0.006214 -0.005036 -0.003747 -0.002653 -0.001942
-0.002942 -0.005275 -0.006599 -0.006832 -0.006184 -0.005029 -0.003758 -0.002676 -0.00197
-0.002906 -0.005215 -0.006534 -0.006779 -0.006154 -0.005022 -0.00377 -0.002698 -0.001998
-0.00287 -0.005155 -0.006469 -0.006727 -0.006123 -0.005015 -0.00378 -0.002721 -0.002026
-0.002835 -0.005097 -0.006406 -0.006675 -0.006093 -0.005007 -0.003791 -0.002742 -0.002054
-0.0028 -0.00504 -0.006343 -0.006624 -0.006063 -0.005 -0.003801 -0.002764 -0.002082
-0.002767 -0.004984 -0.006282 -0.006573 -0.006033 -0.004991 -0.003811 -0.002785 -0.002109
-0.002734 -0.004928 -0.006222 -0.006524 -0.006003 -0.004983 -0.00382 -0.002806 -0.002136
-0.002701 -0.004874 -0.006162 -0.006474 -0.005973 -0.004974 -0.003829 -0.002827 -0.002163
-0.00267 -0.004821 -0.006104 -0.006425 -0.005943 -0.004966 -0.003838 -0.002847 -0.002189
-0.002639 -0.004769 -0.006046 -0.006377 -0.005913 -0.004957 -0.003846 -0.002867 -0.002216
-0.002608 -0.004718 -0.005989 -0.006329 -0.005884 -0.004947 -0.003854 -0.002887 -0.002242
-0.002578 -0.004668 -0.005933 -0.006282 -0.005854 -0.004938 -0.003861 -0.002906 -0.002267
-0.002549 -0.004618 -0.005878 -0.006236 -0.005825 -0.004928 -0.003869 -0.002925 -0.002293
-0.00252 -0.00457 -0.005824 -0.00619 -0.005796 -0.004919 -0.003876 -0.002944 -0.002318
-0.002492 -0.004522 -0.005771 -0.006144 -0.005767 -0.004909 -0.003882 -0.002962 -0.002343
-0.002464 -0.004475 -0.005718 -0.006099 -0.005738 -0.004898 -0.003889 -0.00298 -0.002368
-0.002437 -0.004429 -0.005667 -0.006054 -0.005709 -0.004888 -0.003895 -0.002998 -0.002393
-0.00241 -0.004384 -0.005616 -0.00601 -0.005681 -0.004878 -0.003901 -0.003015 -0.002417
-0.002384 -0.004339 -0.005565 -0.005967 -0.005652 -0.004867 -0.003906 -0.003033 -0.002441
-0.002358 -0.004295 -0.005516 -0.005924 -0.005624 -0.004857 -0.003911 -0.00305 -0.002464
-0.002333 -0.004252 -0.005467 -0.005882 -0.005596 -0.004846 -0.003916 -0.003066 -0.002488
-0.002308 -0.00421 -0.005419 -0.00584 -0.005568 -0.004835 -0.003921 -0.003082 -0.002511
-0.002284 -0.004169 -0.005372 -0.005798 -0.00554 -0.004824 -0.003926 -0.003098 -0.002534
-0.00226 -0.004128 -0.005325 -0.005757 -0.005513 -0.004814 -0.00393 -0.003114 -0.002556
-0.002236 -0.004087 -0.005279 -0.005716 -0.005486 -0.004802 -0.003934 -0.00313 -0.002579
-0.002213 -0.004048 -0.005234 -0.005676 -0.005459 -0.004791 -0.003938 -0.003145 -0.002601
-0.002191 -0.004009 -0.005189 -0.005637 -0.005432 -0.00478 -0.003942 -0.00316 -0.002622
-0.002168 -0.003971 -0.005145 -0.005598 -0.005405 -0.004769 -0.003945 -0.003174 -0.002644
-0.002146 -0.003933 -0.005102 -0.005559 -0.005378 -0.004758 -0.003948 -0.003189 -0.002665
-0.002125 -0.003896 -0.005059 -0.005521 -0.005352 -0.004746 -0.003951 -0.003203 -0.002686
-0.002104 -0.00386 -0.005017 -0.005483 -0.005326 -0.004735 -0.003954 -0.003217 -0.002707
-0.002083 -0.003824 -0.004976 -0.005445 -0.0053 -0.004724 -0.003957 -0.003231 -0.002727
-0.002062 -0.003788 -0.004935 -0.005409 -0.005274 -0.004712 -0.003959 -0.003244 -0.002747
-0.002042 -0.003753 -0.004894 -0.005372 -0.005248 -0.004701 -0.003961 -0.003257 -0.002767
-0.002022 -0.003719 -0.004855 -0.005336 -0.005223 -0.004689 -0.003963 -0.00327 -0.002787
-0.002003 -0.003686 -0.004815 -0.0053 -0.005198 -0.004678 -0.003965 -0.003283 -0.002806
-0.001983 -0.003652 -0.004777 -0.005265 -0.005172 -0.004666 -0.003967 -0.003295 -0.002825
-0.001965 -0.00362 -0.004738 -0.00523 -0.005148 -0.004654 -0.003969 -0.003307 -0.002844
-0.001946 -0.003588 -0.004701 -0.005196 -0.005123 -0.004643 -0.00397 -0.003319 -0.002862
-0.001928 -0.003556 -0.004664 -0.005162 -0.005099 -0.004631 -0.003971 -0.003331 -0.002881
-0.00191 -0.003525 -0.004627 -0.005128 -0.005074 -0.00462 -0.003973 -0.003342 -0.002899
-0.001892 -0.003494 -0.004591 -0.005095 -0.00505 -0.004608 -0.003974 -0.003353 -0.002916
-0.001875 -0.003464 -0.004555 -0.005062 -0.005026 -0.004597 -0.003974 -0.003364 -0.002934
-0.001858 -0.003434 -0.00452 -0.005029 -0.005003 -0.004585 -0.003975 -0.003375 -0.002951
-0.001841 -0.003404 -0.004486 -0.004997 -0.004979 -0.004574 -0.003976 -0.003386 -0.002968
-0.001824 -0.003375 -0.004451 -0.004965 -0.004956 -0.004562 -0.003976 -0.003396 -0.002985
-0.001808 -0.003347 -0.004418 -0.004934 -0.004933 -0.004551 -0.003977 -0.003406 -0.003001
-0.001792 -0.003319 -0.004384 -0.004903 -0.00491 -0.004539 -0.003977 -0.003416 -0.003017
-0.001776 -0.003291 -0.004352 -0.004872 -0.004887 -0.004528 -0.003977 -0.003426 -0.003033
-0.00176 -0.003264 -0.004319 -0.004842 -0.004865 -0.004516 -0.003977 -0.003435 -0.003049
-0.001745 -0.003237 -0.004287 -0.004812 -0.004843 -0.004505 -0.003977 -0.003445 -0.003065
-0.00173 -0.00321 -0.004256 -0.004782 -0.00482 -0.004493 -0.003976 -0.003454 -0.00308
-0.001715 -0.003184 -0.004225 -0.004753 -0.004799 -0.004482 -0.003976 -0.003463 -0.003095
-0.0017 -0.003159 -0.004194 -0.004724 -0.004777 -0.004471 -0.003976 -0.003471 -0.003109
-0.001685 -0.003133 -0.004164 -0.004695 -0.004755 -0.004459 -0.003975 -0.00348 -0.003124
-0.001671 -0.003108 -0.004134 -0.004667 -0.004734 -0.004448 -0.003974 -0.003488 -0.003138
-0.001657 -0.003083 -0.004104 -0.004639 -0.004713 -0.004437 -0.003974 -0.003496 -0.003152
-0.001643 -0.003059 -0.004075 -0.004611 -0.004692 -0.004425 -0.003973 -0.003504 -0.003166
-0.00163 -0.003035 -0.004046 -0.004584 -0.004671 -0.004414 -0.003972 -0.003512 -0.003179
-0.001616 -0.003011 -0.004018 -0.004557 -0.00465 -0.004403 -0.003971 -0.00352 -0.003193
-0.001603 -0.002988 -0.00399 -0.00453 -0.00463 -0.004392 -0.00397 -0.003527 -0.003206
-0.00159 -0.002965 -0.003962 -0.004504 -0.00461 -0.004381 -0.003969 -0.003535 -0.003219
-0.001577 -0.002942 -0.003935 -0.004478 -0.00459 -0.00437 -0.003967 -0.003542 -0.003231
-0.001564 -0.00292 -0.003908 -0.004452 -0.00457 -0.004359 -0.003966 -0.003549 -0.003244
-0.001552 -0.002898 -0.003881 -0.004426 -0.00455 -0.004348 -0.003964 -0.003555 -0.003256
-0.001539 -0.002876 -0.003855 -0.004401 -0.00453 -0.004337 -0.003963 -0.003562 -0.003268
-0.001527 -0.002855 -0.003829 -0.004376 -0.004511 -0.004326 -0.003961 -0.003568 -0.00328
-0.001515 -0.002833 -0.003803 -0.004351 -0.004492 -0.004315 -0.00396 -0.003575 -0.003291
-0.001503 -0.002813 -0.003778 -0.004327 -0.004473 -0.004305 -0.003958 -0.003581 -0.003303
-0.001492 -0.002792 -0.003753 -0.004303 -0.004454 -0.004294 -0.003956 -0.003587 -0.003314
-0.00148 -0.002772 -0.003728 -0.004279 -0.004435 -0.004283 -0.003954 -0.003592 -0.003325
-0.001469 -0.002752 -0.003704 -0.004255 -0.004417 -0.004273 -0.003952 -0.003598 -0.003335
-0.001458 -0.002732 -0.00368 -0.004232 -0.004398 -0.004262 -0.00395 -0.003604 -0.003346
-0.001447 -0.002712 -0.003656 -0.004209 -0.00438 -0.004252 -0.003948 -0.003609 -0.003356
-0.001436 -0.002693 -0.003633 -0.004186 -0.004362 -0.004241 -0.003946 -0.003614 -0.003366
-0.001425 -0.002674 -0.003609 -0.004163 -0.004344 -0.004231 -0.003944 -0.003619 -0.003376
-0.001415 -0.002655 -0.003587 -0.004141 -0.004327 -0.00422 -0.003941 -0.003624 -0.003386
-0.001404 -0.002637 -0.003564 -0.004119 -0.004309 -0.00421 -0.003939 -0.003629 -0.003395
-0.001394 -0.002618 -0.003542 -0.004097 -0.004292 -0.0042 -0.003937 -0.003633 -0.003405
-0.001384 -0.0026 -0.003519 -0.004075 -0.004275 -0.004189 -0.003934 -0.003638 -0.003414
-0.001374 -0.002582 -0.003498 -0.004054 -0.004257 -0.004179 -0.003932 -0.003642 -0.003423
-0.001364 -0.002565 -0.003476 -0.004033 -0.004241 -0.004169 -0.003929 -0.003646 -0.003432
-0.001354 -0.002547 -0.003455 -0.004012 -0.004224 -0.004159 -0.003927 -0.003651 -0.00344
-0.001344 -0.00253 -0.003434 -0.003991 -0.004207 -0.004149 -0.003924 -0.003654 -0.003449
-0.001335 -0.002513 -0.003413 -0.003971 -0.004191 -0.004139 -0.003921 -0.003658 -0.003457

You might also like