Re: New Root Seeking Algorithms Hans de Moor Unregistered Posts: 9 Threads: 3 Joined: Jun 2007 11-13-2007, 11:33 PM Namir, Your comparison with Newton's method is wrong. Newton's method requires that you use f '(x) explicitly. If you don't, what you are doing is really measuring the convergence of fixed point iteration using a function that approximates f '(x). Simply stated, you are comparing your method to the accuracy of the approximation of f '. The last example you give is for a polynomial with 6 unique roots. If you use the derivative explicitly, an initial guess of 7 will give the root to 10 significant digits in 6 iterations: six calls of f and six calls of f '. If you do one more iteration and you are using sufficient significant digits, you will have the root to 20 significant digits. Your example fails because the inaccuracy of f ' slows the convergence rate. One of shortcomings of Newton's method is that you have to know and call the derivative. If you don't know the derivative and only want to call the function itself, use the secant method. This is closer to what you are doing but the error in f ' is reduced with every iteration. I suggest you look at the fixed point iteration theorem and calculate theoretical convergence. This will give you a better picture of how your method behaves. Hans hugh steers Unregistered Posts: 536 Threads: 56 Joined: Jul 2005 11-14-2007, 02:38 PM hi, i understand secant dominates when f' is numerical. is this true? « Next Oldest | Next Newest »

 Possibly Related Threads… Thread Author Replies Views Last Post [HP Prime] Using n-root symbol and exponent problem uklo 7 2,787 11-11-2013, 01:39 AM Last Post: Alberto Candel OT--TI-36X Algorithms Matt Agajanian 48 10,861 09-01-2013, 08:13 PM Last Post: robert rozee Cubic root (-8) = 2 ? Gilles Carpentier 37 9,881 08-12-2013, 10:26 PM Last Post: jep2276 HELP WANTED ON ALGORITHMS Joerg Woerner 19 4,744 04-27-2013, 12:56 PM Last Post: Eric Smith Square Root Simplifier for HP39gII Mic 4 1,909 03-11-2013, 08:25 AM Last Post: Eddie W. Shore Seeking DIY4X update David Griffith 10 2,743 11-15-2012, 12:41 AM Last Post: Eric Smith Cube root on standard calculator Thomas Klemm 22 6,132 11-09-2012, 06:11 AM Last Post: Pierre Seeking new resources for TI NSpire Programming Namir 6 2,233 07-29-2012, 03:57 AM Last Post: fhub ROOT bug? HP 48S/48G Eddie W. Shore 8 3,727 07-13-2012, 07:05 PM Last Post: Eddie W. Shore OT: primitive mult/div algorithms Egan Ford 9 2,627 05-27-2012, 09:27 PM Last Post: Egan Ford

Forum Jump: