% Test 1, calculo numerico (521230), S2-2011
% Tema 6: Sistemas de ecuaciones lineales
% Pregunta 2

B = [6 -1 0 0 0 0;-1 6 -1 0 0 0;0 -1 6 -1 0 0;0 0 -1 6 -1 0;0 0 0 -1 6 -1;0 0 0 0 -1 6];
D = diag(diag(B));

A = [B .5*D zeros(6,6);.5*D B .5*D;zeros(6,6) (2/5)*D B];

% gradiente conjugado no es conveniente, A no es simetrica
A-A'

% descomposicion de A
D=diag(diag(A));
E=-(tril(A)-D);
F=-(triu(A)-D);

% matriz de iteracion de Jacobi
J = inv(D)*(E+F);
max(abs(eig(J)))

% matriz de iteracion de GS
GS = inv(D-E)*F;
max(abs(eig(GS)))

