next up previous
Next: Verification, Validation and Certification Up: Fixed point theorem and Previous: Application to solution of

Application to Newton-Raphson method

Consider the problem of finding the roots of an equation

\begin{displaymath}
f(x) = 0
\end{displaymath}

If an approximation to the root $x_n$ is already available then the next approximation

\begin{displaymath}
x_{n+1}= x_n + \Delta x_n
\end{displaymath}

is calculated by requiring that

\begin{displaymath}
f(x_n + \Delta x_n) = 0
\end{displaymath}

Using Taylors formula and assuming that $f$ is differentiable we can approximate the left hand side,

\begin{displaymath}
f(x_n) + \Delta x_n f^\prime(x_n) = 0    \Longrightarrow   \
\Delta x_n = - \frac{f(x_n)}{ f^\prime(x_n) }
\end{displaymath}

so that the new approximation is given by

\begin{displaymath}
x_{n+1} = x_n - \frac{f(x_n)}{f^\prime(x_n) }
\end{displaymath}

This defines the iterative scheme for Newton-Raphson method and the solution if it exists, is the fixed point of

\begin{displaymath}
F(x) = x - \frac{ f(x) }{ f^\prime(x) }
\end{displaymath}

If we assume that
\begin{displaymath}
\vert F^\prime(x) \vert \le \alpha < 1
\end{displaymath} (11)

in some neighbourhood of the solution then

\begin{displaymath}
F(y) - F(x) = \int^y_x F^\prime(z) \textrm{d}z    \Longrightarrow    \vert F(y) -
F(x)\vert \le \alpha \vert y-x\vert
\end{displaymath}

Thus condition (11) guarantees that $F$ is contractive and the Newton iterations will converge.

Praveen. C, last updated on 18-February-2005