I am having a rather unusual output from linsolve when trying to solve a system of linear equation over algebraic numbers. The system that I have is defined by a matrix and I want to solve the kernel of this matrix. The matrix is rather long (so it is not a good idea to put it in this forum) but I have provided a pastebin link to it: https://pastebin.com/ErJsDxGq
What i do is I defined this matrix as m in giacpy and then I do the following
Code : Tout sélectionner
b = [0]*9 #the answer vector should be 0
print linsolve(m,b) #see m in pastebin link
Code : Tout sélectionner
[-C_1,C_0,C_1,0,-3*C_0,-2*C_1,C_0,C_1,0,C_1]
I don't get such behaviour in XCas. Any ideas? Thanks in advance!Unable to build a single algebraic extension for simplifying.
Trying rational simplification only. This might return a wrong answer if simplifying 0/0!
Jose