Numerical Methods With Applications
Numerical Methods With Applications
APPLICATIONS
(MEC500)
x3 − 0.165x 2 + 3.993 × 10 −4 = 0
Use the bisection method of finding roots of equations to find the depth to
which the ball is submerged under water.
Conduct three iterations to estimate the root of the above equation.
Find the absolute relative approximate error at the end of each iteration,
and the number of significant digits at least correct at the end of each
iteration.
From http://numericalmethods.eng.usf.edu
Engineering Application
Bisection Method
From the physics of the problem, the ball would be submerged between
and, x = 0 and x = 2R, where R is the radius of the ball, that is:
0 ≤ x ≤ 2R
0 ≤ x ≤ 2(0.055)
0 ≤ x ≤ 0.11
Let us assume
xl = 0, xu = 0.11
€
Check if the function changes sign between xl and xu.
f(xl ) = f (0) = (0) 3 − 0.165(0) 2 + 3.993 × 10 −4 = 3.993 × 10 −4
€
f (xu ) = f (0.11) = (0.11) 3 − 0.165(0.11) 2 + 3.993 × 10 −4 = −2.662 × 10 −4
Hence
f(xl ) ⋅ f (xu ) = (3.993 × 10 −4 ) ⋅(−2.662 × 10 −4 ) < 0
€
Engineering Application
Bisection Method
1st Iteration: The estimate of the root is:
xl + xu 0 + 0.11
xr = = = 0.055
2 2
f(xr ) = f (0.055) = (0.055) 3 − 0.165(0.055) 2 + 3.993 × 10 −4 = 6.655 × 10 −5
f (xl ) ⋅ f (xr ) = f (0) ⋅ f (0.055) = (3.993 × 10 −4 ) ⋅(6.655 × 10 −5 ) > 0
€
Hence the root is bracketed between xr and xu, that is, between 0.055
and 0.11. So, the lower and upper limit of the new bracket is:
xl = 0.055, xu = 0.11
At this point, the absolute relative approximate error cannot be
calculated as we do not have a previous approximation. ∈a
€
Engineering Application
Bisection Method
2nd Iteration: The estimate of the root is:
xl + xu 0.055 + 0.11
xr = = = 0.0825
2 2
f(xr ) = f (0.0825) = (0.0825) 3 − 0.165(0.0825) 2 + 3.993 × 10 −4 = −1.622 × 10 −4
f (xl ) ⋅ f (xr ) = f (0.055) ⋅ f (0.0825) = (6.655 × 10 −5 ) ⋅(−1.622 × 10 −4 ) < 0
Hence the root is bracketed between xl and xr, that is, between 0.055 and
0.0.0825. So, the lower and upper limit of the new bracket is:
xl = 0.055, xu = 0.0825
The absolute relative approximate, ∈a is 33.33%. None of the significant
digits are at least correct in the estimated root of xr = 0.0825 because
the absolute relative approximate error is greater than 5%.
€
Engineering Application
Bisection Method
3rd Iteration: The estimate of the root is:
xl + xu 0.055 + 0.0825
xr = = = 0.06875
2 2
f(xr ) = f (0.06875) = (0.06875) 3 − 0.165(0.06875) 2 + 3.993 × 10 −4 = −5.563 × 10 −5
f (xl ) ⋅ f (xr ) = f (0.055) ⋅ f (0.06875) = (6.655 × 10 −5 ) ⋅(−5.563 × 10 −5 ) < 0
Hence the root is bracketed between xl and xr, that is, between 0.055
and 0.0.06875. So, the lower and upper limit of the new bracket is:
xl = 0.055, xu = 0.06875
The absolute relative approximate, ∈a is 20%. Still none of the
significant digits are at least correct in the estimated root of the
€ equation as the absolute relative approximate error is greater than 5%.
Engineering Application
Bisection Method
Seven more iterations were conducted and these iterations are
shown in Table below:
Hence the number of significant digits at least correct is given by the largest
value of for m which
€ ε a ≤0.5 × 10 2−m
0.1721 ≤ 0.5 × 10 2−m
0.3442 ≤ 10 2−m
log(0.3442) ≤ 2 − m
m ≤ 2 − log(0.3442) = 2.463