code: def randorth(p=10):
      """
      A random orthogonal matrix.
      """
      A = symnormal(p)
      return npl.eig(A)[1]
  
imports: from itertools import combinations
  import numpy
  import numpy.linalg
  import pytest
  from numpy.testing import assert_almost_equal, assert_array_equal
  from ..None import intvol
label_desc: Use image node_id 2 for calling a global function with example usage: randorth(p) and returns: npl