Testing the output

> evalm(L&*C&*transpose(L));

> equal(simplify(%),A);