gauss_seidel.m

来自「学习matlab的朋友」· M 代码 · 共 18 行

function[x,i]=Gauss_Seidel(A,b)
D=diag(diag(A));
U=triu(-A,1);
L=tril(-A,-1);
x=zeros(length(A),1);
i=0;
while  true
  x0=x;  
  x=inv(D-L)*U*x+inv(D-L)*b;
  i=i+1;
  format long
    if (norm(x0-x,inf)<1e-4) && i<1000
        break;
    end
    x
end