A = ones(5) + diag(5:9);
b = [-10 -6 0 8 18]';
[x1 x2 elteres] = Gauss_Jacobi_it(A,b,20);
