f(x) = x5 + 5x -1 Z0 0

Z1 0.5 First trial

z Zk f(zk)
0 0 1
1 1 -2 0.5
2 0.5 -0.375 0 0.25
3 0.25 0.390625 0.375 0.625 0
4 0.625 -0.787109375 0.4375 0 0
5 0.4375 -0.1740722656 0 0.34375 0
6 0.34375 0.1167907715 0.390625 0.484375 0
7 0.671875 -0.9432868958 0.5078125 0 0
8 0.507813 -0.4004206657 0 0.425781 0
9 0.425781 -0.1369521022 0 0 0.384766
10 0.384766 -0.0086579546 0 0 0
11 0.364258 0.054447961 0.3745117 0.39502 0.436035
12 0.682129 -0.9774630278 0.5231934 0 0
13 0.523193 -0.4506349075 0 0.443726 0
z f(z) 14 0.443726 -0.1938697695 0 0 0.403992
0 -1 15 0.403992 -0.0684667849 0 0 0
0.2 0.00032 16 0.384125 -0.0066748579 0 0 0
0.2 0.00032
0.3 0.50243
0.4 1.01024
0.5 1.53125
0.6 2.07776
0.7 2.66807
0.8 3.32768
0.9 4.09049
1 5
1.1 6.11051
 Zk+1

0.3125
0 0.21875
0.421875 0.671875 0
0.460938 0 0 0.335938
0 0.378906 0 0 0.253906
0 0 0.337891 0 0 0.212891
0.364258 0 0 0.317383 0 0 0.192383
0.518066 0 0.400879 0.494629 0 0.432129 0.682129 0
0 0 0.512939 0 0 0.466064 0 0 0.341064
0 0 0 0.433472 0 0 0.386597 0 0 0.261597
0 0 0 0 0.393738 0 0 0.346863 0 0
0.384125 0 0 0 0 0.373871 0 0 0.326996 0
0 0.374191 0 0 0 0 0.363937 0 0 0.317062
0.221863
0 0.201996
0 0 0.192062
 f(x) = x5 + 5x -1 Z0 0
Z1 0.5 First trial

z Zk f(zk)
0 -1 433
1 0 -9 -0.020362
2 -0.02036 -4.4963809278 0 -0.03043
3 -0.03043 -2.266891367 0 0 -0.03548
4 -0.03548 -1.1480711912 0 0 0
5 -0.03803 -0.5827707397 0 0 0
6 -0.03932 -0.2961607508 0 0 0
7 -0.03998 -0.1505952314 0 0 0
8 -0.04031 -0.0765991444 0 0 0
9 -0.04048 -0.0389674676 0 0 0
10 -0.04057 -0.0198250271 0 0 0
11 -0.04061 -0.0100865434 0 0 0
12 -0.04064 -0.0051319162 0 0 0
13 -0.04065 -0.0026110858 0 0 0
z f(z) 14 -0.04065 -0.0013285104 0 0 0
0 -1 15 -0.04066 -0.0006759428 0 0 0
0.2 0.00032 16 -0.04066 -0.0003439184 0 0 0
0.2 0.00032 17 -0.04066 -0.0001749852 0 0 0
0.3 0.50243 18 -0.04066 -8.903223E-05 0 0 0
0.4 1.01024 19 -0.04066 -4.529949E-05 0 0 0
0.5 1.53125
0.6 2.07776
0.7 2.66807
0.8 3.32768
0.9 4.09049
1 5
1.1 6.11051
 Z2 = f1Z0 + f0Z1/(f1 - f0)

Zk+1

-0.03803
0 -0.039323
0 0 -0.03998
0 0 0 -0.040314
0 0 0 0 -0.040484
0 0 0 0 0 -0.04057
0 0 0 0 0 0 -0.040614
0 0 0 0 0 0 0 -0.040636
0 0 0 0 0 0 0 0 -0.040648
0 0 0 0 0 0 0 0 0 -0.040653
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
-0.040656
0 -0.040658
0 0 -0.040658
0 0 0 -0.040659
0 0 0 0 -0.040659
0 0 0 0 0
 Zk+1 = Zk - fk/f[Zk-1,Zk]
f[ Zk-1,Zk] = (fk - fk-1)/(Zk - Zk-1)

Z Zk f(Zk) f[ Zk-1,Zk] Zk+1

