Root-finding algorithms

If N is an approximation to  , a better approximation can be found by using the Taylor series of the square root function:

The radical or root may be represented by the infinite series:

with <span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/2f9498d60b2319a4ae7c5607794b537c559a976d" data-alt="|x| . This expression can be derived from the binomial series.

• J.M. McNamee: «Numerical Methods for Roots of Polynomials — Part I», Elsevier (2007).
• J.M. McNamee and Victor Pan: «Numerical Methods for Roots of Polynomials — Part II», Elsevier (2013).

An unknown Babylonian mathematician somehow correctly calculated the square root of 2 to three sexagesimal «digits» after the 1, but it is not known exactly how. The Babylonians knew how to approximate a hypotenuse using

The denominator in the fraction corresponds to the nth root. In the case above the denominator is 2, hence the equation specifies that the square root is to be found. The same identity is used when computing square roots with logarithm tables or slide rules.

the Lucas sequence of the first kind U_n(P,Q) is defined by the recurrence relations:

and the characteristic equation of it is:

it has the discriminant   and the roots:

so when we want  , we can choose   and  , and then calculate   using   and  for large value of  .
The most effective way to calculate   and  is:

then when  :

The initialization step of this method is

while the iterative steps read

Then,   (while  ).

The convergence of  , and therefore also of  , is quadratic.

The proof of the method is rather easy. First, rewrite the iterative definition of   as

 .

Then it is straightforward to prove by induction that

This can be used to construct a rational approximation to the square root by beginning with an integer. If   is an integer chosen so   is close to  , and   is the difference whose absolute value is minimized, then the first iteration can be written as:

Using the same example as given with the Babylonian method, let   Then, the first iteration gives

Likewise the second iteration gives

It was once conjectured that all polynomial equations could be solved algebraically (that is, that all roots of a polynomial could be expressed in terms of a finite number of radicals and elementary operations). However, while this is true for third degree polynomials (cubics) and fourth degree polynomials (quartics), the Abel–Ruffini theorem (1824) shows that this is not true in general when the degree is 5 or greater. For example, the solutions of the equation

cannot be expressed in terms of radicals. (cf. quintic equation)

Brent’s method is a combination of the bisection method, the secant method and inverse quadratic interpolation. At every iteration, Brent’s method decides which method out of these three is likely to do best, and proceeds by doing a step according to that method. This gives a robust and fast method, which therefore enjoys considerable popularity.

Ridders’ method is a hybrid method that uses the value of function at the midpoint of the interval to perform an exponential interpolation to the root. This gives a fast convergence with a guaranteed convergence of at most twice the number of iterations as the bisection method.

Many root-finding processes work by interpolation. This consists in using the last computed approximate values of the root for approximating the function by a polynomial of low degree, which takes the same values at these approximate roots. Then the root of the polynomial is computed and used as a new approximate value of the root of the function, and the process is iterated.

Two values allow interpolating a function by a polynomial of degree one (that is approximating the graph of the function by a line). This is the basis of the secant method. Three values define a quadratic function, which approximates the graph of the function by a parabola. This is Muller’s method.

Regula falsi is also an interpolation method, which differs from the secant method by using, for interpolating by a line, two points that are not necessarily the last two computed points.

• Applying Newton’s method to the equation   produces a method that converges quadratically using three multiplications per step:

   

• Another iteration is obtained by Halley’s method, which is the Householder’s method of order two. This converges cubically, but involves five multiplications per iteration:^[]

   , and

   .

• If doing fixed-point arithmetic, the multiplication by 3 and division by 8 can implemented using shifts and adds. If using floating-point, Halley’s method can be reduced to four multiplications per iteration by precomputing   and adjusting all the other constants to compensate:

   , and

   .

The first way of writing Goldschmidt’s algorithm begins

 

  (typically using a table lookup)

 

 

until   is sufficiently close to 1, or a fixed number of iterations. The iterations converge to

 , and

 .

Note that it is possible to omit either   and   from the computation, and if both are desired then   may be used at the end rather than computing it through in each iteration.

A second form, using fused multiply-add operations, begins

  (typically using a table lookup)

 

 

