簡單迭代法

A=input('請輸入待求解的係數矩陣：\n');
b=input('輸入b：');
nmax=('最大迭代次數是幾？');
eps=('誤差是幾？');
x0=('初始解是幾？');
n=size(A);
for k = 1:100
    for i=1:n
        s=0;
        for j=1:i-1
            s=s+A(i,j)*x(j);
        end
        for j =i+1:n
            s=s+A(i,j)*x(j);
        end
        y(i)=(b(i)-s)/A(i,i);
     end
     for i=1:n
         if abs(y(i)-x(i))<=eps*abs(x(i))
             y;
         else
             disp('nmax次迭代後不收斂.')
     x=y;
     end
     end
end
x