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

   1 double
   2 quadratic (double x, void *params)
   3 {
   4   struct quadratic_params *p
   5     = (struct quadratic_params *) params;
   6
   7   double a = p->a;
   8   double b = p->b;
   9   double c = p->c;
  10
  11   return (a * x + b) * x + c;
  12 }
  13
  14 double
  15 quadratic_deriv (double x, void *params)
  16 {
  17   struct quadratic_params *p
  18     = (struct quadratic_params *) params;
  19
  20   double a = p->a;
  21   double b = p->b;
  22   double c = p->c;
  23
  24   return 2.0 * a * x + b;
  25 }
  26
  27 void
  28 quadratic_fdf (double x, void *params,
  29                double *y, double *dy)
  30 {
  31   struct quadratic_params *p
  32     = (struct quadratic_params *) params;
  33
  34   double a = p->a;
  35   double b = p->b;
  36   double c = p->c;
  37
  38   *y = (a * x + b) * x + c;
  39   *dy = 2.0 * a * x + b;
  40 }