There are $n$ equations.
I need answer for the case $n=3$.
$$ \frac{1}{x_1}(1-x_1)^2+\frac{1}{x_2}(1-x_2)^2+\cdots+\frac{1}{x_n}(1-x_n)^2=0, $$ and $$ \frac{1}{x_1}(1-x_1^2)^j+\frac{1}{x_2}(1-x_2^2)^j+\cdots+\frac{1}{x_n}(1-x_n^2)^j=0,\; 2\leq j\leq n,$$
where $x_j\geq 1$ or $x_j \leq -1$.
If $n=2,3,4$, computer told me that the only real solution is $(1,1,\cdots,1)$. I did not tell computer the condition $|x_j|\geq 1$.
I do not know if some inequality can be applied like chebyshev's inequality and rearrangement inequality. It seems also a little like the Vandermonde determinant.
Comments: : The following three answers perfectly solve the case $n=3$, two of them drop the condition $|x_j|\geq 1$. The case $n=3$ is enough for me, however, for more interested you can try with larger $n\geq 4$. I am also very curious about that.
Solution for $n=3$
For $|x|,|y|,|z|\ge1$, we aim to solve the system \begin{align}\frac{(1-x)^2}x+\frac{(1-y)^2}y+\frac{(1-z)^2}z&=0\\\frac{(1-x^2)^2}x+\frac{(1-y^2)^2}y+\frac{(1-z^2)^2}z&=0\\\frac{(1-x^2)^3}x+\frac{(1-y^2)^3}y+\frac{(1-z^2)^3}z&=0.\end{align} After expanding each equation and subtracting appropriately, we find an interesting equivalence \begin{align}x+y+z+\frac1x+\frac1y+\frac1z&=6\\x^3+y^3+z^3+3\left(\frac1x+\frac1y+\frac1z\right)&=12\\x^5+y^5+z^5+5\left(\frac1x+\frac1y+\frac1z\right)&=18\end{align} which we can attack using the elementary symmetric polynomials $e_1=\sum_{\rm cyc}x,e_2=\sum_{\rm cyc}xy,e_3=xyz$. Using Newton's identities, this is equivalent to \begin{align}e_1+\frac{e_2}{e_3}&=6\\e_1^3-3e_1e_2+3e_3+\frac{3e_2}{e_3}&=12\\e_1^5-5e_1^3e_2+5e_1^2e_3+5e_1e_2^2-5e_2e_3+\frac{5e_2}{e_3}&=18.\end{align}
The rest of my original solution which is quite tedious can be found here, and I will illustrate a much simpler solution by @MykolaPochekai:
Letting $e_2=re_3$, we have $e_1=6-r$ and $e_3=\dfrac{(12-3r)-(6-r)^3}{3(r^2-6r+1)}$. Substituting into the third equation yields $$\frac{(r-3)^4(r^3-24r^2+189r-498)}{9(r^2-6r+1)}=0.$$ Note that $|r|\le3$ from $|x|,|y|,|z|\ge1$, and it is easy to show the cubic has no real roots over this interval. Therefore, $e_1=3$, $e_2=3$ and $e_3=1$ and Vieta tell us each of $x,y,z$ solves $t^3-3t^2+3t-1=0$ so $x=y=z=1$ is the unique solution. Note that for completeness, the denominator is zero when $r=3-2\sqrt2$ but this is irrelevant since $12-3r-(6-r)^3\ne0$ and $e_3$ is finite. $\square$