1 files changed, 30 insertions, 0 deletions
diff --git a/gsl-1.9/doc/examples/ieeeround.c b/gsl-1.9/doc/examples/ieeeround.c
new file mode 100644
index 0000000..04f837d
--- /dev/null
+++ b/gsl-1.9/doc/examples/ieeeround.c
@@ -0,0 +1,30 @@
+#include <stdio.h>
+#include <gsl/gsl_math.h>
+#include <gsl/gsl_ieee_utils.h>
+
+int
+main (void)
+{
+  double x = 1, oldsum = 0, sum = 0; 
+  int i = 0;
+
+  gsl_ieee_env_setup (); /* read GSL_IEEE_MODE */
+
+  do 
+    {
+      i++;
+      
+      oldsum = sum;
+      sum += x;
+      x = x / i;
+      
+      printf ("i=%2d sum=%.18f error=%g\n",
+              i, sum, sum - M_E);
+
+      if (i > 30)
+         break;
+    }  
+  while (sum != oldsum);
+
+  return 0;
+}