Unitary complex metric

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

> S:=SCHMIDT(B); Returns an ortogonal (unitary metric) basis in Col(B)

Testing the output