Article View: sci.math.num-analysis
Article #1308Re: Root finding - Bug in Num Recipes
From: williamroeder
Date: Tue, 24 Nov 1998 00:00
Date: Tue, 24 Nov 1998 00:00
26 lines
1761 bytes
1761 bytes
In article <3653374B.950EFD0F@trc.tellabs.com>, Ken Laberteaux <klaberte@trc.tellabs.com> wrote: > 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)); In Numerical Recipes in Fortran (first edition) the line was h = g - 2.0 * f / b So the 2 has been in there from the beginning. Why don't you post a bug report to their web site: http://www.nr.com/bugstop.html > > Can anyone shed some light into this inconsistency? > > Ken Laberteaux > -----------== Posted via Deja News, The Discussion Network ==---------- http://www.dejanews.com/ Search, Read, Discuss, or Start Your Own
Technical Details
Message-ID:
<73ed79$l4e$1@nnrp1.dejanews.com>Path:
arc1-ca.pugleaf.net!unknown!not-for-mailReferences:
<3653374B.950EFD0F@trc.tellabs.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: 1025b4,9646b0a747ad2fab X-Google-Attributes: gid1025b4,public X-Google-Thread: 1014db,9646b0a747ad2fab X-Google-Attributes: gid1014db,public From: williamroeder@my-dejanews.com Subject: Re: Root finding - Bug in Num Recipes Date: 1998/11/24 Message-ID: <73ed79$l4e$1@nnrp1.dejanews.com> X-Deja-AN: 415147235 References: <3653374B.950EFD0F@trc.tellabs.com> X-Http-Proxy: 1.0 x11.dejanews.com:80 (Squid/1.1.22) for client 159.62.105.224 Organization: Deja News - The Leader in Internet Discussion X-Article-Creation-Date: Tue Nov 24 13:44:10 1998 GMT Newsgroups: sci.math.num-analysis,comp.lang.c,gnu.misc.discuss X-Http-User-Agent: Mozilla/4.5 [en] (X11; I; IRIX 6.3 IP32)