2 Newton Ralphson Method
2 Newton Ralphson Method
2 Newton Ralphson Method
Autar Kaw
After reading this chapter, you should be able to:
1.
2.
3.
4.
Introduction
Methods such as the bisection method and the false position method of finding roots of
a nonlinear equation f ( x) = 0 require bracketing of the root by two guesses. Such
methods are called bracketing methods. These methods are always convergent since
they are based on reducing the interval between the two guesses so as to zero in on the
root of the equation.
In the Newton-Raphson method, the root is not bracketed. In fact, only one initial
guess of the root is needed to get the iterative process started to find the root of an
equation. The method hence falls in the category of open methods. Convergence in
open methods is not guaranteed but if the method does converge, it does so much
faster than the bracketing methods.
Derivation
The Newton-Raphson method is based on the principle that if the initial guess of the
root of f ( x) = 0 is at xi , then if one draws the tangent to the curve at f ( xi ) , the point
xi +1 where the tangent crosses the x -axis is an improved estimate of the root (Figure 1).
Using the definition of the slope of a function, at x = xi
f ( xi ) = tan
=
f ( xi ) 0
,
xi xi +1
Saylor.org
Page 1 of 15
which gives
xi +1 = xi
f ( xi )
f ( xi )
(1)
Equation (1) is called the Newton-Raphson formula for solving nonlinear equations of
the form f ( x ) = 0 . So starting with an initial guess, xi , one can find the next guess, xi +1 ,
by using Equation (1). One can repeat this process until one finds the root within a
desirable tolerance.
Algorithm
The steps of the Newton-Raphson method to find the root of an equation f ( x ) = 0 are
1. Evaluate f ( x ) symbolically
2. Use an initial guess of the root, xi , to estimate the new value of the root, xi +1 ,
as
f (xi )
xi +1 = xi
f ( xi )
3. Find the absolute relative approximate error a as
a =
xi +1 xi
100
xi +1
Saylor.org
Page 2 of 15
f (x)
f (xi)
[xi, f (xi)]
f (xi+1)
xi+2
xi+1
xi
Example 1
You are working for DOWN THE TOILET COMPANY that makes floats for ABC
commodes. The floating ball has a specific gravity of 0.6 and has a radius of 5.5 cm.
You are asked to find the depth to which the ball is submerged when floating in water.
Saylor.org
Page 3 of 15
The equation that gives the depth x in meters to which the ball is submerged under
water is given by
x 3 0.165 x 2 + 3.993 10 4 = 0
x1 = x0
f ( x0 )
f ( x0 )
= 0.05
= 0.05
1.118 10 4
9 10 3
= 0.05 ( 0.01242)
= 0.06242
Saylor.org
Page 4 of 15
a =
x1 x 0
100
x1
0.06242 0.05
100
0.06242
= 19.90%
The number of significant digits at least correct is 0, as you need an absolute relative
approximate error of 5% or less for at least one significant digit to be correct in your
result.
Iteration 2
The estimate of the root is
x2 = x1
f (x1 )
f ( x1 )
3
2
(
0.06242 ) 0.165(0.06242 ) + 3.993 10 4
= 0.06242
2
3(0.06242 ) 0.33(0.06242 )
= 0.06242
3.97781 10 7
8.90973 10 3
= 0.06242 4.4646 10 5
= 0.06238
a =
x2 x1
100
x2
0.06238 0.06242
100
0.06238
= 0.0716%
Saylor.org
Page 5 of 15
The maximum value of m for which a 0.5 102 m is 2.844. Hence, the number of
significant digits at least correct in the answer is 2.
Iteration 3
The estimate of the root is
x3 = x 2
f (x2 )
f ( x 2 )
= 0.06238
= 0.06238
4.44 10 11
8.91171 10 3
= 0.06238 4.9822 10 9
= 0.06238
0.06238 0.06238
100
0.06238
=0
The number of significant digits at least correct is 4, as only 4 significant digits are
carried through in all the calculations.
Saylor.org
Page 6 of 15
f ( x ) = ( x 1) + 0.512 = 0
xi +1 = xi
( xi 1)3 + 0.512
3( xi 1) 2
Starting with an initial guess of x0 = 5.0 , Table 1 shows the iterated values of the root of
the equation. As you can observe, the root starts to diverge at Iteration 6 because the
previous estimate of 0.92589 is close to the inflection point of x = 1 (the value of f ' ( x )
is zero at the inflection point). Eventually, after 12 more iterations the root converges to
the exact value of x = 0.2 .
Table 1 Divergence near inflection point.
Iteration
xi
Number
0
5.0000
3.6560
2.7465
2.1084
1.6000
0.92589
30.119
19.746
12.831
8.2217
10
5.1498
11
3.1044
Saylor.org
Page 7 of 15
12
1.7464
13
0.85356
14
0.28538
15
0.039784
16
0.17475
17
0.19924
18
0.2
2. Division by zero
For the equation
f (x ) = x 3 0.03 x 2 + 2.4 10 6 = 0
xi +1 = xi
Saylor.org
Page 8 of 15
For x0 = 0 or x0 = 0.02 , division by zero occurs (Figure 4). For an initial guess close to
0.02 such as x0 = 0.01999 , one may avoid division by zero, but then the denominator in
the formula is a small number. For this case, as given in Table 2, even after 9
iterations, the Newton-Raphson method does not converge.
f ( xi )
0.019990
1.60000 10-6
2.6480
18.778
100.75
1.7620
5.5638
50.282
1.1714
1.6485
50.422
0.77765
0.48842
50.632
0.51518
0.14470
50.946
0.34025
0.042862
51.413
0.22369
0.012692
52.107
0.14608
0.0037553
53.127
0.094490
0.0011091
54.602
a %
Number
Saylor.org
Page 9 of 15
1.00E-05
f(x)
7.50E-06
5.00E-06
2.50E-06
0.00E+00
-0.03
-0.02
-0.01
0
-2.50E-06
0.01
0.02
0.03
0.04
0.02
-5.00E-06
-7.50E-06
-1.00E-05
Saylor.org
Page 10 of 15
f(x)
3
2
0
-2
-1.75
-1
-0.3040
0.5
3.142
-1
f ( xi )
a %
Number
Saylor.org
Page 11 of 15
1.0000
3.00
0.5
2.25
1.75
5.063
0.30357 2.092
3.1423
11.874
1.2529
3.570
0.17166 2.029
5.7395
34.942
2.6955
9.266
0.97678
2.954
300.00
128.571
476.47
109.66
150.80
829.88
102.99
112.93
175.96
4. Root jumping
In some case where the function f (x) is oscillating and has a number of roots, one may
choose an initial guess close to a root. However, the guesses may jump and converge
to some other root. For example for solving the equation sin x = 0 if you choose
x0 = 2.4 = (7.539822 ) as an initial guess, it converges to the root of x = 0 as shown in
Table 4 and Figure 6. However, one may have chosen this as an initial guess to
converge to x = 2 = 6.2831853 .
f ( xi )
a %
Number
Saylor.org
Page 12 of 15
f(x)
7.539822
0.951
4.462
0.969
0.5499
0.5226
711.44
0.06307
0.06303
971.91
8.376 10 4
8.375 10 5
7.54 10 4
1.95861 10 13
1.95861 10 13
4.28 1010
68.973
1.5
0.5
0
-2
-0.06307
0.5499
4.461
7.539822
10
-0.5
-1
-1.5
down-to-up is called its inflection point. For example, for the function f ( x ) = ( x 1) , the
concavity changes at x = 1 (see Figure 3), and hence (1,0) is an inflection point.
Saylor.org
Page 13 of 15
An inflection points MAY exist at a point where f ( x) = 0 and where f ' ' ( x ) does
not exist. The reason we say that it MAY exist is because if f ( x) = 0 , it only makes it a
possible inflection point. For example, for f ( x ) = x 4 16 , f (0) = 0 , but the concavity
does not change at x = 0 . Hence the point (0, 16) is not an inflection point of
f ( x ) = x 4 16 .
3
f" (xi )
(xi +1 xi )2 + L
2!
As an approximation, taking only the first two terms of the right hand side,
f ( xi +1 ) f ( xi ) + f ( xi )( xi +1 xi )
which gives
xi +1 = xi
f ( xi )
f' ( xi )
This is the same Newton-Raphson method formula series as derived previously using
the geometric method.
Source URL: http://numericalmethods.eng.usf.edu/
Saylor URL: http://www.saylor.org/courses/me205/
Saylor.org
Page 14 of 15
Saylor.org
Page 15 of 15