Vai al contenuto
Hub Nexus
Aggiornato

AutoreAncora senza autorePrendila in carico

Vedi qualcosa da migliorare? Proponi una modifica.

Sostegno

Questa pagina è pubblica.

Theory

Let x be the number of which we want to obtain the square root y = \sqrt{x}; Let's consider its decimal expansion, for example: x = 456.7890123. Let's separate its digits into groups of two around the decimal point in the following way

x = 456.7890123 = 4 ⋅ 56 ⋅ 78 ⋅ 90 ⋅ 12 ⋅ 30 ⋅ 00 ⋅ 00 ⋯

or, in other words, let's define the sequence of integers {α}_{i}

{α}_{1} = 4 , {α}_{2} = 56 , {α}_{3} = 78 , {α}_{4} = 90 , {α}_{5} = 12 , {α}_{6} = 30 , {α}_{7} = 00 , ⋯

and let's build the sequence {x}_{i} recursively from {x}_{0} = 0

{x}_{i} = 100 {x}_{i − 1} + {α}_{i}

and let {y}_{i} be the integer part of the square root of {x}_{i}

{y}_{i} = ⌊ {\sqrt{x}}_{i} ⌋

i.e.  {y}_{i} is the largest integer whose square is not greater than {x}_{i}. Finally, let us call remainders to the differences

{r}_{i} = {x}_{i} − {y}_{i}^{2} ≥ 0

For our example we have:

i{α}_{i}{x}_{i}{y}_{i}{r}_{i}
0000
14420
2564562115
37845678213309
490456789021371121
5124567890122137226628
⋯

Let's see that, by construction, {x}_{i} grows as {}1 {0}^{2 i} (two more digits in each step), in fact the sequence {x}_{i} ⋅ 1 {0}^{4 − 2 i}: (0, 400, 456, 456.78, 456.7890, etc.) tends to x ( {x}_{i} ⋅ 1 {0}^{4 − 2 i} → x ) or x = ( {lim}_{i → ∞} {x}_{i} ⋅ 1 {0}^{4 − 2 i} ). By comparison, {y}_{i} , as the integer part of the square root of {x}_{i}, grows only as {}1 {0}^{i} (one digit more each step). As {y}_{i} is the largest integer whose square is not greater than {x}_{i} we have {r}_{i} = {x}_{i} − {y}_{i}^{2} ≥ 0 as above but

{x}_{i} − ( {y}_{i} + 1 {)}^{2} = {x}_{i} − {y}_{i}^{2} − 2 {y}_{i} − 1 < 0

by definition of {y}_{i}, or

{r}_{i} = {x}_{i} − {y}_{i}^{2} < 2 {y}_{1} + 1

multiplying by {}1 {0}^{4 − 2 i}

{x}_{i} ⋅ 1 {0}^{4 − 2 i} − ( {y}_{i} ⋅ 1 {0}^{2 − i} {)}^{2} < ( 2 {y}_{i} + 1 ) ⋅ 1 {0}^{4 − 2 i}

but as  {y}_{i} grows only as {}1 {0}^{i}, the second term tends to zero as {}1 {0}^{− i}.  With this

{x}_{i} ⋅ 1 {0}^{4 − 2 i} − ( {y}_{i} ⋅ 1 {0}^{2 − i} {)}^{2} → 0

and {x}_{i} ⋅ 1 {0}^{4 − 2 i} → x so that we have

{y}_{i} ⋅ 1 {0}^{2 − i} → y = \sqrt{x}

For other numbers, the above factors are: {}1 {0}^{2 k − 2 i} and {}1 {0}^{k − i}, where k is the number of two-digit groups to the left of the decimal point, negative if it is followed by 00 groups (ex. k = − 1 for x = 0.00456, k = − 2 for x = 0.00000456, etc.).

This is the basis for traditional manual square root methods.

Procedure

We start with i = 0, {x}_{0} = 0, {y}_{0} = 0, {r}_{0} = 0.

First digit

b{b}^{2}
11
24
39
416
525
636
749
864
981

