Combine matrices to block diagonal structure
blockdiag(x, ..., pad = 0)
x
Matrix
...
Additional matrices
pad
Vyalue outside block-diagonal
Klaus K. Holst
A <- diag(3)+1 blockdiag(A,A,A,pad=NA)