htswanalysis/MACS/lib/gsl/gsl-1.11/cdf/rat_eval.h

   1 static double
   2 rat_eval (const double a[], const size_t na,
   3           const double b[], const size_t nb, const double x)
   4 {
   5   size_t i, j;
   6   double u, v, r;
   7
   8   u = a[na - 1];
   9
  10   for (i = na - 1; i > 0; i--)
  11     {
  12       u = x * u + a[i - 1];
  13     }
  14
  15   v = b[nb - 1];
  16
  17   for (j = nb - 1; j > 0; j--)
  18     {
  19       v = x * v + b[j - 1];
  20     }
  21
  22   r = u / v;
  23
  24   return r;
  25 }