雅可比迭代法与高斯-塞德尔迭代法(2)

简要介绍资料的主要内容,以获得更多的关注

x(k 1) x(k)

10 6时，迭代终止，要求输出对应的近似解。

（2）输出对应的迭代的次数。

（3）输出一个矩阵，每行表示每次迭代的向量。

1．程序：

function [x,i,G]=kb(A,b,max,eps)

A=[10 2 3;2 10 1;3 1 10];

b=[14 11 20]';

max=100;

eps=1e-006

D=diag(diag(A));

ID=inv(D);

J=ID*(-A+D);

f=ID*b;

x0=zeros(rank(A),1);

x=x0;

G(1,:)=x0';

for i=1:max

k=x;

x=J*x+f;

G(i+1,:)=x';

n=norm((x-k),inf);

if n<=eps

break;