function [omega_ott]=SOROmegaOttimale(A)
%-------------------------------------------
% Inputs:
%  A=matrice del sistema
% Output:
%  omega_ott: il parametro ottimale
%--------------------------------------------
% costruisco la matrice H(omega)
D=diag(diag(A));
B=D-tril(A); C=D-triu(A);

N=500;
w=linspace(.1,1.9,N);
for i=1:N
  H=inv(D-w(i)*B)*((1-w(i))*D+w(i)*C);
  p(i)=max(abs(eig(H)));
end;

plot(w,p); xlabel(' \omega '); ylabel(' p(H_\omega) ');

title('Raggio spettrale della matrice del metodo SOR');
[mp,imp]=min(p); omega_ott=w(imp);