Square root

From Wikipedia, the free encyclopedia.

In mathematics, the principal square root of a non-negative real number

is denoted

and represents the non-negative real number whose square (the result of multiplying the number by itself) is

For example, since 3² = 3 × 3 = 9.

This example suggests how square roots can arise when solving quadratic equations such as or, more generally,

There are two solutions to the square root of a non-zero number. For a positive real number, the two square roots are the principle square root and the negative square root. For negative real numbers, the concept of imaginary and complex numbers has been developed to provide a mathematical framework to deal with the results.

Square roots of positive integers are often irrational numbers, i.e., numbers not expressible as a quotient of two integers. For example, cannot be written exactly as m/n, where n and m are integers. Nonetheless, it is exactly the length of the diagonal of a square with side length 1.

The discovery that is irrational is attributed to Hippasus, a disciple of Pythagoras.

The square root symbol (√) was first used during the 16th century. It has been suggested that it originated as an altered form of lowercase r, representing the Latin radix (meaning "root").

Table of contents

1 Properties
2 Common errors involving equations with the principal square root function
3 Computing square roots

3.1 Calculators
3.2 Babylonian method
3.3 Simple approximation

3.3.1 Construction

3.4 Bakhshali approximation

3.4.2 Bakhshali approximation example

3.5 An exact "long-division like" algorithm
3.6 Rough estimation
3.7 More accurate rough estimation
3.8 Inverse square root approximation

3.8.3 Steps for Inverse square root approximation
3.8.4 Properties

3.9 Square roots using Newton iteration

3.9.5 First method
3.9.6 Second method
3.9.7 Example
3.9.8 Comparison

3.10 Pell's equation
3.11 Finding square roots using mental arithmetic
3.12 Finding square roots using basic arithmetic
3.13 Continued fraction methods

3.13.9 Steps for finding continued fractions

3.14 Quadratic equation method

3.14.10 Example using the quadratic equation method

4 Square roots of complex numbers
5 Square roots of matrices and operators
6 Infinitely nested square roots
7 Square roots of the first 20 positive integers
8 See also
9 External links

Properties

The principal square root function is a function which maps the non-negative real domain R⁺∪{0} into the non-negative real codomain R⁺∪{0}.
The principal square root function always returns a single unique value.
There are only two solutions to the equation The solution set is { 0,1 }.
To obtain both roots of a positive number, take the value given by the principal square root function as the first root (root₁) and obtain the second root (root₂) by subtracting the first root from zero (ie root₂ = 0 - root₁).
The following important properties of the square root functions are valid for all positive real numbers xx and yy:

for every real number xx (see absolute value)

The square root function generally maps rational numbers to algebraic numbers; is rational if and only if xx is a rational number which, after cancelling, is a ratio of two perfect squares. In particular, is irrational.
In geometrical terms, the square root function maps the area of a square to its side length.
Suppose that xx and aa are reals, and that , and we want to find xx. A common mistake is to "take the square root" and deduce that . This is incorrect, because the principal square root of is not xx, but the absolute value , one of our above rules. Thus, all we can conclude is that , or equivalently .
In calculus, for instance when proving that the square root function is continuous or differentiable or when computing certain limitss, the following identity often comes handy:

It is valid for all non-negative numbers xx and yy which are not both zero.

The function has the following graph, made up of half a parabola lying on its side:

The function is continuous for all non-negative xx, and differentiable for all positive xx (it is not differentiable for x = 0x=0 since the slope of the tangent there is ∞). Its derivative is given by

Its Taylor series about x = 1x=1 can be found using the binomial theorem:

for

Common errors involving equations with the principal square root function

Not recognising as the principal square root of xx in a squaring step followed by principal square root step.
In details we have:
This is because and so
The error comes in taking the principal square root of a square of −1.
Not taking into account false root(s) when squaring equations with the principal square root function
squaring both sides which introduces the false root x=0.38197
move x from the right hand side to the left hand side
Using the quadratic equation formula, we get the solutions
or
But when we substitute back
Error!

Computing square roots

