Maple编写高斯消元程序报错：非法使用形式参数

错误翻译 #

Error, (in Gauss) illegal use of a formal parameter
中文翻译：错误，(在Gauss过程中)非法使用形式参数

错误原因及修正步骤 #

• 非法修改形式参数：
  过程定义中n是形式参数（n::posint），但代码里执行了n := linalg[rowdim](A);，Maple不允许在过程内部给形式参数重新赋值。解决方法是改用新的局部变量存储矩阵行数，比如把局部变量声明改为local a, i, k, j, p, m;，然后用m := linalg[rowdim](A);，后续循环里的n替换成m，或者直接移除形式参数n，从矩阵A直接获取行数即可。

• 赋值符号使用错误：
  Maple中变量赋值必须用:=，而=是等式判断符号。代码里的p = a[i,k]/a[k,k];和a[i,j] = a[i,j]-p*a[k,j]都要改成:=才能完成赋值操作。

• 矩阵行交换未生效：
  swaprow函数不会直接修改原矩阵，需将返回结果重新赋值给a，即改为a := swaprow(a, k, i);。

• 主元判断逻辑错误：
  代码中if a[i,i] = 0的判断逻辑不符合高斯消元规则，应该检查当前主元位置a[k,k]是否为0，若为0才需要交换行，正确判断应为if a[k,k] = 0 then ...，同时要遍历下方行找到第一个非零元素的行来交换。

• 线性代数包混用冲突：
  with(MTM)和linalg包混用可能导致函数命名冲突，建议统一使用Maple推荐的LinearAlgebra包，替换相关函数：rowdim改为RowDimension，det改为Determinant，swaprow改为SwapRows，矩阵拼接用<A | c>更直观。

修正后的代码 #

Gauss := proc(A::Matrix, c::Vector) 
    local a, i, k, j, p, m; 
    uses LinearAlgebra;
    m := RowDimension(A);
    # 先检查矩阵是否为方阵
    if ColumnDimension(A) <> m then
        print("矩阵必须是方阵");
        return;
    end if;
    if Determinant(A) = 0 then
        print("矩阵不可逆");
    else
        # 构造增广矩阵
        a := <A | c>;
        for k from 1 to m-1 do
            # 主元为0时，找下方非零行交换
            if a[k,k] = 0 then
                for i from k+1 to m do
                    if a[i,k] <> 0 then
                        a := SwapRows(a, k, i);
                        break;
                    end if;
                end do;
            end if;
            # 执行消元操作
            for i from k+1 to m do
                p := a[i,k]/a[k,k];
                for j from k to m+1 do
                    a[i,j] := a[i,j] - p*a[k,j];
                end do;
            end do;
        end do;
        # 返回上三角形式的增广矩阵
        return a;
    end if;
end proc;

补充说明 #

• 移除了冗余的形式参数n，直接从矩阵A获取行列数，先做方阵检查避免后续行列式计算报错。
• 使用uses LinearAlgebra;替代with，减少全局命名冲突。
• 优化了主元判断与行交换逻辑，确保消元过程正确执行。
• 最终返回处理后的上三角增广矩阵，方便后续回代求解方程组。

内容的提问来源于stack exchange，提问作者Chafi Aya