Unitary complex metric

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

> T:=Schmidt(B); Returns an ortonormal (unitary metric) basis in Col(B)

Testing the output