Hilbert metric

> F:=linalg[hilbert](3);

> S:=SCHMIDT(A,F); Returns an ortogonal (Hilbert metric) basis in Col(A)

Testing the output