function [T,V]=SYMQRSTEP1(T,V,f,l)
% [T,V]=SYMQRSTEP1(A,V,f,l)
% T=SYMQRSTEP1(A,f,l)
%
% QR-Schritt mit Wilkinson-Shift für symmetrische (Hermitesche)
% Tridiagonalmatrizen. Dabei wird eine (l-f+1)x(l-f+1)-Matrix Q aus dem
% symmetrischen QR-Schritt für T(f:l,f:l) bestimmt. Mit Hilfe dieser Matrix
% wird T mit seinem neuen Wert T_neu=W'*T*W überschrieben. Dabei gilt
% W(l:f,l:f)=Q und sonst ist W überall die Identität.
% Bei Bedarf wird V mit V_neu=V*W aufdatiert.
% 
% -----------------------------------------------------------------------------
%
% Input
% -----
% l,f           1<=f<=l<=n. Block T(f:l,f:l) von T der durch QR-Schritt
%               bearbeitet wird 
% 
% -----------------------------------------------------------------------------
%
% Input/Output
% ------------
% T             nxn Hermitesche Tridiagonalmatrix
% V             optional, Transformationsmatrix


% Bei drei Parametern werden nur die ersten drei (T,V,f) übergeben. Wir
% interpretieren sie aber als (T,f,l)
if nargin==3
   l=f;
   f=V;
   clear V;
end % if

n=size(T,1);
v=f:l;


