On a simultaneous method of Newton-Weierstrass’ type for finding all zeros of a polynomial

Combining a suitable two-point iterative method for solving nonlinear equations and Weierstrass’ correction, a new iterative method for simultaneous finding all zeros of a polynomial is derived. It is proved that the proposed method possesses a cubic convergence locally. Numerical examples demonstrate a good convergence behavior of this method in a global sense. It is shown that its computational efficiency is higher than the existing derivative-free methods.     

Efficient methods for the inclusion of polynomial zeros

Using a suitable zero-relation and the inclusion isotonicity property, new interval iterative methods for the simultaneous inclusion of simple complex zeros of a polynomial are derived. These methods produce disks in the complex plane that contain the polynomial zeros in each iteration, providing in this manner an information about upper error bounds of approximations. Starting from the basic method of the fourth order, two accelerated methods with Newton’s and Halley’s corrections, having the order of convergence five and six respectively, are constructed. This increase of the convergence rate is obtained without any additional operations, which means that the methods with corrections are very efficient. The convergence analysis of the basic method and the methods with corrections is performed under computationally verifiable initial conditions, which is of practical importance. Two numerical examples are presented to demonstrate the convergence behavior of the proposed interval methods.     

Relationships between different types of initial conditions for simultaneous root finding methods

The construction of initial conditions of an iterative method is one of the most important problems in solving nonlinear equations. In this paper, we obtain relationships between different types of initial conditions that guarantee the convergence of iterative methods for simultaneously finding all zeros of a polynomial. In particular, we show that any local convergence theorem for a simultaneous method can be converted into a convergence theorem with computationally verifiable initial conditions which is of practical importance. Thus, we propose a new approach for obtaining semilocal convergence results for simultaneous methods via local convergence results.     

Weierstrass method for quaternionic polynomial root‐finding

Quaternions, introduced by Hamilton in 1843 as a generalization of complex numbers, have found, in more recent years, a wealth of applications in a number of different areas that motivated the design of efficient methods for numerically approximating the zeros of quaternionic polynomials. In fact, one can find in the literature recent contributions to this subject based on the use of complex techniques, but numerical methods relying on quaternion arithmetic remain scarce. In this paper, we propose a Weierstrass‐like method for finding simultaneously all the zeros of unilateral quaternionic polynomials. The convergence analysis and several numerical examples illustrating the performance of the method are also presented.     

Safe convergence of simultaneous methods for polynomial zeros

The theory of point estimation treating the initial conditions for the safe convergence of iterative processes for the simultaneous determination of polynomial zeros is considered. A general approach which makes use of corrections appearing in iterative formulas is given and demonstrated in the case of three well known methods without derivatives and based on Weierstrass’ corrections. The established convergence conditions are of practical importance since they depend only on available data: coefficients of a polynomial and initial approximations to the zeros. This revised version was published online in June 2006 with corrections to the Cover Date.     

A semilocal convergence result for Newton’s method under generalized conditions of Kantorovich

From Kantorovich’s theory we establish a general semilocal convergence result for Newton’s method based fundamentally on a generalization required to the second derivative of the operator involved. As a consequence, we obtain a modification of the domain of starting points for Newton’s method and improve the a priori error estimates. Finally, we illustrate our study with an application to a special case of conservative problems.     

The self-validated method for polynomial zeros of high efficiency

The improved iterative method of Newton’s type for the simultaneous inclusion of all simple complex zeros of a polynomial is proposed. The presented convergence analysis, which uses the concept of the R-order of convergence of mutually dependent sequences, shows that the convergence rate of the basic third order method is increased from 3 to 6 using Ostrowski’s corrections. The new inclusion method with Ostrowski’s corrections is more efficient compared to all existing methods belonging to the same class. To demonstrate the convergence properties of the proposed method, two numerical examples are given.     

A unified semilocal convergence analysis of a family of iterative algorithms for computing all zeros of a polynomial simultaneously

In this paper, we first present a family of iterative algorithms for simultaneous determination of all zeros of a polynomial. This family contains two well-known algorithms: Dochev-Byrnev’s method and Ehrlich’s method. Second, using Proinov’s approach to studying convergence of iterative methods for polynomial zeros, we provide a semilocal convergence theorem that unifies the results of Proinov (Appl. Math. Comput. 284: 102–114, 2016) for Dochev-Byrnev’s and Ehrlich’s methods.     

Concerning the semilocal convergence of Newton’s method and convex majorants

We approximate a locally unique solution of an equation on a Banach space setting using Newton’smethod.Motivated by the work by Ferreira and Svaiter [5] but using more precise majorization sequences, and under the same computational cost we provide: a larger convergence region; finer error bounds on the distances involved, and an at least as precise information on the location of the solution than in [5]. The results can also compare favorably to the corresponding ones given byWang in [10]. Finally we complete the study with two concrete applications.      

A new semilocal convergence theorem for a fast iterative method with nondifferentiable operators

