Regula Falsi Method
Regula Falsi Method
EQUATIONS 57
Similarly by performing subsequent iterations. the successive approximations to the root are given
by
x=0.7137, xs = 0.7343, x,= 0.7421. X, = 0.7382, x = 0.7402, x, =0.7392, xo =0.7388, x 0.7385
From 10th and 11th iteration, it is clear that there is no change in the succeasive approximations to
root upto first three decimal places. the
So a real root of given equation is given by x = 0.738 (corect to three decimal
places).
2. THE METHOD OF INTERPOLATION (OR THE METHOD OF FALSE
POSITION OR
REGULA-FALSI METHOD)
The method of false position is the oldest method and dates back to the ancient Egyptians. This
method is also known as the method of linear interpolation or the Regula-falsi method. It remains an
efiective altermative to the bisection method for solving nonlinear equation of the form f (%) =0for areal
root in the interval (a, b), given that fis continuous on the interval [a, b] and f(a) and f (b) have opposite
signs.
Procedure
Without loss of generality, suppose that a continuous function f is negative at a and positive at b. So
there exists atleast one real root between aand b. (We may also take fas positive at a and negative at b).
Now we connect the two points (a,f(a) and (6,f(b) on the sketch of the curve y *f(*) by a straight
line segment (see fig. 1).
The equation of the chord joining the points (a,f(a) and (6,f(6) is given by
y-f(a) = f(6)-f(a) (*-a)
b-a
The point of intersection of this chord with the x-axis is given by putting y = 0 in the above equation.
i.e. 0-f(a) f(6)-f(0)x-a)
b-a
(6-a)
or x-a =-f(@)
f(b)- f(a)
(b-a)
f(a()
f(b)-
af(b) - bf(a)
or
f(6)-f(a)
This value ofx is taken as the approximate value of the root.
Hence the first approximate value of the root say x, is given by
af(b)- bf(a) .()
f(b)- f(a)
As in the bisection method, there are three possibilities :
) f(x) =0, in which case x, is the root.
(i) f(x ) <0, in which case the root lies in the interval (x,,b)
(iüi) f(x,) >0, in which case the root lies inthe interval (a, x )
BRILLIANI NUMERICAT
occurs.
thc Same
or
cAsc
(i)
pt.
()
case
leminates,
if
ciher
the
s u c c e s s i v e
points where
which
Now if case (0) occurs, the proces d e s i r e d a c c u r a c y
I in
0s repcated until the root is obtaincdto the in fig.
g r a p h i c a l l y
Fig. 1
=0 by (ro,
Xi) instead of (a, b,
General Formula for Regula-Falsi Method equationf(x)
of the
Denoting the initial interval containing the root
next approximation to the root say x is given by
*o1(a)-3fo)
f(x)-fz)
If e-'k and ek+1 are the errors in the approximations x,1, X, and xa respectively then
e-k-| -5, 4 =X-§ and ek+1k|-5
ie. x- $+ e- »X *+ e and x+ =+ e+|
Substituting these values in fornula (4), we have
(6+e;)-(F+ e-)
sG+e4)-fG+e-1)
..(5)
Expanding fE+ e) and f ( ep-) using Taylor's series about the point and substituting
f(E) =0 in (5), we have
,f()+;f")+...
= e
S') +;(e tey-1)s"9) +...
2 f)
60
higher
orders of errors
having
terms
. 2 / " 4
orders of erTors.
higher
having
erms
2 f()
SO
terms
having higher orders of en
2 f(9) (neglecting
the
in called lerror
equation.
(,())
lies above the
The equation (6) in terms of errors and with k.
the points
(*f())
point
varies If the
Assuming the straight line joining fixed and
other
straight line
joining the points (Z:f
by the
approximated
e7+l
2 f(5)
or
constant.
where A is asymptotic error
2 f(9)
i.e. orde
clear that the regula-falsi method is linear order convergent
From cquation (7) it is
convergence of this method is 1.
method. Secant method resembles Regula-falsi method except that no attempt is made to ensure tna
root is enclosed between two approximations.
sOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS 61
Starting with two initial approximations X and x, to the root E, the further approximations
Is, Xj. are computed using the following formula :
y=f)
>X
Fig. 1.
Order of Convergence
Let us assume that is a simple root of equationf(*) = 0. The general iterative formula for secant
method is
...(3)
If e- e and e, + are the errors in the approximations x_,X and xrespectively then
point. ¿ and
about
the ub.
serics
Taylor's
)using
1, 2/"(5)
e
second
term on R.H.S.
.byf (&
denominator
of
numerator and
(By dividing
T*
em
+ terms having higher orders of
errors
e,e,- +terms having higher orders of
2 f()
I.e
So the order of convergence of secant method is 1.618, which sometimes referred to as super linear
convergence.
Note, The speed of convergence of secant method is faster than that of Bisection and Regula-Falsi method
ILLUSTRATIVE EXAMPLES]
Example 1. Find a real root of the equation * +x- | =0 correct to three decimal places by using false
position method.
Sol. The given equation is x+x-l =0.
Let fe) =+*-1
Now f(0.5) =(0.5)' +0.5-|= -0.375 <0
and f() ='+1-| =| >0
So a real root of given equation lies in the interval (0.5, 1)
Iteration I. Taking a= 0.5 and b= Iso that f(a) =-0.375 and f(b) =|
The first approximation to the root is given by
af (b) - bf(a) 0.5(0) - l(-0.375)
0.6364
S(b)- fa) |-(-0.375)
-1.
a n d / ( l )
1058
64
((a)
- 0.
b-|
sothat
Iteration 2. Taking a-0.6364
a and
=0.6712.
given by -I(-0.1058)
1-(-0.1058)
bla)
af(b)-
f(b)-fla) 00264
<0
|=-
+0.6712
-|
-
Now S(0.6712) - (0.6712)
f ) =|>0 (0.6712,
1) =|
and interval and/(b)
f(b)-bf(a) I-(-0.0264)
f(b)-f(a) <0
=-0.0065
-|
Now f(0.6796) =(0.6796)' +0.6796-
and f(1) =1>0 (0.6796, 1)
lies in the
interval
andf(6) =1
So arealroot of the given equation =-0.0065
fla)
Iteration 4. Taking a =0.6796 and b =1 so that
The fourth approximation to the root is given by
0.6796() -I(-0.0065) =0.6817
g6)-bf(a)
1-(-0.0065)
f(b)-f(a)
=-0.0015 <0
Now J0.6817) -(0.6817)°+0.6817- 1
and S) =1>0
the interval (0.68l7, 1)
So a real root of the given eauation lies in
Takinga= 0.6817 and b=1 so that f(a) =-0.0015 and/ (b) = 1
Iteration 5.
The fifth approximation to the root is given by
af(b) -bf(a)0.6817)-I(0.0015) -0 6822
f(b)-f(a) I-(0.0015)
and/(b)-6
66 (aj-17f
f
oIhat
IterationI. Taking a-25 andh-l
given by
The first ayroximation tothe root n 26800
a/(b) h fta) 6175)
9--044712 )
4 (2 6800)
Now (26800) (2 6800)
and 3)60 3)
interval (26800, and/(b)- 6
So arcal root of given cquation lies inthe 0.4712
fla) -
Iteration 2. Takinga 26800 and h 3so that
The second approximation tothe root is given by
2.7033
26800(6) - 3(-0.4712)
/(6)-bf(a)
6-(-04712)
fb)- f(a)
0,0579 <0
Now f(2.7033) =(.7033) - 4 (2.7033) - -9=-
-(1.40279)- I
Now
f(1.40279) =(1.40279) -(1.40279)
and
S(1.5) = 1.95312 >0
So a real root of given equation lies in the interval (1.402/7, 12
=-0.01274 and f(6) = 1.953
Iteration 2. Taking a= 1.40279 and b= 1.5 so that f(a)
The second approximation to the root is given by
(-0.0 1274) = 140342
a f(6)-b f(a) 1.40279 (|I.953 12) -1.5
1.95312-(-0.01274)
f(6)-f(a)
Now
J(1.40342) (1.40342)-(1.40342) -(1.40342) - l=-0.00286 <0
and
S(1.5) =1.95312 >0
So a real root of given equation lies in the interval (1.40342, 1.5)
lteration 3. Taking a-1.40342 and b-1.5 sothat f(a)=-0.00286 and f(b)= 1.953 12
The third approximation to the root is given by
a f(b)-bf(a) 1.40342 (1.953 12) -1.5 (-0.00286) = 1.40356
f(b)- f(a) 1.953 12-(-0.00286)
Now
S(140356) = (1.40356)° -(1.40356) -(1l.40356)|=-0.00066 <0
and (1.5) = 1.953 12>0
So a real root of given equation lies in the interval
(1.40356, 1.5)
Iteration 4. Taking a = 1.40356 and b=1.5 so that
f(a) =-0.00066 and f(b) = 1.95312
The fourth approximation to the root is
given by
a f(b)-b f(a) 1.40356 (1.953 12) -1.5
f(b)- f(a) (-0.00066) =1.40359
From 3rd and 4th iteration, we see 1.95312--(-0.00066)
upto first four decimal places. that there is no change in the successive
approximations to the ro
So a real root of given equation is
given by x= |4035 (correct to four
decimal places)
SOLUTION OFF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS 69
Eample 6. Find the fourth root of 32 correct lo three decimal places using the method of false position.
Sol. Let =(32)/4
r-32 0
Let f(o) = x -32
Now f(2.25) =(2.25) -32 =-6.3 <0
and f(2.5) =(2.5) -32 =7.0625 >0
So a real root of equation (1) lies in the interval (2.25, 2.5)
Iteration 1. Taking a= 2.25 and b=2.5 so that f(a) =-6.3 and f(b) =7.0625
The first approximation to the root is given by
a f(b)-b f(a) 2.25 (7.0625)-2.5(-0.) -2.3574
f(b)-f(a) 7.0625-(-6.3)
Now f(2.3574) =(2.3574-32 =- .|160 <0
and f(2.5) = 7.0625 >0
So a real root of the equation (1) lies in the interval (2.3574, 2.5)
Iteration 2. Taking a=2.3574 and b=2.5 so that f(a) =- 1.1160 and f(b)=7.062)
The second approximation to the root is given by
af(b)-bf(a)
= 2.3574(7.0625)- 2.5(-1.1160) 3768
f(b)- f(a) 7.0625 (-1.1160)
Now f(2.3768) =(2.3768)'-32 =-0.0868 <0
and f(2.5) = 7.0625 > 0
So a real root of the equation (1) lies in the interval (2.3768, 2.5)
Iteration 3. Taking a=2.3768 and b =2.5 so that f(a) =-0.0868 and f(b) =7.0625
The third approximation to the root is given by
af(b)-bf(a) 2.3768 (7.0625) -2.5(-0.0868) = 2.3783
f(b)-f(a) 7.0625-(-0.0868)
Now f(2.3783) = (2.3783y' -32 =-0.0061 <0
and S(2.5) =7.0625 >0
CorectL
position
lalse
method
of
the
70 L2 by
log,o r
equation
r
Example 7. Find areal root ofthe
decimal places. -1.2
0
Sol. Rewriting the given cquation as x logyo
2.7- 1.2
Now f(2.7) = 2.7 logjoo >0
0.23136
3-1.2 = 0.23136
and f(3) = 3log,o (2.7. 3) 6)=0.
interval and
inthe =-0.035318
f(6)-f(a) -0.00080
<0
2.73973- 1.1.2 =
Now f(2.73973) =2.73973 logo
and S3) =0.23 136 > 0 (2.73973. 3) =0.23136
the interval I and S6)
So areal root of the equation lies in =-0.00080
that f(a)
Taking a =2.73973 and b= -3 so
Iteration 2. T
The second approximation to the root is given by (0.23136) -3
(-0.00080)
- 2.74063
2.73973
a f(b)-b f(a) 0 . 2 3 1 3 6 - ( - 0 . 0 0 0 8 0 )
f()-f(a)
2.74063- 1.2=-0.00001 <0
Now J2.74063) =2.74063 log,o
and f(3) =0.23136 > 0
(2.74063. 3)
So a real root of the equation lies in the interval =0.23136
b=3sothat f(a) =-0.00001 and f(6)
Iteration 3. Taking a= 2.74063 and
The third approximation to the root is given by
a f(b) -b fa) 2.74063 (0.23136) -3(-0.00001) =2.74064
0.23136-(-0.00001)
f(b)-f(a)
From second and third iteration, it is clear that there is no change in the successive approximation:
the root upto first four decimal places.
So areal root of the given equation is given by x=2.7406 (correctto 4decimal places)
Example 8. Find a real root of theequation e -r =0corect to four decimal places using the metho
false position.
Sol. The given equation is e-r=0.
Let f) = e-r
Now fU.75) = e-(0.75)' =0.39523 >0
and 061094<0
So a real root of given equation lies in the interval (1.75, 2)
coLUTION OF ALGEBRAIC AND TRANSCENDENTAL 71
EQUATIONS
Iteration . Taking a= 1.75 and b-2 so that f(a) 0.39523 and ) 0.61094
The first approximation to the root is
given by
a f(b)-b f(a) 1.75(-0.61094) -2 (0.39523)
I.84820
f(b)-f(a) 0.61094-0.39523
4
a n d / ( 6 )
-6375
72
t h a t / ( u )
-2 No
Iteration I. Taking - 1.5 and b
"(0)-6Ila)given
54-%-6437L8084
by
The first approximation tothe Toot is
4-(64375)
(b)-f(a) 4
that/(a)
8084so
Iteration2. Taking a 2 and b | 10-11|35
L8084 -
and f(b) - (1.8084)"-
given by
The second approxiNimationtothe root is 24-LI|35)-1.8084(4) 1.8501
aJ6)-bfta) -L.|135-4
1.|135
f(b)- f(a) =-
f(a)
that
1.8501 so
Iteration 3. Taking a= 1.8084 and h =
0.1341
and f(b) =(1.8501 ) - 18501 - 10= -
The third approximation to the root is given by
(6-bf(a) L 8 0 8 4 ( - 0 . 1 3 41 ) - 1 . 8 5 0 1 ( - 1 . I 1 3 5 ) = 1.8558
-0.1341-(-I.135)
f(b)-f(a)
=- 0.1341
that/(a)
Iteration 4. Taking a=1.8501 and b= 1.8558 so
f(6)-f(a) 0.25-(-3)