A=matrix(4, 3, [2, 3, 7, 2, 5, 4, 1, 2, -3, 3, 3, 7])
B=A.transpose()
C=B*A
print "det(A'A)="
print (B*A).det()
D=matrix(4, 1, [1, 0, 1, 0])
print
print "x="
print C.inverse()*B*D
|
det(A'A)=
3703
x=
[ 2/7]
[ 5/161]
[-12/161]
det(A'A)=
3703
x=
[ 2/7]
[ 5/161]
[-12/161]
|