0 6 -6
1 7 22 28 6.2142857143
2 6.214285714 -1.71975218659 30.18877551 6.2712523238
3 6.271252324 -0.43529960729 22.547464031 6.2905582469
4 6.290558247 0.014395415692 23.293111654 6.2899402352
5 6.289940235 -0.0001137813 23.477220301 6.2899450817
6 6.289945082 -2.93317584E-08 23.47116809 6.2899450829
7 6.289945083 0 23.471175079 6.2899450829
8 6.289945083 0 #DIV/0! #DIV/0!
9 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
10 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
11 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
12 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
13 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
14 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
15 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
16 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
17 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
18 #DIV/0! #DIV/0! #DIV/0! #DIV/0!
(fk - fk-1)/(Zk - Zk-1)

x f
-15.5 87030.94
-15 77292
-14.5 68392.19
-14 60282
-13.5 52913.44
-13 46240
-12.5 40216.69
-12 34800
-11.5 29947.94
-11 25620
-10.5 21777.19
-10 18382
-9.5 15398.44
-9 12792
-8.5 10529.69
-8 8580
-7.5 6912.938
-7 5500
-6.5 4314.188
-6 3330
-5.5 2523.438
-5 1872
-4.5 1354.688
-4 952
-3.5 645.9375
-3 420
-2.5 259.1875
-2 150
-1.5 80.4375
-1 40
-0.5 19.6875
0 12
0.5 10.9375
1 12
1.5 12.1875
2 10
2.5 5.4375
3 0
 3.5 -3.3125
4 0
4.5 15.9375
5 52
5.5 117.1875
6 222
6.5 378.4375
7 600
7.5 901.6875
8 1300
8.5 1812.938
9 2460
9.5 3262.188
10 4242
10.5 5423.438
11 6832
11.5 8494.688
12 10440
12.5 12697.94
13 15300
13.5 18279.19
14 21670
14.5 25508.44
15 29832
15.5 34679.69
16 40092
16.5 46110.94
17 52780
17.5 60144.19
18 68250
18.5 77145.44
19 86880
19.5 97504.69
20 109072
20.5 121635.9
21 135252
21.5 149977.2
22 165870
22.5 182990.4
23 201400
23.5 221161.7
24 242340
24.5 265000.9
25 289212
25.5 315042.2
 26 342562
26.5 371843.4
27 402960
27.5 435986.7
28 471000
28.5 508077.9
29 547300
29.5 588747.2
30 632502
30.5 678648.4
31 727272
31.5 778459.7
32 832300
32.5 888882.9
33 948300
33.5 1010644
34 1076010
34.5 1144493
35 1216192
35.5 1291205
36 1369632
36.5 1451576
37 1537140
37.5 1626429
38 1719550
38.5 1816610
39 1917720
39.5 2022990
40 2132532
40.5 2246461
41 2364892
41.5 2487942
42 2615730
42.5 2748375
43 2886000
43.5 3028727
44 3176680
44.5 3329986
45 3488772
45.5 3653167
46 3823302
46.5 3999308
47 4181320
47.5 4369472
48 4563900
48.5 4764743
49 4972140
49.5 5186232
50 5407162
50.5 5635073
51 5870112
51.5 6112425
 f(x) = x5 + 5x -1
