🚀 go-pugleaf

>Many of the examples from Numerical recipies don't work either.
>The explanation I got for it is that a lot of numerical work in C
>makes undue assumptions about the treatment of floats and doubles
>in functions. You could try #define double float to get a rough idea
>of what's going on. I gave up and hand coded.

I've done a number of serious numerical projects on the Amiga
(some orbital mechanics stuff for modeling planetary motion, an
ODE solver, finite element stuff for PDE's, etc) using the FFP,
IEEE double precision, and the '882 IEEE libraries.  My results
were rock-solid when compared against the same code on other 
machines (VAXen).  

I'm running Lattice 5.04.

It's true that the C version of Numerical Recipes has some funny
code in it; best to develop an understanding of the specific algorithm,
then make sure their code is doing what it says it is...

Regards,
 Art Baker			| "Whatever happened to Fay Wrey --
 (usual disclaimers)		|   that delicate, satin-draped frame..."
=========================================================================