Metrica dada por una matriz hermitiana definida positiva

> F:=matrix(4,4,[1,1+I,2,3-I,1-I,4,2-I,3*I,2,2+I,8,3-I,3+I,-3*I,3+I,50]):

> rango_signatura(F);

> B := matrix([[1+I, 2-I], [4*I, 2], [5-I, 3+2*I],[0,2-3*I]]);

> S:=SCHMIDT(B,F); Da una base ortogonal de Col(B) en la metrica unitaria definida por F

Comprobacion