function [x,n]=jacobi(A,b,x0,eps)
if nargin==3
eps= 1.0e-6;
M = 1000;
elseif nargin<3
disp ('error');
return
elseif nargin ==5
M = 1000;
end
D=diag(diag(A)); %求A的对角矩阵
L=-tril(A,-1); %求A的下三角阵
U=-triu(A,1); %求A的上三角阵
B=D\(L+U);
f=D\b;
x=B*x0+f;
n=0; %迭代次数
while norm(x-x0)>=eps
x0=x;
x=B*x0+f;
n=n+1;
if(n>=M)
disp('Warning: 迭代次数太多，可能不收敛！');
return;
end
end
