Rod calculus


Rod calculus or rod calculation was the mechanical method of algorithmic computation with counting rods in China from the Warring States to Ming dynasty before the counting rods were replaced by the more convenient and faster abacus. Rod calculus played a key role in the development of Chinese mathematics to its height in Song Dynasty and Yuan Dynasty, culminating in the invention
of polynomial equations of up to four unknowns in the work of Zhu Shijie.

Hardware

The basic equipment for carrying out rod calculus is a bundle of counting rods and a counting board. The counting rods are usually made of bamboo sticks, about 12 cm- 15 cm in length, 2mm to 4 mm diameter, sometimes from animal bones, or ivory and jade. A counting board could be a table top, a wooden board with or without grid, on the floor or on sand.
In 1971 Chinese archaeologists unearthed a bundle of well-preserved animal bone counting rods stored in a silk pouch from a tomb in Qian Yang county in Shanxi province, dated back to the first half of Han dynasty. In 1975 a bundle of bamboo counting rods was unearthed.
The use of counting rods for rod calculus flourished in the Warring States, although no archaeological artefacts were found earlier than the Western Han Dynasty ; since the rod calculus software must have gone along with rod calculus hardware, there is no doubt that rod calculus was already flourishing during the Warring States more than 2,200 years ago.

Software

The key software required for rod calculus was a simple 45 phrase positional decimal multiplication table used in China since antiquity, called the nine-nine table, which were learned by heart by pupils, merchants, government officials and mathematicians alike.

Rod numerals

Displaying numbers

Rod numerals is the only numeric system that uses different placement combination of a single symbol to convey any number or fraction in the Decimal System. For numbers in the units place, every vertical rod represent 1. Two vertical rods represent 2, and so on, until 5 vertical rods, which represents 5. For number between 6 and 9, a biquinary system is used, in which a horizontal bar on top of the vertical bars represent 5. The first row are the number 1 to 9 in rod numerals, and the second row is the same numbers in horizontal form.
For numbers larger than 9, a decimal system is used. Rods placed one place to the left of the units place represent 10 times that number. For the hundreds place, another set of rods is placed to the left which represents 100 times of that number, and so on. As shown in the adjacent image, the number 231 is represented in rod numerals in the top row, with one rod in the units place representing 1, three rods in the tens place representing 30, and two rods in the hundreds place representing 200, with a sum of 231.
When doing calculation, usually there was no grid on the surface. If rod numerals two, three, and one is placed consecutively in the vertical form, there's a possibility of it being mistaken for 51 or 24, as shown in the second and third row of the adjacent image. To avoid confusion, number in consecutive places are placed in alternating vertical and horizontal form, with the units place in vertical form, as shown in the bottom row on the right.

Displaying zeroes

In Rod numerals, zeroes are represented by a space, which serves both as a number and a place holder value. Unlike in Arabic numerals, there is no specific symbol to represent zero. In the adjacent image, the number zero is merely represented with a space.

Negative and positive numbers

mathematicians used red to represent positive numbers and black for negative numbers. However, another way is to add a slash to the last place to show that the number is negative.

Decimal fraction

The Mathematical Treatise of Sunzi used decimal fraction metrology. The unit of length was 1 chi,
1 chi = 10 cun, 1 cun = 10 fen, 1 fen = 10 li, 1 li = 10 hao, 10 hao = 1 shi, 1 shi = 10 hu.
1 chi 2 cun 3 fen 4 li 5 hao 6 shi 7 hu is laid out on counting board as
where is the unit measurement chi.
Southern Song dynasty mathematician Qin Jiushao extended the use of decimal fraction beyond metrology. In his book Mathematical Treatise in Nine Sections, he formally expressed 1.1446154 day as
He marked the unit with a word “日” underneath it.

Addition

Rod calculus works on the principle of addition. Unlike Arabic numerals, digits represented by counting rods have additive properties. The process of addition involves mechanically moving the rods without the need of memorising an addition table. This is the biggest difference with Arabic numerals, as one cannot mechanically put 1 and 2 together to form 3, or 2 and 3 together to form 5.
The adjacent image presents the steps in adding 3748 to 289:
  1. Place the augend 3748 in the first row, and the addend 289 in the second.
  2. Calculate from LEFT to RIGHT, from the 2 of 289 first.
  3. Take away two rod from the bottom add to 7 on top to make 9.
  4. Move 2 rods from top to bottom 8, carry one to forward to 9, which becomes zero and carries to 3 to make 4, remove 8 from bottom row.
  5. Move one rod from 8 on top row to 9 on bottom to form a carry one to next rank and add one rod to 2 rods on top row to make 3 rods, top row left 7.
  6. Result 3748+289=4037
