Why do I keep getting an infinity loop when I try to implement the Jacobi method?

1 次查看（过去 30 天）

显示更早的评论

Ruten9 2015-12-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/259055-why-do-i-keep-getting-an-infinity-loop-when-i-try-to-implement-the-jacobi-method

关闭： MATLAB Answer Bot 2021-8-20

在 MATLAB Online 中打开

I derived B from this linear system:

1a - 1b + 2c= 8 
1a + 1b + 2c = 10
1a + 0b + 1c = 2
clear all;clc;
c=[8 10 2]';
B= [0 1 -2;-1 0 -2;-1 0 0];
tol = 10^(-10);
x0=[1,1,1]'; x(:,1)= x0;
k=1;
while 1
    xnew(:,k+1) = (B*x(:,k)) + c;
    if norm(xnew(:,k+1)-x(:,1)) < tol
        break;
    end
x=xnew;
k=k+1;
end
k