Finding the roots of polynomials has occupied mathematicians for many centuries. For equations up to fourth order, there are algebraic expressions for the roots. For higher order equations, many excellent numerical methods are available, but the results are not always reliable.

James Wilkinson (1963) examined the behaviour of a high-order polynomial

whose roots for are at the first twenty positive integers. He found that a tiny perturbation of the coefficients may have a major impact on the nature and location of the roots.

For simplicity, we will consider a tenth-order polynomial

Expanded into powers of , this becomes

The leading coefficient is . The coefficient of is .

In the figure above, we plot the polynomial and also the polynomial

with . That is, the coefficient of is changed by one part in a million. Surely, this could not have much influence on the roots? But it does!

As increases, the roots move closer to each other in pairs. At some point, the seventh and eighth roots coalesce and then become a pair of complex conjugate roots. This is shown in the figure below. We can also see how the ninth and tenth roots are approaching each other. For a slightly larger value of , these will morph into a complex conjugate pair.

The above example shows how the extraction of polynomial roots can be an ill-conditioned problem, even for polynomials with well-separated zeros. Wilkinson (1984) described his own personal reaction to this discovery:

*“I regard it as the most traumatic experience in my career as a numerical analyst.”*

** Sources **

Wilkinson, James H., 1963: *Rounding Errors in Algebraic Processes*. Prentice Hall, Englewood Cliffs, New Jersey

Wilkinson, James H., 1984: *The perfidious polynomial*. Studies in Numerical Analysis, ed. by G. H. Golub, pp. 1–28. (Studies in Mathematics, vol. 24). Math. Assoc. America. ISBN 978-0-88385-126-5. PDF

Wikipedia article *Wilkinson’s Polynomial*: http://www.wikipedia.org/