Euclidean metric
> A := matrix([[1, 2], [4, 2], [5, 3]]);
> S:=SCHMIDT(A); Returns an ortogonal (usual metric) basis in Col(A)
Testing the output