function [RA,RB,n,X]=gaus(A,b)
%用高斯消元法解线性方程组
B=[A b];n=length(b);RA=rank(A);RB=rank(B);zhica=RB-RA;
if zhica>0
disp('RA~=RB,无解');
end
if RA==RB
if RA==n
disp('有唯一解');
X=zeros(n,1);
for p=1:n-1 %n阶矩阵的1到n-1列
for k=p+1:n % n阶矩阵的p+1到n行
m=B(k,p)/B(p,p);
B(k,p:n+1)=B(k,p:n+1)-m*B(p,p:n+1);
end
end
b=B(1:n,n+1);A=B(1:n,1:n);X(n)=b(n)/A(n,n);
for q=n-1:-1:1
X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q);
end
else
disp('无穷多解');
end
end
%用高斯消元法解线性方程组
B=[A b];n=length(b);RA=rank(A);RB=rank(B);zhica=RB-RA;
if zhica>0
disp('RA~=RB,无解');
end
if RA==RB
if RA==n
disp('有唯一解');
X=zeros(n,1);
for p=1:n-1 %n阶矩阵的1到n-1列
for k=p+1:n % n阶矩阵的p+1到n行
m=B(k,p)/B(p,p);
B(k,p:n+1)=B(k,p:n+1)-m*B(p,p:n+1);
end
end
b=B(1:n,n+1);A=B(1:n,1:n);X(n)=b(n)/A(n,n);
for q=n-1:-1:1
X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q);
end
else
disp('无穷多解');
end
end