Bernoulli method

A method for finding the real root of algebraic equations of the type

(*)

with the largest modulus (absolute value). The method was proposed by D. Bernoulli [1] and is based on the following principle. Let be random numbers and let the values of be calculated by the following difference equation:

In general, as , the expression tends to the value of the root of equation (*) with the largest modulus.

References

[1]	D. Bernoulli, Comm. Acad. Sci. Imper. Petropolitanae , 3 , Petropolis (1732) pp. 62–69
[2]	Sir Edmund Whittaker, G. Robinson, "The calculus of observations. A treatise on numerical mathematics" (4th ed.) Blackie & Son (1954) p.99. Zbl 0058.33603

Comments

The expression tends to the root of largest modulus only if this root is simple and if no other roots of the same maximum modulus occur. If not, then one can proceed as follows (cf. [a1]). Let there be one (not necessarily real) root of largest modulus, with multiplicity . Then the solution of the linear difference equation for is

as . This expression is the solution of a linear difference equation with characteristic polynomial , so that satisfies the linear difference equation

(a1)

as . If , then this equation can be solved explicitly in , otherwise (a1) is rewritten with replaced by and the resulting equations are solved in .

If two roots and (of multiplicity and ) of the same maximum modulus occur, then satisfies, as , a linear difference equation with as characteristic polynomial, and (a1) must be replaced by the corresponding difference equation of order .

References

[a1]	F.B. Hildebrand, "Introduction to numerical analysis" , McGraw-Hill (1956)
[a2]	A.S. Householder, "The numerical treatment of a single nonlinear equation" , McGraw-Hill (1970)