Return a block diagonal matrix formed from the input list of matrices
bdiag(x)
x
a list of matrices
Returned is a matrix of the form diag(x1, ..., xn) where n is the length of the list.