[Maxima] Exponential Integrals - Complex Bigfloat algorithm

Barton Willis willisb at unk.edu
Sun Jul 20 11:45:36 CDT 2008


-----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


More information about the Maxima mailing list