Comprobando

> M:=l[1]:

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

> modulariza(%,2);