SOR法
请高手帮忙看看,哪错了。。。一运行就死机。。、
function Sor(A,b,w)
D=diag(diag(A));
L=-tril(A,-1);
U=-triu(A,1);
Lw=(D-w*L)\((1-w)*D+w*U);
f=w*(D-w*L)\b;
x=zeros(size(b));
while (norm(b-A*x,inf))>=1e-5
x=Lw*x+f;
end
disp(x)
回复 楼主 cdrs001 的帖子
死机? 我会猜while回圈有问题!
页:
[1]