TMC: Nonlinear Equation

Presenter: Dr. Ha Viet Uyen Synh.

 The Mind map of Engineering Mathematics
4

6
3 5

2
1
9 7

8
Motivation
•Why?
− b  b 2
− 4ac
ax 2 + bx + c = 0  x =
2a

•But
ax5 + bx 4 + cx 3 + dx 2 + ex + f = 0  x = ?
sin x + x = 0  x = ?

3
 Nonlinear Equation
Solvers

Bracketing Graphical Open Methods

Bisection Newton Raphson

False Position Secant

All Iterative
4
PART A

BRACKETING METHODS
Bracketing Methods
(Or, two point methods for finding roots)

Two initial guesses for the root are required.

These guesses must “bracket” or be on either side
of the root.

If one root of a real and continuous function,

f(x)=0, is bounded by values x=xl, x=xu then
f(xl) . f(xu) <0.

6
Parts (a) and (c) indicate that if
both f (x ) and f (x ) have the
l u

same sign, either there will be

no roots or there will be an even
number of roots within the
interval.

Parts (b) and (d ) indicate that if

the function has different signs
at the end points, there will be
an odd number of roots in the
interval

7
1. The Bisection Method
For the arbitrary equation of one variable, f(x)=0

1. Pick xl and xu such that they bound the root of interest,

check if f(xl).f(xu) <0.

2. Estimate the root by evaluating f[(xl+xu)/2].

3. Find the pair

