#17468: Strange integration error/hang with log(abs(sin(x)))
-------------------------------------------------+-------------------------
Reporter: kcrisman | Owner:
Type: defect | Status: new
Priority: major | Milestone: sage-6.5
Component: calculus | Resolution:
Keywords: | Merged in:
Authors: | Reviewers:
Report Upstream: Fixed upstream, but not in a | Work issues:
stable release. | Commit:
Branch: | Stopgaps:
Dependencies: |
-------------------------------------------------+-------------------------
Comment (by pbruin):
Also directly in Maxima (with a very long time between "encountered
log(0)" and "encountered a Lisp error"):
{{{
Maxima 5.35.1 http://maxima.sourceforge.net
using Lisp ECL 13.5.1
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) integrate(log(abs(2*sin(u))),u,0,%pi/3);
sqrt(3) %i + 1 sqrt(3) %i + 1
(%o1) (18 %i li (--------------) + 18 %i li (- --------------) - 3 %pi
log(3)
2 2 2 2
2
2 %pi log(3) %i
%pi
+ 2 %i %pi )/18 + ---------- -
-------
6
12
(%i2) load(abs_integrate);
(%o2)
/home/bruinpj/src/sage/local/share/maxima/5.35.1/share/contrib/integration/abs_integrate.mac
(%i3) integrate(log(abs(2*sin(u))),u,0,%pi/3);
log: encountered log(0).
Maxima encountered a Lisp error:
In form
(PROGN (REQUIRE :DEFSYSTEM) (RUN))
CALL-ARGUMENTS-LIMIT exceeded
Automatically continuing.
To enable the Lisp debugger set *debugger-hook* to nil.
}}}
--
Ticket URL: <http://trac.sagemath.org/ticket/17468#comment:7>
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 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-trac.
For more options, visit https://groups.google.com/d/optout.