powered by
Y = BLKDIAG(A,B,...) produces diag(A,B,...)
blkdiag(x)
a list of matrix
E - Block diagonal concatenation of matrix