until   is sufficiently close to 0, or a fixed number of iterations. This converges to

 , and

 .

Bracketing methods determine successively smaller intervals (brackets) that contain a root. When the interval is small enough, then a root has been found. They generally use the intermediate value theorem, which asserts that if a continuous function has values of opposite signs at the end points of an interval, then the function has at least one root in the interval. Therefore, they require to start with an interval such that the function takes opposite signs at the end points of the interval. However, in the case of polynomials there are other methods (Descartes’ rule of signs, Budan’s theorem and Sturm’s theorem) for getting information on the number of roots in an interval. They lead to efficient algorithms for real-root isolation of polynomials, which ensure finding all real roots with a guaranteed accuracy.

False position (regula falsi)

The false position method, also called the regula falsi method, is similar to the bisection method, but instead of using bisection search’s middle of the interval it uses the -intercept of the line that connects the plotted function values at the endpoints of the interval, that is

False position is similar to the secant method, except that, instead of retaining the last two points, it makes sure to keep one point on either side of the root. The false position method can be faster than the bisection method and will never diverge like the secant method; however, it may fail to converge in some naive implementations due to roundoff errors that may lead to a wrong sign for ; typically, this may occur if the rate of variation of is large in the neighborhood of the root.

«Heron’s method» redirects here. For the formula used to find the area of a triangle, see Heron’s formula.

More precisely, if is our initial guess of   and is the error in our estimate such that S = (x+ ε)², then we can expand the binomial

and solve for the error term

  since  .

Therefore, we can compensate for the error and update our old estimate as

1. Begin with an arbitrary positive starting value (the closer to the actual square root of , the better).
2. Let x_{n + 1} be the average of and (using the arithmetic mean to approximate the geometric mean).
3. Repeat step 2 until the desired accuracy is achieved.

It can also be represented as:

This algorithm works equally well in the -adic numbers, but cannot be used to identify real square roots with -adic square roots; one can, for example, construct a sequence of rational numbers by this method that converges to +3 in the reals, but to −3 in the 2-adics.

To calculate , where = 125348, to six significant figures, use the rough estimation method above to get

Therefore, ≈ 354.045.

Suppose that x₀ > 0 and S > 0. Then for any natural number n, x_n > 0. Let the relative error in x_n be defined by

-1}» data-class=»mwe-math-fallback-image-inline»> 

Then it can be shown that

And thus that

and consequently that convergence is assured, and quadratic.

Worst case for convergence

