 function  Z=rot(k,l,n,a)
% ======================================================================
% Rotation
%	S.F.
% ----------------------------------------------------------------------
Z=eye(n)    ;
cos(a)
Z(k,k)=cos(a);
Z(k,l)=sin(a);
Z(l,k)=-sin(a);
Z(l,l)=cos(a);

return;















