Exercise 3.  Prove that Newton's method always works for polynomials of degree  1  (functions of the form  [Graphics:Images/JuliaMandelbrotModHome_gr_392.gif],  where  [Graphics:Images/JuliaMandelbrotModHome_gr_393.gif].  

How many iterations are necessary before Newton's method produces the solution  [Graphics:Images/JuliaMandelbrotModHome_gr_394.gif]  to  [Graphics:Images/JuliaMandelbrotModHome_gr_395.gif] ?  

Solution 3.

See text and/or instructor's solution manual.

        For   [Graphics:../Images/JuliaMandelbrotModHome_gr_396.gif],   we have   [Graphics:../Images/JuliaMandelbrotModHome_gr_397.gif],   and

if our initial guess is  [Graphics:../Images/JuliaMandelbrotModHome_gr_398.gif],   then  

                    [Graphics:../Images/JuliaMandelbrotModHome_gr_399.gif].  

But this is the solution to the equation  [Graphics:../Images/JuliaMandelbrotModHome_gr_400.gif],  so our iteration either stops here or with  [Graphics:../Images/JuliaMandelbrotModHome_gr_401.gif]  if by chance we had set  [Graphics:../Images/JuliaMandelbrotModHome_gr_402.gif].  

We are done.   

Caveat.   It seems that Newton's method will also work for some polynomials of degree greater than  1.  

Aside.  We can use Mathematica to explore this situation.  

Extra Example. (For a Polynomial)  Iteration using a polynomial of degree greater than  1.

Consider using Newton's method to find the complex roots of the polynomial   

                    [Graphics:../Images/JuliaMandelbrotModHome_gr_403.gif].

Let us explore what happens when we use the starting value  [Graphics:../Images/JuliaMandelbrotModHome_gr_404.gif].

[Graphics:../Images/JuliaMandelbrotModHome_gr_405.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_406.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_407.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_408.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_409.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_410.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_411.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_412.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_413.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_414.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_415.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_416.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_417.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_418.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_419.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_420.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_421.gif]

Newton's iteration is converging to the complex root   [Graphics:../Images/JuliaMandelbrotModHome_gr_422.gif].

                    [Graphics:../Images/JuliaMandelbrotModHome_gr_423.gif]

                    Some iterates for Newton's method for  [Graphics:../Images/JuliaMandelbrotModHome_gr_424.gif],

                    starting with [Graphics:../Images/JuliaMandelbrotModHome_gr_425.gif]  we have  
          
                    [Graphics:../Images/JuliaMandelbrotModHome_gr_426.gif]
          
                    This sequence is converging to the complex root   [Graphics:../Images/JuliaMandelbrotModHome_gr_427.gif].

We are really done.   

Consider using Newton's method to find the complex roots of the polynomial   [Graphics:../Images/JuliaMandelbrotModHome_gr_428.gif].

Let us explore what happens when we use the starting value  [Graphics:../Images/JuliaMandelbrotModHome_gr_429.gif].

[Graphics:../Images/JuliaMandelbrotModHome_gr_430.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_431.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_432.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_433.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_434.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_435.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_436.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_437.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_438.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_439.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_440.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_441.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_442.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_443.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_444.gif]

Newton's iteration is converging to the complex root   [Graphics:../Images/JuliaMandelbrotModHome_gr_445.gif].

                    [Graphics:../Images/JuliaMandelbrotModHome_gr_446.gif]

                    Some iterates for Newton's method for  [Graphics:../Images/JuliaMandelbrotModHome_gr_447.gif],

                    starting with [Graphics:../Images/JuliaMandelbrotModHome_gr_448.gif]  we have
          
                    [Graphics:../Images/JuliaMandelbrotModHome_gr_449.gif]

                    This sequence is converging to the complex root   [Graphics:../Images/JuliaMandelbrotModHome_gr_450.gif].

We are really really done.   

        Let us explore what happens when we use the starting value  [Graphics:../Images/JuliaMandelbrotModHome_gr_451.gif].

[Graphics:../Images/JuliaMandelbrotModHome_gr_452.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_453.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_454.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_455.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_456.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_457.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_458.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_459.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_460.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_461.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_462.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_463.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_464.gif]


[Graphics:../Images/JuliaMandelbrotModHome_gr_465.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_466.gif]

In this case we have chosen an unlucky starting value.  The iteration produces an oscillating sequence that is divergent.

                    [Graphics:../Images/JuliaMandelbrotModHome_gr_467.gif]

                    Some iterates for Newton's method for  [Graphics:../Images/JuliaMandelbrotModHome_gr_468.gif],

                     starting with   [Graphics:../Images/JuliaMandelbrotModHome_gr_469.gif]   we have

                              [Graphics:../Images/JuliaMandelbrotModHome_gr_470.gif]  

                    This oscillating sequence is divergent.

We are really really really done.   

Aside.  We can use Mathematica to find the roots.

[Graphics:../Images/JuliaMandelbrotModHome_gr_471.gif]

[Graphics:../Images/JuliaMandelbrotModHome_gr_472.gif]

We are really really really  really done.   

        It is tempting to mimic the proof of Newton's method for real numbers and try to make a proof

for complex numbers which would show that the Newton iteration will construct a convergent sequence.  

But the proof for real numbers uses the Mean Value Theorem for derivatives.  In Exercise 13 in Section 3.1

we asked you to show that the Mean Value Theorem for derivatives does not extend to complex functions.  

So a proof that  [Graphics:../Images/JuliaMandelbrotModHome_gr_473.gif]  converges is not straight forward and is beyond the scope of this course.

Caveat.  

Newton's method can sometimes be successful in locating the real solution to  [Graphics:../Images/JuliaMandelbrotModHome_gr_474.gif],  

even though the starting value  [Graphics:../Images/JuliaMandelbrotModHome_gr_475.gif]  is a complex number,

but proving why it will work or when it will work is not an easy task.    

 

Section 3.1, Exercise 13.

 

 


























 

This solution is complements of the authors.

 

 


























 


























 

(c) 2008 John H. Mathews, Russell W. Howell