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)
|