1 | initial version |
If you want the Gram-Schmidt basis you could compute it and then normalize by hand. This is suboptimal though.
sage: M = Matrix([[1,2],[3,4]])
sage: [v.normalized() for v in M.gram_schmidt()[0].rows()]
[(1/5*sqrt(5), 2/5*sqrt(5)), (2*sqrt(1/5), -sqrt(1/5))]