Comprobando

> evalm(inverse(P)&*A&*P):

> modulariza(%,5);