blkdiag function

Block diagonal matrix