blodiag function

Block diagonal matrix