高斯赛德尔
function[x] = G_S(A,b,x0,e,M)
n = size(A,1);
x = zeros(n,1);
for m = 1 : M
for p=1 : n
s = 0;
t = 0;
for q=1 : p-1
s = s + A(p,q) * x(q);
end
for q = p + 1 : n
t = t + A(p,q) * x0(q);
end
x§ = -(s + t - b§)/A(p,p);
end
if norm(x - x0, Inf) < e
break;
el
原创
2021-06-07 10:21:09 ·
440 阅读 ·
0 评论