Testing the output

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

> modulariza(%,3); We have worked on Z/3Z.