Jacobi迭代算法matlab版

function Jacobi(A,b,n,x0,E,N)
k=1
x=zeros(n,1);
 fprintf('/n**********start**************/n',k);
for k=1:N
    fprintf('%d/t',k);
    for i=1:n
        x(i)=b(i)/A(i,i);
        for j=1:n
            if j~=i
                x(i)=x(i)-A(i,j)*x0(j)/A(i,i);
            end
        end
        fprintf('%f/t',x(i));
    end
    if norm(x-x0)<E
            break;
    else
        k=k+1;
        x0=x;
    end
    fprintf('/n');
end
 fprintf('/n***********end*************/n',k);
x0
x