function s=shift(A)
%
%  s=shift(A)
%
%     Copyright: Eduardo Casas (Octubre, 1999)
%
%  Esta función calcula el shift de Wilkinson correspondiente a una matriz.
%  La salida s es el valor propio de la matriz A, 2x2, más próximo a A(2,2).
%
% Esta función necesita las funciones eig2.m y root2.m.
%
l=eig2(A);
t=abs(A(2,2)-l);
if t(1) < t(2)
   s=l(1);
else
   s=l(2);
end