#11233: failing calculation of an integral
-------------------------+--------------------------------------------------
Reporter: casamayou | Owner: burcin
Type: defect | Status: new
Priority: major | Milestone: sage-4.7
Component: calculus | Keywords: integrate
Author: | Upstream: Not yet reported upstream; Will do
shortly.
Reviewer: | Merged:
Work_issues: |
-------------------------+--------------------------------------------------
Changes (by kcrisman):
* upstream: N/A => Not yet reported upstream; Will do shortly.
Comment:
I can confirm this in Maxima 5.24.0:
{{{
Maxima 5.24.0 http://maxima.sourceforge.net
using Lisp SBCL 1.0.24
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) assume(a>0);
(%o1) [a > 0]
(%i2) assume(x>0);
(%o2) [x > 0]
(%i3) integrate(log(1+a/(x*t)^2),t,1,inf);
2
x + a sqrt(a)
sqrt(a) x log(------) - 2 a atan(-------)
2 x
x
(%o3) -----------------------------------------
sqrt(a) x
}}}
Looks like a different branch of log is being used. Also for comparison,
here is what GSL says:
{{{
sage: numerical_integral(log(1+7/(x^2)),1,oo)
(4.3202562566855764, 2.2384237933839469e-07)
}}}
We should file something with the Maxima tracker.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11233#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.
