Engineering
Computation
Lecture 3
E. T. S. I. Caminos, Canales y Puertos
1
Roots of Equations
Objective:
Solve for x, given that f(x) = 0
-orEquivalently, solve for x such that
g(x) = h(x) ==> f(x) = g(x) h(x) = 0
E. T. S. I. Caminos, Canales y Puertos
2
Roots of Equations
Chemical Engineering (C&C 8.1, p. 187):
van der Waals equation; v = V/n (= volume/# moles)
Find the molal volume v such that
f(v) = ( p +
a
v
)(v - b) - RT = 0
p = pressure,
T = temperature,
R = universal gas constant,
a & b = empirical constants
E. T. S. I. Caminos, Canales y Puertos
3
Roots of Equations
Civil Engineering (C&C Prob. 8.17, p. 205):
Find the horizontal component of tension, H, in
a cable that passes through (0,y0) and (x,y)
H
wx
f(H) = cosh
- 1 + y 0 - y = 0
w
H
w = weight per unit length of cable
E. T. S. I. Caminos, Canales y Puertos
4
Roots of Equations
Electrical Engineering (C&C 8.3, p. 194):
Find the resistance, R, of a circuit such that the
charge reaches q at specified time t
f(R) = e
-Rt 2L
cos
2
1 R
-
LC 2L
t - q = 0
q0
L = inductance,
C = capacitance,
q0 = initial charge
E. T. S. I. Caminos, Canales y Puertos
5
Roots of Equations
Mechanical Engineering (C&C 8.4, p. 196):
Find the value of stiffness k of a vibrating mechanical
system such that the displacement x(t) becomes zero at
t= 0.5 s. The initial displacement is x0 and the initial
velocity is zero. The mass m and damping c are known,
and = c/(2m).
x(t) = x0e-t cos t + sin t = 0
in which
k c2
=
m 4m 2
E. T. S. I. Caminos, Canales y Puertos
6
Roots of Equations
Determine real roots of :
Algebraic equations (including polynomials)
Transcendental equations such as f(x) = sin(x) + e-x
Combinations thereof
E. T. S. I. Caminos, Canales y Puertos
7
Roots of Equations
Engineering Economics Example:
A municipal bond has an annual payout of $1,000
for 20 years. It costs $7,500 to purchase now.
What is the implicit interest rate, i ?
1 1 i n
Solution: Present-value, PV, is: PV A
i
in which:
PV = present value or purchase price = $7,500
A = annual payment = $1,000/yr
n = number of years = 20 yrs
i
= interest rate = ? (as a fraction, e.g., 0.05
= 5%)
E. T. S. I. Caminos, Canales y Puertos
8
Roots of Equations
Engineering Economics Example (cont.):
We need to solve the equation for i:
1 1 i 20
7,500 1, 000
i
Equivalently, find the root of:
1 1 i 20
0
f (i) 7,500 1,000
i
E. T. S. I. Caminos, Canales y Puertos
9
Roots of Equations
Engineering Economics Example
$10,000
f(i) = PV - g(i)
$5,000
$0
0%
20%
40%
60%
80%
100%
-$5,000
-$10,000
-$15,000
Excel
Interest Rate, i
E. T. S. I. Caminos, Canales y Puertos
10
Roots of Equations
Graphical methods:
Determine the friction coefficient c necessary for a
parachutist of mass 68.1 kg to have a speed of 40
m/seg at 10 seconds.
Reorganizing.
E. T. S. I. Caminos, Canales y Puertos
11
Roots of Equations
Two Fundamental Approaches
1. Bracketing or Closed Methods
- Bisection Method
- False-position Method (Regula falsi).
2. Open Methods
- Newton-Raphson Method
- Secant Method
- Fixed point Methods
E. T. S. I. Caminos, Canales y Puertos
12
Bracketing Methods
f(x)
a)
f(x)
b)
f(x)
c)
f(x)
xl
d)
xu
In Figure a) we have the case of f(xl)
and f(xu) with the same sign, and there
is no root in the interval (xl,xu).
In Figure b) we have the case of f(xl)
and f(xu) With different sign, and there
is a root in the interval (xl,xu).
In Figure c) we have the case of f(xl)
and f(xu) with the same sign, and there
are two roots.
In Figure d) we have the case of f(xl)
and f(xu) with different sign, and there is
an odd number of roots.
E. T. S. I. Caminos, Canales y Puertos
13
Bracketing Methods
Though the cases above are generally valid, there are
cases in which they do not hold.
In Figure a) we have the case of f(xl)
and f(xu) with different sign, but there is
a double root.
f(x)
f(x)
f(x)
xl
a)
b)
c)
xu
In Figure b) We have the case of f(xl)
and f(xu) With different sign, but there
are two discontinuities.
In Figure c) we have the case of f(xl)
and f(xu) with the same sign, but there
is a multiple root.
E. T. S. I. Caminos, Canales y Puertos
14
Bracketing Methods (Bisection method)
Bisection Method
f(x)
xl xu
xr
2
f(xu)
f(xu)
f(x1) f(xr) > 0
(x1)
f(xr)
f(x1)
f(x)
(xu) x
xr => x1
f(xr)
(x1)
(xr)(xu) x
f(x1)
E. T. S. I. Caminos, Canales y Puertos
15
Bracketing Methods (Bisection method)
Bisection Method
Advantages:
1. Simple
2. Estimate of maximum error:
3. Convergence guaranteed
E max
xl xu
1
Eimax
0.5 Eimax
Disadvantages:
1. Slow
2. Requires two good initial estimates which define
an interval around root:
use graph of function,
incremental search, or
trial & error
E. T. S. I. Caminos, Canales y Puertos
16
Bracketing Methods (False-position Method)
False-position Method
f(x)
f(x)
xr xu
f (x u )(x1 x u )
f (x1 ) f (x u )
(xr)
f(xu)
f(x1) f(xr) > 0
(x1)
(xu)
f(xr)
f(x1)
x1 = x r
(x1)
f(x1)
f(xu)
(xu)
f(xr)
E. T. S. I. Caminos, Canales y Puertos
17
Bracketing Methods (False-position Method)
There are some cases in which the false position
method is very slow, and the bisection method gives
a faster solution.
E. T. S. I. Caminos, Canales y Puertos
18
Bracketing Methods (False-position Method)
Summary of False-Position Method:
Advantages:
1. Simple
2. Brackets the Root
Disadvantages:
1. Can be VERY slow
2. Like Bisection, need an initial interval around the
root.
E. T. S. I. Caminos, Canales y Puertos
19
Open Methods
Roots of Equations - Open Methods
Characteristics:
1. Initial estimates need not bracket the root
2. Generally converge faster
3. NOT guaranteed to converge
Open Methods Considered:
- Fixed-point Methods
- Newton-Raphson Iteration
- Secant Method
E. T. S. I. Caminos, Canales y Puertos
20
Roots of Equations
Two Fundamental Approaches
1. Bracketing or Closed Methods
- Bisection Method
- False-position Method
2. Open Methods
- One Point Iteration
- Newton-Raphson Iteration
- Secant Method
E. T. S. I. Caminos, Canales y Puertos
21
Open Methods (Newton-Raphson Method)
Newton-Raphson Method:
Geometrical Derivation:
Slope of tangent at xi is
f (x i ) 0
f '(x i )
x i x i 1
Solve for xi+1:
f (x i )
x i1 x i
f '(x i )
[Note that this is the same form as the generalized onepoint iteration, xi+1 = g(xi)]
E. T. S. I. Caminos, Canales y Puertos
22
Open Methods (Newton-Raphson Method)
Newton-Raphson Method
f(x)
Tangent w/slope=f '(xi )
f (x i ) 0
f '(x i )
x i x i 1
f '(x i )
f(xi)
f (x i ) 0
x i x i 1
f(xi)
f(xi+1)
xi+1 xi
f(x)
f(xi+1)
x
x i 1 x i
f (x i )
f '(x i )
(xi)
xi+1
xi = xi+1
E. T. S. I. Caminos, Canales y Puertos
23
Open Methods
a) Inflection point in the
neighboor of a root.
b) Oscilation in the neighboor
of a maximum or minimum.
c) Jumps in functions with
several roots.
d) Existence of a null
derivative.
E. T. S. I. Caminos, Canales y Puertos
24
Open Methods (Newton-Raphson Method)
Bond Example:
To apply Newton-Raphson method to:
1 1 i 20
0
f (i) 7,500 1, 000
i
We need the derivative of the function:
1, 000
f '(i)
1 1 i 20
21
20 1 i
i
E. T. S. I. Caminos, Canales y Puertos
25
Open Methods (Secant Method)
Secant Method
Approx. f '(x) with backward
FDD:
Substitute this into the N-R
equation:
f (x i1 ) f (x i )
f '(x)
x i1 x i
x i1 x i
f (x i )
f '(x i )
f (x i )(x i 1 x i )
x i1 x i
f (x i1 ) f (x i )
to obtain the iterative
expression:
E. T. S. I. Caminos, Canales y Puertos
26
Open Methods (Secant Method)
Secant Method
f(x)
f(x)
f '(x i )
f(x i-1 ) - f(x i )
x i-1 - x i
f '(x i )
f(x i-1 ) - f(x i )
x i-1 - x i
f(xi-1)
f(xi)
f(xi)
x
xi+1 xi xi-1
x i 1 x i
xi+1 xi
f(xi-1)
xi-1
f (x i )(x i 1 x i )
f (x i 1 ) f (x i )
xi = xi+1
E. T. S. I. Caminos, Canales y Puertos
27
