#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:          
----------------------------+-----------------------------------------------
Changes (by kcrisman):

 * cc: benjaminfjones (added)


Old description:

> See extensive discussion at
>
> Reported by Ben Jones:
> {{{
> sage: integrate(ln(1+4/5*sin(x)), x, -3.1415, 3.1415)
>

> ;;;
> ;;; Binding stack overflow.
> ;;; Jumping to the outermost toplevel prompt
> ;;;
>
> ...
>
> ;;;
> ;;; Binding stack overflow.
> ;;; Jumping to the outermost toplevel prompt
> ;;;
>

> /Users/jonesbe/sage/sage-5.1.beta2/spkg/bin/sage: line 335: 86594 Illegal
> instruction     sage-ipython "$@" -i
> }}}
>
> There don't seem to be any problems with Maxima.  Although it's been
> traced down (thanks to Dan Drake) to 5.0.beta3, the only obvious
> candidate for a change in our interface is when we added `abs_integrate`,
> which doesn't seem to impact things.
> {{{
> (%i2) load(abs_integrate);
> (%o2) /Users/.../sage-5.0/local/share/maxima/5.26.0/s\
> hare/contrib/integration/abs_integrate.mac
> (%i4) keepfloat:true;
> (%o4)                                true
> (%i7) integrate(ln(1+4/5*sin(x)), x, -3.1415, 3.1415);
>                           3.1415
>                          /
>                          [            4 sin(x)
> (%o7)                    I         ln(-------- + 1) dx
>                          ]               5
>                          /
>                           - 3.1415
> }}}
> And we have
> {{{
> sage: maxima_calculus(" integrate(ln(1+4/5*sin(x)), x, -3.1415,
> 3.1415);")
> 'integrate(ln(4*sin(x)/5+1),x,-3.1415,3.1415)
> }}}
> so something must indeed be wrong in how we're sending this to Maxima.

New description:

 See extensive discussion at
 [https://groups.google.com/forum/?fromgroups#!topic/sage-
 support/dZFPiWRWSoc this sage-support thread].

 Reported by Ben Jones:
 {{{
 sage: integrate(ln(1+4/5*sin(x)), x, -3.1415, 3.1415)


 ;;;
 ;;; Binding stack overflow.
 ;;; Jumping to the outermost toplevel prompt
 ;;;

 ...

 ;;;
 ;;; Binding stack overflow.
 ;;; Jumping to the outermost toplevel prompt
 ;;;


 /Users/jonesbe/sage/sage-5.1.beta2/spkg/bin/sage: line 335: 86594 Illegal
 instruction     sage-ipython "$@" -i
 }}}

 There don't seem to be any problems with Maxima.  Although it's been
 traced down (thanks to Dan Drake) to 5.0.beta3, the only obvious candidate
 for a change in our interface is when we added `abs_integrate`, which
 doesn't seem to impact things.
 {{{
 (%i2) load(abs_integrate);
 (%o2) /Users/.../sage-5.0/local/share/maxima/5.26.0/s\
 hare/contrib/integration/abs_integrate.mac
 (%i4) keepfloat:true;
 (%o4)                                true
 (%i7) integrate(ln(1+4/5*sin(x)), x, -3.1415, 3.1415);
                           3.1415
                          /
                          [            4 sin(x)
 (%o7)                    I         ln(-------- + 1) dx
                          ]               5
                          /
                           - 3.1415
 }}}
 And we have
 {{{
 sage: maxima_calculus(" integrate(ln(1+4/5*sin(x)), x, -3.1415, 3.1415);")
 'integrate(ln(4*sin(x)/5+1),x,-3.1415,3.1415)
 }}}
 so something must indeed be wrong in how we're sending this to Maxima.

--

-- 
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/13097#comment:1>
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.

Reply via email to