九章算术方程术

九章算术方程术

今天看到一篇关于九章算术方程术的文章，看得瞠目结舌，叹为观止！中国古人果真是太厉害了，现博文记录。
方程：今有上禾三秉,中禾二秉,下禾一秉,实三十九斗；上禾二秉,中禾三秉,下禾一秉,实三十四斗；上禾一秉,中禾二秉,下禾三秉,实二十六斗；问上、中、下禾实一秉各几何?

答曰：上禾一秉，九斗、四分斗之一，中禾一秉，四斗、四分斗之一，下禾一秉，二斗、四分斗之三。

术曰：置上禾三秉，中禾二秉，下禾一秉，實三十九斗，於右方。中、左禾列如右方。以右行上禾遍乘中行而以直除。又乘其次，亦以直除。然以中行中禾不盡者遍乘左行而以直除。左方下禾不盡者，上為法，下為實。實即下禾之實。求中禾，以法乘中行下實，而除下禾之實。餘如中禾秉數而一，即中禾之實。求上禾亦以法乘右行下實，而除下禾、中禾之實。餘如上禾秉數而一，即上禾之實。實皆如法，各得一斗。

翻译：今有上等稻3捆、中等稻2捆、下等稻1捆，共打出39斗米；有上等稻2捆、中等稻3捆、下等稻1捆，共打出34斗米；有上等稻1捆、中等稻2捆、下等稻3捆，共打出26斗米。问上等稻、中等稻、下等稻各1捆能打出多少斗米？现在我们用x,y,z分别代替上等稻、中等稻、下等稻各1捆能打出的斗米数。有

$$\left\{\begin{matrix} 3x+2y+z=39\\ 2x+3y+z=34\\ x+2y+3z=26 \end{matrix}\right.$$
九章算术解法：将方程组的系数和值从右往左排列：【置上禾三秉，中禾二秉，下禾一秉，實三十九斗，於右方。中、左禾列如右方。】得到
$$\begin{matrix} 1 & 2 & 3\\ 2 & 3 & 2\\ 3 & 1 & 1\\ 26 & 34 &39 \end{matrix}$$
【以右行上禾遍乘中行而以直除。又乘其次，亦以直除。】:中行各元素乘以右行的上禾，然后中行各元素减掉上行对应的元素，直到中行上禾为0；左行各元素乘以右行的上禾，然后左行各元素减掉右行对应的元素，直到左行上禾为0。

123262313432139→3697869310232