The rods in the augend change throughout the addition, while the rods in the addend at the bottom "disappear".

Subtraction

Without borrowing

In situation in which no borrowing is needed, one only needs to take the number of rods in the subtrahend from the minuend. The result of the calculation is the difference. The adjacent image shows the steps in subtracting 23 from 54.

Borrowing

In situations in which borrowing is needed such as 4231–789, one need use a more complicated procedure. The steps for this example are shown on the left.
  1. Place the minuend 4231 on top, the subtrahend 789 on the bottom. Calculate from the left to the right.
  2. Borrow 1 from the thousands place for a ten in the hundreds place, minus 7 from the row below, the difference 3 is added to the 2 on top to form 5. The 7 on the bottom is subtracted, shown by the space.
  3. Borrow 1 from the hundreds place, which leaves 4. The 10 in the tens place minus the 8 below results in 2, which is added to the 3 above to form 5. The top row now is 3451, the bottom 9.
  4. Borrow 1 from the 5 in the tens place on top, which leaves 4. The 1 borrowed from the tens is 10 in the units place, subtracting 9 which results in 1, which are added to the top to form 2. With all rods in the bottom row subtracted, the 3442 in the top row is then, the result of the calculation

    Multiplication

Sunzi Suanjing described in detail the algorithm of multiplication. On the left are the steps to calculate 38×76:
  1. Place the multiplicand on top, the multiplier on bottom. Line up the units place of the multiplier with the highest place of the multiplicand. Leave room in the middle for recording.
  2. Start calculating from the highest place of the multiplicand. Using the multiplication table 3 times 7 is 21. Place 21 in rods in the middle, with 1 aligned with the tens place of the multiplier. Then, 3 times 6 equals 18, place 18 as it is shown in the image. With the 3 in the multiplicand multiplied totally, take the rods off.
  3. Move the multiplier one place to the right. Change 7 to horizontal form, 6 to vertical.
  4. 8×7 = 56, place 56 in the second row in the middle, with the units place aligned with the digits multiplied in the multiplier. Take 7 out of the multiplier since it has been multiplied.
  5. 8×6 = 48, 4 added to the 6 of the last step makes 10, carry 1 over. Take off 8 of the units place in the multiplicand, and take off 6 in the units place of the multiplier.
  6. Sum the 2380 and 508 in the middle, which results in 2888: the product.

    Division

The animation on the left shows the steps for calculating.
  1. Place the dividend, 309, in the middle row and the divisor, 7, in the bottom row. Leave space for the top row.
  2. Move the divisor, 7, one place to the left, changing it to horizontal form.
  3. Using the Chinese multiplication table and division, 30÷7 equals 4 remainder 2. Place the quotient, 4, in the top row and the remainder, 2, in the middle row.
  4. Move the divisor one place to the right, changing it to vertical form. 29÷7 equals 4 remainder 1. Place the quotient, 4, on top, leaving the divisor in place. Place the remainder in the middle row in place of the dividend in this step. The result is the quotient is 44 with a remainder of 1
The Sunzi algorithm for division was transmitted in toto by al Khwarizmi to Islamic country from Indian sources in 825AD. Al Khwarizmi's book was translated into Latin in the 13th century, The Sunzi division algorithm later evolved into Galley division in Europe. The division algorithm in Abu'l-Hasan al-Uqlidisi's 925AD book Kitab al-Fusul fi al-Hisab al-Hindi and in 11th century Kushyar ibn Labban's Principles of Hindu Reckoning were identical to Sunzu's division algorithm.

Fractions

If there is a remainder in a place value decimal rod calculus division, both the remainder and the divisor must be left in place with one on top of another. In Liu Hui's notes to Jiuzhang suanshu, the number on top is called "shi", while the one at bottom is called "fa". In Sunzi Suanjing, the number on top is called "zi" or "fenzi", and the one on the bottom is called "mu" or "fenmu". Fenzi and Fenmu are also the modern Chinese name for numerator and denominator, respectively. As shown on the right, 1 is the numerator remainder, 7 is the denominator divisor, formed a fraction. The quotient of the division is 44 +.
Liu Hui's used a lot of calculations with fraction in Haidao Suanjing.
This form of fraction with numerator on top and denominator at bottom without a horizontal bar in between, was transmitted to Arabic country in an 825AD book by al Khwarizmi via India, and in use by 10th century Abu'l-Hasan al-Uqlidisi and 15th century Jamshīd al-Kāshī's work "Arithematic Key".

