x, y = var('x, y')
h(x, y) = [x-y,x+y]
T = linear_transformation(QQ^2, QQ^2, h)
x1=vector([1, 1])
x2=vector([-1, 0])
y1=vector([1, 1])
y2=vector([-1, 0])
B=column_matrix([y1, y2, T(x1), T(x2)])
print B
print
C=B.echelon_form()
print C
print
A=C.submatrix(0, 2, 2, 2)
print A
|
[ 1 -1 0 -1]
[ 1 0 2 -1]
[ 1 0 2 -1]
[ 0 1 2 0]
[ 2 -1]
[ 2 0]
[ 1 -1 0 -1]
[ 1 0 2 -1]
[ 1 0 2 -1]
[ 0 1 2 0]
[ 2 -1]
[ 2 0]
|