Affiliation:
1. Univ. of Toronto, Toronto, Ont., Canada
Abstract
Formulas developed originally by Weierstrass have been used since the 1960s by many others for the simultaneous determination of all the roots of a polynomial. Convergence to simple roots is quadratic, but individual approximations to a multiple root converge only linearly. However, it is shown here that the mean of such individual approximations converges quadratically to that root. This result, along with some detail about the behavior of such approximations in the neighborhood of the multiple root, suggests a new approach to the design of polynomial rootfinders. It should also be noted that the technique is well suited to take advantage of a parallel environment. This article first provides the relevant mathematical results: a short derivation of the formulas, convergence proofs, an indication of the behavior near a multiple root, and some error bounds. It then provides the outline of an algorithm based on these results, along with some graphical and numerical results to illustrate the major theoretical points. Finally, a new program based on this algorithm, but with a more efficient way of choosing starting values, is described and then compared with corresponding programs from IMSL and NAG with good results. This program is available from Mathon (combin@cs.utoronto.ca).
Publisher
Association for Computing Machinery (ACM)
Subject
Applied Mathematics,Software
Reference19 articles.
1. Iteration methods for finding all zeros of a polynomial simultaneously;ABERTH O.;Math. Comput.,1973
2. ABERTH O. 1988. Precise Numerical Analysis. William C. Brown Dubuque Iowa. ABERTH O. 1988. Precise Numerical Analysis. William C. Brown Dubuque Iowa.
3. A stopping criterion for polynomial root finding
4. A modified Newton method for the simultaneous approximate calculation of all the roots of a given algebraic equation;HE;Phys. Math. J. Bulgarian Acad. Sci.,1962
5. DURANI) E. 1960. Solutions Numdriques des t~quations Algdbriques. Vol. 1 t~quations du Type F(x) = O Racines d'un Polyn6rne. Masson Paris. DURANI) E. 1960. Solutions Numdriques des t~quations Algdbriques. Vol. 1 t~quations du Type F(x) = O Racines d'un Polyn6rne. Masson Paris.
Cited by
8 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献