function [s1,s2]=shifte(A)
%
%  [s1,s2]=shifte(A)
%
%     Copyright: Eduardo Casas (Octubre, 1999)
%
%  Esta función calcula el valor de dos shifts excepcionales reemplazando al shift de 
%  Wilkinson correspondiente a una matriz A, nxn, con n > 2. 
%
%  La salida [s1,s2] son los dos shifts excepcionales a usar en las dos siguientes
%  iteraciones del método QR.
%
n=length(A);
p=abs(A(n-1,n-2))+abs(A(n,n-1));
q=p^2;
p=1.5*p;
s1=(p+sqrt(p^2-4*q))/2;
s2=conj(s1);