Evaluando el cambio de base
> A:=matrix(3,3,[0,2+I,3-I,2-I,0,4-2*I,3+I,4+2*I,0]):
> rango_signatura(A,P);
> op(P);
Comprobando