function x= Gauss_C(A,b)
% Gauss列主元消元法求方程组Ax=b
% 此处显示详细说明
% A为系数矩阵
Ab = [A b]; % 增广矩阵A|b
m = size(A,1); % m为A的行数
n = size(A,2); % n为A的列数
x = zeros(m,1);
for i=1:m-1
[~,r]=max(abs(Ab(i:m,i))); %找到第i到m行,第i列最大的元素k及其所在的行,
Ab([i,r+i-1],:)=Ab([r+i-1,i],:); %做行交换
for j =i:m-1
Ab(j+1,:) = Ab(j+1,:)-Ab(j+1,i)/Ab(i,i)*Ab(i,:); % 得到上三角矩阵
end
end
A = Ab(:,1:n);
b = Ab(:,n+1); %将A,b分离出来
x(m) = b(m)/A(m,m); %先求出x(m)
for h = m-1:-1:1 % 回代求x(m-1):x(1)
temp=0;
for l = h+1:m
temp = temp+A(h,l)*x(l);
end
x(h) = (b(h)-temp)/A(h,h);
end
end