function[x,XA]=GaussXQByOrder(A,b)%高斯順序消元法N=size(A);n=N(1);fori=1:(n-1)forj=(i+1):nif(A(i,i)==0)disp("對角元素為0！");%防止對角元素為0return;endl=A(j,i);m=A(i,i);A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m;%消元方程b(j)=b(j)-l*b(i)/m;endendx=SolveUpTriangle(A,b);%通用的求上三角係數矩陣線性方程組的函式XA=A;%消元后的係數矩陣functionx=SolveUpTriangle(A,b)N=size(A);n=N(1);fori=n:-1:1if(i<n)s=A(i,(i+1):n)*x((i+1):n,1);elses=0;endx(i,1)=(b(i)-s)/A(i,i);end程式如上，自己算吧。這個東西應該嘗試自己程式設計 參考資料： matlab常用演算法
function[x,XA]=GaussXQByOrder(A,b)%高斯順序消元法N=size(A);n=N(1);fori=1:(n-1)forj=(i+1):nif(A(i,i)==0)disp("對角元素為0！");%防止對角元素為0return;endl=A(j,i);m=A(i,i);A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m;%消元方程b(j)=b(j)-l*b(i)/m;endendx=SolveUpTriangle(A,b);%通用的求上三角係數矩陣線性方程組的函式XA=A;%消元后的係數矩陣functionx=SolveUpTriangle(A,b)N=size(A);n=N(1);fori=n:-1:1if(i<n)s=A(i,(i+1):n)*x((i+1):n,1);elses=0;endx(i,1)=(b(i)-s)/A(i,i);end程式如上，自己算吧。這個東西應該嘗試自己程式設計 參考資料： matlab常用演算法