function [x,n]=jacobi(A,b,x0,eps,varargin)
if nargin==3
    eps= 1.0e-6;
    M  = 200;
elseif nargin<3
    error
    return
elseif nargin ==5
    M  = varargin{1};
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=1;                  %迭代次數

while norm(x-x0)>=eps
    x0=x;
     x=B*x0+f;
    n=n+1;
    if(n>=M)
        disp('Warning: 迭代次數太多，可能不收斂！');
        return;
    end
end