[fpc-pascal]Double-precision recursion

[Rich Pasco]

My Pascal implementation of recursive floating-point algorithm fails,
but a line-by-line C translation of the same program works fine.

I stripped the program down to "just the bug" and posted the minimized
code at http://home.earthlink.net/~richpasco/pasbug.html

I would deeply appreciate anyone telling me just why the Pascal version
fails.  Trivial suggestions like "avoid floating point" or "restructure
it to avoid recursion" or "use the C version instead" are not helpful.

     - Rich