Scribd
Upload
168 views

Method of Successive Substitution (MOSS) : I FX X FX

The document provides examples of numerical methods for finding roots of equations, including the method of successive substitution, incremental search method, half-interval method, regula-falsi method, secant method, and Newton-Raphson method. For each method, the document shows the formulation, initialization, iterative process, and convergence criteria using sample problems.

Uploaded by

Naci John Trance
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
168 views

Method of Successive Substitution (MOSS) : I FX X FX

The document provides examples of numerical methods for finding roots of equations, including the method of successive substitution, incremental search method, half-interval method, regula-falsi method, secant method, and Newton-Raphson method. For each method, the document shows the formulation, initialization, iterative process, and convergence criteria using sample problems.

Uploaded by

Naci John Trance
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

Find the other root of x 12 e 2 0, d 0.001.

Problem:

Method of Successive Substitution (MOSS)


i

xi

F(xi)

x0

F(xo)

IF Remarks = continue
F(xo)

F(x1)

x 12 e 2 0

F xo xo
F xo

Remarks
IF error < d, "stop"
else, "continue"

F x1 x1
F x1

IF error < d, "stop"


else, "continue"

x 12 e 2
i

x 12 e 2 0

Error (di)

xi

F(xi)

Error (di)

Remarks

0.824360635

0.213061319

continue

1
2
3
4
5
6

0.824360635
0.755053355
0.729336179
0.720017983
0.716671151
0.715472865

0.755053355
0.729336179
0.720017983
0.716671151
0.715472865
0.715044323

0.091791235
0.035261074
0.012941615
0.004669968
0.001674817
0.000599322

continue
continue
continue
continue
continue
stop

F(xi)

Error (di)

Remarks

x ln 4 x 2
i

xi

50

9.210340372

4.428681024

continue

1
2
3
4
5
6
7
8
9

9.210340372
5.826947974
4.91128108
4.569364002
4.425042415
4.360854089
4.331630221
4.318182292
4.311963457

5.826947974
4.91128108
4.569364002
4.425042415
4.360854089
4.331630221
4.318182292
4.311963457
4.309081079

0.580645719
0.186441557
0.074828155
0.032614735
0.01471921
0.006746621
0.003114257
0.001442228
0.000668908

continue
continue
continue
continue
continue
continue
continue
continue
stop

Incremental Search Method (ISM)


xo = a, f xo 0
i

xi

f(xi)

Error

depend

x0

f(xo)

depend

Remarks
IF f(xo) > 0, "revert"
IF error < d, "stop"
IF error > d, "continue"

xo =
error =

0
0.00001

x =

f(xo) =

-0.5

OK

remarks

xi

f(xi)

|f(xi)|
error

0
1
1
2
3
4
5
6
7

1
1
0.1
0.1
0.1
0.1
0.1
0.1
0.1

0
1
0.1
0.2
0.3
0.4
0.5
0.6
0.7

-0.5
0.175639365
-0.425635548
-0.352585459
-0.280917121
-0.210701379
-0.142012708
-0.074929404
-0.009533774

0.5
0.175639365
0.425635548
0.352585459
0.280917121
0.210701379
0.142012708
0.074929404
0.009533774

(-) continue
revert
(-) continue
(-) continue
(-) continue
(-) continue
(-) continue
(-) continue
(-) continue

8
8
9
9
10
11
12

0.1
0.01
0.01
0.001
0.001
0.001
0.001

0.8
0.71
0.72
0.711
0.712
0.713
0.714

0.054087651
-0.003090327
0.003335293
-0.002446961
-0.001803774
-0.001160765
-0.000517935

0.054087651
0.003090327
0.003335293
0.002446961
0.001803774
0.001160765
0.000517935

revert
(-) continue
revert
(-) continue
(-) continue
(-) continue
STOP

xo =
error =
x =
i

4
0.001
1
x

f(xo) =

xi

1
1

-0.30547

OK

f(xi)

|f(xi)|
error

remarks

-0.30547

0.305472

(-) continue

1
0.1

5
4.1

1.091247
-0.21605

1.091247
0.216049

revert
(-) continue

2
3
4
4

0.1
0.1
0.1
0.01

4.2
4.3
4.4
4.31

-0.11692
-0.00757
0.112507
0.003945

0.116915
0.007571
0.112507
0.003945

(-) continue
(-) continue
revert
revert

4
5

0.001
0.001

4.301
4.302

-0.00642
-0.00528

0.006424
0.005276

(-) continue
(-) continue

6
7

0.001
0.001

4.303
4.304

-0.00413
-0.00298

0.004127
0.002977

