Chapter 6
Nonlinear Equations
6.1 The Problem of Nonlinear Root-finding
In this module we consider the problem of using numerical techniques to find the roots of nonlinear
x = 0 equations, f . Initially we examine the case where the nonlinear equations are a scalar
function of a single independent variable, x. Later, we shall consider the more difficult problem
n f x= 0
where we have a system of n nonlinear equations in independent variables, .
x = 0 Since the roots of a nonlinear equation, f , cannot in general be expressed in closed form, we must use approximate methods to solve the problem. Usually, iterative techniques are
used to start from an initial approximation to a root, and produce a sequence of approximations
0 1 2
x ;x ;x ;::: ;
which converge toward a root. Convergence to a root is usually possible provided that the function
x f is sufficiently smooth and the initial approximation is close enough to the root.
6.2 Rate of Convergence
Let
0 1 2
x ;x ;x ;:::
c Department of Engineering Science, University of Auckland, New Zealand. All rights reserved. July 11, 2002
70 NONLINEAR EQUATIONS
k k
= x p be a sequence which converges to a root ,andlet . If there exists a number and
a non-zero constant c such that
k +1
lim = c;
p (6.1)
k
k !1
j j
p = 1; 2; 3 then p is called the order of convergence of the sequence. For the order is said to be linear, quadratic and cubic respectively.
6.3 Termination Criteria
Since numerical computations are performed using floating point arithmetic, there will always be
a finite precision to which a function can be evaluated. Let denote the limiting accuracy of the nonlinear function near a root. This limiting accuracy of the function imposes a limit on the
accuracy, , of the root. The accuracy to which the root may be found is given by:
: =
(6.2)
0
jf j
This is the best error bound for any root finding method. Note that is large when the first derivative
0
f j of the function at the root, j , is small. In this case the problem of finding the root is ill-
conditioned. This is shown graphically in Figure 6.1.
f x
f x
2 2
1 2
2
1
2
2
x
0 1
large 2 small
f 1 ill-conditioned well-conditioned
FIGURE 6.1:
When the sequence
0 1 2
x ;x ;x ;::: ;
k k 1 k
x x converges to the root, , then the differences x will decrease until .
c Department of Engineering Science, University of Auckland, New Zealand. All rights reserved. July 11, 2002 6.4 BISECTION METHOD 71
With further iterations, rounding errors will dominate and the differences will vary irregularly.
k
The iterations should be terminated and x be accepted as the estimate for the root when the
following two conditions are satisfied simultaneously:
k +1 k k k 1
x x x
1. x and (6.3)
k k 1
x x
:
2. < (6.4)
k
1+jx j
k
x
is some coarse tolerance to prevent the iterations from being terminated before is close to ,
k k 1