附matlab源码如下：

function  x = GuassFull( A,B )
%说明A是方阵，b是列向量
%example:
%A=[1,2,3;4,8,9;10,11,12];
%B = [2,3,5];
%x =GuassFull(A,B');
% Created By wangxiaole 14.9.18
% This function use Guass methods to calculate Linear Equations
% Choose the maximum element in all Matrix A;
[m,n] = size(A);
if rank(A) == mif  m == 1 %只有一个方程if A ~=0x = B/A;elsex = inf;endelseMaxEle = max(max(abs(A)));[m1,n1] = find(abs(A) == MaxEle);Row = A(:,n1(1));A(:,n1(1))=[];A = [Row,A];Column  =A(m1(1),:);A(m1(1),:) =[];A = [Column;A];B_Column = B(m1(1));B(m1(1)) = [];B =[B_Column ;B];for i = 2:mtemp  = A(i,1);A(i,:) = A(i,:) + -1*temp/A(1,1)*A(1,:);B(i)= -1*temp/A(1,1)*B(1) +B(i);endA1 = A;B1 = B;A1(1,:)=[];A1(:,1)=[];B1(1)= [];x = GuassFull(A1,B1);tempx = (sum(-1*A(1,2:end).*x') +B(1))/A(1,1);if(x == inf)x = inf;elseif Row == 1x = [tempx;x];elsex = [x(1:n1(1)-1);tempx;x(n1(1):end)];endendend
elsex = inf;
endend