(-) continue
(-) continue

8
9

0.001
0.001

4.305
4.306

-0.00183
-0.00067

0.001826
0.000674

(-) continue
STOP

Half-Interval Method (HIM)


Given f(xi-) < 0. Then, IF f(xi-)f(xi+1) < 0, change xi+ to xi+1, ELSE change xi- to xi+1.
i

xi -

xi +

xi

f(xi)

Error

Remarks

xo-

xo+

(xo- + xo+)/2

f(xo)

depend

IF f(xo) > 0, "change xo+"


IFf(xo) < 0, "change xo-"
IF error < d, "stop"

xi -

f(xo-) =
f(xo+) =
xi +

-0.5 OK
0.175639365 OK
xi

xo- = 0
xo+ = 1
i

Proceed
f(xi)

0.5

-0.142012708

1
2
3
4
5
6
7
8
9
10

0.5
0.5
0.625
0.6875
0.6875
0.703125
0.7109375
0.7109375
0.712890625
0.713867188

1
0.75
0.75
0.75
0.71875
0.71875
0.71875
0.71484375
0.71484375
0.71484375

0.75
0.625
0.6875
0.71875
0.703125
0.7109375
0.71484375
0.712890625
0.713867188
0.714355469

0.022504293
-0.058418971
-0.017613017
0.002533068
-0.007518287
-0.002487167
2.43142E-05
-0.001231086
-0.000603301
-0.000289472

xo- = 4

f(xo-) =

-0.30547

OK

xo+ = 5
i

f(xo+) =
xi +

1.091247
xi

OK
f(xi)

xi -

Error (di)

remarks
Change x-

0.333333333
0.2
0.090909091
0.043478261
0.022222222
0.010989011
0.005464481
0.002739726
0.001367989
0.000683527

Change x+
Change xChange xChange x+
Change xChange xChange x+
Change xChange xSTOP

Proceed
Error (di)

Remarks

0
1
2
3
4

4
4
4.25
4.25
4.25

5
4.5
4.5
4.375
4.3125

4.5
4.25
4.375
4.3125
4.28125

0.243868
-0.06355
0.081451
0.006841
-0.02887

0.058824
0.028571
0.014493
0.007299

Change x+
Change xChange x+
Change x+
Change x-

5
6
7

4.28125
4.296875
4.3046875

4.3125
4.3125
4.3125

4.296875
4.304688
4.308594

-0.01115
-0.00219
0.002319

0.003636
0.001815
0.000907

Change xChange xSTOP

Regula-Falsi
Given f(xi) < 0. Then, IF f(xi)f(xi+2) < 0, change xi+1 to xi+2, ELSE change xi to xi+2.
i

xi

xi+1

x1

f(xi)

x2

f(x1)

x3 x1

Formula:

f(x1) =

-0.5

x2 = 1
i

f(x2) =
xi

xi+1

0.17563936
f(xi)

1
2
3
4

0
0
0
0

1
0.74003977
0.71687319
0.71497415

1
2
3
4
5

f(x1) =
f(x2) =

f(x2)

Formula
below

OK
OK
f(xi+1)

-0.5
-0.5
-0.5
-0.5

0.17563936
0.01615806
0.00132804
0.00010811

OK
OK
f(xi+1)
1.09124698
1.09124698
1.09124698
1.09124698
1.09124698

xi

xi+1

-0.30547195
1.09124698
f(xi)

4
4.2187068
4.2826381
4.3001527
4.3048638

5
5
5
5
5

-0.30547195
-0.09725193
-0.02730998
-0.00739568
-0.00198302

Secant Method
Newton-Raphson Method (NRM)

xi+2

f(xi+2)

f(x3)

error

Remarks

depend

Given f(x1) < 0.


IF f(x3) > 0, "change x2+"
IF f(x3) < 0, "change x1-"
IF error < d, "stop"

x2 x1 f x1
f x2 f x1

x1 = 0

x1 = 4
x2 = 5
i

f(xi+1)

Proceed
xi+2

f(xi+2)

Error, di

Remarks

0.74003977
0.71687319
0.71497415
0.71481959

0.01615806
0.00132804
0.00010811
8.7938E-06

0.03231613
0.00265609
0.00021622

Change x2
Change x2
Change x2
STOP

xi+2

f(xi+2)

Error, di

Remarks

4.218706816
4.28263811
4.300152767
4.304863898
4.306124816

-0.09725193
-0.02730998
-0.00739568
-0.00198302
-0.00053029

0.01492801
0.00407303
0.00109437
0.00029282

Change xChange xChange xChange xSTOP

Proceed

You might also like