Newton's Method for Root Finding - Infinite Loop

Click For Summary

Homework Help Overview

The original poster attempts to construct a function f(x) for Newton's method that results in an infinite cycle, specifically where xn = (-1)n x0, regardless of the starting value x0. The discussion revolves around the properties of Newton's method and the behavior of the function involved.

Discussion Character

  • Exploratory, Conceptual clarification, Mathematical reasoning

Approaches and Questions Raised

  • Participants explore the relationship between f(x) and its derivative, with attempts to derive conditions under which Newton's method leads to the desired infinite cycle. Questions arise regarding the formulation of f(x) and its implications for different domains of x.

Discussion Status

Some participants have suggested potential forms for f(x), including f(x) = e∫(1/2x) and f(x) = sqrt(|x|). There is acknowledgment of different approaches and a recognition of the need to consider how f(x) behaves for negative values of x. Multiple interpretations of the function's definition are being explored.

Contextual Notes

Participants note the importance of defining f(x) appropriately for x < 0 and discuss the implications of different choices for the function's formulation.

Scootertaj
Messages
97
Reaction score
0
1. Construct a function f (x) so that Newton's method gets 'hanging' in an infinite cycle xn = (-1)n x0 , no matter how the
start value x0 is chosen.




2. Homework Equations :
xn+1 = xn - f(xn) / f'(xn)



The Attempt at a Solution


xn+1 = xn - f(xn) / f'(xn) = (-1)n+1x0 = (-1)nx0 - f(xn) / f'(xn) [itex]\Rightarrow[/itex] f(xn) / f'(xn) = 2(-1)nx0
But, I don't know if that's what I want to do or what to do with it.

Any ideas?
 
Physics news on Phys.org
If you apply Newtons iteration to any value of x, you want to wind up with -x. Try solving the differential equation -x=x-f(x)/f'(x).
 
Ah, so you would get f(x) / f'(x) = 1/2x [itex]\Rightarrow[/itex] f'(x) = 1/2x * f(x) [itex]\Rightarrow[/itex] f(x) = e[itex]\int[/itex] 1/2x[itex]\Rightarrow[/itex] f(x) = [itex]\sqrt{x}[/itex] ?
 
Last edited:
Dick, you're right that worked perfectly. I had it wrong at first since I had f(x) = e-[itex]\int 1/2x[/itex] but we get the double negative and thus get f(x) = e[itex]\int 1/2x[/itex]

Thank you!
 
Scootertaj said:
Dick, you're right that worked perfectly. I had it wrong at first since I had f(x) = e-[itex]\int 1/2x[/itex] but we get the double negative and thus get f(x) = e[itex]\int 1/2x[/itex]

Thank you!

Good! You'll also need to think a little about how to define f(x) for x<0. But drawing a picture should make that pretty clear.
 
Would sqrt(|x|) work?
 
Scootertaj said:
Would sqrt(|x|) work?

Try it. I had been thinking about a different choice, but that works too.
 
Hmm, if you don't mind, what were you thinking of?
 
Scootertaj said:
Hmm, if you don't mind, what were you thinking of?

sqrt(x) for x>=0, -sqrt(-x) for x<0. The graph looks a little 'smoother'. But I don't think there's anything wrong with sqrt(|x|) either.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 6 ·
Replies
6
Views
5K
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 3 ·
Replies
3
Views
5K
  • · Replies 12 ·
Replies
12
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 12 ·
Replies
12
Views
6K
Replies
7
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K