Ch-7-Prob-22-New-김강섭

1858 days ago by matrix

p=matrix(3,6,[1,2,1,1,1,1,2,3,3,3,4,2,3,5,7,2,6,4]) a=p.rref() b=matrix(3,3,[5,4,0,-5/3,-5/3,1/3,-2/3,1/3,1/3]) c=b.inverse() d=vector([1,1,1]) e=vector([3,2,1]) print a print print b print print c print print c*d print print b*e 
       
[   1    0    0    5    4    0]
[   0    1    0 -5/3 -5/3  1/3]
[   0    0    1 -2/3  1/3  1/3]

[   5    4    0]
[-5/3 -5/3  1/3]
[-2/3  1/3  1/3]

[ 1/3  2/3 -2/3]
[-1/6 -5/6  5/6]
[ 5/6 13/6  5/6]

(1/3, -1/6, 23/6)

(23, -8, -1)
[   1    0    0    5    4    0]
[   0    1    0 -5/3 -5/3  1/3]
[   0    0    1 -2/3  1/3  1/3]

[   5    4    0]
[-5/3 -5/3  1/3]
[-2/3  1/3  1/3]

[ 1/3  2/3 -2/3]
[-1/6 -5/6  5/6]
[ 5/6 13/6  5/6]

(1/3, -1/6, 23/6)

(23, -8, -1)