简单迭代法

简单迭代法

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