%通過Gauss-Seidel迭代求解線性方程組
function x=G_S(A,b,x0,delta,max1)

%Input  - A要求解的線性方程組的系數矩陣
%       - b是線性方程組標準形式Ax=b中的b列向量
%       - x0是初始迭代值
%       - delta可以是接受的誤差
%       - max1為最大迭代次數
%Output - x為所求的解

n=length(b);
x=zeros(n,1);
for i=1:max1
    for j=1:n
        x(j)=(b(j)-A(j,1:j-1)*x(1:j-1)-A(j,j+1:n)*x0(j+1:n))/A(j,j);
    end
    err=norm(x-x0);
    relerr=err/norm(x);
    x0=x;
    if (err<delta)|(relerr<delta)
        break;
    end
end