jacobi.m

matlab下雅克比在数值计算中有很好的应用!

function x=jacobi(A,b,p,delta,max1)
%p是初值,delta是误差界,max1是给定的迭代最高次数
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