用集算器解多元一次方程组

给出一个多元一次方程组，对于n个未知数，方程组中需要包含n个方程，如：

3x+4y+5z=26

5x+6y+10z=47

4x+8y+7z=41

那么，如何用集算器来求解这个多元一次方程组呢？首先，需要将方程组整理为标准形式，并把参数填入一个文本文件中，同一行的参数之间用Tab分隔，如下：

如何用编程的方式求解多元一次方程组？常用的方法就是加减消元法，如下面的方程组：

a₁₁x₁+a₁₂x₂+a₁₃x₃+…+a_1nx_n=a₁₀

a₂₁x₁+a₂₂x₂+a₂₃x₃+…+a_2nx_n=a₂₀

a₃₁x₁+a₃₂x₂+a₃₃x₃+…+a_3nx_n=a₃₀

……

a_n1x₁+a_n2x₂+a_n3x₃+…+a_nnx_n=a_n0

当a₁₁不为0时，将第1个方程各项同时乘以a₂₁/a₁₁，并用第2个方程减去，就可以把第2个方程中x₁的系数消为0。如果a₁₁等于0，则需要找到1个x₁的系数不为0的方程，将其移到第一个。用同样的方法，可以消去第1个之后所有方程式中，x₁的系数，如下：

a₁₁x₁+a₁₂x₂+a₁₃x₃+…+a_1nx_n=a₁₀

0+b₂₂x₂+b₂₃x₃+…+b_2nx_n=b₂₀

0+b₃₂x₂+b₃₃x₃+…+b_3nx_n=b₃₀

……

0+b_n2x₂+b_n3x₃+…+b_nnx_n=b_n0

此时，继续再来考察b₂₂的情况，用类似的方法，继续消去第2个之后的所有方程式中，x₂的系数……这样执行下去，直到消去第n个方程式中x_n-1的系数为止。此时的方程组已经变为下面的形式：

a₁₁x₁+a₁₂x₂+a₁₃x₃+…+a_1nx_n=a₁₀

0+b₂₂x₂+b₂₃x₃+…+b_2nx_n=b₂₀

0+0+c₃₃x₃+…+c_3nx_n=c₃₀

……

0+0+0+…+m_nnx_n=m_n0

现在，就可以从第n个方程开始，一步一步求解了：

x_n=m_n0/m_nn

_……

x₃=(c₃₀-…-c_3nx_n)/c₃₃

x₂=(c₃₀-b₂₃x₃-…-b_2nx_n)/b₂₂

x₁=(c₃₀-a₁₂x₂-a₁₃x₃-…-a_1nx_n)/a₁₁

如果在计算中，a₁₁，b₂₂，c₃₃等参数中出现0，则说明方程组无解或有无穷解。

用集算器求解多元一次方程组的解法如下：

在B1中，从文本文件中读取数据，并在A2中将方程组参数转换为序列组：

在C1中求出方程组的方程数，在A3中循环，用加减消元法逐步消除x_n的系数。首先，在B4的循环中，从第k个方程式起，寻找x_k的系数不为0的第1个方程式，如果未找到，则说明该项无法消元。在第7行的代码中，将x_k的系数不为0的方程式移到第k行。在第8、9行，将k+1行起的方程式中，x_k的系数消去。

各步消元执行完毕后，在A10中，就可以从最后一个方程起，逐步求解。其中，用序列result来存储解。如果出现无解或无穷解的情况，返回Error。考虑到双精度数在计算中可能存在误差，求解时保留4位小数。

计算完成后，在A13中，可以查看到结果：

如果需要求解新的方程组：

7x₁+2x₂+9x₃-x₄=0

2x₁+9x₂-x₄=0

9x₁+11x₃-x₄=0

x₁+x₂+x₃ =1

只需要修改文件EquInput.txt：

用同样的程序就可以计算，计算后，可以在A13中读取结果：

另外，参数列表也可以不用文件读入，而是直接填写在集算器中，计算过程是类似的。