-----maxima-bounces at math.utexas.edu wrote: -----
>I suppose the problems arise because we need the constants BIGFLOAT%E and
>BIGFLOAT%GAMMA in the calculation.
That's a good start, I think. Try something like
(let ((bigfloat%e ($bfloat '$%e)) ...)
Also, if your code changes the value of fpprec, look at the macro
bind-fpprec (defined in nummod.lisp).
Barton