1 of 65

NUMERICAL METHODS

Solutions of Transcendental Equations

US03CMTH21

By: Dr. Snehal Shukla

2 of 65

BASIC CONCEPT OF � NUMERICAL METHODS:

  • Numerical methods are concerned with all aspects of the numerical solution of a problem, from the theoretical development and understanding of numerical methods to their practical implementation. Most numerical analysts specialize in small subfields, but they share some common concerns, perspectives, and mathematical methods of analysis.

What are numerical methods?

3 of 65

is a root finding method which repeatedly bisects an interval and then selects a subinterval in which a root must lie for further processing.

It is very simple and robust method.

Draw Back of the method:

It is a relatively slow and lengthy process

The Bi-section method in

numerical methods

4 of 65

Exa. Find the missing term of the following data.

(1)

x

0

1

2

3

4

y

3

6

11

?

27

S0l. Here we have four values. So, we need to go

up-to ∆4 in Binomial expansion .

5 of 65

(2)

x

0

1

2

3

4

y

1

3

9

?

81

(3)

x

0

1

2

3

4

5

y

0

-

8

15

81

35

S0l. Here we have five values of y. So, the fourth difference vanish, i.e. ∆4=0 for all values of x.

6 of 65

7 of 65

Bi-section Method

8 of 65

(1) Consider the transcendental equation

f(x)=0 --------(1)

which has zero in interval [a, b] and f(a)*f(b)<0.

Suppose f(a)<0 and f(b)>0

(2) Find first approximate

Explanation of Method

x1=a+b/2

Calculate f(x1) and examine its sign.

9 of 65

(2.1) If f(x1)<0 implies root lies between x1 and b,

than second approximate is given by

x2= x1 +b/2

(2.2) If f(x1)>0 implies root lies between x1 and a,

than second approximate is given by

x2= x1 +a/2

Calculate f(x2) and repeat steps (2.1) and (2.2).

Continuous this process until one can achieve the desire root of the equation up to required decimal places.

10 of 65

Using Bi-section method

find a root of the equation

correct up to three decimal places.

11 of 65

Sol.- First of all we have to find two nearest values for which function keeps opposite sign for respectively value.

Here f(0)=-1 and f(1)=1. (opposite sign)

Therefore, root of the equation lies in (0,1)

Compare [0,1] with [a,b], we have a=0 and b=1.

12 of 65

Iteration no.

a

b

f(a)

f(b)

f(xi)

update

1

0

1

0.5

-1<0

1>0

-0.625<0

a=

2

0.5

1

0.75

-0.625<0

1>0

-0.015625<0

a=

3

0.75

1

0.875

-0.015625<0

1>0

0.435546875>0

4

0.75

0.875

0.8125

-0.015625<0

0.435546875>0

0.196533203125>0

5

0.75

0.8125

0.78125

-0.015625<0

0.196533203125>0

0.0871887207031>0

13 of 65

7

0.75

0.78125

0.765625

-0.015625<0

0.0871887207031>0

0.0349769592285>0

8

0.75

0.765625

0.7578125

-0.015625<0

0.0349769592285>0

0.00947618484497

>0

9

0.75

0.7578125

0.75390625

-0.015625<0

0.00947618484497

>0

-0.0031241774559<0

10

0.75390625

0.7578125

0.755859375

-0.0031241774559<0

0.00947618484497

>0

0.0031635388732>0

14 of 65

11

0.75390625

0.755859375

0.7548828125

-0.0031241774559<0

0.00947618484497

>0

1.65672972798E-005>0

12

0.75390625

0.7548828125

0.75439453125

From the table it clear that has

same values after three decimal places.

Our desire root is 0.754

15 of 65

FALSE POSITION

METHOD

Object: To find real root of an equation.

16 of 65

Let f(x)=0 ---------------(1)

be a continuous function on [a,b], and f(a) and f(b) have opposite signs.

Therefore, there must be a real root of the equation in [a,b].

a

b

B(b, f(b))

A(a, f(a))

C(c,0)

17 of 65

A better approximation is find to point (c,0) on X-axis.

Find the slope of =

Similarly, find the slope of

=

In this way we have

=

18 of 65

c can be treated as an approximation of actual real root.

Now, we have

(1) f(a) and f(c) have opposite sign than root lies in [a,c].

(2) f(c) and f(b) have opposite sign than root lies in [c,b].

Once more we obtain next approximation x-coordinate of the point and joining the chord and do as above.

19 of 65

Continuing similarly we get a sequence of approximations x1,x2,---------

In this way one can find the desire root of the equation of (1) up to given (m) decimal places.

20 of 65

Example:

Find the real root of the equation

By method of False position or Regular falsi method

Up to three decimal places.

Sol.- First of all we have to find two nearest values for which function keeps opposite sign for respectively value.

Here f(2)=-9 and f(3)=6. (opposite sign)

Therefore, root of the equation lies in (2,3)

Compare [2,3] with [a,b], we have a=2 and b=3.

21 of 65

Iteration No.

A

a

B

f(a)

C

b

D

f(b)

E

B-D

F

AD

G

BC

H

G-F

I

H/E

f(Xi)

UPDATE

1

2

-9<0

3

6>0

-15

12

-27

-39

2.6

-1.824<0

2

2.6

-1.824<0

3

6>0

-7.824

15.6

-5.472

-21.072

2.69325153374

-0.237226510807<0

3

2.69325153374

-0.237226510807<0

3

6>0

-6.237226510807

16.1595092024

-0.7116795324

-16.8711887348

2.70491839693

-0.0289121838676<0

22 of 65

Iteration No.

2.70633348696

2.70650446856

23 of 65

From the table it clear that and

has same values after three decimal places.

Our desire root is 2.70 6

