网站设计与网页制作教程广州市新闻发布

求解线性方程组的G-S迭代法

// 运行不成功啊function [x,k,index]= Gau_Seid(A,b,ep,it_max)
% 求解线性方程组的G-S迭代法，其中
% A为方程组的系数矩阵
% b为方程组的右端项
% ep为精度要求，省缺为1e-5
% it_max为最大迭代次数，省缺为100
% x为方程组的解
% k为迭代次数
if nargin <4 it_max = 100; end
if nargin <3 ep = 1e-5; end
n = length(A);
k = 0;
x = zeros(n,1);
y = zeros(n,1);
index = 1;
while 1y = x;for i = 1:nz = b(i);for j = 1:nif j~=iz = z-A(i,j)*x(j);endendif abs(A(i,j))<1e-10 | k==it_maxindex = 0;return;endz = z/A(i,j);x(i) = z;endif norm(y-x,inf)<epbreak;endk = k+1;
end

matlab中调用上述函数结果显示：

哪里出问题了啊？