#定义矩阵,默认定义在实数域
A = Matrix([[1,2,3,5],[3,2,1,2],[1,1,1,0],[3,7,2,2]])
A^-1
#定义在其他域上的矩阵,如有限域
A = Matrix(GF(13),[[1,2,3,5],[3,2,1,2],[1,1,1,0],[3,7,2,2]])
A^-1
#可以看到两个逆矩阵不一样
#定义向量,定义在有限域,默认定义在实数域
w = vector(GF(13),[1,1,4,3])
Y=A*w;Y
Z=w*A;Z
#解线性方程组
X = A.solve_right(Y);X
#也可以使用符号\
A\Y
X = A.solve_left(Z);X
#格基约减
A = Matrix([[1,2,3,5],[3,2,1,2],[1,1,1,0],[3,7,2,2]])
A.LLL()