Ch.8 Problem.15

1300 days ago by g1.kjiwon

E=matrix(2, 2, [sqrt(3), 0, 0, 1]) F=E.inverse() y=matrix(2, 1, [0, 0]) print F.augment(y).transpose() 
        
[1/3*sqrt(3)           0]
[          0           1]
[          0           0]
[1/3*sqrt(3)           0]
[          0           1]
[          0           0]
V=matrix(3, 3, [1/6*sqrt(6), 2/6*sqrt(6), 1/6*sqrt(6), 1/2*sqrt(2), 0, -1/2*sqrt(2), 1/3*sqrt(3), -1/3*sqrt(3), 1/3*sqrt(3)]) U=matrix(2, 2, [1/2*sqrt(2), -1/2*sqrt(2), 1/2*sqrt(2), 1/2*sqrt(2)]) print V.transpose()*F.augment(y).transpose()*U.transpose() 
        
[1/36*sqrt(2)*sqrt(3)*sqrt(6) - 1/2 1/36*sqrt(2)*sqrt(3)*sqrt(6) + 1/2]
[      1/18*sqrt(2)*sqrt(3)*sqrt(6)       1/18*sqrt(2)*sqrt(3)*sqrt(6)]
[1/36*sqrt(2)*sqrt(3)*sqrt(6) + 1/2 1/36*sqrt(2)*sqrt(3)*sqrt(6) - 1/2]
[1/36*sqrt(2)*sqrt(3)*sqrt(6) - 1/2 1/36*sqrt(2)*sqrt(3)*sqrt(6) + 1/2]
[      1/18*sqrt(2)*sqrt(3)*sqrt(6)       1/18*sqrt(2)*sqrt(3)*sqrt(6)]
[1/36*sqrt(2)*sqrt(3)*sqrt(6) + 1/2 1/36*sqrt(2)*sqrt(3)*sqrt(6) - 1/2]