Ticket #8321: trac_8321_numeric_int_mpmath.patch

sage/symbolic/integration/integral.py

@@ -176 +176 @@
                 pass
         return None
 
-    def _evalf_(self, f, x, a, b, parent=None):
+    def _evalf_(self, f, x, a, b, parent = None): 
         """
         Returns numerical approximation of the integral
         
@@ -186 +186 @@
             sage: h = definite_integral(sin(x)/x^2, x, 1, 2); h
             integrate(sin(x)/x^2, x, 1, 2)
             sage: h.n() # indirect doctest
-            0.4723991772689525...
+            0.472399177268953
 
         TESTS:
 
@@ -194 +194 @@
 
             sage: integrate(x^2.7 * e^(-2.4*x), x, 0, 3).n()
             0.154572952320790
+
+        #Check if #8321 is fixed:
+
+            sage: d = definite_integral(sin(x)/x^2, x, 1, 2)
+            sage: d.n(77)
+            0.4723991772689525199904
         """
-        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]
+
+        import sage.libs.mpmath.all as mpmath
+
+        try:
+            precision = parent.prec()
+            mpmath.mp.prec = precision
+        except AttributeError:
+            precision = mpmath.mp.prec
+
+        mp_f = lambda z: \
+               f(x = mpmath.mpmath_to_sage(z,precision))
+
+        return mpmath.call(mpmath.quad,mp_f,[a,b])
 
     def _tderivative_(self, f, x, a, b, diff_param=None):
         """