I must be doing something wrong but it appears that the linsolve() function in SymPy returns the incorrect solution in the following code:
A = Matrix([ [ 0, 0
