Testing the output

> evalm(htranspose(S)&*S); The columns of S are ortogonal