Ruffini's rule

In mathematics, Ruffini's rule is a practical way for paper-and-pencil computation of the Euclidean division of a polynomial by a binomial of the form xr. It was described by Paolo Ruffini in 1804.[1] The rule is a special case of synthetic division in which the divisor is a linear factor.

Algorithm

The rule establishes a method for dividing the polynomial

by the binomial

to obtain the quotient polynomial

;

The algorithm is in fact the long division of P(x) by Q(x).

To divide P(x) by Q(x):

  1. Take the coefficients of P(x) and write them down in order. Then, write r at the bottom-left edge just over the line:
  2. Pass the leftmost coefficient (an) to the bottom just under the line.
  3. Multiply the rightmost number under the line by r, and write it over the line and one position to the right.
  4. Add the two values just placed in the same column.
  5. Repeat steps 3 and 4 until no numbers remain.

The b values are the coefficients of the result (R(x)) polynomial, the degree of which is one less than that of P(x). The final value obtained, s, is the remainder. The polynomial remainder theorem asserts that the remainder is equal to P(r), and the value of the polynomial at r.

Example

Here is a worked example of polynomial division as described above.

Let:

P(x) will be divided by Q(x) using Ruffini's rule. The main problem is that Q(x) is not a binomial of the form x r, but rather x + r. Q(x) must be rewritten as

Now the algorithm is applied:

1. Write down the coefficients and r. Note that, as P(x) didn't contain a coefficient for x, 0 is written:

    |     2     3     0     -4
    |                                    
 -1 |                                    
----|----------------------------
    |                                    
    |

2. Pass the first coefficient down:

    |     2     3     0     -4
    |                                    
 -1 |                                    
----|----------------------------
    |     2                              
    |

3. Multiply the last obtained value by r:

    |     2     3     0     -4
    |                                    
 -1 |          -2                         
----|----------------------------
    |     2                              
    |

4. Add the values:

    |     2     3     0     -4
    |
 -1 |          -2
----|----------------------------
    |     2     1
    |

5. Repeat steps 3 and 4 until it's finished:

    |     2     3     0     -4
    |
 -1 |          -2    -1      1
----|----------------------------
    |     2     1    -1     -3
    |{result coefficients}{remainder}


So, if original number = divisor × quotient + remainder, then

, where
and

Uses

Ruffini's rule has many practical applications, most of which rely on simple division (as demonstrated below) or the common extensions given still further below.

Polynomial root-finding

The rational root theorem states that for a polynomial f(x) = anxn + an1xn1 + ... + a1x + a0 all of whose coefficients (an through a0) are integers, the real rational roots are always of the form p/q, where p is an integer divisor of a0 and q is an integer divisor of an. Thus if our polynomial is

the possible rational roots are all the integer divisors of a0 (2):

(The example is simple because the polynomial is monic (an = 1). Fornon-monic polynomials the set of possible roots will include some fractions but only a finite number of them since an and a0 have only a finite number of integer divisors each.) In any case, for monic polynomials, every rational root is an integer and so every integer root is just a divisor of the constant term (a0). It can be shown that to remain true for non-monic polynomials: to find the integer roots of any polynomials with integer coefficients, checking the divisors of the constant term suffices.

Therefore, setting r equal to each of the possible roots in turn, the polynomial is divided by (x  r). If the resulting quotient has no remainder, a root was found.

