#14274: Numerical approximation of a divergent integral
-------------------------------+--------------------------------------------
Reporter: eviatarbach | Owner: burcin
Type: defect | Status: new
Priority: major | Milestone: sage-5.9
Component: calculus | Resolution:
Keywords: | Work issues:
Report Upstream: N/A | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
-------------------------------+--------------------------------------------
Comment (by kcrisman):
Regarding the first one; basically, when we return a noun form
{{{
sage: integrate(x^3/sqrt(x^7+1), x, 1, oo)
integrate(x^3/sqrt(x^7 + 1), x, 1, +Infinity)
}}}
and call `n` we do
{{{
sage: A._convert(RR)
-2.0585298599985333
}}}
because we do
{{{
sage: sage.symbolic.integration.integral.DefiniteIntegral._evalf_??
from sage.gsl.integration import numerical_integral
# The gsl routine returns a tuple, which also contains the error.
# We only return the result.
return numerical_integral(f, a, b)[0]
}}}
So these are both manifestations of the same thing.
----
So... is it user error to numerically integrate a divergent integral? I
certainly don't know that we should be checking every numerical integral
for divergence, particularly since Maxima apparently can't (yet) do the
first one in any case!
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/14274#comment:2>
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?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.
