format long
  
A1=[3 0 4;7 4 2;-1 1 2];
D1=[3 0 0;0 4 0;0 0 2];

A2=[-3 3 -6;-4 7 -8;5 7 -9];
D2=[-3 0 0;0 7 0;0 0 -9];

A3=[4 1 1;2 -9 0;0 -8 -6];
D3=[4 0 0;0 -9 0;0 0 -6];

b=[1 2 3];
x0=[1 1 1];
x=x0';
b=b';

A=A1;
D=D1;
eplot(1)=norm(A*x-b);
for j=1:100
 x=inv(D)*b+(eye(3,3)-inv(D)*A)*x;
 eplot(j+1)=norm(A*x-b);
end
plot(eplot)