Sobre Z/5Z(a,b,c,d)

> A:=matrix(2,2,[a/c,5*b/d,5*b/d,a/c]);

> DIAGONALIZA(A,P,5);

> op(P);

Comprobando