<span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/97b965d9c243f4fa62fe825c21a5c171a8e5d8c2" data-alt="{\begin{aligned}S&=1;&x_{0}&=2;&x_{1}&=1.250;&\varepsilon _{1}&=0.250.\\S&=10;&x_{0}&=2;&x_{1}&=3.500;&\varepsilon _{1}&<0.107.\\S&=10;&x_{0}&=6;&x_{1}&=3.833;&\varepsilon _{1}&<0.213.\\S&=100;&x_{0}&=6;&x_{1}&=11.333;&\varepsilon _{1}& 

Thus in any case,

 

<span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/5421d029f479921e0c154a095df9a6314eab837a" data-alt="\varepsilon _{2}<2^{-5} 

<span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/edfb1c400b6a50ee05f9788500fdb114d6201ec8" data-alt="\varepsilon _{3}<2^{-11} 

<span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/3bc199c55d6ae1f04351b7b7aae217d0fb9a5b21" data-alt="\varepsilon _{4}<2^{-23} 

<span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/b3425039f9411fb9109b4174cb04d837dc38fe81" data-alt="\varepsilon _{5}<2^{-47} 

<span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/d73f882414b129be7dc41538e4ce0d60c51891d7" data-alt="\varepsilon _{6}<2^{-95} 

<span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/bb09dcdf99c0c344075f8f6097c36fa976be0469" data-alt="\varepsilon _{7}<2^{-191} 

<span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/ca951a28fe668eb7abc2f3e5af2b968c1689f5c5" data-alt="\varepsilon _{8}<2^{-383} 

Rounding errors will slow the convergence. It is recommended to keep at least one extra digit beyond the desired accuracy of the being calculated to minimize round off error.

Solving an equation f(x) = g(x) is the same as finding the roots of the function h(x) = f(x) – g(x). Thus root-finding algorithms allow solving any equation defined by continuous functions. However, most root-finding algorithms do not guarantee that they will find all the roots; in particular, if such an algorithm does not find any root, that does not mean that no root exists.

Most numerical root-finding methods use iteration, producing a sequence of numbers that hopefully converges towards the root as its limit. They require one or more initial guesses of the root as starting values, then each iteration of the algorithm produces a successively more accurate approximation to the root. Since the iteration must be stopped at some point, these methods produce an approximation to the root, not an exact solution. Many methods compute subsequent values by evaluating an auxiliary function on the preceding values. The limit is thus a fixed point of the auxiliary function, which is chosen for having the roots of the original equation as fixed points, and for converging rapidly to these fixed points.

The behavior of general root-finding algorithms is studied in numerical analysis. However, for polynomials, root-finding study belongs generally to computer algebra, since algebraic properties of polynomials are fundamental for the most efficient algorithms. The efficiency of an algorithm may depend dramatically on the characteristics of the given functions. For example, many algorithms use the derivative of the input function, while others work on every continuous function. In general, numerical algorithms are not guaranteed to find all the roots of a function, so failing to find a root does not prove that there is no root. However, for polynomials, there are specific algorithms that use algebraic properties for certifying that no root is missed, and locating the roots in separate intervals (or disks for complex roots) that are small enough to ensure the convergence of numerical methods (typically Newton’s method) to the unique root so located.

Methods of computing square roots are numerical analysis algorithms for approximating the principal, or non-negative, square root (usually denoted
,
, or
) of a real number. Arithmetically, it means given
, a procedure for finding a number which when multiplied by itself, yields
; algebraically, it means a procedure for finding the non-negative root of the equation
; geometrically, it means given two line segments, a procedure for constructing their geometric mean.

The continued fraction representation of a real number can be used instead of its decimal or binary expansion and this representation has the property that the square root of any rational number (which is not already a perfect square) has a periodic, repeating expansion, similar to how rational numbers have repeating expansions in the decimal notation system.

The method employed depends on what the result is to be used for (i.e. how accurate it has to be), how much effort one is willing to put into the procedure, and what tools are at hand. The methods may be roughly classified as those suitable for mental calculation, those usually requiring at least paper and pencil, and those which are implemented as programs to be executed on a digital electronic computer or other computing device. Algorithms may take into account convergence (how many iterations are required to achieve a specified precision), computational complexity of individual operations (i.e. division) or iterations, and error propagation (the accuracy of the final result).

Procedures for finding square roots (particularly the square root of 2) have been known since at least the period of ancient Babylon in the 17th century BCE. Heron’s method from first century Egypt was the first ascertainable algorithm for computing square root. Modern analytic methods began to be developed after introduction of the Arabic numeral system to western Europe in the early Renaissance. Today, nearly all computing devices have a fast and accurate square root function, either as a programming language construct, a compiler intrinsic or library function, or as a hardware operator, based on one of the described procedures.

Every complex number other than 0 has n different nth roots.

The two square roots of a complex number are always negatives of each other. For example, the square roots of are and , and the square roots of are

If we express a complex number in polar form, then the square root can be obtained by taking the square root of the radius and halving the angle:

A principal root of a complex number may be chosen in various ways, for example

which introduces a branch cut in the complex plane along the positive real axis with the condition , or along the negative real axis with .

Using the first(last) branch cut the principal square root   maps   to the half plane with non-negative imaginary(real) part. The last branch cut is presupposed in mathematical software like Matlab or Scilab.

Roots of unity

The number 1 has n different nth roots in the complex plane, namely

These roots are evenly spaced around the unit circle in the complex plane, at angles which are multiples of  . For example, the square roots of unity are 1 and −1, and the fourth roots of unity are 1,  , −1, and  .

Every complex number has n different nth roots in the complex plane. These are

In polar form, a single nth root may be found by the formula

Here r is the magnitude (the modulus, also called the absolute value) of the number whose root is to be taken; if the number can be written as a+bi then  . Also,   is the angle formed as one pivots on the origin counterclockwise from the positive horizontal axis to a ray going from the origin to the number; it has the properties that     and  

Thus finding nth roots in the complex plane can be segmented into two steps. First, the magnitude of all the nth roots is the nth root of the magnitude of the original number. Second, the angle between the positive horizontal axis and a ray from the origin to one of the nth roots is  , where   is the angle defined in the same way for the number whose root is being taken. Furthermore, all n of the nth roots are at equally spaced angles from each other.

If n is even, a complex number’s nth roots, of which there are an even number, come in additive inverse pairs, so that if a number r₁ is one of the nth roots then r₂ = –r₁ is another. This is because raising the latter’s coefficient –1 to the nth power for even n yields 1: that is, (–r₁)ⁿ = (–1)ⁿ × r₁ⁿ = r₁ⁿ.

As with square roots, the formula above does not define a continuous function over the entire complex plane, but instead has a branch cut at points where θ / n is discontinuous.

In mathematics, an nth root of a number x is a number r which, when raised to the power n, yields x:

where n is a positive integer, sometimes called the degree of the root. A root of degree 2 is called a square root and a root of degree 3, a cube root. Roots of higher degree are referred by using ordinal numbers, as in fourth root, twentieth root, etc. The computation of an th root is a root extraction.

For example, 3 is a square root of 9, since 3² = 9, and −3 is also a square root of 9, since (−3)² = 9.

Any non-zero number considered as a complex number has different complex th roots, including the real ones (at most two). The th root of 0 is zero for all positive integers , since 0ⁿ = 0. In particular, if is even and is a positive real number, one of its th roots is real and positive, one is negative, and the others (when ) are non-real complex numbers; if is even and is a negative real number, none of the th roots is real. If is odd and is real, one th root is real and has the same sign as , while the other () roots are not real. Finally, if is not real, then none of its th roots are real.

Roots of real numbers are usually written using the radical symbol or radix
, with
denoting the positive square root of if is positive; for higher roots,
denotes the real th root if is odd, and the positive nth root if is even and is positive. In the other cases, the symbol is not commonly used as being ambiguous. In the expression
, the integer n is called the index and is called the radicand.

When complex th roots are considered, it is often useful to choose one of the roots, called principal root, as a principal value. The common choice is to choose the principal th root of as the th root with the greatest real part, and when there are two (for real and negative), the one with a positive imaginary part. This makes the th root a function that is real and positive for real and positive, and is continuous in the whole complex plane, except for values of that are real and negative.

A difficulty with this choice is that, for a negative real number and an odd index, the principal th root is not the real one. For example,
has three cube roots,
,
and
The real cube root is
and the principal cube root is

Roots can also be defined as special cases of exponentiation, where the exponent is a fraction:

Roots are used for determining the radius of convergence of a power series with the root test. The th roots of 1 are called roots of unity and play a fundamental role in various areas of mathematics, such as number theory, theory of equations, and Fourier transform.

Although all root-finding algorithms proceed by iteration, an iterative root-finding method generally uses a specific type of iteration, consisting of defining an auxiliary function, which is applied to the last computed approximations of a root for getting a new approximation. The iteration stops when a fixed point (up to the desired precision) of the auxiliary function is reached, that is when the new computed value is sufficiently close to the preceding ones.

Newton’s method (and similar derivative-based methods)

Newton’s method assumes the function f to have a continuous derivative. Newton’s method may not converge if started too far away from a root. However, when it does converge, it is faster than the bisection method, and is usually quadratic. Newton’s method is also important because it readily generalizes to higher-dimensional problems. Newton-like methods with higher orders of convergence are the Householder’s methods. The first one after Newton’s method is Halley’s method with cubic order of convergence.

Replacing the derivative in Newton’s method with a finite difference, we get the secant method. This method does not require the computation (nor the existence) of a derivative, but the price is slower convergence (the order is approximately 1.6 (golden ratio)). A generalization of the secant method in higher dimensions is Broyden’s method.

If we use a polynomial fit to remove the quadratic part of the finite difference used in the Secant method, so that it better approximates the derivative, we obtain Steffensen’s method, which has quadratic convergence, and whose behavior (both good and bad) is essentially the same as Newton’s method but does not require a derivative.

Fixed point iteration method

We can use the fixed-point iteration to find the root of a function. Given a function   which we have set to zero to find the root ( ), we rewrite the equation in terms of   so that   becomes   (note, there are often many   functions for each   function). Next, we relabel each side of the equation as   so that we can perform the iteration. Next, we pick a value for   and perform the iteration until it converges towards a root of the function. If the iteration converges, it will converge to a root. The iteration will only converge if <span data-src="https://wikimedia.org/api/rest_v1/media/math/render/svg/e72ed01b48be0f78bea75e50a39c992d0a02b94e" data-alt="{\displaystyle |g'(root)| .

 ,

 ,

 ,

 , or

 .

The appearance of complex values in interpolation methods can be avoided by interpolating the inverse of f, resulting in the inverse quadratic interpolation method. Again, convergence is asymptotically faster than the secant method, but inverse quadratic interpolation often behaves poorly when the iterates are not close to the root.

Quadratic irrationals (numbers of the form  , where a, b and c are integers), and in particular, square roots of integers, have periodic continued fractions. Sometimes what is desired is finding not the numerical value of a square root, but rather its continued fraction expansion, and hence its rational approximation. Let S be the positive number for which we are required to find the square root. Then assuming a to be a number that serves as an initial guess and r to be the remainder term, we can write   Since we have  , we can express the square root of S as

By applying this expression for   to the denominator term of the fraction, we have

Proceeding this way, we get a generalized continued fraction for the square root as
 

Step 2 is to reduce the continued fraction from the bottom up, one denominator at a time, to yield a rational fraction whose numerator and denominator are integers. The reduction proceeds thus (taking the first three denominators):

Finally (step 3), divide the numerator by the denominator of the rational fraction to obtain the approximate value of the root:

  rounded to three digits of precision.

The actual value of is 1.41 to three significant digits. The relative error is 0.17%, so the rational fraction is good to almost three digits of precision. Taking more denominators gives successively better approximations: four denominators yields the fraction  , good to almost 4 digits of precision, etc.

In general, the larger the denominator of a rational fraction, the better the approximation. It can also be shown that truncating a continued fraction yields a rational fraction that is the best approximation to the root of any fraction with denominator less than or equal to the denominator of that fraction — e.g., no fraction with a denominator less than or equal to 70 is as good an approximation to as 99/70.

Using Newton’s method

The th root of a number can be computed with Newton’s method, which starts with an initial guess and then iterates using the recurrence relation

until the desired precision is reached. For computational efficiency, the recurrence relation is commonly rewritten

This allows to have only one exponentiation, and to compute once for all the first factor of each term.

The approximation is accurate to 25 decimal places and is good for 51.

Newton’s method can be modified to produce various generalized continued fractions for the nth root. For example,

Digit-by-digit calculation of principal roots of decimal (base 10) numbers

Write the original number in decimal form. The numbers are written similar to the long division algorithm, and, as in long division, the root will be written on the line above. Now separate the digits into groups of digits equating to the root being taken, starting from the decimal point and going both left and right. The decimal point of the root will be above the decimal point of the radicand. One digit of the root will appear above each group of digits of the original number.

1. Starting on the left, bring down the most significant (leftmost) group of digits not yet used (if all the digits have been used, write «0» the number of times required to make a group) and write them to the right of the remainder from the previous step (on the first step, there will be no remainder). In other words, multiply the remainder by   and add the digits from the next group. This will be the current value c.
2. Find p and x, as follows:
   • Let   be the part of the root found so far, ignoring any decimal point. (For the first step,  ).
   • Determine the greatest digit   such that  .
   • Place the digit   as the next digit of the root, i.e., above the group of digits you just brought down. Thus the next p will be the old p times 10 plus x.
3. Subtract   from   to form a new remainder.
4. If the remainder is zero and there are no more digits to bring down, then the algorithm has terminated. Otherwise go back to step 1 for another iteration.

Find the square root of 152.2756.

  1 2. 3 4    / \/ 01 52.27 56

 01 100·1·00·12 + 101·2·01·11 ≤ 1 < 100·1·00·22 + 101·2·01·21 x = 1  01  y = 100·1·00·12 + 101·2·01·11 = 1 + 0 = 1 00 52 100·1·10·22 + 101·2·11·21 ≤ 52 < 100·1·10·32 + 101·2·11·31 x = 2  00 44  y = 100·1·10·22 + 101·2·11·21 = 4 + 40 = 44 08 27 100·1·120·32 + 101·2·121·31 ≤ 827 < 100·1·120·42 + 101·2·121·41 x = 3  07 29  y = 100·1·120·32 + 101·2·121·31 = 9 + 720 = 729 98 56 100·1·1230·42 + 101·2·1231·41 ≤ 9856 < 100·1·1230·52 + 101·2·1231·51 x = 4  98 56  y = 100·1·1230·42 + 101·2·1231·41 = 16 + 9840 = 9856 00 00 Algorithm terminates: Answer is 12.34

Find the cube root of 4192 to the nearest hundredth.

  1 6. 1 2 4 3 / \/ 004 192.000 000 000

 004 100·1·00·13 + 101·3·01·12 + 102·3·02·11 ≤ 4 < 100·1·00·23 + 101·3·01·22 + 102·3·02·21 x = 1  001  y = 100·1·00·13 + 101·3·01·12 + 102·3·02·11 = 1 + 0 + 0 = 1 003 192 100·1·10·63 + 101·3·11·62 + 102·3·12·61 ≤ 3192 < 100·1·10·73 + 101·3·11·72 + 102·3·12·71 x = 6  003 096  y = 100·1·10·63 + 101·3·11·62 + 102·3·12·61 = 216 + 1,080 + 1,800 = 3,096 096 000 100·1·160·13 + 101·3·161·12 + 102·3·162·11 ≤ 96000 < 100·1·160·23 + 101·3·161·22 + 102·3·162·21 x = 1  077 281  y = 100·1·160·13 + 101·3·161·12 + 102·3·162·11 = 1 + 480 + 76,800 = 77,281 018 719 000 100·1·1610·23 + 101·3·1611·22 + 102·3·1612·21 ≤ 18719000 < 100·1·1610·33 + 101·3·1611·32 + 102·3·1612·31 x = 2  015 571 928  y = 100·1·1610·23 + 101·3·1611·22 + 102·3·1612·21 = 8 + 19,320 + 15,552,600 = 15,571,928 003 147 072 000 100·1·16120·43 + 101·3·16121·42 + 102·3·16122·41 ≤ 3147072000 < 100·1·16120·53 + 101·3·16121·52 + 102·3·16122·51 x = 4 The desired precision is achieved: The cube root of 4192 is about 16.12

The principal nth root of a positive number can be computed using logarithms. Starting from the equation that defines r as an nth root of x, namely   with x positive and therefore its principal root r also positive, one takes logarithms of both sides (any base of the logarithm will do) to obtain

The root r is recovered from this by taking the antilog:

(Note: That formula shows b raised to the power of the result of the division, not b multiplied by the result of the division.)

Assume that   is rational. That is, it can be reduced to a fraction  , where and are integers without a common factor.

This means that  .

Since   and  ,  .

This means that   and thus,  . This implies that   is an integer. Since x is not a perfect nth power, this is impossible. Thus   is irrational.

If S < 0, then its principal square root is

If S = a+bi where a and b are real and b ≠ 0, then its principal square root is

is the modulus of S. The principal square root of a complex number is defined to be the root with the non-negative real part.

1. There is no factor of the radicand that can be written as a power greater than or equal to the index.
2. There are no fractions under the radical sign.
3. There are no radicals in the denominator.

Simplifying radical expressions involving nested radicals can be quite difficult. It is not obvious for instance that:

The above can be derived through:

Let  , with and coprime and positive integers. Then   is rational if and only if both   and   are integers, which means that both and are nth powers of some integer.

Expressing the degree of an nth root in its exponent form, as in  , makes it easier to manipulate powers and roots. If   is a non-negative real number,

Every non-negative number has exactly one non-negative real nth root, and so the rules for operations with surds involving non-negative radicands   and   are straightforward within the real numbers:

Subtleties can occur when taking the nth roots of negative or complex numbers. For instance:

  but, rather,  

Since the rule   strictly holds for non-negative real radicands only, its application leads to the inequality in the first step above.

The Poincaré–Miranda theorem gives a criterion for the existence of a root in a rectangle, but it is hard to verify, since it requires to evaluate the function on the entire boundary of the triangle.

A number is represented in a floating point format as   which is also called scientific notation. Its square root is   and similar formulae would apply for cube roots and logarithms. On the face of it, this is no improvement in simplicity, but suppose that only an approximation is required: then just   is good to an order of magnitude. Next, recognise that some powers, , will be odd, thus for 3141.59 = 3.1415910 rather than deal with fractional powers of the base, multiply the mantissa by the base and subtract one from the power to make it even. The adjusted representation will become the equivalent of 31.415910 so that the square root will be 10.

A table with only three entries could be enlarged by incorporating additional bits of the mantissa. However, with computers, rather than calculate an interpolation into a table, it is often better to find some simpler calculation giving equivalent results. Everything now depends on the exact details of the format of the representation, plus what operations are available to access and manipulate the parts of the number. For example, Fortran offers an EXPONENT(x) function to obtain the power. Effort expended in devising a good initial approximation is to be recouped by thereby avoiding the additional iterations of the refinement process that would have been needed for a poor approximation. Since these are few (one iteration requires a divide, an add, and a halving) the constraint is severe.

So for a 32-bit single precision floating point number in IEEE format (where notably, the power has a bias of 127 added for the represented form) you can get the approximate logarithm by interpreting its binary representation as a 32-bit integer, scaling it by  , and removing a bias of 127, i.e.

For example, 1.0 is represented by a hexadecimal number 0x3F800000, which would represent   if taken as an integer. Using the formula above you get  , as expected from  . In a similar fashion you get 0.5 from 1.5 (0x3FC00000).

/* Assumes that float is in the IEEE 754 single precision floating point format */            /* Convert type, preserving bit pattern */ * ((((val.i / 2^m) - b) / 2) + b) * 2^m = ((val.i - 2^m) / 2) + ((b + 1) / 2) * 2^m) * b = exponent bias * m = number of mantissa bits    /* Subtract 2^m. */  /* Divide by 2. */    /* Add ((b + 1) / 2) * 2^m. */ /* Interpret again as float */