Addition

3 × 5
The algorithm for finding the highest common factor of two numbers and reduction of
fraction was laid out in Jiuzhang suanshu.
The highest common factor is found by successive division with remainders until
the last two remainders are identical.
The animation on the right illustrates the algorithm for finding the highest common factor of and reduction of a fraction.
In this case the hcf is 25.
Divide the numerator and denominator by 25. The reduced fraction is.

Interpolation

Calendarist and mathematician He Chengtian used fraction interpolation method, called "harmonisation of the divisor of the day" to obtain a better approximate value than the old one by iteratively adding the numerators and denominators a "weaker" fraction with a "stronger fraction". Zu Chongzhi's legendary could be obtained with He Chengtian's method

System of linear equations

Chapter Eight ''Rectangular Arrays of Jiuzhang suanshu provided an algorithm for solving System of linear equations by method of elimination:
Problem 8-1: Suppose we have 3 bundles of top quality cereals, 2 bundles of medium quality cereals, and a bundle of low quality cereal with accumulative weight of 39 dou. We also have 2, 3 and 1 bundles of respective cereals amounting to 34 dou; we also have 1,2 and 3 bundles of respective cereals, totaling 26 dou.
Find the quantity of top, medium, and poor quality cereals.
In algebra, this problem can be expressed in three system equations with three unknowns.
This problem was solved in Jiuzhang suanshu with counting rods laid out on a counting board in a tabular format similar to a 3x4 matrix:
Algorithm:
The amount of on bundle of low quality cereal =
From which the amount of one bundle of top and medium quality cereals can be found easily:
One bundle of top quality cereals=9 dou
One bundle of medium cereal=4 dou >

Extraction of Square root

Algorithm for extraction of square root was described in Jiuzhang suanshu and with minor difference in terminology in Sunzi Suanjing.
The animation shows the algorithm for rod calculus extraction of an approximation of the square root from the algorithm in chap 2 problem 19 of Sunzi Suanjing:
The algorithm is as follows:
North Song dynasty mathematician Jia Xian developed an additive multiplicative algorithm for square root extraction, in which he replaced the traditional "doubling" of "fang fa" by adding
shang digit to fang fa digit, with same effect.

Extraction of cubic root

vol iv "shaoguang" provided algorithm for extraction of cubic root.
〔一九〕今有積一百八十六萬八百六十七尺。問為立方幾何？答曰：一百二十三尺。
problem 19: We have a 1860867 cubic chi, what is the length of a side ? Answer:123 chi.
North Song dynasty mathematician Jia Xian invented a method similar to simplified form of Horner scheme for extraction of cubic root.
The animation at right shows Jia Xian's algorithm for solving problem 19 in Jiuzhang suanshu vol 4.

Polynomial equation

North Song dynasty mathematician Jia Xian invented Horner scheme for solving simple 4th order equation of the form
South Song dynasty mathematician Qin Jiushao improved Jia Xian's Horner method to solve polynomial equation up to 10th order.
The following is algorithm for solving
This equation was arranged bottom up with counting rods on counting board in tabular form
Algorithm:
  1. Arrange the coefficients in tabular form, constant at shi, coeffienct of x at shang lian, the coeffiecnt of X^4 at yi yu;align the numbers at unit rank.
  2. Advance shang lian two ranks
  3. Advance yi yu three ranks
  4. Estimate shang=20
  5. let xia lian =shang * yi yu
  6. let fu lian=shang *yi yu
  7. merge fu lian with shang lian
  8. let fang=shang * shang lian
  9. subtract shang*fang from shi
  10. add shang * yi yu to xia lian
  11. retract xia lian 3 ranks, retract yi yu 4 ranks
  12. The second digit of shang is 0
  13. merge shang lian into fang
  14. merge yi yu into xia lian
  15. Add yi yu to fu lian, subtract the result from fang, let the result be denominator
  16. find the highest common factor =25 and simplies the fraction
  17. solution

    Tian Yuan shu

Yuan dynasty mathematician Li Zhi developed rod calculus into Tian yuan shu
Example Li Zhi Ceyuan haijing vol II, problem 14 equation of one unknown:

Polynomial equations of four unknowns

Mathematician Zhu Shijie further developed rod calculus to include polynomial equations of 2 to four unknowns.
For example, polynomials of three unknowns:
Equation 1:
Equation 2:
Equation 3:
After successive elimination of two unknowns, the polynomial equations of three unknowns was reduced to
a polynomial equation of one unknown:
Solved x=5;