X <- array(c(rnorm(120)),c(6,5,4))dim(X)# matricized arrayXa <- supermat(X)$Xa
# matricized X with the A-mode entities in its rowsdim(Xa)# matricized X with the B-mode entities in its rowsXb <- permnew(Xa,6,5,4)dim(Xb)# matricized X with the C-mode entities in its rowsXc <- permnew(Xb,5,4,6)dim(Xc)