Any of the following three methods can be chosen since all of them yield the same results, with the exception that only through the second method and the third method (when applying Ruffini's rule to obtain a factorization) can discover if a given root is repeated. (Neither method will discover irrational or complex roots.)

Method 1

Division P(x) by the binomial (x each possible root). If the remainder is 0, the selected number is a root (and vice versa):

    |    +1    +2    -1     -2                      |    +1    +2    -1    -2
    |                                               |
 +1 |          +1    +3     +2                   -1 |          -1    -1    +2
----|----------------------------               ----|---------------------------
    |    +1    +3    +2      0                      |    +1    +1    -2     0

    |    +1    +2    -1     -2                      |    +1    +2    -1    -2
    |                                               |
 +2 |          +2    +8    +14                   -2 |          -2     0    +2
----|----------------------------               ----|---------------------------
    |    +1    +4    +7    +12                      |    +1     0    -1     0

In the example, P(x) is a degree three polynomial. By the fundamental theorem of algebra, it can have no more than three complex solutions. Hence, the polynomial is factored as follows:

Method 2

Start just as in Method 1 until a valid root is found. Then, instead of restarting the process with the other possible roots, continue testing the possible roots against the result of the Ruffini on the valid root found currently until only a coefficient remains (remember that roots can be repeated: if you get stuck, try each valid root twice):

    |    +1    +2    -1    -2                      |    +1    +2    -1    -2
    |                                              |
 -1 |          -1    -1    +2                   -1 |          -1    -1    +2
----|---------------------------               ----|---------------------------
    |    +1    +1    -2   | 0                      |    +1    +1    -2   | 0
    |                                              |
 +2 |          +2    +6                         +1 |          +1    +2
-------------------------                      -------------------------
    |    +1    +3   |+4                            |    +1    +2   | 0
                                                   |
                                                -2 |          -2
                                               -------------------
                                                   |    +1   | 0

Method 3

  • Determine the set of the possible integer or rational roots of the polynomial according to the rational root theorem.
  • For each possible root r, instead of performing the division P(x)/(xr), apply the polynomial remainder theorem, which states that the remainder of the division is P(r), the polynomial evaluated for x = r.

Thus, for each r in our set, r is actually a root of the polynomial if and only if P(r)=0

That shows that finding integer and rational roots of a polynomial no division or the application of Ruffini's rule.

However, once a valid root has been found, call it r1: Ruffini's rule can be applied to determine

Q(x)=P(x)/(xr1).

That allows a partial factorization of the polynomial as

P(x)=(xr1Q(x)

Any additional (rational) root of the polynomial is also a root of Q(x) and, of course, is still to be found among the possible roots determined earlier, which have not yet been checked (any value already determined not to be a root of P(x) is not a root of Q(x) either; more formally, P(r)≠0 → Q(r)≠0 ).

Thus, you can proceed evaluating Q(r) instead of P(r), and (as long as you can find another root, r2) dividing Q(r) by (xr2).

Even if you search only for roots, that allows you to evaluate polynomials of successively-smaller degree as the factorization proceeds.

If, as is often the case, you are also factorizing a polynomial of degree n:

  • if you have found p=n rational solutions you end up with a complete factorization (see below) into p=n linear factors;
  • if you have found p<n rational solutions you end up with a partial factorization (see below) into p linear factors and another non-linear factor of degree np, which, in turn, may have irrational or complex roots.
Finding roots without applying Ruffini's Rule
P(x)=x³+2x²–x–2

Possible roots = {1, –1, 2, -2}

  • P(1) = 0 → x1 = 1
  • P(-1) = 0 → x2 = -1
  • P(2) = 12 → 2 is not a root of the polynomial

and the remainder of (x³+2x²-x-2)/(x-2) is 12

  • P(-2) = 0 → x3 = -2
Finding roots applying Ruffini's Rule and obtaining a (complete) factorization
P(x) = x³+2x²-x-2

Possible roots = {1, -1, 2, -2}

  • P(1) = 0 → x1 = 1

Then, applying Ruffini's Rule:

(x³+2x²-x-2)/(x-1) = (x²+3x+2)
x³+2x²-x-2 = (x-1)(x²+3x+2)

Here, r1=-1 and Q(x) = x²+3x+2

  • Q(-1) = 0 → x2 = -1

Again, applying Ruffini's Rule:

(x²+3x+2)/(x+1) = (x+2)
x³+2x²-x-2 = (x-1)(x²+3x+2) = (x-1)(x+1)(x+2)

As it was possible to completely factorize the polynomial, it's clear that the last root is -2 (the previous procedure would have given the same result, with a final quotient of 1).

Polynomial factoring

Having used the "p/q" result above (or, to be fair, any other means) to find all the real rational roots of a particular polynomial, it is but a trivial step further to partially factor that polynomial using those roots. As is well-known, each linear factor (x  r), which divides a given polynomial corresponds with a root r and vice versa.

So if

 is our polynomial; and
are the found roots, then consider the product

By the fundamental theorem of algebra, R(x) should be equal to P(x), if all the roots of P(x) are rational. However, since the method finds only rational roots, it is very likely that R(x) is not equal to P(x); it is very likely that P(x) has some irrational or complex roots not in R. So consider

, which can be calculated using polynomial long division.

If S(x) = 1, then it is known R(x) = P(x) and procedure is done. Otherwise, S(x) will itself be a polynomial, which is another factor of P(x), which has no real rational roots. Therefore, write out the right-hand-side of the following equation in full:

One can call it a complete factorization of P(x) over Q (the rationals) if S(x) = 1. Otherwise, there is only a partial factorization of P(x) over Q, which may or may not be further factorable over the rationals but will certainly be further factorable over the reals or at worst the complex plane. (Note that a "complete factorization" of P(x) over Q means a factorization as a product of polynomials with rational coefficients, such that each factor is irreducible over Q, ans "irreducible over Q" means that the factor cannot be written as the product of two non-constant polynomials with rational coefficients and smaller degree.)

Example 1: no remainder

Let

Using the methods described above, the rational roots of P(x) are:

Then, the product of (x each root) is

And P(x)/R(x):

Hence the factored polynomial is P(x) = R(x) · 1 = R(x):

Example 2: with remainder

Let

Using the methods described above, the rational roots of P(x) are:

Then, the product of (x each root) is

And P(x)/R(x)

As , the factored polynomial is P(x) = R(x) · S(x):

Factoring over the complexes

To completely factor a given polynomial over C, the complex numbers, all of its roots must be known (and that could include irrational and/or complex numbers). For example, consider the polynomial above:

Extracting its rational roots and factoring it, yields:

But that is not completely factored over C. If factorization of a polynomial must finish to a product of linear factors, quadratic factor must be dealt with:

The easiest way is to use quadratic formula, which yields

and the solutions

So the completely factored polynomial over C will be:

However, one cannot, in every case expect things to be so easy; the quadratic formula's analogue for fourth-order polynomials is very convoluted and no such analogue exists for fifth-order or higher polynomials. See Galois theory for a theoretical explanation of why that is so, and see numerical analysis for ways to approximate roots of polynomials numerically.

Limitations

It is entirely possible that in looking for a given polynomial's roots, an intricate higher-order polynomial for S(x) is obtained, which is further factorable over the rationals even before irrational or complex factorings are considered such as for the polynomial x5  3x4 + 3x3  9x2 + 2x  6. Using Ruffini's method, only one root is found (x = 3), factoring it out to P(x) = (x4 + 3x2 + 2)(x  3).

As explained above, if stated assignment was to "factor into irreducibles over C", it would be necessary to find some way to dissect the quartic and look for its irrational and/or complex roots. But if the assignment were "factor into irreducibles over Q", one might think it is already done, but it is important to realize that might not be the case.

In that instance, the quartic is actually factorable as the product of two quadratics (x2 + 1)(x2 + 2). They, at last, are irreducible over the rationals (and the reals as well in that example), which concludes the method; P(x) = (x2 + 1)(x2 + 2)(x  3). In that instance, it is in fact easy to factor the quartic by treating it as a biquadratic equation; but finding such factorings of a higher degree polynomial can be very difficult.

History

The method was invented by Paolo Ruffini, eho took part in a competition organised by the Italian Scientific Society (of Forty). The question to be answered was a method to find the roots of any polynomial. Five submissions were received. In 1804, Ruffini's was awarded the first place, and his method was published. Ruffini published refinements of his method in 1807 and 1813.

Horner's method was published in 1819 and in a refined version in 1845.

See also

References

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.