Testing the output

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

> map(simplify,%);