On Thu, Dec 29, 2011 at 2:16 PM, Edwin Woollett <woollett at charter.net>wrote:
> attempted evaluation of bfloat of an expression
> involving gamma_incomplete's leads to an error message:
> ------------------------------**----------------------------
> Maxima 5.25.1 http://maxima.sourceforge.net
> using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL)
>
>
> (%i1) r:integrate(integrate(cos(x*y***log(2)),y,1,100),x,1,100);
>
> (%o1) -(%i*gamma_incomplete(0,10000***%i*log(2))
> -2*%i*gamma_incomplete(0,100*%**i*log(2))
> +%i*gamma_incomplete(0,%i*log(**2))-%i*gamma_incomplete(0,-%i***
> log(2))
> +2*%i*gamma_incomplete(0,-100***%i*log(2))
> -%i*gamma_incomplete(0,-10000***%i*log(2))) /(2*log(2))
>
> (%i2) block([fpprec:20],bfloat(r));
>
> expintegral_e: series failed.
> -- an error. To debug this try: debugmode(true);
>
>
Could you file a bug report on this? I know about some of these issues and
am working on them (slowly). The basic parts work, but the corner cases
need more work (especically since my approach fails in some cases where the
original gave correct answers).
Thanks,
Ray