function [x,n]=Songchidiedai(A,b,w,eps)
% w為松弛迭代因子，目前對于一般的矩陣A
% 還沒有最好的w，只是對于一些特殊的矩陣才有
% w 為0～2
if nargin==3
    eps=1e-5;
end
D=diag(diag(A));
L=-tril(A,-1);
U=-triu(A,1);
Sw=inv(D-w*L)*((1-w)*D+w*U);
f=w*inv(D-w*L)*b;
if max(abs(eig(Sw)))<1
    x0=zeros(length(A),1);
    x1=Sw*x0+f;
    k=0;
    while norm(x1-x0,inf)>eps
        x0=x1;
        x1=Sw*x0+f;
        k=k+1;
    end
    x=x1;
    n=k;
else
    disp('不滿足jacobi迭代條件！');
end
end