Example 1.  Use Newton's method and Muller's method to find numerical approximations to the multiple root  [Graphics:Images/MullersMethodMod_gr_6.gif]  of the function  [Graphics:Images/MullersMethodMod_gr_7.gif].  
Show details of the computations for the starting value  [Graphics:Images/MullersMethodMod_gr_8.gif].  Compare the number of iterations for the two methods.

Solution 1.

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


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

Graph the function.

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


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

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

Starting with  [Graphics:../Images/MullersMethodMod_gr_14.gif], use the Newton-Raphson method to find a numerical approximation to the root.

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



[Graphics:../Images/MullersMethodMod_gr_16.gif]
[Graphics:../Images/MullersMethodMod_gr_17.gif]
[Graphics:../Images/MullersMethodMod_gr_18.gif]
[Graphics:../Images/MullersMethodMod_gr_19.gif]
[Graphics:../Images/MullersMethodMod_gr_20.gif]
[Graphics:../Images/MullersMethodMod_gr_21.gif]
[Graphics:../Images/MullersMethodMod_gr_22.gif]
[Graphics:../Images/MullersMethodMod_gr_23.gif]
[Graphics:../Images/MullersMethodMod_gr_24.gif]

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

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

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


We can use Mathematica's Solve procedure to determine some of the roots.

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


[Graphics:../Images/MullersMethodMod_gr_29.gif]
[Graphics:../Images/MullersMethodMod_gr_30.gif]

For Newton's method, how far away is the eighth iteration  [Graphics:../Images/MullersMethodMod_gr_31.gif]  from the root  [Graphics:../Images/MullersMethodMod_gr_32.gif] ?
Note. The last iteration is actually stored in  [Graphics:../Images/MullersMethodMod_gr_33.gif].

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


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

Starting with  [Graphics:../Images/MullersMethodMod_gr_36.gif], [Graphics:../Images/MullersMethodMod_gr_37.gif], and [Graphics:../Images/MullersMethodMod_gr_38.gif], use Muller's method to find a numerical approximation to the root.

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



[Graphics:../Images/MullersMethodMod_gr_40.gif]
[Graphics:../Images/MullersMethodMod_gr_41.gif]
[Graphics:../Images/MullersMethodMod_gr_42.gif]
[Graphics:../Images/MullersMethodMod_gr_43.gif]
[Graphics:../Images/MullersMethodMod_gr_44.gif]
[Graphics:../Images/MullersMethodMod_gr_45.gif]
[Graphics:../Images/MullersMethodMod_gr_46.gif]
[Graphics:../Images/MullersMethodMod_gr_47.gif]
[Graphics:../Images/MullersMethodMod_gr_48.gif]
[Graphics:../Images/MullersMethodMod_gr_49.gif]
[Graphics:../Images/MullersMethodMod_gr_50.gif]
[Graphics:../Images/MullersMethodMod_gr_51.gif]
[Graphics:../Images/MullersMethodMod_gr_52.gif]

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

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

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

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


For Muller's method, how far away is the eighth iteration from the root  [Graphics:../Images/MullersMethodMod_gr_57.gif] ?
Note. The last iteration is actually stored in  [Graphics:../Images/MullersMethodMod_gr_58.gif].

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


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

This is closer than  [Graphics:../Images/MullersMethodMod_gr_61.gif]  which was obtained with Newton's method.

We are done.

Aside.  Compare with Mathematica's built in routine.

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


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

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

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

Mathematica's answer is not so good, need to adjust the number of iterations and the working precision.

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


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

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

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


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


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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004