cdrs001 发表于 2009-5-10 21:14

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)

ChaChing 发表于 2009-5-10 23:42

回复 楼主 cdrs001 的帖子

死机? 我会猜while回圈有问题!
页: [1]
查看完整版本: SOR法