shiyan56.m

用Matlab软件以及雅克比迭代和高斯－赛德尔迭代解方程组Ax=b

r=1;R=6;v=18;n=10;
b1=sparse(1,1,v,n,1);
b=full(b1);
a1=triu(r*ones(n,n));
a2=diag(R*ones(1,n));
a3=-tril(R*ones(n,n),-1)+tril(R*ones(n,n),-2);
a=a1+a2+a3;
I=a\b
I0=sum(I)

r=1;R=6;v=18;n=20;
b1=sparse(1,1,v,n,1);
b=full(b1);
a1=triu(r*ones(n,n));
a2=diag(R*ones(1,n));
a3=-tril(R*ones(n,n),-1)+tril(R*ones(n,n),-2);
a=a1+a2+a3;
I1=a\b
I2=sum(I1)