数值计算方法雅可比迭代matlab迭代

clear;
A=[1,0.4,0.4;0.4,1,0.8;0.4,0.8,1];
b=[1;1;1];
[m,n]=size(A);
L=-tril(A);
U=-triu(A);
D=diag(diag(A));
B=eye(m,n)-D\A;
f=D\b;
e=1e-3;
x0=zeros(m,1);
x1=B*x0+f;
while(max(abs(x0-x1))>e)
    x0=x1;
    x1=B*x0+f;
end
disp(x0)