To ﬁnd all the real and complex roots of a polynomial with real coeﬃcients, i.e. calculate the zeros of
The user can supply error bounds on the coeﬃcients 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.