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

   1 #include <stdio.h>
   2 #include <gsl/gsl_math.h>
   3 #include <gsl/gsl_deriv.h>
   4
   5 double f (double x, void * params)
   6 {
   7   return pow (x, 1.5);
   8 }
   9
  10 int
  11 main (void)
  12 {
  13   gsl_function F;
  14   double result, abserr;
  15
  16   F.function = &f;
  17   F.params = 0;
  18
  19   printf ("f(x) = x^(3/2)\n");
  20
  21   gsl_deriv_central (&F, 2.0, 1e-8, &result, &abserr);
  22   printf ("x = 2.0\n");
  23   printf ("f'(x) = %.10f +/- %.10f\n", result, abserr);
  24   printf ("exact = %.10f\n\n", 1.5 * sqrt(2.0));
  25
  26   gsl_deriv_forward (&F, 0.0, 1e-8, &result, &abserr);
  27   printf ("x = 0.0\n");
  28   printf ("f'(x) = %.10f +/- %.10f\n", result, abserr);
  29   printf ("exact = %.10f\n", 0.0);
  30
  31   return 0;
  32 }