Numerical Analysis - MTH603 Handouts Lecture 5

Numerical Analysis –MTH603 VU
Example
Using Regula-Falsi method, find the real root of the following equation correct, to three
decimal places: x log 10 x =1.2
Solution:
Let f (x) = x log10 x − 1.2
f (2) = 2 log10 2 − 1.2 = − 0.5979,
f (3) = 3 log10 3 − 1.2 = 0.2314.
Since f (2) and f (3) are of opposite signs, the real root lies betweenx1 = 2, x2 = 3.
The first approximation is obtained from
x2 − x1 3− 2
x3 = x2 − f ( x2 ) = 3 − (0.2314)
f ( x2 ) − f ( x1 ) 0.2314 + 0.5979
0.2314
= 3− = 2.72097
0.8293
f ( x3 ) = Let f (x) = 2.72097log10 2.72097 − 1.2 = − 0.01713.
Since f (x2) and f (x3) are of opposite signs, the root of f (x) = 0 lies between x2 and
x3. Now, the second approximation is given by
x3 − x2 2.72097 − 3
x4 = x3 − f ( x3 ) = 2.72097 − (−0.1713) = 2.7402
f ( x3 ) − f ( x2 ) −0.1713 − 0.2314
f ( x4 ) = 2.7402 log10 2.7402 − 1.2 − 3.8905 ×10−4
Thus, the root of the given equation correct to three decimal places is 2.740
NOTE: Here if TOL is given then we can simply find the value of TOL by subtracting
both the consecutive roots and write it in the exponential notation if the required TOL is
obtained then we stop.
Method of Iteration
Method of iterations can be applied to find a real root of the equation f (x) = 0 by
rewriting the same in the form.
x = φ ( x)
Let x = x0 be the initial approximation to the actual root, say, α of the equation .then the
first approximation is x1 = φ ( x0 ) and the successive approximation are x2 = φ ( x1 )
x3 = φ ( x2 ), x4 = φ ( x3 ),..., xn = φ ( xn −1 ) if sequence of the approximate roots, x1 , x2 , x3 ,...xn
converges to α it is taken as the root of the equation f(x)=0.
For convergence purpose the initial approximation is to be done carefully.the choice
of the x0 is done according to the theorem.
© Copyright Virtual University of Pakistan 1

Theorem
If α be a root of f(x) =0 which is equivalent to x = φ ( x) I be any interval containing the
point x= α and | φ '( x) |< 1 ∀ xε I , then the sequence of approximations x1 , x2 , x3 ,...xn will
converge to the root α provided that the initial approximation x0 is chosen in I
Example,
f (x) = cos x - 2x + 3 = 0.
It can be
1
Rewritten as x = (cos x + 3) = φ ( x)
2
1
φ ( x) = (cos x + 3)
2
f (x) = cos x - 2x + 3 = 0.
f(1)=cos 1 - 2(1) + 3 =1.54030>0
f(2)=cos 1 - 2(2) + 3 =-0.041614-4+3=-1.41614<0
so root lies between 1and 2
1
φ '( x) = − (sin x)
2
both φ '(1)andφ '(2) < 1 so the method of iterations can be applied
let x0 = 1.5
1 1
x1 = (cos x0 + 3) = (cos(1.5) + 3) = 1.999825
2 2
1 1
x2 = (cos x1 + 3) = (cos(1.999825) + 3) = 1.999695
2 2
1 1
x3 = (cos x2 + 3) = (cos(1.999625) + 3) = 1.999695
2 2
So this is the required root correct up to 5 places of decimal.
Example
Find the real root of the equation x3 + x 2 − 1 = 0 by method of iterations
Solution
let f ( x) = x 3 + x 2 − 1
now f (0) = 03 + 02 − 1 = −1 < 0
f (1) = 13 + 12 − 1 = 1 > 0
hence a real root lies between 0 and 1

