jacobi.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 20 行

function X=jacobi(A,b,P,delta,max1)
%  A是n维非奇异阵。
%  b是n维向量。
%  P是初值。
%  delta是误差界。
%  max1是给定的迭代最高次数。
%  X为所求的方程组AX=b的近似解。
N=length(b);
for k=1:max1
for j=1:N
  X(j)=(b(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j);
    end
err=abs(norm(X'-P));
P=X';
    if(err<delta)
           break
  end
end
X=X';k,err;