绍兴网站制作公司,ipv6网站建设东莞,网络营销专业属于什么类别,网站做二级域名干什么用求解线性方程组的G-S迭代法
// 运行不成功啊function [x,k,index] Gau_Seid(A,b,ep,it_max)
% 求解线性方程组的G-S迭代法#xff0c;其中
% A为方程组的系数矩阵
% b为方程组的右端项
% ep为精度要求#xff0c;省缺为1e-5
% it_max为最大迭代次数#xff0c;省缺为100
% …求解线性方程组的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 | kit_maxindex 0;return;endz z/A(i,j);x(i) z;endif norm(y-x,inf)epbreak;endk k1;
endmatlab中调用上述函数结果显示 哪里出问题了啊
