function s=shifts(A)
%
%  s=shifts(A)
%
%     Copyright: Eduardo Casas (Octubre, 1999)
%
%  Esta función calcula el shift de Wilkinson correspondiente a una matriz
%  simétrica. La salida s es el valor propio de la matriz A, 2x2, más próximo
%  a A(2,2).
%
t=abs(A(1,1)-A(2,2))/2;
r=A(2,1)^2;
if A(1,1) < A(2,2)
   s=A(2,2)+(r/(t+sqrt(t^2+r)));
else
   s=A(2,2)-(r/(t+sqrt(t^2+r)));
end