clear;
A=[1,0.4,0.4;0.4,1,0.8;0.4,0.8,1];
b=[1;2;3];
[m,n]=size(A);
L=-tril(A,-1);
U=-triu(A,1);
D=diag(diag(A));
B=(D-L)\U;
f=(D-L)\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)
数值计算方法高斯赛德尔迭代matlab实现
最新推荐文章于 2023-10-11 10:05:29 发布