function [ Q ] = randQ(n)
%randQ - Generate a random rotation matrix
%   returns a random element of SO(n), uniformly distributed according
%   to Haar measure

A = randn(n,n);
[Q,~] = qr(A);
if det(Q) < 0
    Q(1,:) = -Q(1,:);
end

end