For i = 1 and {x}_{1} = {α}_{1}, it is trivial to find {y}_{1} such that its square does not exceed {x}_{1} through the use of the following table of squares that is easily retained in memory as it is only a subset of the multiplication table. In the case of the example we find {y}_{1} = 2.

Rest of digits

For i > 1, we have {x}_{i} = 100 {x}_{i − 1} + {α}_{i} as defined above and we try to build {y}_{i} in the way:

{y}_{i} = 10 {y}_{i − 1} + b

where b is a one-digit integer ranging from 0 to 9. To obtain b we have to choose the greatest digit from 0 to 9 such that:

{x}_{i} − {y}_{i}^{2} = {x}_{i} − ( 10 {y}_{i − 1} + b {)}^{2} ≤ 0

or

{x}_{i} − ( a + b {)}^{2} ≤ 0

if we write a = 10 {y}_{i − 1}. Expanding the binomial we have

{x}_{i} − {a}^{2} − 2 a b − {b}^{2} = 100 {x}_{i − 1} + {α}_{i} − ( 10 {y}_{i − 1} {)}^{2} − 2 a b − {b}^{2} ≤ 0

or

100 {r}_{i − 1} + {α}_{i} ≥ 2 a b + {b}^{2} = ( 2 a + b ) b = ( 20 {y}_{i − 1} + b ) b

The left side of the above expression may be seen as simply the previous remainder with the next two-digits group appended to it, and the parenthesis of the last term as twice the previous root with digit b appended to it. In our example, for i = 2 we have 56 on the left and the above expression is

56 ≥ ( 40 + b ) b

which holds only for b = 0 or b = 1 so that 1 is our next root digit but, how can we proceed in the general case without having to systematically explore every possibility (b = 0 , 1 , ⋯ , 9)?

Here Knott[1] distinguishes two different approaches:

  • Preparing the divisor
  • Preparing the dividend.

Preparing the divisor

This corresponds with the above expression

100 {r}_{i − 1} + {α}_{i} ≥ ( 20 {y}_{i − 1} + b ) b

And it is the strategy usually used with paper and pencil and can also be implemented, of course, on the abacus. In the above expression, if we see the left part as the dividend and the parenthesis as the divisor, b is the first digit of the division

b = ( 100 {r}_{i − 1} + {α}_{i} ) / ( 20 {y}_{i − 1} + b )

but since we don't know b yet, we approximate it using only the main part of the divisor

b = ( 100 {r}_{i − 1} + {α}_{i} ) / ( 20 {y}_{i − 1} )

This gives us a guess as to what the value of b might be, but we need:

  • Verify that the value thus obtained is correct, or, where appropriate, correct it up or down as needed.
  • Obtain the next remainder to prepare the calculation of the next digit of the root.

Both steps require subtracting ( 20 {y}_{i − 1} + b ) b or {}20 {y}_{i − 1} b and {b}^{2} from {}100 {r}_{i − 1} + {α}_{i}; checking that we are not going to negative values and that what remains is less than {}2 {y}_{i} + 1 (otherwise we would have to revise b up or down). If we do this correctly, what we are left with is the new remainder {r}_{i}. It should be noted that, as we proceed in the calculations (i increasing) b is progressively a smaller and smaller contribution to the divisor ( 20 {y}_{i − 1} + b ); so the process indicated above will look more and more like a mere division.