A new semilocal convergence theorem for a fast iterative method in Banach spaces is provided for approximating a solution of a nondifferentiable operator equation. A condition for divided differences of order one is considered in this paper, which generalizes the usual ones, i.e., Lipschitz continuous or Hölder continuous conditions. Note that no conditions of divided differences of order two are used. Therefore our results are of theoretical and practical interest. Finally, a numerical example is provided to show that the new iterative method compares favorably with earlier ones.     

A note on some improvements of the simultaneous methods for determination of polynomial zeros

Applying Gauss-Seidel approach to the improvements of two simultaneous methods for finding polynomial zeros, presented in [9], two iterative methods with faster convergence are obtained. The lower bounds of the R-order of convergence for the accelerated methods are given. The improved methods and their accelerated modifications are discussed in view of the convergence order and the number of numerical operations. The considered methods are illustrated numerically in the example of an algebraic equation.     

Finding roots of a real polynomial simultaneously by means of Bairstow's method

Aberth's method for finding the roots of a polynomial was shown to be robust. However, complex arithmetic is needed in this method even if the polynomial is real, because it starts with complex initial approximations. A novel method is proposed for real polynomials that does not require any complex arithmetic within iterations. It is based on the observation that Aberth's method is a systematic use of Newton's method. The analogous technique is then applied to Bairstow's procedure in the proposed method. As a result, the method needs half the computations per iteration than Aberth's method. Numerical experiments showed that the new method exhibited a competitive overall performance for the test polynomials.     

A higher order family for the simultaneous inclusion of multiple zeros of polynomials

Starting from a suitable fixed point relation, a new family of iterative methods for the simultaneous inclusion of multiple complex zeros in circular complex arithmetic is constructed. The order of convergence of the basic family is four. Using Newtons and Halleys corrections, we obtain families with improved convergence. Faster convergence of accelerated methods is attained with only few additional numerical operations, which provides a high computational efficiency of these methods. Convergence analysis of the presented methods and numerical results are given. AMS subject classification 65H05, 65G20, 30C15     

Local convergence of a secant type method for solving least squares problems

Local convergence of a secant type iterative method for approximating a solution of nonlinear least squares problems is investigated in this paper. The radius of convergence is determined as well as usable error estimates. Numerical examples are also provided.     

Recurrence relations for semilocal convergence of a Newton-like method in Banach spaces

The aim of this paper is to establish the semilocal convergence of a multipoint third order Newton-like method for solving F(x)=0 in Banach spaces by using recurrence relations. The convergence of this method is studied under the assumption that the second Fréchet derivative of F satisfies Hölder continuity condition. This continuity condition is milder than the usual Lipschitz continuity condition. A new family of recurrence relations are defined based on the two new constants which depend on the operator F. These recurrence relations give a priori error bounds for the method. Two numerical examples are worked out to demonstrate the applicability of the method in cases where the Lipschitz continuity condition over second derivative of F fails but Hölder continuity condition holds.     

A new semi-local convergence theorem for the inexact Newton methods

We present a new semi-local convergence theorem for the inexact Newton methods in the assumption that the derivative satisfies some kind of weak Lipschitz conditions. As special cases of our main result we re-obtain some well-known convergence theorems for Newton methods.     

A new simultaneous method of fourth order for finding complex zeros in circular interval arithmetic

Starting from disjoint disks which contain polynomial complex zeros, the new iterative interval method for simultaneous finding of inclusive disks for complex zeros is formulated. The convergence theorem and the conditions for convergence are considered, and the convergence is shown to be fourth. Numerical examples are included.     

A computational comparison of the first nine members of a determinantal family of root-finding methods

For each natural number m greater than one, and each natural number k less than or equal to m, there exists a root-finding iteration function, B_m^(k) defined as the ratio of two determinants that depend on the first m−k derivatives of the given function. This infinite family is derived in Kalantari (J. Comput. Appl. Math. 126 (2000) 287–318) and its order of convergence is analyzed in Kalantari (BIT 39 (1999) 96–109). In this paper we give a computational study of the first nine root-finding methods. These include Newton, secant, and Halley methods. Our computational results with polynomials of degree up to 30 reveal that for small degree polynomials B_m^(k−1) is more efficient than B_m^(k), but as the degree increases, B_m^(k) becomes more efficient than B_m^(k−1). The most efficient of the nine methods is B₄⁽⁴⁾, having theoretical order of convergence equal to 1.927. Newton's method which is often viewed as the method of choice is in fact the least efficient method.     

On the simultaneous determination of the zeros of an analytic function inside a simple smooth closed contour in the complex plane

A family of higher-order iterative methods for the simultaneous determination of all simple or multiple zeros of an analytic function (inside a simple smooth closed contour) is obtained using earlier results of the author. With the help of circular arithmetic, the interval variant of this family is proposed. Many parallel iterative methods of the literature are special cases of this family.