Evaluando el cambio de base

> A := matrix([[-5, 99, -50, 31], [99, -61, -12, -26], [-50, -12, -18, -62], [31, -26, -62, 1]]):

> rango_signatura(A,P);

> op(P);

Comprobando