CircularUnitaryMatrixDistribution
CircularUnitaryMatrixDistribution[n] represents a circular unitary matrix distribution with matrix dimensions {n, n}.
Examples
Sample a random unitary matrix:
m = RandomVariate[CircularUnitaryMatrixDistribution[3]];
Chop[m . ConjugateTranspose[m]]
(* IdentityMatrix[3] *)Please visit the official Wolfram Language Reference for more details.