Comprobando
>
evalm(inverse(P)&*A&*P);