System of multiple nonlinear equations vs. system of linear equations with nonlinear constraints in Python

Question

I am trying to solve in Python a system of 22 nonlinear equations with 22 variables x_1, ..., x_{22}. The system is linear on each variable, i.e., I only have products of difere