Calculators

Pocket calculatorss typically implement good routines to compute the exponential function and the natural logarithm, and then compute the square root of xx using the identity

The same identity is exploited when computing square roots with logarithm tables or slide rules.

Babylonian method

A commonly used algorithm for approximating

is known as the "Babylonian method" and is based on Newton's method. It proceeds as follows:

start with an arbitrary positive start value xx (the closer to the root the better)
replace xx by the average of xx and r/x
go to 2

This is a quadratically convergent algorithm, which means that the number of correct digits of rr roughly doubles with each step.

This could be represented as

where

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

Simple approximation

The simple approximation is rather simple and can be highly inaccurate. The amount of inaccuracy for this approximation is dependent on the value of the expression d/2N , the larger the value of this expression, the more inaccurate the value of the approximated result.

Construction

If N > 0 and d > 0 then

So the solution for

must be between

and

Thus

Bakhshali approximation

The Bakhshali square root approximation is a mathematical method for finding an approximation to a square root which was described in an ancient manuscript by the name of "The Bakhshali Manuscript" because it was discovered in 1881 near the village of Bakhshali (or Bakhshalai) in the Yusufzai subdivision of the Peshawar district (now part of Pakistan). The manuscript was made out of the leaves of birch bark tree and was written in Sarada script.

The Bakhshali square root approximation is just the simple approximation applied twice.

Let

Let A = N + PA=N+P

When expanded, the equation becomes

Bakhshali approximation example

Find

Using

An exact "long-division like" algorithm

This method, while much slower than the Babylonian method, has the advantage that it is exact: if the given number has a square root whose decimal representation terminates, then the algorithm terminates and produces the correct square root after finitely many steps. It can thus be used to check whether a given integer is a square number.

Write the number in decimal and divide it into pairs of digits starting from the decimal point. The numbers are laid out similar to the long division algorithm and the final square root will appear above the original number.

For each iteration:

Bring down the most significant pair of digits not yet used and append them to any remainder. This is the current value referred to in steps 2 and 3.
If ss denotes the part of the result found so far, determine the greatest digit xx that does not make exceed the current value. Place the new digit xx on the quotient line.
Subtract yy from the current value to form a new remainder.
If the remainder is zero and there are no more digits to bring down the algorithm has terminated. Otherwise continue with step 1.

Example: What is the square root of 152.2756?

1 2. 3 4 | 01 52.27 56 1 The digit 1 is a solution digit x 01 1(20*0+1)=1 +1 00 52 22 The digit 2 is a solution digit 2x 00 44 2(20*1+2)=44 + 2 08 27 243 The digit 3 is a solution digit 24x 07 29 3(20*12+3)=729 + 3 98 56 2464 The digit 4 is a solution digit 246x 98 56 4(20*123+4)=9856 4 00 00 Algorithm terminates: answer is 12.34

Although demonstrated here for base 10 numbers, the procedure works for any base, including base 2. In the description above, 20 means double the number base used, in the case of binary this would really be 100. The algorithm is in fact much easier to perform in base 2, as in every step only the two digits 0 and 1 have to be tested. See Shifting nth-root algorithm.

Rough estimation

Many of the methods used for finding the square root of a number requires an initial seed value which should ideally be close to the actual value of the square root. One way of obtaining a very rough estimate of the value of the square root is as follows.

Provided that r » 1

Steps

Take the integer part of the number r. Z = int(r)
Count the number of digits in Z. Let D be the number of digits.
Calculate the value of 3^D.
The rough estimate is half the value obtained in step 3. E = (3^D) / 2

Example

r	Z	D	3 ^ D	Estimate (E)	Actual Value of √r
723.47	723	3	27	13.5	26.89
5396.37	5396	4	81	40.5	73.45
24956.41	24956	5	243	121.5	157.96
789345.464	789345	6	729	364.5	888.45

More accurate rough estimation

Similar to the rough estimation procedure above. This version is slightly more accurate because it also uses the additional information provided by the first digit of the number r. The steps are almost the same as the one for previous rough estimate. This estimate can be taken one extra step further by applying the Bakhshali square root approximation.