• If f(xl). f[(xl+xu)/2]<0, root lies in the lower interval, then
xu=(xl+xu)/2 and go to step 2.
• If f(xl). f[(xl+xu)/2]>0, root lies in the upper interval, then
xl= [(xl+xu)/2, go to step 2.
• If f(xl). f[(xl+xu)/2]=0, then root is (xl+xu)/2 and terminate.

4. Compare 𝜀𝑠 with 𝜀𝑎

5. If 𝜀𝑎 < 𝜀𝑠, stop. Otherwise repeat the process.

8
Termination Criteria and Error Estimates
An approximate percent relative error εa can be
calculated
𝑥𝑟𝑛𝑒𝑤 − 𝑥𝑟𝑜𝑙𝑑
𝜀𝑎 = 𝑛𝑒𝑤
100%
𝑥𝑟
where
xrnew is the root for the present iteration and
xrold is the root from the previous iteration.

When εa becomes less than a pre-specified stopping

criterion εs, the computation is terminated.
Example #1
Determine the root of the following equation
667.38
𝑓 𝑐 = 1 − 𝑒 −0.146843𝑐 − 40
𝑐
Example #1 (cont.)
The first step in bisection is to guess two values of the unknown that give values for f
(c) with different signs. From the above figure, we can see that the function changes
sign between values of 12 and 16. Therefore, the initial estimate of the root xr lies at
the midpoint of the interval
12 + 16
𝑥𝑟 = = 14
2
Next we compute the product of the function value at the lower bound and at the
midpoint:
f(12) f(14) = 6.067(1.569) = 9.517
Consequently, the root must be located between 14 and 16. The method can be
repeated until the result is accurate enough to satisfy your needs.

Thus, after six iterations εa finally falls below εs = 0.5%, and the computation can be
terminated.
The False-Position Method
(Regula-Falsi)
If a real root is bounded by xl and xu of f(x)=0, then we can
approximate the solution by doing a linear interpolation
between the points [xl, f(xl)] and [xu, f(xu)] to find the xr value
such that l(xr)=0, l(x) is the linear approximation of f(x).

12
Procedure
1. Find a pair of values of x, xl and xu such that fl=f(xl) <0
and fu=f(xu) >0.
2. Estimate the value of the root from the following
formula
xl f u − xu f l
xr =
fu − fl
and evaluate f(xr).
3. Use the new point to replace one of the original
points, keeping the two points on opposite sides of
the x axis.
If f(xr)<0 then xl=xr == > fl=f(xr)
If f(xr)>0 then xu=xr == > fu=f(xr)
If f(xr)=0 then you have found the root and need go no further!

4. See if the new xl and xu are close enough for

convergence to be declared. If they are not, go back
to step 2.

13
Example #2
Use the false-position method to determine the root of the same equation
investigated in Example #1

As in Example #1, initiate the computation with guesses of xl =12 and xu =16.
PART B

OPEN METHODS
Open Methods
Open methods are based on formulas that require only a
single starting value of x or two starting values that do not
necessarily bracket the root.

16
Simple Fixed-point Iteration
Rearrange the function so that x is on the left side of
the equation:

f ( x) = 0  g ( x) = x
xk = g ( xk −1 ) xo given , k = 1, 2, ...

•Bracketing methods are “convergent”.

•Fixed-point methods may sometime “diverge”,
depending on the stating point (initial guess) and how
the function behaves.

17
Example :
f ( x) = x 2 − x − 2 x0
g ( x) = x 2 − 2
or
g ( x) = x + 2
or
2
g ( x) = 1 +
x


18
Example #2
Use simple fixed-point iteration to locate the root of
f (x) = e−x − x.
The function can be separated as
xi+1 = e−xi
Starting with an initial guess of x0 = 0, this iterative
equation can be applied to compute

Thus, each iteration brings the estimate closer to the true

value of the root: 0.56714329.
Convergence
x=g(x) can be expressed as a pair of equations:
y1=x
y2=g(x) (component equations)
Plot them separately.

20
Conclusion
Fixed-point iteration converges if
g ( x)  1 (slope of the line f(x) = x)

When the method converges, the error is roughly

proportional to or less than the error of the previous step,
therefore it is called “linearly convergent.”

21
1. Newton-Raphson Method
Most widely used method.

Based on Taylor series expansion:

x 2
f ( xi +1 ) = f ( xi ) + f ( xi )x + f ( xi ) + Ox 3
2!
The root is the value of x i +1 when f(x i +1 ) = 0
Rearrangin g,
0 = f(x i ) + f (x i )( xi +1 − xi )
f ( xi )
xi +1 = xi −
f ( xi )

22
A convenient method for functions whose derivatives can
be evaluated analytically. It may not be convenient for
functions whose derivatives cannot be evaluated
analytically.

23
Example #3
Use the Newton-Raphson method to estimate the root of f (x) =e−x − x, employing an
initial guess of x0 = 0.
Four cases where the Newton-
Raphson method exhibits poor
convergence

25
2. Secant Method
A slight variation of Newton’s method for functions whose
derivatives are difficult to evaluate. For these cases the
derivative can be approximated by a backward finite
divided difference.

f ( xi ) − f ( xi −1 )
f ( xi ) 
xi − xi −1
xi − xi −1
xi +1 = xi − f ( xi ) i = 1,2,3, 
f ( xi ) − f ( xi −1 )

26
•Requires two initial estimates of x , e.g, xo, x1.
However, because f(x) is not required to change signs
between estimates, it is not classified as a “bracketing”
method.
•The scant method has the same properties as
Newton’s method. Convergence is not guaranteed for all
xo, f(x).

27
28
Example #4
Use the secant method to estimate the root of f(x)=e−x − x . Start with initial
estimates of x−1 = 0 and x0 = 1.0.
Multiple Roots
“Multiple root” corresponds to a point where a function is
tangent to the x axis.

Difficulties
Function does not change sign at the multiple root, therefore,
cannot use bracketing methods.
Both f(x) and f′(x)=0, division by zero with Newton’s and Secant
methods.

30
Multiple Roots

None of the methods deal with multiple roots efficiently,

however, one way to deal with problems is as follows:
This function has roots at
f ( xi )
Set u ( xi ) = all the same locations as
f ( xi ) the original function
u ( xi )
Then find x i +1 = xi −
u ( xi )
f ' ( x) f ' ( x) − f ( x) f " ( x)
u ' ( x) =
 f ' ( x)2
f ( xi ) f ' ( xi )
= x i +1 = xi −
 f ' ( x)2 − f ( x) f " ( x)

31
Example #4
Use both the standard and modified Newton-Raphson methods to evaluate the multiple root
of f(x) = (x − 3)(x − 1)(x − 1), with an initial guess of x0 = 0.

The first derivative of f(x) is f’(x) = 3x2 − 10x + 7, and therefore, the standard Newton-Raphson
method for this problem is
Systems of Linear Equations
To this point, we have focused on the determination of
the roots of a single equation. A related problem is to
locate the roots of a set of simultaneous equations,

f1 ( x1 , x2 , x3 , , xn ) = 0
f 2 ( x1 , x2 , x3 , , xn ) = 0

f n ( x1 , x2 , x3 , , xn ) = 0

33
Fixed-Point Iteration
Use fixed-point iteration to determine the roots of
u(x, y) = x2 + xy − 10 = 0
v(x, y) = y + 3xy2 − 57 = 0
Note that a correct pair of roots is x = 2 and y = 3. Initiate the computation with
guesses of x = 1.5 and y = 3.5.

We have
𝑥= 10 − 𝑥𝑦
57−𝑦
𝑦=
3𝑥

Now the results are more satisfactory:

𝑥= 10 − 1.5(3.5) = 2.17945
57−3.5
𝑦= = 2.86051
3(2.17945)

𝑥= 10 − 2.17945(2.86051) = 1.94053
57−2.86051
𝑦= = 3.04955
3(1.94053)

Thus, the approach is converging on the true values of x = 2 and y = 3.

Newton-Raphson
Taylor series expansion of a function of more than one
variable
ui u
ui +1 = ui + ( xi +1 − xi ) + i ( yi +1 − yi )
x y
vi vi
vi +1 = vi + ( xi +1 − xi ) + ( yi +1 − yi )
x y
The root of the equation occurs at the value of x and
y where ui+1 and vi+1 equal to zero.
ui u u u
xi +1 + i yi +1 = −ui + xi i + yi i
x y x y
vi vi vi vi
xi +1 + yi +1 = −vi + xi + yi
x y x y
A set of two linear equations with two unknowns that
can be solved for.
35
 vi ui
ui − vi
y y
xi +1 = xi −
ui vi ui vi
−
x y y x
vi ui
ui − vi
yi +1 = yi − x x
ui vi ui vi
−
x y y x

Determinant of the
Jacobian of the
system.
36
Example #5
Use the multiple-equation Newton-Raphson method to determine roots of
u(x, y) = x2 + xy − 10 = 0
v(x, y) = y + 3xy2 − 57 = 0
Note that a correct pair of roots is x = 2 and y = 3. Initiate the computation with guesses of x =
1.5 and y = 3.5.
PART C

ROOTS OF POLYNOMIALS
Roots of Polynomials
The roots of polynomials such as

f n ( x) = ao + a1 x + a2 x 2 +  + an x n

Follow these rules:

1.For an nth order equation, there are n real or complex
roots.
2.If n is odd, there is at least one real root.
3.If complex root exist in conjugate pairs (that is, l+mi
and l-mi), where i=sqrt(-1).

40
Conventional Methods
The efficacy of bracketing and open methods depends
on whether the problem being solved involves complex
roots. If only real roots exist, these methods could be
used.

However, finding good initial guesses complicates both

the open and bracketing methods, also the open
methods could be susceptible to divergence.

A special method has been developed to find the real

and complex roots of polynomials – Müller method.

41
Müller Method
Müller’s method obtains a root estimate by
projecting a parabola to the x axis through three
function values.

42
 Müller Method
The method consists of deriving the coefficients of
parabola that goes through the three points:

1. Write the equation in a convenient form:

f 2 ( x) = a( x − x2 ) 2 + b( x − x2 ) + c

A comparison of two related

approaches for locating roots:
(a) the secant method and
(b) Müller’s method. 43
2. The parabola should intersect the three points
[xo, f(xo)], [x1, f(x1)], [x2, f(x2)]. The coefficients
of the polynomial can be estimated by
substituting three points to give
f ( xo ) = a( xo − x2 ) 2 + b( xo − x2 ) + c
f ( x1 ) = a( x1 − x2 ) 2 + b( x1 − x2 ) + c
f ( x 2 ) = a ( x2 − x 2 ) 2 + b ( x 2 − x 2 ) + c
3. Three equations can be solved for three
unknowns, a, b, c. Since two of the terms in the
3rd equation are zero, it can be immediately
solved for c=f(x2).

f ( xo ) − f ( x2 ) = a ( xo − x2 ) 2 + b( xo − x2 )
f ( x1 ) − f ( x2 ) = a ( x1 − x2 ) 2 + b( x1 − x2 )
44
If

h o = x1 - x o h1 = x 2 - x1

f ( x1 ) − f ( x o ) f ( x 2 ) − f ( x1 )
o = 1 =
x1 − x o x 2 − x1

We have

x 2 − x 0 = h1 + h 0

 0 h 0 = f ( x1) − f ( x 0 )

 1h1 = f ( x 2 ) − f ( x1)

 f ( x 2 ) − f ( x 0 ) =  0 h 0 +  1h1

So ,
2
( ho + h1 ) b − ( ho + h1 ) a = ho  o + h11
2
h1b − h1 a = h11

The result can be summarized as

1 −  o
a= b = ah1 + 1 c = f ( x2 )
h1 + ho 45
Roots can be found by applying an alternative form of
quadratic formula:

− 2c
x3 = x2 +
b  b 2 − 4ac
The error can be calculated as

x3 − x2
a = 100%
x3
±term yields two roots, the sign is chosen to agree with
b. This will result in a largest denominator, and will give
root estimate that is closest to x2.

46
 Once x3 is determined, the process is repeated
using the following guidelines:
1. If only real roots are being located, choose the
two original points that are nearest the new
root estimate, x3.
2. If both real and complex roots are estimated,
employ a sequential approach just like in
secant method, x1, x2, and x3 to replace xo,
x1, and x2.

47
Example #6
Use Müller’s method with guesses of x0, x1, and x2 = 4.5, 5.5, and 5, respectively, to
determine a root of the equation f(x) = x3 − 13x − 12
Note that the roots of this equation are −3, −1, and 4.
Polynomial order reduction
It is based on dividing a polynomial by a factor (x-t)
Any Questions?

 [email protected]