TOPICS
Search

Newton's Iteration


Newton's iteration is an algorithm for computing the square root sqrt(n) of a positive real number n via the recurrence equation

 x_(k+1)=1/2(x_k+n/(x_k)),
(1)

where x_0=1. For n!=1, the recurrence equation converges to sqrt(n) with quadratic convergence order. For n=1, every iterate equals 1.

Newton's iteration is simply an application of Newton's method for solving the equation

 x^2-n=0.
(2)

For example, when applied numerically, the first few iterations to Pythagoras's constant sqrt(2)=1.4142... are 1, 1.5, 1.41667, 1.41422, 1.41421, ....

The first few approximants x_0, x_1, ... to sqrt(n) are given by

 1,1/2(1+n),(1+6n+n^2)/(4(n+1)),(1+28n+70n^2+28n^3+n^4)/(8(1+n)(1+6n+n^2)),....
(3)

For nonnegative integer k, these are given by the analytic formula

x_k=sqrt(n)((1+sqrt(n))^(2^k)+(1-sqrt(n))^(2^k))/((1+sqrt(n))^(2^k)-(1-sqrt(n))^(2^k)).
(4)

For n!=1, this also equals sqrt(n)coth(2^ktanh^(-1)sqrt(n)), using the principal value of the inverse hyperbolic tangent. The formula can be derived by noting that, for k>=1, the recurrence equation can be written as

 (x_k-sqrt(n))/(x_k+sqrt(n))=((x_(k-1)-sqrt(n))/(x_(k-1)+sqrt(n)))^2,
(5)

which has the closed-form solution

 (x_k-sqrt(n))/(x_k+sqrt(n))=((1-sqrt(n))/(1+sqrt(n)))^(2^k).
(6)

Solving for x_k then gives the solution derived above.

The following table summarizes the first few iterates for small positive integer n

nOEISx_0, x_1, ...
11, 1, 1, 1, 1, 1, 1, 1, ...
2A001601/A0510091, 3/2, 17/12, 577/408, 665857/470832, ...
3A002812/A0715791, 2, 7/4, 97/56, 18817/10864, 708158977/408855776, ...

See also

Newton's Method, Square Root, Square Root Algorithms, Wolfram's Iteration

Explore with Wolfram|Alpha

References

Sloane, N. J. A. Sequences A001601/M3042, A002812/M1817, A051009, A071579 in "The On-Line Encyclopedia of Integer Sequences."Wolfram, S. A New Kind of Science. Champaign, IL: Wolfram Media, p. 913, 2002.

Referenced on Wolfram|Alpha

Newton's Iteration

Cite this as:

Weisstein, Eric W. "Newton's Iteration." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/NewtonsIteration.html

Subject classifications