Provided that r » 1

Steps

Take the integer part of the number r. Z = int(r)
Count the number of digits in Z. Let D be the number of digits.
Calculate the value of 3.16^D. Note: 3.16 ≈ √10
The estimate is the value obtain in step 3 multiple by an adjustment factor based on the first (left most) digit of the number n. 'E\' = ADJ * 3.16^D
Applying the Bakhshali square root approximation to the estimate E.

A = E + PA=E+P

Table of ADJ

First Digit	1	2	3	4	5	6	7	8	9
ADJ	0.32	0.45	0.55	0.63	0.71	0.78	0.84	0.90	0.95

Example

r	Z	D	3.16 ^ D	ADJ	Estimate (E)	Bakhshali (B)	Actual Value of √r
723.47	723	3	31.55	0.84	26.50	26.89	26.89
5396.37	5396	4	99.71	0.71	70.79	73.46	73.45
24956.41	24956	5	315.09	0.45	141.79	157.98	157.96
789345.464	789345	6	995.7	0.84	836.37	888.45	888.45

Inverse square root approximation

There is an algorithmn for finding the inverse square root

very quickly.

Because , if we know the value of the inverse square root, we can easily get the value for the square root.

Steps for Inverse square root approximation

Step 1 Start the initial value

an approximate value of the inverse square root which is accurate to 2 significant digits.

Step 2 Calculate the next iteration with the following algorithmn.

where

Step 3 Finally obtain the value of

with

Because

is the inverse square root of r.

Properties

The greatest weakness of this algorithmn is that it is not stable. The initial value must be near the actual value for the algorithmn to converge correctly. If the initial value isn't closer enough, the algorithmn will diverge away from the actual value.
However when the algorithmn converges, it converges very very quickly. You can get an result with 20 digits of accuracy with just 6 iterations.

Square roots using Newton iteration

Basic Newton iteration finds a single root of a function given a sufficiently precise approximation to the root. The nature of which root will be given based on an approximation is dependent on the Newton fractal which we will not discuss here any further. The basic iteration is given by:

There are two widely used functions

and

used to find the square root of a number denoted by "z".

First method

The first method finds the square root of "z"

Note that both

and

are roots of the function

. ie

The first derivative of is

Thus iteration for is derived where:

and




	.

Second method

The second method finds the reciprocal of the square root of "z".

The two roots to

are

and

The derivative of is .

Thus iteration for is derived where:

and






	.

Example

Find the

using both methods.

Because we are looking for the square root of 7

Comparison

The iteration for

involves a division which is more time consuming than a multiplication in computer integer arithmetic. The iteration for

involves no division and is thus recommended for large integers z.

This iteration using "g" involves only a squaring and two multiplications, as opposed to a division in the case of f. In practical implementations of large integer square roots, the iteration involving g is faster for large integers z since division is at best , a constant times the time function of multiplication. The constant term is almost always 3 or more, meaning that a single division can almost never be faster than 3 multiplications.

Pell's equation

Pell's equation yields a method for finding rational approximations of square roots of integers.

Finding square roots using mental arithmetic

Based on Pell's equation there is a method to calculate square roots simply by subtracting odd numbers.

Ex: For we start with the following sequence:

27 - 1 = 26
26 - 3 = 23
23 - 5 = 18
18 - 7 = 11
11 - 9 = 2

Five steps have been taken and thus the integer part of the square root of 27 is 5. We do not proceed any further because the sixth subtraction would give a negative answer.

Now take the rest (2) and multiply it by 100 to get the starting number for the next step. Take the answer we have already got (5) and multiply it by 20, then add 1 to get the first odd number we subtract by. This gives us 2 × 100 = 200 as the starting number and 5 × 20 + 1 = 101 as the first odd number. Subtract the successive odd numbers, 103, 105, etc. until we get to a step where the next subtraction would result in a negative number. So,

200 - 101 = 99

and 99 is less than 103 so this is the place to stop and since we only took one step we get that the next digit is 1.

