Comprobando

> M:=l[1]:

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

> map(simplify,%);