Newton&'s Method
If
are continuous near a root
,
then this extra information regarding the nature of
can be used to develop algorithms that will produce sequences
that converge faster to
than either the bisection or false position method. The Newton-Raphson (or
simply Newton's) method is one of the most useful and best known algorithms that
relies on the continuity of
.The
method is attributed to
Sir Isaac Newton (1643-1727) and
Joseph Raphson (1648-1715).
Theorem ( Newton-Raphson Theorem ).Assume that
and there exists a number
,
where
.If,
then there exists a
such that the sequence
defined by the iteration
for
will converge to
for any initial approximation.
Algorithm ( Newton-Raphson
Iteration ).To
find a root ofgiven
an initial approximationusing
the iteration
for.
Mathematica Subroutine (Newton-Raphson Iteration).
Example.Use Newton's method to find
the three roots of the cubic polynomial.
Determine the Newton-Raphson iteration formulathat
is used.Show details of the computations for the starting value.
Solution.
Definition (Order of a Root)Assume
thatf(x)and its derivativesare
defined and continuous on an interval aboutx = p.We say thatf(x) = 0has
a root of ordermatx = pif and only
if
.
A root of orderm
= 1is often called a
simple root, and ifm
> 1it is called a
multiple root.A root of orderm
= 2is sometimes called a
double root, and so on.The next
result will illuminate these concepts.
Definition (Order of
Convergence)Assume that
converges top,and set.If
two positive constantsexist,
and
then the sequence is said to converge topwith
order of convergence R.The
numberAis called the
asymptotic error constant.The casesare
given specialconsideration.
(i)If,
the convergence ofis
called linear.
(ii)If,
the convergence ofis
called quadratic.
Theorem (Convergence Rate for
Newton-Raphson Iteration)Assume that
Newton-Raphson iteration produces a sequence
that converges to the rootpof the function.
Ifpis a simple root, then convergence is
quadratic andforksufficiently
large.
Ifpis a multiple root of orderm,then
convergence is linear andforksufficiently
large.
Reduce the volume of printout.
After you have debugged you program and it is working properly,
delete the unnecessary print statements
and
Concise Program
for the Newton-Raphson Method
Now test this subroutine using the function in Example 1.
Error Checking in the Newton-Raphson Method
Error checking can be added to the Newton-Raphson method.Here we have added
a third parameterto
the subroutine which estimate the accuracy of the numerical solution.
The following subroutine call uses a maximum of 20 iterations, just to make
sure enough iterations are performed.However, it will terminate when the
difference between consecutive iterations is less than.By
interrogatingkafterward we can see how many iterations were actually
performed.
Various
Scenarios
|