#3863: numerical integration of x^2.7 * e^(-2.4*x) fails
----------------------+-----------------------------------------------------
Reporter: ddrake | Owner: gfurnish
Type: defect | Status: new
Priority: major | Milestone: sage-4.0.1
Component: calculus | Keywords: integration integral calculus symbolic
----------------------+-----------------------------------------------------
Comment(by kcrisman):
Oh, totally, and of course we have numerical_integrate and nintegral (one
is Maxima, one is GSL I think). So I guess what I'm saying is that when
someone has time (which won't be me until June due to the (US) holiday and
then a conference) it should be relatively straightforward to at least do
something in the generic code for .n() on whatever kind of object
(symbolic?) an expression like integrate(f(x),x,a,b) is like
{{{
try:
usual code
except TypeError:
if "integrate" is in the string:
either reraise the TypeError with a message suggesting the use of
numerical_integrate
or actually try to string-magic replace "integrate" with
"numerical_integrate", which should be possible (you'd have to take [0] of
that result, of course)
else:
do whatever used to happen with a TypeError
}}}
I think that it is reasonable for a user to expect that .n() numerically
evaluates an expression as best as possible, even one Maxima can't
evaluate! And all these methods are much more sophisticated than Simpson,
certainly, so that's not the issue.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/3863#comment:4>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
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
-~----------~----~----~----~------~----~------~--~---
