Sunday, March 16, 2008

Polynomials

A polynomial is a mathematical expression involving a sum of Powers in one or more variables multiplied by coefficient. A polynomial in one variable with constant coefficient is given by

 a_nx^n+...+a_2x^2+a_1x+a_0.
(1)

The individual summands with the coefficient (usually) included are called monomials (Becker and Weispfenning 1993, p. 191), whereas the products of the form x_1^(a_1)...x_n^(a_n) in the multivariate case, i.e., with the coefficients omitted, are called terms (Becker and Weispfenning 1993, p. 188). The highest power in a univariate polynomial is called its order, or sometimes its degree.

Any polynomial P(x) with P(0)!=0 can be expressed as

 P(x)=P(0)product_(rho)(1-x/rho),
(2)

where the product runs over the roots rho of P(rho)=0 and it is understood that multiple roots are counted with multiplicity.

A polynomial in two variables (i.e., a bivariate polynomial) with constant coefficient is given by

 a_(nm)x^ny^m+...+a_(22)x^2y^2+a_(21)x^2y+a_(12)xy^2+a_(11)xy+a_(10)x+a_(01)y+a_(00).
(3)

The sum of two polynomials is obtained by adding together the coefficient sharing the same powers of variables (i.e., the same terms) so, for example,

 (a_2x^2+a_1x+a_0)+(b_1x+b_0)=a_2x^2+(a_1+b_1)x+(a_0+b_0)
(4)

and has order less than (in the case of cancellation of leading terms) or equal to the maximum order of the original two polynomials. Similarly, the product of two polynomials is obtained by multiplying term by term and combining the results, for example

(a_2x^2+a_1x+a_0)(b_1x+b_0)=a_2x^2(b_1x+b_0)+a_1x(b_1x+b_0)+a_0(b_1x+b_0)
(5)
=a_2b_1x^3+(a_2b_0+a_1b_1)x^2+(a_1b_0+a_0b_1)x+a_0b_0,
(6)

and has order equal to the sum of the orders of the two original polynomials.

A polynomial quotient

 R(z)=(P(z))/(Q(z))
(7)

of two polynomials P(z) and Q(z) is known as a rational funtion. The process of performing such a division is called longdivision, with synthetic division being a simplified method of recording the division.

For any polynomial P(x), P(x)-x divides P(P(x))-x, meaning that the polynomial quotient is a rational polynomial or, in the case of an integer polynomial, another integer polynomial (pers. comm., N. Sato, Nov. 23, 2004).

Exchanging the coefficient of a univariate polynomial end-to-end produces a polynomial

 a_0x^n+a_1x^(n-1)+...+a_(n-1)x+a_n=0
(8)

whose roots are reciprocal1/x_i of the original rootsx_i.

horner's root provides a computationally efficient method of forming a polynomial from a list of its coefficients, and can be implemented in Mathematica as follows.

  Polynomial[l_List, x_] := Fold[x #1 + #2&, 0, l]

The following table gives special names given to polynomials of low orders.

polynomial orderpolynomial name
2quadratic polynomial
3cubic polynomial
4quartic
5quintic
6sextic

Polynomials of fourth degree may be computed using three multiplications and five additions if a few quantities are calculated first (Press et al. 1989):

 a_0+a_1x+a_2x^2+a_3x^3+a_4x^4=[(Ax+B)^2+Ax+C][(Ax+B)^2+D]+E,
(9)

where

A=(a_4)^(1/4)
(10)
B=(a_3-A^3)/(4A^3)
(11)
D=3B^2+8B^3+(a_1A-2a_2B)/(A^2)
(12)
C=(a_2)/(A^2)-2B-6B^2-D
(13)
E=a_0-B^4-B^2(C+D)-CD.
(14)

Similarly, a polynomial of fifth degree may be computed with four multiplications and five additions, and a polynomial of sixth degree may be computed with four multiplications and seven additions.

Polynomials of orders one to four are solvable using only rational operations and finite root extractions. A first-order equation is trivially solvable. A second-order equation is soluble using the quadratic equation. A third-order equation is solvable using the cubic equation. A fourth-order equation is solvable using the quartic equation. It was proved by Abel and Galois using group theory that general equations of fifth and higher order cannot be solved rationally with finite root extractions (Abel's impossibility theorem).

However, solutions of the general quintic equation may be given in terms of Jacobi theta functions or hypergeometric functions in one variable. Hermite and Kronecker proved that higher order polynomials are not soluble in the same manner. Klein showed that the work of Hermite was implicit in the group properties of the icosahedron. Klein's method of solving the quintic in terms of hypergeometric functions in one variable can be extended to the sextic, but for higher order polynomials, either hypergeometric functions in several variables or "Siegel functions" must be used (Belardinelli 1960, King 1996, Chow 1999). In the 1880s, Poincaré created functions which give the solution to the nth order polynomial equation in finite form. These functions turned out to be "natural" generalizations of the elliptic functions.

No comments:

test