Euclidean metric in R^n

> A := matrix([[1, 2], [4, 2], [5, 3]]); Returns a free family in R^3

> S:=GRAM_SCHMIDT(A); Returns an ortogonal basis in R^3

Testing the output