For the next step the starting number will be 99 × 100 = 9900 and the answer we have already got is 51 so the first odd number will be 51 × 20 + 1 = 1021

9900 - 1021 = 8879
8879 - 1023 = 7856
7856 - 1025 = 6831
6831 - 1027 = 5804
5804 - 1029 = 4775
4775 - 1031 = 3744
3744 - 1033 = 2711
2711 - 1035 = 1676
1676 - 1037 = 639

We took nine steps so the next digit is 9.

Continuing with this procedure, the next starting number will be 639 × 100 = 63900 and the first odd number will be 519 × 20 + 1 = 10381

63900 - 10381 = 53519
53519 - 10383 = 43136
43136 - 10385 = 32751
32751 - 10387 = 22364
22364 - 10389 = 11975
11975 - 10391 = 1584

We took six steps so the next digit is 6.

The result gives us 5.196 as an approximation of the square root of 27.

Finding square roots using basic arithmetic

If you have access to a calculator capable of addition, subtraction, multiplication and division then the following method can be used to find the square root of a number.

Find

The starting values are

Now find new values for each variable

We found a digit when the value of “

” is equal to the value of “”

After a new digit has been found.

New

The step are shown in the table below

Find
and
nn			Remarks	SOLN
n=0	99	55
n=1	55	55	STOP	55
and
n=0	99	22
n=1	22	11
n=2	11	11	STOP	5151
and
n=0	99	99	STOP	519519
and
n=0	99	66
n=1	66	66	STOP	51965196
and
n=0	99	11
n=1	11	11	STOP	5196151961

Continued fraction methods

Quadratic irrationals, that is numbers involving square roots in the form (a + √b)/c, have periodic continued fractions. This makes them easy to calculate recursively given the period. For example, to calculate √2, we make use of the fact that √2 − 1 = [0; 2, 2, 2, 2, 2, ...], and use the recurrence relation

a_{n + 1} = 1/(2 + a_n) with a₀ = 0

to obtain √2 − 1 to some specific precision specified through n levels of recurrence, and add 1 to the result to obtain √2.

Steps for finding continued fractions

Find using continued fractions

Let

thus

Step 0. We shall assume that r is not a perfect square. In other words:

Step 1. Find

using some other method. The best method is

using some other algorithmn to determine the integer square root.

Step 2. Find the highest lower bound (L) and lowest upper bound (U) for

where both (L) and (U) are integers.

Hence

Step 3. Write

in terms of

and

Step 4. Substitute

with

after substitution

Since

is less than one, we can determine

from the numeric values of

and

because

(ie.

is an integer). Determine the numeric value of

Step 5. Once we know the value of , we can rework the equation for

and

Step 6. Write

in terms of

and

Step 7. Substitute

with

after substitution

Since

is less than one, we can determine

from the numeric values of

and

because

(ie.

is an integer). Determine the numeric value of

Step 8. Once we know the value of , we can rework the equation for

and

Step 9. Repeat step 6 , 7 and 8 .

Quadratic equation method

The solution to a properly setup quadratic equation can be used to find where 1 < r < 100

For square root of values greater than 100, use the following identity:

For square root of values less than 1, use the following indentity:

Using the equation where and

Solve for XX using quadratic equation formula, choose the solution that satisfy the restriction

Final solution is :

The obvious problem is that we cannot evaluate the solutions to the quadratic equation without the usage of the square root function. However we can make the snake bite its own tail.

Let

Thus

And

So the quadratic equation becomes :

Solving for X as far as possible brings

Reciprocal

So the final solution becomes

Lets go deeper by substituting √r on the right hand side with its own definition.

Renormalise

Deeper still

And so on and so forth

Example using the quadratic equation method

Find

Using identity

First find

Hence

because

Using the quadratic equation formula, we get the two solutions. soln1 = - 0.165 or soln2 = 25.7519

Choose soln2 as it satisfy the restriction .

Hence

Alternatively

And the final solution is

Square roots of complex numbers