24 of 65

ITERATION

METHOD

Object: To find real root of an equation.

25 of 65

Statement:

Let be a root of f(x)=0 and let I be an

interval containing . Let and

Be continuous in I, where is defined by

which is equivalent to f(x)=0. Then if

For all x € I, the sequence of approximations

defined by

Converges to , provided that the initial approximation is chosen I.

26 of 65

Proof:

f(x)=0 is a real valued function, and be a root.

If we express the equation equivalently as

Then we have

Now, let be the initial approximation chosen in some interval I.

Suppose,

Apply Langrage’s mean value therefore

27 of 65

--------(1)

For some between and

Using we obtain the sequence

And corresponding inequalities given below

For some between and

Same we apply for the next we have

28 of 65

Continuous the process we have

Now, we have assume for some k

Then above inequalities(1), (2), (n+1) give

This implies that if for all i , then all

Successive approximations remain in I provided the initial approximation is chosen in I

29 of 65

Finally, we show that the sequence of successive approximations converges to .

Multiplying the inequalities(1), (2), (n+1) and then simplifying we get,

Since

Also since k<1, the R.H.S. tends to zero as

Therefore the sequence of successive approximations will converges to

If the condition in the theorem are satisfied then using iteration formula

30 of 65

We can find successive approximations .

31 of 65

Find a real root of

by Iteration method correct up to three decimal places.

32 of 65

Sol:

For f(x)=sinx-10(x-1) ----------(1)

f(1)=sin(1)-10(1-1)

= 0.841470984808>0

f(2) = sin(2)-10(2-1)

=-9.09070257317<0

There fore, root lies between 1 and 2.

We obtain initial approximation as a mid-value

33 of 65

By equ. (1)

According to theorem

Now working with the formula

34 of 65

As the condition

Which satisfies for all x between 1 and 2.

35 of 65

Functional

value

Iteration

value

1.09974949866

1.08910937057

1.08862146598

1.08859885204

From the table it is clear and has

Same value after three decimal places.

Our desire root is 1.088

36 of 65

AITKEN’S

-PROCESS

Object: To find real root of an equation.

37 of 65

Let be a real root of f(x)=0. So if we express the equation equivalently as

then we have

Now, let be the initial approximation chosen in some interval I.

Suppose, using Iteration method we have three successive approximations

to

Then for some k such that we have

PROCESS:

38 of 65

Dividing we can eliminate k and obtain

Simplifying

39 of 65

On the right hand side adding and subtracting

40 of 65

In terms of forward difference notation we have

Therefore above formula can be expressed as

As it gives next approximation to we have the

Aitken’s process approximation formula as given

41 of 65

EXAMPLE:

Solve by Aitken’s Process correct up to three decimal places.

Sol.

Root of the equation lies between 0 and 1.

Now, we obtain initial approximation as mid-value of an interval,

42 of 65

Now, let us express the given equation as subject

to x we

43 of 65

We get following successive approximation.

Finding first three values we can apply Iteration method.

44 of 65

Formulae

In (*) i=2

45 of 65

In this way

46 of 65

For i=3

47 of 65

In this way,

48 of 65

NEWTON-RAPHSON METHOD

N-R METHOD

Object: To find real root of an equation.

49 of 65

Let f(x)=0 be an equation and be initial

approximation to be a root, say , of the

Equation.

Suppose for some h.

As is a root of the equation , we have

Therefore

50 of 65

Initial x0

First

Approximation

51 of 65

If f satisfies all the conditions for Taylor’s

Series expansion then

Neglecting the terms with second and higher order derivatives in above expansion we have

52 of 65

If f satisfies all the conditions for Taylor’s

Series expansion then

Neglecting the terms with second and higher order derivatives in above expansion we have

53 of 65

Substituting in we get

Using similar arguments we can obtain a

generalized formula called Newton-Raphson approximation formula as

54 of 65

EXAMPLE:

Find a real root of , correct up to three decimal places, by using Newton-Raphson method.

Solution:

The initial approximation is

55 of 65

A

Initial value =xi

B

F(xi)

C

F’(xi)

D=B/C

F(xi)/f’(xi)

A-D

UP-DATE

VALUE

X1=0.5

0.625/-4.25

=-0.147058824

X2=0.5+0.147058824

0.647058824

X2=0.647058824

VALUE OF B

-3.74394463

VALUE OF C

-0.00951397136

VALUE OF D

0.656572795

VALUE OF X3

56 of 65

APPLICATION OF

N-R METHOD

APP-1: Formula for finding cube root of a equation

If is an approximation to the root of an equation f(x)=0, where f is differentiable on an interval containing , then by N-R formula

Now, suppose x is a real cube root of N.

57 of 65

Then we have

Now, to approximate the cube root, let us define

Here,

Therefore, if is an approximation to x then using N-R approximation formula

58 of 65

We get the next approximation by

59 of 65

Thus, the iterative formula for finding successive

Successive approximations of cube root of a

number N is

Exa.: Find the cube root of 12

60 of 65

APP-I1: Formula for finding square root of a equation

If is an approximation to the root of an equation f(x)=0, where f is differentiable on an interval containing , then by N-R formula

Now, suppose x is a real square root of N.

61 of 65

We get the next approximation by

62 of 65

Thus, the iterative formula for finding successive

approximations of square root of a number N is

Exa.: Find the square root of 8

63 of 65

APP-I1: Formula for finding square root of a equation

If is an approximation to the root of an equation f(x)=0, where f is differentiable on an interval containing , then by N-R formula

Now, suppose x is a real square root of N.

64 of 65

We get the next approximation by

65 of 65

Thus, the iterative formula for finding successive

approximations of square root of a number N is

Exa.: Find the square root of 8