b.diag function

Create a Block Diagonal Matrix