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