here
x3 + x 2 − 1 = 0
x 2 ( x + 1) = 1
1 1
x2 = ⇒x= = φ ( x)
( x + 1) x +1
3
here φ '( x) = −1/ 2[1/( x + 1) 2 ]
5
φ '(0) = 1/ 2 < 1 and φ '(1) = 1/ 2 2 < 1
so φ '( x) < 1 for all the values in the int erval
let x0 = 0.65
1 1
x1 = φ ( x0 ) = = = 0.7784989
x0 + 1 1.65
1 1
x2 = φ ( x1 ) = = = 0.7498479
x1 + 1 1.7784989
1 1
x3 = φ ( x2 ) = = = 0.7559617
x2 + 1 1.7498479
1 1
x4 = φ ( x3 ) = = = 0.7546446
x3 + 1 1.7559617
1 1
x5 = φ ( x4 ) = = = 0.7549278
x4 + 1 1.7546446
1 1
x6 = φ ( x5 ) = = = 0.7548668
x5 + 1 1.7549278
1 1
x7 = φ ( x6 ) = = = 0.7548799
x6 + 1 1.7548668
1 1
x8 = φ ( x7 ) = = = 0.7548771
x7 + 1 1.7548799
1 1
x9 = φ ( x8 ) = = = 0.7548777
x8 + 1 1.7548771
1 1
x10 = φ ( x9 ) = = = 0.7548776
x9 + 1 1.7548777
1 1
x11 = φ ( x10 ) = = = 0.7548776
x10 + 1 1.7548776
hence root is 0.7548776
Note: In this question the accuracy up to 7 places is acquires or here the TOL is 10−7

Example
Find a real root of the equation cos x = 3 x − 1 correct to seven places of decimal.
Solution
Here it is a transcendental function and all the calculation must be done in the radians
mode and value of pi should be 3.14
f ( x) = cos x − 3 x + 1
f (0) = cos 0 − 3(0) + 1 = 1 > 0
f (π / 2) = cos(1.57) − 3(1.57) + 1 = 0.0007963 − 4.71 + 1 = −3.7092037 < 0
so a real root lies between 0 and π / 2
1
here φ ( x) = (cos x + 1)
3
1
we have φ '( x) = − sin x
3
it is clearly less than 1 as sin is a bounded function and it ' s values lies between − 1 and 1
hence iteration method can be applied
let x0 = 0.5 be the inital approximation then
1
x1 = φ ( x0 ) = [cos(0.5) + 1] = 0.6258608
3
1
x2 = φ ( x1 ) = [cos(0.6258608) + 1] = 0.6034863
3
1
x3 = φ ( x2 ) = [cos(0.6034863) + 1] = 0.6077873
3
1
x4 = φ ( x3 ) = [cos(0.6077873) + 1] = 0.6069711
3
1
x5 = φ ( x4 ) = [cos(0.6069711) + 1] = 0.6071264
3
1
x6 = φ ( x5 ) = [cos(0.6071264) + 1] = 0.6070969
3

1
x7 = φ ( x6 ) = [cos(0.6070969) + 1] = 0.6071025
3
1
x8 = φ ( x7 ) = [cos(0.6071025) + 1] = 0.6071014
3
1
x9 = φ ( x8 ) = [cos(0.6071014) + 1] = 0.6071016
3
1
x10 = φ ( x9 ) = [cos(0.6071016) + 1] = 0.6071016
3

Numerical Analysis - MTH603 Handouts Lecture 5

Uploaded by

Copyright:

Available Formats

Numerical Analysis - MTH603 Handouts Lecture 5

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Numerical Analysis - MTH603 Handouts Lecture 5

Uploaded by

Copyright:

Available Formats

Numerical Analysis –MTH603 VU

The first approximation is obtained from

© Copyright Virtual University of Pakistan 1

So this is the required root correct up to 5 places of decimal.

Find the real root of the equation x3 + x 2 − 1 = 0 by method of iterations

© Copyright Virtual University of Pakistan 2

© Copyright Virtual University of Pakistan 3

© Copyright Virtual University of Pakistan 4

© Copyright Virtual University of Pakistan 5

You might also like