Add an example for math.fsum() and elaborate on the accurary note.

diff --git a/Doc/library/math.rst b/Doc/library/math.rst
index b33c597d7bdeda6e584f439f19614d0a016c3a79..4ecd14efd2f10c3320b5afe68bf65b907ffdb26b 100644 (file)
--- a/Doc/library/math.rst
+++ b/Doc/library/math.rst
@@ -87,14 +87,18 @@ Number-theoretic and representation functions
 .. function:: fsum(iterable)
 
    Return an accurate floating point sum of values in the iterable.  Avoids
-   loss of precision by tracking multiple intermediate partial sums.  The
-   algorithm's accuracy depends on IEEE-754 arithmetic guarantees and the
-   typical case where the rounding mode is half-even.
-
-   .. note::
-
-      The accuracy of fsum() may be impaired on builds that use
-      extended precision addition and then double-round the results.
+   loss of precision by tracking multiple intermediate partial sums::
+
+        >>> sum([.1, .1, .1, .1, .1, .1, .1, .1, .1, .1])
+        0.99999999999999989
+        >>> fsum([.1, .1, .1, .1, .1, .1, .1, .1, .1, .1])
+        1.0
+
+   The algorithm's accuracy depends on IEEE-754 arithmetic guarantees and the
+   typical case where the rounding mode is half-even.  On some non-Windows
+   builds, the underlying C library uses extended precision addition and may
+   occasionally double-round an intermediate sum causing it to be off in its
+   least significant bit.
 
    .. versionadded:: 2.6