Euclidean metric in R^n
> S:=Gram_schmidt(A); It returns an ortonormal basis in R^3
Testing the output