Newton's method and complex roots

Click For Summary
Newton's method can be applied to find complex roots by substituting the variable with a complex number, z = x + yi, and using the function in complex form. The function should be expressed as p(z) = z^8 - 7z^7 + 14z^6 - 14z^5 + 27z^4 - 14z^3 + 14z^2 - 7z + 1, and the Newton iteration formula is z_{k+1} = z_k - p(z_k)/p'(z_k). To successfully converge to a non-real root, it is essential to start with a non-real initial value, as certain starting points can lead to chaotic limit cycles instead of stable roots. Additional resources and articles are available online to further explore Newton's method for complex roots, including practical implementations. Understanding these concepts is crucial for effectively applying Newton's method to complex functions.
Delta31415
Messages
90
Reaction score
8
1) the problem
I understand Newton's method and I was able to find all the real roots of the function.However, I don't understand how to find the complex roots. I know that z=x+yi, and that I can plug in z for the formula. However I, don't know how to change the function ( x^8-7x^7+14x^6-14x^5+27x^4-14x^3+14x^2-7x+1) to into complex form. because of that I don't know what the initial value would.

2) Relevant equations
Newton's method and z=x+yi)
3) attempts
found real roots

btw the function I wrote is the one I am solving,however that doesn't matter I just need to know how to find the complex roots for any function given in that form.



 
Physics news on Phys.org
You are not very verbose with your problem statement :rolleyes: but I assume you want to solve ##f({\bf z}) = 0## for ##{\bf z} \in {\mathbb C}## ?

That's in fact two equations (e.g ## \ \operatorname { Re } (f( {\bf z} ) ) = 0 \ ## and ## \ \operatorname {Im} (f( {\bf z} ) )= 0 ## ) and you'll have to make those explicit before embarking on Newton's method for finding zeroes of a function.
 
Delta31415 said:
1) the problem
I understand Newton's method and I was able to find all the real roots of the function.However, I don't understand how to find the complex roots. I know that z=x+yi, and that I can plug in z for the formula. However I, don't know how to change the function ( x^8-7x^7+14x^6-14x^5+27x^4-14x^3+14x^2-7x+1) to into complex form. because of that I don't know what the initial value would.

2) Relevant equations
Newton's method and z=x+yi)
3) attempts
found real roots

btw the function I wrote is the one I am solving,however that doesn't matter I just need to know how to find the complex roots for any function given in that form.



Just replace ##p(x) = x^8-7 x^7+14 x^6 -14 x^5 +27 x^4 -14 x^3 +14 x^2 -7 x +1## by ##p(z) = z^8-7 z^7+14 z^6 -14 z^5 +27 z^4 -14 z^3 +14 z^2 -7 z +1## and use the usual Newton iteration
$$z_{k+1} = z_k - \frac{p(z_k)}{p'(z_k)}.$$
However, to converge to a non-real root you need to start at a non-real ##z##, and even then there are some starting points that can yield chaotic "limit cycles" rather than single points. Google "Newtons method for complex roots" to find numerous articles on the issue.
 
  • Like
Likes BvU and Delta31415
Ray Vickson said:
Just replace ##p(x) = x^8-7 x^7+14 x^6 -14 x^5 +27 x^4 -14 x^3 +14 x^2 -7 x +1## by ##p(z) = z^8-7 z^7+14 z^6 -14 z^5 +27 z^4 -14 z^3 +14 z^2 -7 z +1## and use the usual Newton iteration
$$z_{k+1} = z_k - \frac{p(z_k)}{p'(z_k)}.$$
However, to converge to a non-real root you need to start at a non-real ##z##, and even then there are some starting points that can yield chaotic "limit cycles" rather than single points. Google "Newtons method for complex roots" to find numerous articles on the issue

Thank, I was able to find the complex roots and I also looked up chaotic limit cycle the graph isn't a circle or a spiral so there isn't any orbit.right?
 
Last edited by a moderator:
Delta31415 said:
Thank, I was able to find the complex roots and I also looked up chaotic limit cycle the graph isn't a circle or a spiral so there isn't any orbit.right?
I don't know if it is right or not; perhaps I should have termed it differently, but the articles on the subject used such terms--possibly more accurately. As I said, there are numerous articles on the topic, and they will tell you everything you need to know. Some of them even include Matlab codes for the problem.
 
Last edited:
Question: A clock's minute hand has length 4 and its hour hand has length 3. What is the distance between the tips at the moment when it is increasing most rapidly?(Putnam Exam Question) Answer: Making assumption that both the hands moves at constant angular velocities, the answer is ## \sqrt{7} .## But don't you think this assumption is somewhat doubtful and wrong?

Similar threads

  • · Replies 22 ·
Replies
22
Views
966
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 6 ·
Replies
6
Views
5K
  • · Replies 7 ·
Replies
7
Views
1K
  • · Replies 17 ·
Replies
17
Views
3K
  • · Replies 6 ·
Replies
6
Views
3K
Replies
16
Views
7K
Replies
4
Views
2K