function varargout=Jacobiliu(varargin)A=[5 2 1;2 8 -3;1 -3 -6];b=[8;21;1];x0=[0;0;0];x=Jacobi(A,b,x0)function x=Jacobi(A,b,x0)n=100;%最大迭代次数ee=0.0001;%精度n1=length(b);x1=zeros(n1,1);for i=1:n
for j=1:n1
s=0;
for k=1:n1
if k~=j
s=s+A(j,k)*x0(k);
e
function [j,c]=jacobi(A,b,x0,eps,N)n=size(A);n=n(1);D=diag(diag(A));H=eye(n)-inv(D)*A;f=inv(D)*b;j=H*x0+f;c=0;while (norm(j-x0)>=eps)&(c<N)
c=c+1;
x0=j;