赛德尔迭代(MatLab)

%计算方法书上p62 例3.5
%赛德尔迭代
%程序：张建波
%2005-11-18 PM23:00
a=[10 -2 -1;-2 10 -1;-1 -2 5;-1 -2 5];
b=[3 15 10];
x=[0;0;0];
N=3;%N 是解得个数，X1,X2,X3
es=0.0001; %精度控制
while(0<1)
    es0=0;
    %开始迭代
    for i=1:N
        s=x(i);
        x(i)=0;
        x(i)=(b(i)-a(i,:)*x)/a(i,i);
        es0=max(abs(s-x(i)),es0);
    end
    x
    if es0<es break;
    end
  end