Revision history [back]

Dimension 1 tells us about the system having too many "degrees of freedom". Since we work in a finite field, it is easy to deal with it by fixing one of the variables, say, x:

variety = sum( (Ideal(polynomials+[x-x0]).variety() for x0 in Fq), [] )

which produces 16 solutions.