This is the method proposed by Takashi Kojima in his second book: Advanced Abacus - Theory and Practice2], and that you can see described in Square roots as solved by Kojima[3] in [Totton heffelfinger’s website, works to which I refer the reader to see explanations and examples. What follows here, for purposes of illustration, is a sketch of how the calculation might be started in our example: x = 456.7890123

Napier's bones arranged to help with the third root digit of the example (Image: Jccsvq, CC BY-SA 4.0)

AbacusComment
ABCDEFGHIJKLM
4567890123Entering radicand starting in CD (first group)
2First root digit in B
-4Subtract square of B from first group
2  567890123Null remainder
4   56 7890123Doubling B. Appending next group (56)to remainder
41 5678901235/4≈1, try 1 as next root digit
-4Continue division by 41, subtract 1✕41 from EF
-1
41 15789012315 as remainder
42 157890123Double second root digit
42 15 78 90123Append next group (78)
423157890123157/42≈3, try 3 as next root digit
-12Continue division by 423, subtract 3✕423 from E-H
-06
-09
423 30990123309 as remainder
426 30990123Double third root digit
426 309 90 123Append next group (90)
etc.

As can be seen, twice the root grows to the left of the abacus to the detriment of the radicand / remainder and the groups of two figures still unused. This is contrary to what happens with the rest of the elementary operations on the abacus, where the result sought replaces the operand (or one of them). For this reason, the traditionally preferred method for obtaining square roots seems to have been the following, where we will see the root appear directly on the abacus, not its double .

Preparing the dividend

Starting again with

100 {r}_{i − 1} + {α}_{i} ≥ ( 20 {y}_{i − 1} + b ) b = 20 {y}_{i − 1} b + {b}^{2}

dividing by 2

( 100 {r}_{i − 1} + {α}_{i} ) / 2 ≥ 10 {y}_{i − 1} b + {b}^{2} / 2

This modified expression will allow us to directly obtain the square root in the abacus following practically the same procedure as above with only keeping half-remainders on our instrument. Here, with square roots, the change is almost trivial, but it will be more important when dealing with cube roots. As can be seen in the above expression, neglecting the term {b}^{2} / 2 we obtain a guess of b by simply dividing the extended half-remainder ( 100 {r}_{i − 1} + {α}_{i} ) / 2 by the previous root {y}_{i − 1} (in fact {}10 {y}_{i − 1}).  After that, we need again:

  • Verify that the value thus obtained is correct, or, where appropriate, correct it up or down as needed.
  • Obtain the next half-remainder to prepare for the next digit of the root.

This is done by subtracting {}10 {y}_{i − 1} b and {b}^{2} / 2 from the half remainder, and this makes it convenient to memorize the following table of semi-squares:

b{b}^{2} / 2
10.5
22
34.5
48
512.5
618
724.5
832
940.5

Fortunately, since 2 is a divisor of our base (10), the decimal fractions in the table have a finite expression; which will not happen when we try to extend this procedure to cube roots and we have to deal with thirds of cubes. According to Knott, this makes cube roots a problem that is not well suited to being treated with abacus.

Examples

Here three examples are presented, for additional examples please see Further reading and specially External resources below.

Square root of 961

In this example we have two groups of two figures: 09 and 61. The first group informs us that the first digit of the root is 3.

There are two ways to start square roots:

  • Aligning the groups to the left of the abacus from column B and using the traditional division to obtain the semi-remainder.
ABCDE

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) | ... |
| | | 0 | 9 | 6 | 1 | |

This is the form used in the old books and also the one used in Murakami's Square root tutor with Kijoho (see below External resources).

AbacusComment
ABCDE
0961Align the radicand with B
30961Enter first root digit in A
-9Subtract the square of first root digit (9)
30061
30305Divide the remainder B-E by 2 (帰除法)
  • Aligning the groups to the left of the abacus from column A and using in-situ division, as explained in chapter: Division by powers of two, to obtain the semi-remainder.
ABCDE

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) |

(Image: Jccsvq, CC BY-SA 4.0) | ... |
| | 0 | 9 | 6 | 1 | | |

This method is somewhat faster.

AbacusComment
ABCDE
0961Align the radicand with A
-9Subtract the square of first root digit (9)
0061
0305Divide in situ the remainder by 2
30305Enter first root digit in A

From here, the state of the abacus coincides and we can continue:

AbacusComment
ABCDE
30305
+1Divide half-remainder B-E by 3. Revising up B
-3
31005
-05Subtract b^2/2 =0.5 from D
31000Half-remainder is 0, Done! Root is 31
31Root is 31

Square root of 998001

