Changeset 7800:5184e5015a4c

Timestamp:

12/06/07 13:58:39 (5 years ago)

Author:

Mike Hansen <mhansen@…>

Branch:

default

Message:

Fixed #1139.

File:

• sage/calculus/calculus.py (modified) (2 diffs)

sage/calculus/calculus.py

@@ -1997 +1997 @@
             There is also a function \code{numerical_integral} that implements
             numerical integration using the GSL C library.  It is potentially
-            much faster and applies to arbitrary user defined functions. 
+            much faster and applies to arbitrary user defined functions.
+
+            Also, there are limits to the precision that Maxima can compute
+            the integral to due to limitations in quadpack.
+
+            sage: f = x
+            sage: f = f.nintegral(x,0,1,1e-14)
+            Traceback (most recent call last):
+            ...
+            ValueError: Maxima (via quadpack) cannot compute the integral to that precision
 
         EXAMPLES:
@@ -2049 +2058 @@
         syntax. 
         """
-        v = self._maxima_().quad_qags(var(x),
+        try:
+            v = self._maxima_().quad_qags(var(x),
                                       a, b, desired_relative_error,
                                       maximum_num_subintervals)
+        except TypeError, err:
+            if "ERROR NUMBER = 6" in str(err):
+                raise ValueError, "Maxima (via quadpack) cannot compute the integral to that precision"
+            else:
+                raise TypeError, err
+            
         return float(v[0]), float(v[1]), Integer(v[2]), Integer(v[3])