Article View: sci.math.num-analysis
Article #1312Re: Root finding - Bug in Num Recipes
From: Ken Laberteaux
Date: Thu, 19 Nov 1998 00:00
Date: Thu, 19 Nov 1998 00:00
35 lines
1801 bytes
1801 bytes
Where did you find your version? Perhaps it is from some freely used source. Please let me know. Ken p.s. thanks for checking. Eric Backus wrote: > > Ken Laberteaux <klaberte@trc.tellabs.com> writes: > > > I am reading the 2nd Edition of _Numerical Recipes in C_. I am > > specifically looking at laguer.c. > > > > On p. 373 the line of code which implements Eq 9.5.7 (p. 372) is > > > > h=Csub(g2,RCmul(2.0,Cdiv(f,b))); > > > > However Eq 9.5.7 (p. 372) makes no reference of the 2, indicating that > > the correct version should be > > > > h=Csub(g2,Cdiv(f,b)); > > > > Can anyone shed some light into this inconsistency? > > I agree that this appears to be inconsistant, and that your proposed > modification looks correct. > > My own code for the Laguerre method does not have the factor of two, > and seems to work well. > > -- > Eric Backus <eric_backus@hp.com> > http://labejb.lsid.hp.com/ > (425) 335-2495
Technical Details
Message-ID:
<3654C552.D6BB5C80@trc.tellabs.com>Path:
arc1-ca.pugleaf.net!unknown!not-for-mailReferences:
<3653374B.950EFD0F@trc.tellabs.com> <tuogq3mpgr.fsf@labejb.lsid.hp.com>Show Full Headers
Complete Article Headers
X-Google-Language: ENGLISH,ASCII-7-bit X-Google-Thread: 107079,9646b0a747ad2fab X-Google-Attributes: gid107079,public X-Google-Thread: 1014db,9646b0a747ad2fab X-Google-Attributes: gid1014db,public X-Google-Thread: 1025b4,9646b0a747ad2fab X-Google-Attributes: gid1025b4,public From: Ken Laberteaux <klaberte@trc.tellabs.com> Subject: Re: Root finding - Bug in Num Recipes Date: 1998/11/19 Message-ID: <3654C552.D6BB5C80@trc.tellabs.com> X-Deja-AN: 413707706 Content-Transfer-Encoding: 7bit References: <3653374B.950EFD0F@trc.tellabs.com> <tuogq3mpgr.fsf@labejb.lsid.hp.com> X-Accept-Language: en Content-Type: text/plain; charset=us-ascii Mime-Version: 1.0 Newsgroups: sci.math.num-analysis,comp.lang.c,gnu.misc.discuss