To every non-zero complex number z there exist precisely two numbers w such that w² = z. The usual definition of √z is as follows: if z = r exp(iφ) is represented in polar coordinates with -π < φ ≤ π, then we set √z = √r exp(iφ/2). Thus defined, the square root function is holomorphic everywhere except on the non-positive real numbers (where it isn't even continuous). The above Taylor series for √(1+x) remains valid for complex numbers x with |x| < 1.

When the number is in rectangular form the following formula can be used:

where the sign of the imaginary part of the root is the same as the sign of the imaginary part of the original number.

Note that because of the discontinuous nature of the square root function in the complex plane, the law √(zw) = √(z)√(w) is in general not true. Wrongly assuming this law underlies several faulty "proofs", for instance the following one showing that -1 = 1:

The third equality cannot be justified. (See invalid proof.)

However the law can only be wrong up to a factor -1, √(zw) = ±√(z)√(w), is true for either ± as + or as - (but not both at the same time). Note that √(c²) = ±c, therefore √(a²b²) = ±ab and therefore √(zw) = ±√(z)√(w), using a = √(z) and b = √(w).

Square roots of matrices and operators

If A is a positive definite matrix or operator, then there exists precisely one positive definite matrix or operator B with B² = A; we then define √A = B.

More generally, to every normal matrix or operator A there exist normal operators B such that B² = A. In general, there are several such operators B for every A and the square root function cannot be defined for normal operators in a satisfactory manner. Positive definite operators are akin to positive real numbers, and normal operators are akin to complex numbers.

Infinitely nested square roots

Under certain condition infinitely nested radicals such as

represent rational numbers. This rational number can be found by realizing that x also appears under the radical sign, which gives the equation

If we solve this equation, we find that x = 2. More generally, we find that

The same procedure also works to get

This method will give a rational xx value for all values of nn such that

Square roots of the first 20 positive integers

√ 1 = 1
√ 2 ≈ 1.4142135623 7309504880 1688724209 6980785696 7187537694 8073176679 7379907324 78462
√ 3 ≈ 1.7320508075 6887729352 7446341505 8723669428 0525381038 0628055806 9794519330 16909
√ 4 = 2
√ 5 ≈ 2.2360679774 9978969640 9173668731 2762354406 1835961152 5724270897 2454105209 25638
√ 6 ≈ 2.4494897427 8317809819 7284074705 8913919659 4748065667 0128432692 5672509603 77457
√ 7 ≈ 2.6457513110 6459059050 1615753639 2604257102 5918308245 0180368334 4592010688 23230
√ 8 ≈ 2.8284271247 4619009760 3377448419 3961571393 4375075389 6146353359 4759814649 56924
√ 9 = 3
√10 ≈ 3.1622776601 6837933199 8893544432 7185337195 5513932521 6826857504 8527925944 38639
√11 ≈ 3.3166247903 5539984911 4932736670 6866839270 8854558935 3597058682 1461164846 42609
√12 ≈ 3.4641016151 3775458705 4892683011 7447338856 1050762076 1256111613 9589038660 33818
√13 ≈ 3.6055512754 6398929311 9221267470 4959462512 9657384524 6212710453 0562271669 48293
√14 ≈ 3.7416573867 7394138558 3748732316 5493017560 1980777872 6946303745 4673200351 56307
√15 ≈ 3.8729833462 0741688517 9265399782 3996108329 2170529159 0826587573 7661134830 91937
√16 = 4
√17 ≈ 4.1231056256 1766054982 1409855974 0770251471 9922537362 0434398633 5730949543 46338
√18 ≈ 4.2426406871 1928514640 5066172629 0942357090 1562613084 4219530039 2139721974 35386
√19 ≈ 4.3588989435 4067355223 6981983859 6156591370 0392523244 4936890344 1381595573 28203
√20 ≈ 4.4721359549 9957939281 8347337462 5524708812 3671922305 1448541794 4908210418 51276

External links

Japanese soroban techniques - Professor Fukutaro Kato's method
Japanese soroban techniques - Takashi Kojima's method