I consider a matrix $M$ which I transform into a system of equations similar as this (but with a different $M$)
sage: M=matrix(3,3,[1,2,3,4,5,6,7,8,9])
sage: P=PolynomialRing(GF(p),M.nrows(),names="x")
sage: (vector(P.gen(i) for i in range(3))*M).list()
[x0 + 4x1 + 7x2, 2x0 + 5x1 + 8x2, 3x0 + 6x1 + 9x2]
(Example taken from another equation where $p$ is a prime).
When I try to create the ideal generated by this system, I get the following error : TypeError: R must be a commutative ring. Any idea how I can fix this?