function [x,XA]=GaussXQByOrder(A,b)
%高斯順序消元法
N = size(A);
n = N(1);
for i=1:(n-1)
for j=(i+1):n
if(A(i,i)==0)
disp("對角元素為0!"); %防止對角元素為0
return;
end
l = 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;
x=SolveUpTriangle(A,b); %通用的求上三角係數矩陣線性方程組的函式
XA = A; %消元后的係數矩陣
function x=SolveUpTriangle(A,b)
N=size(A);
n=N(1);
for i=n:-1:1
if(i
function [x,XA]=GaussXQByOrder(A,b)
%高斯順序消元法
N = size(A);
n = N(1);
for i=1:(n-1)
for j=(i+1):n
if(A(i,i)==0)
disp("對角元素為0!"); %防止對角元素為0
return;
end
l = 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;
end
end
x=SolveUpTriangle(A,b); %通用的求上三角係數矩陣線性方程組的函式
XA = A; %消元后的係數矩陣
function x=SolveUpTriangle(A,b)
N=size(A);
n=N(1);
for i=n:-1:1
if(i