x,y = var('x, y')
h(x,y) = [-5*x+6*y, -3*x+4*y]
T = linear_transformation(QQ^2, QQ^2, h)
x1 = vector([1,0]);x2 = vector([0,1])
y1 = vector([1,1]);y2 = vector([2,1])
B=column_matrix([x1,x2,y1,y2])
C=B.echelon_form()
P=C.submatrix(0,2,2,2)
A = T.matrix(side='right')
print A
print
print P
print
print P.inverse()*A*P
|
|
[-5 6]
[-3 4]
[1 2]
[1 1]
[ 1 0]
[ 0 -2]
[-5 6]
[-3 4]
[1 2]
[1 1]
[ 1 0]
[ 0 -2]
|