AbacusComment
ABCDEFG
998001Enter the radicand
-81Subtract 9^2=81 from first group
188001
940005Halve the remainder in situ
9940005Enter first root digit into A
9 9 30005B: Rule 9/9>9+9
-405Subtract 9^2/2=40.5 from D
9989505
99 8 7505C: Rule 8/9>8+8
-72Subtract CxB=72 from DE
998T305Revise up C
+1
-99
9990405
-405Subtract 9^2/2=40.5 from F
9990000Remainder is 0. Done!
999Root is 999

Root of 456.7890123

Our example above...

AbacusComment
ABCDEFGHIJKL
04567890123Enter x aligning digit pairs from AB, CD, etc.
-4Subtract 2^2 from first group
567890123
2839450615Halve remainder and rest of digits pairs
2 2839450615Enter first root digit in A
+1Divide BCD by A (revise up B)
-2
-05Subtract B^2/2=0.5 from D
21 789450615
+3Divide CDEF by AB (revise up C three times)
-6
-3
-45Subtract C^2/2=4.5 from F
213154950615
213554950615Divide DEFGH by ABC. D: Rule 1/2>5+0
-5Subtract DxB=5 from EF
-15Subtract DxC=15 from FG
213548450615
+2Revise up D twice
-426
213705850615
-245Subtract 7^2/2=24.5 From H
21370560F615Root so far: 21.37
etc.etc.

The root 2137… (21.37…) is appearing on the left. See chapter: Abbreviated operations to see how to quickly approximate the next 4 digits.

Conclusion

The method explained as: Preparing the dividend is known as 半九九法 (Hankukuhou in Japanese, Bàn jiǔjiǔ fǎ in Chinese) that we freely translate here as the Half-remainder method and is by far the most convenient to use, at least for two reasons:

  • The root, and not its double, replaces the operand (radicand) as in the rest of operations on the abacus.
  • (The most important) Dividing by numbers that start with 1 is awkward. Think of the first group of two digits, its value is between 1 and 99 and determines the first figure of the square root. For values ​​of the first pair between 25 and 99 (75% of the cases) the first digit of the root is between 5 and 9 and its double begins with one! Therefore, if we use the method preparing the divisor, we will be dividing by numbers that begin with 1 in 75% of the cases. On the contrary, if we use the method preparing the dividend, only in the case that the first group is 1, 2 or 3 (3% of the cases) will we have to divide by numbers that start with one.

The superiority of the half-remainder or preparing the dividend method is undeniable.

References

Further reading

External resources

  • Murakami's Square root tutor with Kijoho, a JavaScript application that you can run directly in your browser or download to your computer from its GitHub repository. You just have to enter the root in the small input box on the left and repeatedly press the "next" button on the screen to see the development of the process step by step. So you can generate as many examples or exercises as you want.

Where this page came from

This page was imported from Wikibooks. From “Traditional Abacus and Bead Arithmetic” on Wikibooks, by its contributors, under CC BY-SA 4.0. Changed here: set as a page; navigation and edit links left out; each image under its own licence, credited in its caption.

Nobody has written it yet — it is the source material at a new address, which is why search engines are asked to skip it and why no one earns from it. It is up for grabs: take it on, and it is yours to rewrite and to earn from.

LingueEnglish

Licenza: CC BY-SA 4.0 · Tratto da en.wikibooks.org

Visualizzazioni

Commenti

Spinner Logo
Versione: 2CC0 1.0 (pubblico dominio)
The runaway star that left the Tarantula Nebula
Versione: 2CC0 1.0 (pubblico dominio)
The Eagle Nebula, seen in the infrared
Versione: 2CC0 1.0 (pubblico dominio)
The Cosmic Cliffs are not cliffs
Versione: 2CC0 1.0 (pubblico dominio)
The recommended targets are real, but the largest health gain is between doing nothing and doing a little. That is the part worth knowing if you are starting from nothing.
Versione: 2CC0 1.0 (pubblico dominio)
Sleep is not one state but a cycle through several, each doing different work. Most sleep disorders are a failure to complete that cycle rather than a shortage of hours.
Versione: 2CC0 1.0 (pubblico dominio)
Four nebulae that revealed the shape of our galaxy