高斯消元法 matlab程序

function[A,b,x]=jvzhen(A,b);
[z,n]=size(A);
%A为方阵，所以也可以为m=size(A,2)
for i=1:z-1
    
    for t=i+1:z
        %从第i+1行始消第i行的，一直消到最后一行：假设i为1，那么就是2，3，4……行全都乘第一行把第一个元素消掉
        m(t,i)=-A(t,i)/A(i,i);
        %后一行首位/前一行首位，得到消元乘数，一行只需要消首位所以只需要一个m，不需要参与列的循环
        for j=i:z
            %从第i列开始乘乘数，消去第i个元素
            A(t,j)=A(t,j)+m(t,i)*A(i,j);
            End
            %完成一行的转换，上循环为t行的第i个元素消元
            b(t)=b(t)+m(t,i)*b(i);
            %该行b也随之变化
        end
        %完成所有行第i列元素的消元，因为最后一行最后一列元素不用消，所以循环到z-1
        End
        %回代求解
        for i=1:z
            y(i)=b(i);
        end
        for i=z:-1:1
            t=z;
            while i<t
                y(i)=y(i)-A(i,t)*x(t);
                t=t-1;
            end
            x(i)=y(i)/A(i,i);
        end
    end
end