Unitary metric in C^n

> B:=matrix([[1, 1+I], [1-I, 2], [1+I, 3]]):

> S:=GRAM_SCHMIDT(B);

Testing the output