Version 1.1.0

Code Download

  • Single
  • Double

PA17 Real coefficients: all roots by the method of Madsen and Reid

To find all the real and complex roots of a polynomial with real coefficients, i.e. calculate the zeros of

a0 + a1x + a2x2 + ... + a nxn = 0

The user can supply error bounds on the coefficients of the polynomial and the routine returns bounds on the moduli of the errors in the roots. The roots are found by the method of Madsen and error bounds by the application of Rouche’s theorem.

Precision: At least 8-byte arithmetic is recommended.