clear all
A=ones(6)+diag(7*ones(6,1));
b=(1:6)';
x1=A\b;
D = diag(diag(A));
T = D-A;
x2 = zeros(size(b));
for k = 1:20
	x2 = D\(T*x2+b);
end
norm(x1-x2)