LA.Ch7.P9.SC.Kim

2045 days ago by seongchan

x1=vector([1,1,0]) x2=vector([0,1,2]) x3=vector([1,2,2]) A=matrix([x1,x2,x3]) [G,mu]=A.gram_schmidt() print G print B=matrix([G.row(i) / G.row(i).norm() for i in range(0, 2)]); B print print B*B.transpose() print print B.transpose()*B 
       
[   1    1    0]
[-1/2  1/2    2]


[1 0]
[0 1]

[ 5/9  4/9 -2/9]
[ 4/9  5/9  2/9]
[-2/9  2/9  8/9]
[   1    1    0]
[-1/2  1/2    2]


[1 0]
[0 1]

[ 5/9  4/9 -2/9]
[ 4/9  5/9  2/9]
[-2/9  2/9  8/9]