Comprobando

> DIAGONAL(A):

> evalm(transpose(P)&*A&*P-%);