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