高斯消元法 Gaussian elimination

function [x]=gauss(A,b)         % 将A转换为上三角矩阵
    n=size(A,1);                % 获取矩阵A的行数
    for k=1:n-1                 % 遍历每一列进行消元
        for i=k+1:n             % 处理当前列下方的每一行
            t=A(i,k)/A(k,k);    % 计算消元系数
            for j=k+1:n         % 子循环：更新当前行的元素
                A(i,j)=A(i,j)-t*A(k,j);
            end
            b(i)=b(i)-t*b(k);   % 更新右端向量b
        end
    end
   
    % 回代过程
    x(n)=b(n)/A(n,n);           % 求解最后一个未知数xn
    for k=n-1:-1:1              % 从倒数第二个未知数开始回代
        s=0;                    % 初始化求和变量
        for l=k+1:n
            s=s+A(k,l)*x(l);    % 计算已知项的线性组合
        end
        x(k)=(b(k)-s)/A(k,k);   % 求解当前未知数
    end
end