I have 10 binary variables x0,..,,x4 and y0,..,y4. I know yi=xi+x(i+1)mod. I want to write x_i in terms of y variables. I know for this toy case we can use Gaussian elimination. But in actual case y_is are nonlinear equations. I tried this idea:
n = 5
RING = PolynomialRing(GF(2),2*n,['x%d'%(i) for i in range(n)] + ['y%d'%(i) for i in range(n)])
RING.inject_variables()
Z = list(RING.gens())
E=[y0==x0+x1, y1==x1+x2, y2==x2+x3, y3==x3+x4, y4==x4+x0]
solve(E, x0,x1,x2,x3,x4)
It is not working.