列主元高斯消去法解线性方程组的程序（MATLAB版）

clear;clc;
%列主元Gauss校区法解线性方程组
A=[3 -5 6 4 -2 -3 8;    
    1 1 -9 15 1 -9 2;
    2 -1 7 5 -1 6 11;   
    -1 1 3 2 7 -1 -2;
    4 3 1 -7 2 1 1;    
    2 9 -8 11 -1 -4 -1;
    7 2 -1 2 7 -1 9];%系数矩阵
b=[11 2 29 9 5 8 25]';%n维向量
y=inv(A)*b   %matlab的计算结果
n=length(b);%方程个数n
x=zeros(n,1);%未知向量
%-------------消去-----------
for k=1:n-1
    Auk=A(k:n,k);
    [m,u]=max(abs(Auk));
    u=u+k-1 %u为最大元所在的列
    %------交换最大的行和当前行的值-------
    for j=k:n
        temp=A(u,j);A(u,j)=A(k,j);A(k,j)=temp;
    end
    temp=b(k);b(k)=b(u);b(u)=temp;

%     if A(k,k)==0;
%         error('Error');
%     end

    for i=k+1:n
%         A(i,k)=A(i,k)/A(k,k);