Generates a random unitary matrix (square)
randomUnitary(size)
size
dimension of matrix
m <- 25 # size of matrix (m x m) set.seed(7) U <- randomUnitary(m) # n <- 5 # First n columns A <- U[,1:n]