Why Does Newton Raphson Method Fail for Some Functions?

  • Thread starter Thread starter relinquished™
  • Start date Start date
  • Tags Tags
    Newton
relinquished™
Messages
79
Reaction score
0
Hello. I've been approached with a problem of explaining why Newton Raphson method fails for some functions. I came across a book in Numerical Analysis (Kellison's book) that the method may fail if

1.) f'(x)=0
2.) The initial value is taken at a maximum or minumum point,
3.) The initial value is taken at a point of inflection,
4.) The initial value is taken near a maximum point and a minimum point,
5.) The initial value is taken near a point of inflection.

Now I can explain (1). Newton's Method will fail since the iteration is given by

<br /> x_{n+1} = x_n - \frac{f(x_n)}{f&#039;(x_n)}<br />

Therefore making the whole thing undefined.

As for (2), since minima and maxima have f'(x)=0, they will fail for the same reasons as (1)

As for (3), ... I'm totally clueless. I don't know how a point of inflection (f''(x)=0) could be related to the Iteration that uses f'

As for (4), I also know that choosing near a minima or maxima might make the method "oscillate", but what I'm looking for are more concrete answers; something that can be related to f'(x) or something that resembles a proof.

As for (5), I have no idea why it may fail for an initial value near a point of inflection..

All help is appreciated,

reli~
 
Last edited:
Physics news on Phys.org
I'm still in the early stages of calc myself but...

for 5 I try to explain it to myself like this: imagine you're using Newton on a function that behaves a bit like x^2 - C as it heads for the x-axis but with a couple of inflexion points near the area you choose as your guess...If for your first guess the concavity of the curve is up then the tangent at the curve will be steeper than the secant between the root you want and the value of your first guess...so as you follow the tangent you are moving in the same direction as the secant but not as fast (along the x-axis)...such that your guess takes you closer to the root, future guesses under favourable conditions such as this would get you closer still.

But now suppose that you're unlucky for the second guess. The function has changed it's concavity between the x value of our first guess and the x value of our second guess. Now it's the secant that is steeper than the tangent (unless the tangent changes it's direction), for any point on this tangent you will be above the curve...and you will hit the x-axis at the other side of your root or further away from it...from this point onwards you might find yourself attracted to a different root or worse still sent on a wild goose chase

http://img132.imageshack.us/img132/24/Newtonspy8.gif
 
Last edited by a moderator:
Also if the derivative is not continuous at the root, then convergence may fail to occur.

Indeed, let f(0) = 0 and f(x) = x + x^2\sin(2/x) elsewhere.
 
For 4/, when the initial value x0 is near the minima, or maxima, then, it may be possible that: f&#039;(x_0) \approx 0. Hence, it will make x1 considerably large: x_1 = x_0 - \frac{f(x_0)}{f&#039;(x_0)}
Not really sure about 3/, and 5/...
 
http://www.karlscalculus.org/NRbox.html makes a pretty good analysis of the method. There you will find answers to your questions.
 
Last edited by a moderator:
Prove $$\int\limits_0^{\sqrt2/4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx = \frac{\pi^2}{8}.$$ Let $$I = \int\limits_0^{\sqrt 2 / 4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx. \tag{1}$$ The representation integral of ##\arcsin## is $$\arcsin u = \int\limits_{0}^{1} \frac{\mathrm dt}{\sqrt{1-t^2}}, \qquad 0 \leqslant u \leqslant 1.$$ Plugging identity above into ##(1)## with ##u...
Back
Top