Root finding - Bug in Num Recipes - Article 1322 - go-pugleaf

Article View: sci.math.num-analysis

Article #1322

Back to sci.math.num-analysis Previous Next

Root finding - Bug in Num Recipes

#1322

From: Ken Laberteaux
Date: Wed, 18 Nov 1998 00:00

15 lines
1106 bytes

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?

Ken Laberteaux

Technical Details

Message-ID: <3653374B.950EFD0F@trc.tellabs.com>

Path: arc1-ca.pugleaf.net!unknown!not-for-mail

Show Full Headers

Complete Article Headers

X-Google-Language: ENGLISH,ASCII-7-bit
X-Google-Thread: 1014db,9646b0a747ad2fab,start
X-Google-Attributes: gid1014db,public
X-Google-Thread: 1025b4,9646b0a747ad2fab,start
X-Google-Attributes: gid1025b4,public
X-Google-Thread: 107079,9646b0a747ad2fab,start
X-Google-Attributes: gid107079,public
From: Ken Laberteaux <klaberte@trc.tellabs.com>
Subject: Root finding - Bug in Num Recipes
Date: 1998/11/18
Message-ID: <3653374B.950EFD0F@trc.tellabs.com>
X-Deja-AN: 413270069
Content-Transfer-Encoding: 7bit
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