The three mathematical operations forming the core of the above function can be expressed in a single line. An additional adjustment can be added to reduce the maximum relative error. So, the three operations, not including the cast, can be rewritten as

Reciprocal of the square root

                            /* The next line can be repeated any number of times to increase accuracy */             

An nth root of a number x, where n is a positive integer, is any of the n real or complex numbers r whose nth power is x:

Every positive real number x has a single positive nth root, called the principal nth root, which is written  . For n equal to 2 this is called the principal square root and the n is omitted. The nth root can also be represented using exponentiation as x^1/n.

For even values of n, positive numbers also have a negative nth root, while negative numbers do not have a real nth root. For odd values of n, every negative number x has a real negative nth root. For example, −2 has a real 5th root,   but −2 does not have any real 6th roots.

Every non-zero number x, real or complex, has n different complex number nth roots. (In the case x is real, this count includes any real nth roots.) The only complex root of 0 is 0.

The nth roots of almost all numbers (all integers except the nth powers, and all rationals except the quotients of two nth powers) are irrational. For example,

All nth roots of rational numbers are algebraic numbers, and all nth roots of integers are algebraic integers.

A square root of a number x is a number r which, when squared, becomes x:

Every positive real number has two square roots, one positive and one negative. For example, the two square roots of 25 are 5 and −5. The positive square root is also known as the principal square root, and is denoted with a radical sign:

Since the square of every real number is nonnegative, negative numbers do not have real square roots. However, for every negative real number there are two imaginary square roots. For example, the square roots of −25 are 5i and −5i, where i represents a number whose square is .

A cube root of a number x is a number r whose cube is x:

Every real number x has exactly one real cube root, written  . For example,

  and  

Every real number has two additional complex cube roots.