版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/qq_38669138/article/details/83272152
简单迭代法
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