🚀 go-pugleaf

Ken Laberteaux <klaberte@trc.tellabs.com> writes:
> 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?

It's my (paranoid, perhaps) hypothesis that Numerical Recipes includes
harmless obfuscations in their code to protect their copyright, much
like Rand McNally introduces harmless small errors in their maps in
order to identify illegal copies.

I don't know about this particular routine, but based on prior
experience with NR, my bet would be that if you follow the logic of
the code, the routine is correct -- but in a way that you may not have
ever written it had you started from scratch.

-- 
- Sean Eddy
- Dept. of Genetics, Washington University in St. Louis
- http://www.genetics.wustl.edu/eddy/