%Jacobian迭代程序
function [x,sp]=jac(a,b,n,x1)
%誤差
e=ones(n,1);
%迭代的解向量
x2=zeros(n,1);
%迭代的次數(shù)
k=0;
%當誤差沒有滿足要求時繼續(xù)迭代
while norm(e,2)>1e-6
   %每隔5步顯示迭代結(jié)果
   if (rem(k,5)==0)
      str=sprintf('X%d= %g %g %g %g\n',k,x1);
      disp(str)
   end
   k=k+1;
   %計算x2（i）
   for i=1:n
      s=0;
      for j=1:n
         s=s+a(i,j)*x1(j);
      end
      x2(i)=(b(i)-s)/a(i,i)+x1(i);
   end
   %計算誤差向量
   e=x2-x1;
   %把x2作為新的迭代初值
   x1=x2;
end
%返回迭代步數(shù)
sp=k;
%返回迭代值
x=x2;