Create block diagonal matrix
Bdiag(A, k)
A
k
Return a block diagonal matrix from the matrix A.
Bdiag(matrix(1:4,2,2), 3)