1 files changed, 21 insertions, 0 deletions
diff --git a/gsl-1.9/doc/examples/specfun_e.c b/gsl-1.9/doc/examples/specfun_e.c
new file mode 100644
index 0000000..4f41f2a
--- /dev/null
+++ b/gsl-1.9/doc/examples/specfun_e.c
@@ -0,0 +1,21 @@
+#include <stdio.h>
+#include <gsl/gsl_errno.h>
+#include <gsl/gsl_sf_bessel.h>
+
+int
+main (void)
+{
+  double x = 5.0;
+  gsl_sf_result result;
+
+  double expected = -0.17759677131433830434739701;
+  
+  int status = gsl_sf_bessel_J0_e (x, &result);
+
+  printf ("status  = %s\n", gsl_strerror(status));
+  printf ("J0(5.0) = %.18f\n"
+          "      +/- % .18f\n", 
+          result.val, result.err);
+  printf ("exact   = %.18f\n", expected);
+  return status;
+}