1234567891011121314151617181920212223242526 |
- static double
- rat_eval (const double a[], const size_t na,
- const double b[], const size_t nb, const double x)
- {
- size_t i, j;
- double u, v, r;
- u = a[na - 1];
- for (i = na - 1; i > 0; i--)
- {
- u = x * u + a[i - 1];
- }
- v = b[nb - 1];
- for (j = nb - 1; j > 0; j--)
- {
- v = x * v + b[j - 1];
- }
- r = u / v;
- return r;
- }
|