Zk+1 = Zk - f(Zk)/f'(Zk)
f'(x) = 5x4 + 5
z Zk f(zk) f'(zk) Zk+1
0 3.2 -1.7984 -7.768 2.968486 1
1 2.968486 0.318956 -10.23734 2.999642 2
2 2.999642 0.003578 -10.00286 3 3
3 3 5.114E-07 -10 3 4
4 3 4.263E-14 -10 3 5
5 3 0 -10 3 6
6 3 0 -10 3 7
7 3 0 -10 3 8
8 3 0 -10 3 9
9 3 0 -10 3 10
10 3 0 -10 3 11
11 3 0 -10 3 12
12 3 0 -10 3 13
13 3 0 -10 3 14
cos(#/3) cos(#/7) sum b^2 4ac b^2-4ac
0.707107 0.92388 2.660116 a 2.660116 256 170.2474 85.75255
b 16
c 16 9.260267

1.266812
Xk yk f(xk,yk) g(xk,yk) fx(xk,yk) fy(xk,yk) gx(xk,yk) gy(xk,yk) DXk
1 2 -79 -17 -96 -108 16 -32 -0.144167
0.855833 1.396667 -17.56786 -5.081382 -43.87889 -75.42 9.562511 -14.6263 0.092653
0.948486 1.109828 -3.130098 -1.107265 -22.46088 -59.93071 8.421253 -7.486959 0.063794
1.01228 1.033691 -0.055044 -0.076606 -16.5534 -55.81929 8.371078 -5.517799 0.007111
1.019392 1.030596 0.003702 -5.88E-05 -16.20752 -55.65216 8.404644 -5.402508 4.191E-05
1.019434 1.03065 -5.88E-05 -1.56E-09 -16.20952 -55.65509 8.405433 -5.403174 -5.72E-07
1.019433 1.030649 9.626E-07 -1.85E-12 -16.20949 -55.65505 8.405421 -5.403162 9.365E-09
1.019433 1.030649 -1.58E-08 0 -16.20949 -55.65505 8.405421 -5.403162 -1.53E-10
1.019433 1.030649 2.579E-10 0 -16.20949 -55.65505 8.405421 -5.403162 2.509E-12
1.019433 1.030649 -4.23E-12 0 -16.20949 -55.65505 8.405421 -5.403162 -4.11E-14
1.019433 1.030649 0 0 -16.20949 -55.65505 8.405421 -5.403162 0
1.019433 1.030649 0 0 -16.20949 -55.65505 8.405421 -5.403162 0
1.019433 1.030649 0 0 -16.20949 -55.65505 8.405421 -5.403162 0
1.019433 1.030649 0 0 -16.20949 -55.65505 8.405421 -5.403162 0
1.019433 1.030649 0 0 -16.20949 -55.65505 8.405421 -5.403162 0
1.019433 1.030649 0 0 -16.20949 -55.65505 8.405421 -5.403162 0
1.019433 1.030649 0 0 -16.20949 -55.65505 8.405421 -5.403162 0
Dyk Xk+1 yk+1
-0.603333 0.855833 1.396667
-0.286839 0.948486 1.109828
-0.076137 1.01228 1.033691
-0.003095 1.019392 1.030596
5.432E-05 1.019434 1.03065
-8.9E-07 1.019433 1.030649
1.457E-08 1.019433 1.030649
-2.38E-10 1.019433 1.030649
3.903E-12 1.019433 1.030649
-6.4E-14 1.019433 1.030649
0 1.019433 1.030649
0 1.019433 1.030649
0 1.019433 1.030649
0 1.019433 1.030649
0 1.019433 1.030649
0 1.019433 1.030649
0 1.019433 1.030649
X f(x)
I = h*[f(x0) + 2*åf(xi) + f(xn)]/2
2 5.54518
2.2 6.93843
2.4 8.8445 Points f(x0) 2*åf(xi) f(xn) n
2.6 9.93732 2 5.54518 0 6.93843 1
2.6 10.46675 3 5.54518 13.87686 8.8445 2
1.4 3.3886 4 5.54518 31.56586 9.93732 3
1.6 3.81 5 5.54518 51.4405 10.46675 4
1.8 4.3511 6 5.54518 72.374 3.3886 5
2 5 7 5.54518 79.1512 3.81 6
2.2 5.7491 8 5.54518 86.7712 4.3511 7
10 22026.47 9 5.54518 95.4734 5 8
11 59874.14 10 5.54518 105.4734 5.7491 9
12 162754.8 11 5.54518 116.9716 22026.47 10
13 442413.4 12 5.54518 44169.9 59874.14 11
14 1202604 13 5.54518 163918.2 162754.8 12
15 3269017 14 5.54518 244713.9 442413.4 13
15 5.54518 1374255 1202604 14
16 5.54518 3779463 3269017 15
 b 2.2
a 2

h I
0.2 1.248361
0.1 1.413327
0.066667 1.5682786667
0.05 1.68631075
0.04 1.6261556
0.033333 1.4751063333
0.028571 1.380964
0.025 1.32523225
0.022222 1.2974186667
0.02 221.48982575
0.018182 945.90536441
0.016667 2722.3210268
0.015385 5285.6370917
0.014286 18406.174163
0.013333 46989.90693
 b
X f(x) a
2 5.54518
2.2 6.93843 Points f(x0) 4*åf(xi) 2*åf(xi) f(xn)
2.4 8.8445 3 5.54518 27.75372 0 8.8445
2.6 9.93732 5 5.54518 67.503 17.689 10.46675
2.6 10.46675 7 5.54518 81.0574 38.6225 3.81
1.4 3.3886 9 5.54518 98.4618 46.2425 5
1.6 3.81 11 5.54518 121.4582 56.2425 22026.47
1.8 4.3511 13 5.54518 239618 44109.17 162754.8
2 5 15 5.54518 2009272 369618.8 1202604
2.2 5.7491
10 22026.47
11 59874.14
12 162754.8
13 442413.4
14 1202604
15 3269017
 2.2
2

n h I
1 0.1 1.40478
2 0.05 1.686732
3 0.033333 1.433723
4 0.025 1.293746
5 0.02 148.0647
6 0.016667 2480.486
7 0.014286 17054.76
X f(x) Points f(x0) 3*åf(xi) 3*åf(xi) 2*åf(xi)
2 5.54518 4 5.54518 20.81529 26.5335 0
2.2 6.93843 7 5.54518 30.98109 57.93375 19.87464
2.4 8.8445 10 5.54518 44.03439 72.93375 27.49464
2.6 9.93732 13 5.54518 179666.5 66152.33 38.99284
2.6 10.46675
1.4 3.3886
1.6 3.81
1.8 4.3511
2 5
2.2 5.7491
10 22026.47
11 59874.14
12 162754.8
13 442413.4
14 1202604
15 3269017
 b 2.2
f(xn) n h I a 2
9.93732 1 0.066667 1.570782
3.81 2 0.033333 1.476808
5.7491 3 0.022222 1.297976
162754.8 4 0.016667 2553.863
Given
4 -1 0 -1 0 0 0 0 0
-1 4 -1 0 -1 0 0 0 0
0 -1 4 0 0 -1 0 0 0
-1 0 0 4 -1 0 -1 0 0
0 -1 0 -1 4 -1 0 -1 0
0 0 -1 0 -1 4 0 0 -1
0 0 0 -1 0 0 4 -1 0
0 0 0 0 -1 0 -1 4 0
0 0 0 0 0 -1 0 -1 4

1st iteration 2 3 4 5 6 7 8 9
X1 6.25 9.375 12.89063 15.16113 16.72668 17.54637 17.94017 18.12646 18.21433
X2 12.5 24.21875 29.83398 33.44116 35.20012 36.01644 36.39835 36.57789 36.66244
X3 37.5 46.67969 52.00195 54.14429 55.03607 55.43888 55.6263 55.71431 55.75575
X4 0 2.34375 5.810547 8.465576 9.985352 10.74424 11.1075 11.27944 11.36057
0 9.765625 16.60156 19.9295 21.48333 22.21435 22.55879 22.72112 22.79764
12.5 28.17383 33.13599 34.94415 35.73909 36.10685 36.27935 36.36056 36.39882
0 0.585938 2.099609 3.285217 3.947258 4.275471 4.432489 4.506814 4.541888
0 2.587891 4.675293 5.80368 6.357646 6.622456 6.74782 6.806985 6.834881
6.25 13.94043 15.70282 16.43696 16.77418 16.93233 17.00679 17.04189 17.05842
 Iterative equations
25
50 X1 = 1/a11*(C1 - a12*X2 - a13*X3 - a14*X4)
150 X2 = 1/a22*(C2 - a21*X1' - a23*X3 - a24*X4)
0 X3 = 1/a33*(C3 - a31*X1' - a32*X2' - a34*X4)
0 X4 = 1/a44*(C4 - a41*X1' - a42*X2' - a43*X3')
50
0
0
25

10 11 12 13 14 15 16 17 18 19
18.25575 18.27528 18.28448 18.28882 18.29086 18.29182 18.29228 18.29249 18.29259 18.29264
36.70228 36.72106 36.72991 36.73409 36.73605 36.73698 36.73742 36.73762 36.73772 36.73776
55.77528 55.78448 55.78882 55.79086 55.79182 55.79228 55.79249 55.79259 55.79264 55.79266
11.39882 11.41685 11.42535 11.42936 11.43124 11.43213 11.43255 11.43275 11.43284 11.43289
22.8337 22.8507 22.85871 22.86249 22.86427 22.86511 22.8655 22.86569 22.86578 22.86582
36.41685 36.42535 36.42936 36.43124 36.43213 36.43255 36.43275 36.43284 36.43289 36.43291
4.558425 4.566221 4.569895 4.571627 4.572443 4.572828 4.573009 4.573095 4.573135 4.573154
6.848032 6.85423 6.857151 6.858528 6.859177 6.859483 6.859628 6.859696 6.859728 6.859743
17.06622 17.06989 17.07163 17.07244 17.07283 17.07301 17.07309 17.07313 17.07315 17.07316
 20 21 22 23 24 25
18.29266 18.29267 18.29268 18.29268 18.29268 18.29268 18.29268 18.29268 18.29268 18.29268
36.73779 36.7378 36.7378 36.7378 36.7378 36.7378 36.7378 36.7378 36.7378 36.7378
55.79267 55.79268 55.79268 55.79268 55.79268 55.79268 55.79268 55.79268 55.79268 55.79268
11.43291 11.43292 11.43292 11.43292 11.43293 11.43293 11.43293 11.43293 11.43293 11.43293
22.86584 22.86585 22.86585 22.86585 22.86585 22.86585 22.86585 22.86585 22.86585 22.86585
36.43292 36.43292 36.43292 36.43293 36.43293 36.43293 36.43293 36.43293 36.43293 36.43293
4.573163 4.573167 4.573169 4.57317 4.57317 4.573171 4.573171 4.573171 4.573171 4.573171
6.85975 6.859753 6.859755 6.859755 6.859756 6.859756 6.859756 6.859756 6.859756 6.859756
17.07317 17.07317 17.07317 17.07317 17.07317 17.07317 17.07317 17.07317 17.07317 17.07317
18.29268
36.7378
55.79268
11.43293
22.86585
36.43293
4.573171
6.859756
17.07317
 assumed l 0.5 p 7.8 for l = 0.5
Dx 0.25 k 0.13
Dt 0.206 c 0.11
Values of x
0.25 0.5 0.75
time steps t comput anal comput anal comput anal
0 0 38.26834 38.26834 70.71068 70.71068 92.38795 92.38795
1 0.206 35.35534 35.43217 65.32815 65.47012 85.35534 85.54083
2 0.412 32.66407 32.8062 60.35534 60.61795 55.76106 79.20116
3 0.618 30.17767 30.37484 44.21257 56.12538 51.5165 73.33135
4 0.824 22.10628 28.12368 40.84709 51.96578 36.04655 67.89656
5 1.03 20.42354 26.03935 29.07642 48.11445 33.30267 62.86456
6 1.236 14.53821 24.10951 26.86311 44.54856 23.54985 58.2055
7 1.442 13.43155 22.32268 19.04403 41.24694 21.75722 53.89173
8 1.648 9.522014 20.66829 17.59439 38.19002 15.40948 49.89766
9 1.854 8.797193 19.13651 12.46574 35.35965 14.2365 46.19961
10 2.06 6.232872 17.71825 11.51685 32.73905 10.08524 42.77563
11 2.266 5.758423 16.4051 8.159057 30.31267 9.317548 39.60541
12 2.472 4.079528 15.18927 7.537985 28.06612 6.600839 36.67015
13 2.678 3.768993 14.06355 5.340183 25.98606 6.09838 33.95242
14 2.884 2.670092 13.02127 4.933686 24.06016 4.320301 31.43612
15 3.09 2.466843 12.05622 3.495197 22.277 3.991438 29.1063
for l = 0.5

1 0.25 0.5 0.75 1

comput anal
100 100
46.19398 92.58873
42.67767 85.72672
27.88053 79.37328
25.75825 73.49071
18.02327 68.04411
16.65133 63.00118
11.77492 58.33199
10.87861 54.00884
7.704738 50.0061
7.118249 46.30001
5.04262 42.86859
4.658774 39.69148
3.300419 36.74984
3.04919 34.02621
2.160151 31.50443
Given x y c
r 0 0 100 1 Xr
s 2 0 200 [M] = 1 Xs
t 0 1 300 1 xt

a1
{a} = a2
a3

Cr
{a} = [M-1] {c} {c} = Cs
Ct
yr 1 0 0
ys [M] = 1 2 0
yt 1 0 1

1 0 0 100
[M-1] -0.5 0.5 0 C 200
-1 0 1 300

a 100
50
200

V(x,y) = a1 + a2X + a3y = (1 x y)[M-1]{C}

N = (1 x y)[M-1]
v(x,y) = 100 + 50x + 200y