雅克比迭代矩阵.txt

本代码是本人在做作业时所编

a=[10 -1 -2;-1 10 -2;-1 -1 5];%雅克比方法计算迭代矩阵
b=[72;83;42];
d1=(diag(a))';
[n,n]=size(a);
I=eye(n,n);
d=diag(d1);
c=inv(d);%D的逆矩阵
f=c*a;
B=I-f;
g=c*b