function GS=GS(A,b,x0,eps)
%該函數文件是通過Gauss-Seidel選代格式求解線性方程組的近似解
%其中,A為方程組的系數矩陣,b為方程組的右端向量,x0為選代初值
if nargin==3
    eps=1.0e-6;
elseif nargin<3
    error
    return
end
D=diag(diag(A));    %求解對角陣
L=-tril(A,-1);    %求解出嚴格下對角陣
U=-triu(A,1);    %求解出嚴格上對角陣
G=(D-L)\U;    %求解出迭代矩陣
f=(D-L)\b;
GS=G*x0+f;
while norm(GS-x0)>=eps
    x0=GS;
    GS=G*x0+f;
end