use bmatrix