Horner scheme

The Horner scheme is an algorithm for the efficient evaluation of polynomial functions, and for dividing polynomials by linear polynomials.

Given a number x and a polynomial p(T) = a₀ + a₁T + ... + a_nT ⁿ, the Horner scheme computes the number

p(x) = a₀ + a₁x + a₂x² + ... + a_n xⁿ

as well as a polynomial q(T) = b₀ + b₁T + ... + b_n-1T ^n-1 such that

p(T) = (T - x) · q(T) + p(x).

The algorithm works as follows:

1. set i := n - 1
2. set b_i := a_n
3. if i < 0, stop; the result p(x) is in b_-1.
4. set i := i - 1
5. set b_i := b_i+1 * x + a_i+1
6. Go to step 3.

This is the method of choice for evaluating polynomials; it is faster and more numerically stable than the "normal" method, which involves computing the powers of x and multiplying them with the coefficients. The Horner scheme is often used to convert between different positional number systems (in which case x is the base of the number system, and the a_i are the digits) and can also be used if x is a matrix, in which case the gain is even larger.