#13097: integral of ln(1+4/5*sin(x)) from -3.1415 to 3.1415 crashes Maxima and
Sage
----------------------------+-----------------------------------------------
Reporter: kcrisman | Owner: burcin
Type: defect | Status: new
Priority: major | Milestone: sage-5.1
Component: calculus | Resolution:
Keywords: | Work issues:
Report Upstream: N/A | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
----------------------------+-----------------------------------------------
Comment (by nbruin):
I don't think maxima understands `ln`. If you do
{{{
maxima_calculus(" integrate(log(1+4/5*sin(x)), x, -3.1415, 3.1415);")
}}}
you see the same thing happening. No problem in vanilla maxima. But when
we load `abs_integrate`:
{{{
(%i1) load(abs_integrate);
(%o1)
/usr/local/sage/5.0/local/share/maxima/5.26.0/share/contrib/integration/abs_integrate.mac
(%i2) integrate(log(1+4/5*sin(x)), x, -3.1415, 3.1415);
...
Condition of type: STACK-OVERFLOW
BINDING-STACK overflow at size 8448. Stack can probably be resized.
}}}
FYI: with
{{{
from sage.interfaces.maxima_lib import *
from sage.libs.ecl import *
args=(ln(1+4/5*sin(x)), x, -3.1415, 3.1415)
I=EclObject(([max_integrate],[sr_to_max(SR(a)) for a in args]))
I
}}}
you see what expression gets passed to maxima_lib. With
{{{
EclObject('#$integrate(log(1+4/5*sin(x)), x, -3.1415, 3.1415)$')
}}}
you see what maxima's parser would produce. With
{{{
maxima_eval(I)
}}}
you can trigger the same explosion. Note that with this trick, you can
quickly see that
{{{
EclObject('#$integrate(ln(1+4/5*sin(x)), x, -3.1415, 3.1415)$')
}}}
creates an essentially different expression.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/13097#comment:4>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
