Comprobando el resultado

> op(L);

> evalm(L&*C&*htranspose(L)-B);