14.04.2017 Views

Numerical Analysis By Shanker Rao

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

22 NUMERICAL ANALYSIS<br />

2.3 METHOD OF BISECTION<br />

Consider the equation<br />

f(x) = 0, refer (1)<br />

where f(x) is continuous on (a, b) and f(a) f(b) < 0. In order to find a root of (1) lying in the interval<br />

(a, b). We shall determine a very small interval (a 0<br />

, b 0<br />

) (by graphical method) in which f(a 0<br />

) f(b 0<br />

)<br />

< 0 and f ′ ( x) maintains the same sign in (a 0<br />

, b 0<br />

), so that there is only one real root of the equation<br />

f(x) = 0.<br />

Divide the interval in half and let<br />

x<br />

1<br />

a<br />

=<br />

+ b<br />

2<br />

0 0<br />

If f(x) = 0 then x 1<br />

is a root of the equation. If f(x 1<br />

) ≠ 0 then either f(a 0<br />

) f(x 1<br />

) < 0 of f(b 0<br />

)<br />

f(x 1<br />

) < 0. If f(a 0<br />

) f(x 1<br />

) < 0 then the root of the equation lies in (a 0<br />

, x 1<br />

) otherwise the root of the<br />

equation lies in (x 1<br />

, b 0<br />

). We rename the interval in which the root lies as (a 1<br />

, b 1<br />

) so that<br />

1<br />

b1<br />

– a1 = bb0 − a0g,<br />

2<br />

now we take<br />

x<br />

2<br />

a<br />

=<br />

+ b<br />

2<br />

1 1<br />

If f(x 2<br />

) = 0 then x 2<br />

is the root of f(x) = 0. If f(x 2<br />

) ≠ 0 and f(x 2<br />

) f(a 1<br />

) < 0, then the root lies<br />

in (a 1<br />

, x 2<br />

). In which case we rename the interval as (a 2<br />

, b 2<br />

), otherwise (x 2<br />

, b 1<br />

) is renamed as (a 2<br />

,<br />

b 2<br />

) where<br />

Proceeding in this manner, we find<br />

a2 − b2 = 1 b a<br />

2<br />

b 0 − 0g.<br />

2<br />

x<br />

n+ 1 =<br />

a<br />

n<br />

+ b<br />

2<br />

n<br />

which gives us the (n + 1)th approximation of the root of f(x) = 0, and the root lies (a n<br />

, b n<br />

) where<br />

1<br />

bn − an = b − a<br />

n<br />

b 0 0g,<br />

2<br />

since the left end points a 1<br />

, a 2<br />

, …, a n<br />

, … form a monotonic non-decreasing bounded sequence, and<br />

the right end points b 1<br />

, b 2<br />

, b 2<br />

, …, b n<br />

, … form a monotonic non-increasing bounded sequence, then<br />

there is a common limit<br />

c=<br />

lim<br />

n →∞<br />

a<br />

n<br />

=<br />

n →∞<br />

such that f(c) = 0 which means that c is a root of equation (1).<br />

The bisection method is well suited to electronic computers. The method may be conveniently<br />

used in rough approximations of the root of the given equation. The bisection method is a simple<br />

but slowly convergent method.<br />

lim<br />

b<br />

n

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!