gauss.asv

来自「Numerical analysis （Gauss、Bisection、jaco」· ASV 代码 · 共 23 行

function s=Gauss(a,b,x0,eps,N,)
% Gauss-Seidel迭代法解线性方程组
% a为系数矩阵，b为方程组ax=b的右端值，x0为初值
if nargin==4
    eps=1.0e-6;
elseif nargin<4
%    error
    return
end
D=diag(diag(a));%求对角矩阵
L=tril(a,-1);%求严格下三角矩阵
U=triu(a,1);%求严格上三角矩阵
C=inv(D+L);
B=-C*U;
f=C*b;
s=B*x0+f;
i=0;
while (norm(s-x0)>=eps&&(i<=N))
    x0=s;
    s=B*x0+f;
    i=i+1;
end
return