htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/specfun_e.c

   1 #include <stdio.h>
   2 #include <gsl/gsl_errno.h>
   3 #include <gsl/gsl_sf_bessel.h>
   4
   5 int
   6 main (void)
   7 {
   8   double x = 5.0;
   9   gsl_sf_result result;
  10
  11   double expected = -0.17759677131433830434739701;
  12
  13   int status = gsl_sf_bessel_J0_e (x, &result);
  14
  15   printf ("status  = %s\n", gsl_strerror(status));
  16   printf ("J0(5.0) = %.18f\n"
  17           "      +/- % .18f\n",
  18           result.val, result.err);
  19   printf ("exact   = %.18f\n", expected);
  20   return status;
  21 }