% A為方程組的增廣矩陣

% e為誤差容量

function x=myJacobi(A,ep)

if nargin<2 

    ep=1.0e-5;

end

[n,m]=size(A);

 x=zeros(n,1);

 y=zeros(n,1);

 k=0;

 while 1

   for i=1:1:n

        s=0.0;

    for j=1:1:n

     if j~=i

     s=s+A(i,j)*x(j);    

     end

     y(i)=(A(i,n+1)-s)/A(i,i);

     end

    end

     xmx=0.0;

     for i=1:1:n

         xmx=max(xmx,abs(x(i)-y(i)));    

     end

     if xmx<=ep

         for i=1:1:n

             x(i)=y(i);

         end

         return;

     end

     for i=1:1:n

      x(i)=y(i);

      y(i)=0.0;

     end

     k=k+1;

     if k>30

         error('failure');

         return;

     end

 end