Verfahren zur (numerischen) Berechnung der betragsgrößten Nullstelle x₁ eines (normierten) Polynoms \begin{eqnarray}{p}_{n}(x)={a}_{1}{x}^{n-1}+\cdots +{a}_{n},\end{eqnarray} nach D. Bernoulli.

Das Verfahren ist i.w. identisch mit dem Vektoriterationsverfahren bei Anwendung auf die BegleitmatrixA von P_n.

Nach Vorgabe von fast beliebigen ζ₁, …, ζ_n mit

ζ_n ≠ 0 werden nacheinander die Zahlen \begin{eqnarray}{\zeta }_{n+k}:=-\sum _{j=0}^{n-1}{a}_{n-j}{\zeta }_{k+1},\,k=1,2,\ldots \end{eqnarray} gebildet. Die Quotienten ζ_n+k/ζ_n+k−1 konvergieren dann gegen die betragsgrößte Nullstelle x₁ von P_n.