🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

Article View: sci.math.num-analysis
Article #1312

Re: Root finding - Bug in Num Recipes

#1312
From: Ken Laberteaux
Date: Thu, 19 Nov 1998 00:00
35 lines
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-mail
References: <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