Nonlinear Seidel Iteration.
To solve the nonlinear fixed-point system X= G(X), given one initial approximation (Po,qo).
Solving for nonlinear system
y2 - 4x2 + 1 = 0 eq. 1 (hyperbola)
x2 + y2 - 2x = 3 eq. 2 (circle)
Rearranging the equations
x = (y2 - 4x2 + 8x + 1)/8 (adding 8x to both sides of the equation)
y = (x2 + y2 - 2x + 4y- 3)/4 (adding 4y to both sides of the equation)
Using Seidel Iteration Formulas.
pk+1 = g1 (pk, qk) qk = g2 (pk+1, qk)
Input:
p0 = 1 q0 = 1
Iteration spreadsheet.
k p q
0 1 1
8
1 0.75 0.265625
2 0.6025695801 -0.677248105
3 0.6033576552 -1.523250569 6
4 0.836373962 -1.936484125
5 1.0803596056 -1.997377017
4
6 1.1204605353 -1.996370595
7 1.1159315737 -1.996636674
8 1.1165996863 -1.9965983 2
9 1.1165028532 -1.996603878
10 1.1165169235 -1.996603068 0
11 1.1165148798 -1.996603186 -3 -2 -1 0 1 2 3 4
12 1.1165151766 -1.996603169
13 1.1165151335 -1.996603171 -2
14 1.1165151398 -1.996603171
15 1.1165151389 -1.996603171 -4
16 1.116515139 -1.996603171
17 1.116515139 -1.996603171
-6
18 1.116515139 -1.996603171
19 1.116515139 -1.996603171
20 1.116515139 -1.996603171 -8
X = 1.1165139 Y = -1.99660317
-2 3.872983
-1.925 3.717862
-1.9 3.666061
-1.85 3.562303
-1.8 3.458323
-1.75 3.354102
-1.7 3.249615
-1.65 3.144837
-1.6 3.039737
-1.55 2.93428
-1.5 2.828427
-1.45 2.722132
-1.4 2.615339
-1.35 2.507987
-1.3 2.4
-1.25 2.291288
-1.2 2.181742
-1.15 2.071232
-1.1 1.959592
-1.05 1.846619
-1 1.732051
-0.95 1.615549
-0.9 1.496663
-0.85 1.374773
-0.8 1.249
-0.75 1.118034
-0.7 0.979796
-0.65 0.830662
-0.6 0.663325
-0.55 0.458258
-0.5 0
-0.45
-0.4
-0.35
-0.3
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.45
0.5 0
0.55 0.458258
0.6 0.663325
0.65 0.830662
0.7 0.979796
0.75 1.118034
0.8 1.249
0.85 1.374773
0.9 1.496663
0.95 1.615549
1 1.732051
1.05 1.846619
1.1 1.959592
1.15 2.071232
1.2 2.181742
1.25 2.291288
1.3 2.4
1.35 2.507987
1.4 2.615339
1.45 2.722132
1.5 2.828427
1.55 2.93428
1.6 3.039737
1.65 3.144837
1.7 3.249615
1.75 3.354102
1.8 3.458323
1.85 3.562303
1.9 3.666061
1.95 3.769615
2 3.872983
2.05 3.976179
2.1 4.079216
2.15 4.182105
2.2 4.284857
2.25 4.387482
2.3 4.489989
2.35 4.592385
2.4 4.694678
2.45 4.796874
2.5 4.898979
2.55 5.001
2.6 5.10294
2.65 5.204805
2.7 5.3066
2.75 5.408327
2.8 5.509991
2.85 5.611595
2.9 5.713143
2.95 5.814637
3 5.91608
-3.872983
-3.717862
-3.666061
-3.562303
-3.458323
-3.354102
-3.249615
-3.144837
-3.039737
-2.93428
-2.828427
-2.722132
-2.615339
-2.507987
-2.4
-2.291288
-2.181742
-2.071232
-1.959592
-1.846619
-1.732051 0 0
-1.615549 0.44441 -0.44441
-1.496663 0.6245 -0.6245
-1.374773 0.759934 -0.759934
-1.249 0.87178 -0.87178
-1.118034 0.968246 -0.968246
-0.979796 1.053565 -1.053565
-0.830662 1.130265 -1.130265
-0.663325 1.2 -1.2
-0.458258 1.263922 -1.263922
0 1.322876 -1.322876
1.377498 -1.377498
1.428286 -1.428286
1.475635 -1.475635
1.519868 -1.519868
1.561249 -1.561249
1.6 -1.6
1.636307 -1.636307
1.670329 -1.670329
1.702204 -1.702204
1.732051 -1.732051
1.759972 -1.759972
1.786057 -1.786057
1.810387 -1.810387
1.83303 -1.83303
1.85405 -1.85405
1.873499 -1.873499
1.891428 -1.891428
1.907878 -1.907878
1.922888 -1.922888
0 1.936492 -1.936492
-0.458258 1.948718 -1.948718
-0.663325 1.959592 -1.959592
-0.830662 1.969137 -1.969137
-0.979796 1.977372 -1.977372
-1.118034 1.984313 -1.984313
-1.249 1.989975 -1.989975
-1.374773 1.994367 -1.994367
-1.496663 1.997498 -1.997498
-1.615549 1.999375 -1.999375
-1.732051 2 -2
-1.846619 1.999375 -1.999375
-1.959592 1.997498 -1.997498
-2.071232 1.994367 -1.994367
-2.181742 1.989975 -1.989975
-2.291288 1.984313 -1.984313
-2.4 1.977372 -1.977372
-2.507987 1.969137 -1.969137
-2.615339 1.959592 -1.959592
-2.722132 1.948718 -1.948718
-2.828427 1.936492 -1.936492
-2.93428 1.922888 -1.922888
-3.039737 1.907878 -1.907878
-3.144837 1.891428 -1.891428
-3.249615 1.873499 -1.873499
-3.354102 1.85405 -1.85405
-3.458323 1.83303 -1.83303
-3.562303 1.810387 -1.810387
-3.666061 1.786057 -1.786057
-3.769615 1.759972 -1.759972
-3.872983 1.732051 -1.732051
-3.976179 1.702204 -1.702204
-4.079216 1.670329 -1.670329
-4.182105 1.636307 -1.636307
-4.284857 1.6 -1.6
-4.387482 1.561249 -1.561249
-4.489989 1.519868 -1.519868
-4.592385 1.475635 -1.475635
-4.694678 1.428286 -1.428286
-4.796874 1.377498 -1.377498
-4.898979 1.322876 -1.322876
-5.001 1.263922 -1.263922
-5.10294 1.2 -1.2
-5.204805 1.130265 -1.130265
-5.3066 1.053565 -1.053565
-5.408327 0.968246 -0.968246
-5.509991 0.87178 -0.87178
-5.611595 0.759934 -0.759934
-5.713143 0.6245 -0.6245
-5.814637 0.44441 -0.44441
-5.91608 0 0
Po= 1
qo= 1
k pk qk
0 1 1
1 0.750000000000 0.265625000000
2 0.602569580078 -0.677248105174
3 0.603357655156 -1.523250568753
4 0.836373962039 -1.936484124877
5 1.080359605588 -1.997377016849
6 1.120460535325 -1.996370594847
7 1.115931573711 -1.996636674409
8 1.116599686308 -1.996598300298
9 1.116502853171 -1.996603878411
10 1.116516923512 -1.996603068223
11 1.116514879773 -1.996603185912
12 1.116515176645 -1.996603168816
13 1.116515133522 -1.996603171300
14 1.116515139786 -1.996603170939
15 1.116515138876 -1.996603170991
15 1.116515139008 -1.996603170984 STOP
15 1.116515138989 -1.996603170985 STOP
15 1.116515138992 -1.996603170985 STOP
15 1.116515138991 -1.996603170985
