ch8.16_noh

328 days ago by kyunganoh

A=matrix(QQ,[[1,1,2],[0,2,3],[3,7,0]]) print A.rank() B=A.transpose() * A print B Pseudo=B^-1*A.transpose() print Pseudo 
       
3
[10 22  2]
[22 54  8]
[ 2  8 13]
[  7/8 -7/12  1/24]
[ -3/8   1/4   1/8]
[  1/4   1/6 -1/12]
3
[10 22  2]
[22 54  8]
[ 2  8 13]
[  7/8 -7/12  1/24]
[ -3/8   1/4   1/8]
[  1/4   1/6 -1/12]