🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

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

Re: Root finding - Bug in Num Recipes

#1308
From: williamroeder
Date: Tue, 24 Nov 1998 00:00
26 lines
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-mail
References: <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)