%順序gauss消去法，gauss函數
function [A,u]=gauss(a,n)
for k=1:n-1
   %消去過程
   for i=k+1:n
      for  j=k+1:n+1
         %如果a(k,k)=0,則不能消去
         if abs(a(k,k))>1e-6
            %計算第k步的增廣矩陣
            a(i,j)=a(i,j)-a(i,k)/a(k,k)*a(k,j);
         else
            %a(k,k)=0,順序gauss消去失敗
            disp('順序gauss消去失敗');
            pause;
            exit;
         end
      end
   end
end
%回代過程
x(n)=a(n,n+1)/a(n,n);
for i=n-1:-1:1
   s=0;
   for j=i+1:n
      s=s+a(i,j)*x(j);
   end
   x(i)=(a(i,n+1)-s)/a(i,i);
end
%返回gauss消去后的增廣矩陣
A=triu(a);
%返回方程組的解
u=x;