De una matriz simetrica real

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

> rango_signatura(A);

Evaluando el cambio de base