    function Q = flip(Y,Z);
                
        [n, toss] =  size(Y);
        N = Y-Z; c = sqrt(N'*N);
        if c > 1.0E-9     
             N = N/c;  % normalize       
             Q = eye(n) - 2*N*N';
        else
             Q = eye(n);
        end