Changeset 7511:38578c1fd58a

Timestamp:

12/02/07 11:12:42 (5 years ago)

Author:

William Stein <wstein@…>

Branch:

default

Message:

Followup patch for trac #1134 -- improve documentation of integration in calculus.py

File:

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

sage/calculus/calculus.py

@@ -1915 +1915 @@
         Return a floating point machine precision numerical
         approximation to the integral of self from a to b, computed
-        using floating point arithmetic and the GSL scientific
-        library.
+        using floating point arithmetic via maxima.
 
         INPUT:
@@ -1973 +1972 @@
         Computing to higher precision we see the truth:
             sage: f.n(200)
-            -0.00000000000074992740280181431112064614366496792309675391526978827185055
+            -0.00000000000074992740280181431112064614366622348652078895136533593355718
             sage: f.n(300)
-            -0.000000000000749927402801814311120646143662663009137292462589621789352092802261939388897590086687280282
+            -0.000000000000749927402801814311120646143662663009137292462589621789352095066181709095575681963967103004
 
         Now numerically integrating, we see why the answer is wrong:
@@ -1989 +1988 @@
             '2.565728500561051482917356396 E-127'
             sage: old_prec = gp.set_real_precision(50)
+            sage: gp.eval('intnum(x=17,42,exp(-x^2)*log(x))')
             '2.5657285005610514829173563961304785900147709554020 E-127'
-
+            sage: gp.set_real_precision(old_prec)
+            57
+            
         Note that the input function above is a string in PARI
         syntax.