On the Bernoulli–Euler–Lagrange–Aitken numerical method for roots of polynomials

The article presents a development of the Euler–Lagrange method for calculation of all roots of an arbitrary polynomial P(z) with complex coefficients based on the calculation of the limits of ratios of determinants (as in the Bernoulli–Aitken–Nikiporets methods) built by means of the Taylor and Laurent series coefficients for the function P′(z) / P(z).