function x=nags(A,b,x0,e,N)
% 用途：用向量形式(滿元素存儲)的Gauss-Seidel迭代解線性方程組Ax=b
% 格式: x=nags(A,b，x0,e,N),A為系數(shù)矩陣，b為右端向量，x返回解向量。
%       x0為初值向量(默認(rèn)原點(diǎn))，e為精度(默認(rèn)1e-4),設(shè)置迭代次數(shù)上限以防發(fā)散(默認(rèn)500)。
n=length(b);
if nargin<5, N=500;end
if nargin<4, e=1e-4;end
if nargin<3, x0=zeros(n,1);end
x=x0;x0=x+2*e; 
k=0; Al=tril(A);iAl=inv(Al);
while norm(x0-x,inf)>e&k<N,
   k=k+1;
   x0=x;x=-iAl*(A-Al)*x0+iAl*b;
   %disp(x')
end
if k==N, warning('已達(dá)迭代次數(shù)上限');end