急寻用matlab编写的解决雅各比迭代法的程序

function X=jacobi(A,b,P,delta,max1)
% A是n维非奇异阵
% B是n维向量
% P是初值
% delta是误差界
% 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 (e