Comprobando

> DIAGONAL(A,complejo):

> evalm(htranspose(P)&*A&*P-%); Se trata de una metrica hermitiana