How to solve an equation of the form $f(x)=f(a)$ for a fixed real a.
Let $$f(x)=\frac{(x^2-x+1)^3}{x^2(x-1)^2}\tag{1}$$
(representative curve in Fig. 1).
Your question can be reformulated in the following way :
$$\text{for a given} \ a, \ \ \text{find all} \ x \ \text{such that} \ \ f(x)=f(a) \tag{2}$$
This answer will have 3 parts : A) Solving the problem itself, B) Understanding a little the associated group, C) Considering function $f$ in a larger scope.
Fig. 1. The red dots correspond to: $f(-3)=f(-1/3)=f(1/4)=f(3/4)=f(4/3)=f(4)=2197/144$ where the values $-3,-1/3,\cdots 4$ are the values $\varphi_k(-3)$ (the first one corresponding to identity function with notations found below). The vertical asymptotes and the asymptotic parabola (with equation $y=x^2-x+3$) are also featured.
A) Solving the problem itself
First of all, the variations of $f$ are easily established once we remark that
$$\tag{3}f'(x)=\underbrace{(x^2-x+1)}_{>0 \ \text{for all} \ x}\dfrac{(x-2)(2x-1)(x+1)}{x^3(x-1)^3}.$$
We deduce from (3) the existence of 3 minima in $(-1,a), (0.5,a)$ and $(2,a)$ where $a=27/4=6.75.$
We can observe (by computations) the following facts, with the introduction of spécific notations
$f(x)=f(\underbrace{1-x}_{\varphi_2(x)})$.
$f(x)=f(\underbrace{1/x}_{\varphi_3(x)})$.
In a natural way the composition of these two invariant transforms, $\varphi_3\circ\varphi_2$ and $\varphi_2\circ\varphi_3$, gives rise to a new invariance property :
$f(x)=f(1-1/x)$ meaning that if $x$ is solution, $a=1-1/x$ is also a solution ; in the same vein :
$f(x)=f(1/(1-x))$ meaning that if $x$ is solution, $a=1/(1-x)$ is also a solution.
A last invariant transformation :
- $f(x)=f(x/(x-1))$ meaning that if $x$ is solution, $a=x/(x-1)$ is also a solution.
We can summarize all this in the following manner:
$$\tag{4}f(x)=f(\varphi_k(x))$$
for all "homographic" functions belonging to the so-called (6 elements) anharmonic group (in the sense of group theory):
$$\phi_1(x)=x, \ \phi_2(x)=1-x, \ \phi_3(x)=\tfrac{1}{x}, \ \phi_4(x)=1-\tfrac{1}{x}, \ \phi_5(x)=\tfrac{1}{1-x}, \ \phi_6(x)=\tfrac{x}{x-1}.$$
Conclusion : Equation (2) has
6 solutions if $a\neq -1,1/2,2$ : $ \ x=a, \ 1-a, \ \tfrac{1}{a}, \ 1-\tfrac{1}{a}, \ \tfrac{1}{1-a}, \ \tfrac{a}{a-1} $
3 solutions if $a=-1,1/2,2$ which are... $x=-1,1/2,2$
Are there other transformations that we haven't considered ? No, because for a given $a$, i.e, for a given $b=f(a)$, there are 5 other points of the curve with a common ordinate $b$, not less not more (set apart the special case of the minimal points). But in fact the fundamental interplay between function $f$ and the group just described can only be understood through "higher algebra" considerations that we will see in part C.
Remark : we could have taken profit of the fact that the curve has a symmetry with respect to a vertical axis situated at $x=\dfrac12$. Indeed, setting $x=X+\dfrac12$, we get a new (even) equation simpler than the initial one :
$$Y=\dfrac{(4X^2+3)^2}{4 (4X^2-1)^2}$$
B) Understanding a little the associated anharmonic group
(see (https://en.wikipedia.org/wiki/Cross-ratio))
It can be considered as a (finite) subgroup of $PGL(2,\mathbb{Z})$ through a classical linear representation (in short : representation by matrices) :
$$\dfrac{az+b}{cz+d} \ \ \ \leftrightarrow \ \ \ \begin{pmatrix}a&b\\c&d\end{pmatrix}$$
It is isomorphic to $S_3$, the group of permutation on 3 objects.
It is generated by $\varphi_2$ and $\varphi_3$.
It has a nice vizualization in terms of "fundamental regions":
Fig. 2 : if $z$ belongs to the green "triangular" region, then $1-z$ belongs to... etc. Please note that straight line $x=\tfrac12$ is an axis of symmetry like for Fig. 1.
C) Considering function $f$ in a larger scope
Why is this working so well ? Because this is the tip of an iceberg. Indeed, function $f$ is the "Klein j-invariant" (see (https://arxiv.org/pdf/1810.08742.pdf)(https://math.stackexchange.com/q/3418400) (http://www.cis.upenn.edu/~jean/gma-v2-chap5.pdf) exercice 5.14) associated with anharmonic group in projective geometry ; it would be too long to explain it here. See for that the Wikipedia article (https://en.wikipedia.org/wiki/Cross-ratio). This invariant is associated with elliptic curves.
Remarks : other connections : How is $E = S(x)$, where $E = k(x), k$ a field, and S= k(I) of all rationals functions of $I = I(x) = \frac{(x^2 -x+1)^3}{x^2(x-1)^2}$? ; about Lüroth theorem https://upload.wikimedia.org/wikipedia/commons/3/3e/Another_elementary_proof_of_Luroth%27s_theorem-06.2004.pdf ; Klein's invariant and negative discriminant} ; an interesting inequality :https://artofproblemsolving.com/community/c6h1447345p8277143
factorizing the given equation and cancel the denominators we get $$(ax-a-x)(ax-x+1)(x-1+a)(-x+a)(ax-1)(ax-a+1)=0$$
$$\frac{(x^2-x+1)^3}{x^2(x-1)^2}=\frac{(a^2-a+1)^3}{a^2(a-1)^2}$$
Now multiply both sides by $x^2(x-1)^2$ :
$$(x^2-x+1)^3=x^2(x-1)^2\frac{(a^2-a+1)^3}{a^2(a-1)^2} \\ \implies (x^2-x+1)^3-x^2(x-1)^2\frac{(a^2-a+1)^3}{a^2(a-1)^2}=0$$
Without expanding, it can be written as :
$$x^6+a_1x^5+a_2x^4+a_3x^3+a_4x^2+a_5x+1=0$$
Since, product of roots of this equation is $1$, each solution's reciprocal will also be a solution to this equation.