for n = [ 1, 10, 25, 50, 100 ]
  U = random_orthogonal_matrix(n);
  actual = norm(U);
  expected = 1;

  msg = sprintf('random_orthogonal_matrix(%d) has norm 1', n);
  unit_test_equals(msg, expected, actual);
end
