| 37 | */ |
| 38 | template <typename T, size_t N> |
| 39 | T polyHorner(const T (&a)[N], T x) { |
| 40 | if (N == 0) |
| 41 | return 0; |
| 42 | |
| 43 | T y = a[N - 1]; |
| 44 | for (size_t n = 1; n < N; n++) { |
| 45 | y = a[N - 1 - n] + y * x; |
| 46 | } |
| 47 | return y; |
| 48 | } |
| 49 | |
| 50 | /** Evaluates a polynomial with coefficients `a[n]` at `x`. |
| 51 | Uses Estrin's method. |