Ten methods to bound multiple roots of polynomials

  • Authors:
  • Siegfried M. Rump

  • Affiliations:
  • Institut fur Informatik III, Technical University Hamburg-Harburg, Schwarzenbergstr. 95, Hamburg 21071, Germany

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2003

Quantified Score

Hi-index 7.31

Visualization

Abstract

Given a univariate polynomial P with a k-fold multiple root or a k-fold root cluster near some z, we discuss nine different methods to compute a disc near z which either contains exactly or contains at least k roots of P. Many of the presented methods are known, some are new. We are especially interested in the behaviour of methods when implemented in a rigorous way, that is when taking into account all possible effects of rounding errors. In other words every result shall be mathematically correct. We display extensive test sets comparing the methods under different circumstances. Based on the results we present a tenth, a hybrid method combining five of the previous methods which, for given z, (i) detects the number k of roots near z and (ii) computes an including disc with in most cases a radius of the order of the numerical sensitivity of the root cluster. Therefore, the resulting discs are numerically nearly optimal.