Question for Tim Wessman Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 08-31-2010, 02:15 PM Tim, What algorithm does the HP-50G use to find the roots of polynomials? Is it the Jenkins-Traub algorithm? The HP-50G seem to have a polynomial roots solver that works pretty well. Namir Tim Wessman Posting Freak Posts: 1,278 Threads: 44 Joined: Jul 2007 08-31-2010, 04:04 PM I assume you are talking about the PROOT command, or the Solve Poly interface? This is what is found in the source there: ** The algorithm is Laguerre iteration with stepsize control ** and is based upon the ZERPOL program contained in the thesis ** entitled "A Zero Finding Algorithm Using Laguerre's Method" ** by Brian T. Smith (directed by W. Kahan). ** The implementation is a modified version of the HP-71 Math ** ROM PROOT runtime execution code implemented by Laurence Grodd. TW Edited: 31 Aug 2010, 4:06 p.m. Gerson W. Barbosa Posting Freak Posts: 2,761 Threads: 100 Joined: Jul 2005 08-31-2010, 05:19 PM More information about the HP-71B Math ROM Polynomial Root Finder can be found in the July 1984 issue of Hewlett-Packard Journal, pages 33 through 36: Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 08-31-2010, 08:49 PM Thanks Gerson! Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 08-31-2010, 08:50 PM Thanks Tim ... gives me a good insight of how the 50G (and other graphing machines) handle getting the roots of polynomials. Namir Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 08-31-2010, 09:35 PM Tim and Gerson, I was looking at the HP-71B Math Pac manual and it does mention the Laguerre method as being used for that ROM (and I am sure for many other HP calculators). I have a Matlab function that implements of the Laguerre method and it works very well ... so the method is good, even though my implementation is bare bone -- not scaling or other programming tricks. Namir Edited: 31 Aug 2010, 10:33 p.m. Didier Lachieze Member Posts: 248 Threads: 5 Joined: Feb 2008 09-03-2010, 03:04 AM Hi Namir, You may be interested by the discussion/comparison of several Polynomial Rootfinder programs as well as the matlab sample program (Appendix A) in this document: Iterative Methods for Roots of Polynomials Regarding PROOT on HP calculators, the following post from Bill Wickes on comp.sys.hp48 back in 1992 shows that it was quite easy to port HP 71B assembler functions to the HP-48: HP 48 Polynomial Rootfinder Given this, I'm a bit surprised that the results for the second test on the HP 71B (+ Math Pac) and the HP 48GX (same CPU, same algorithm, same assembler code?) or the 49G+ in the following post on this forum are somewhat different (The Turtle (HP-71B) and the Hare (HP49G+) [LONG]): ```0.999999999944, 1.312E-12 for the HP 71B 0.999999994032, 2.066E-12 for the 48GX 0.999999994031, 1.441E-12 for the 49G+ ``` ... the HP 71B while being the oldest machine is the closest to the exact result (1,0)! Edited: 5 Sept 2010, 6:03 a.m. « Next Oldest | Next Newest »

 Possibly Related Threads... Thread Author Replies Views Last Post Thanks Tim for blue HP50g's keyboard colors Carey 4 732 11-07-2011, 02:33 PM Last Post: Jan Tim, Help! 15C LE emulator activation not working M. Joury 16 1,843 09-12-2011, 11:42 PM Last Post: M. Joury OT ->Tim Wessman: TI-Nspire current draw Joerg Woerner 1 513 05-08-2011, 12:11 PM Last Post: Tim Wessman Tim Wessman e-mail Hugh Evans 21 2,198 02-25-2006, 04:01 AM Last Post: Marcus von Cube, Germany

Forum Jump: