On 2014-10-15, William Stein <[email protected]> wrote:

> If I try the following integrals, the first three work, but the last
> returns an error (RuntimeError: ECL says: CEXPT only defined for
> non-negative integral exponents.).
>
> integral(1/x^(1/2),x,1,10)
> integral(1/x^.5,x,1,10)
> integral(1/x^(3/2),x,1,10)
> integral(1/x^1.5,x,1,10)

The immediate cause of the error is that Sage enables the keepfloat
flag, otherwise Maxima converts floats to rationals when trying to
compute integrals. Some (maybe a lot) of the code for working with
polynomials and rational functions expects only integers and rationals,
so the presence of floats causes trouble.

I suppose that we have to decide between fixing up Maxima to handle
floats everywhere that keepfloat might introduce them, or do away
with keepfloat. I'm in favor of the latter.

Probably would be helpful if someone would file a bug report.
http://sourceforge.net/p/maxima/bugs

best,

Robert Dodier

-- 
You received this message because you are subscribed to the Google Groups 
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/d/optout.

Reply via email to