Putting $y=f(x)$, your equation becomes:
$y^4 - 2xy^2 + x^2 - x - y = 0$
This is hard to solve for $y$, but we have the suggestion that $y$ has the form $y=\frac{a+\sqrt{cx+d}}{b}$. This can be substituted into our polynomial and then simplified:
$\left(\frac{a+\sqrt{cx+d}}{b}\right)^4 - 2x\left(\frac{a+\sqrt{cx+d}}{b}\right)^2 + x^2 - x - \frac{a+\sqrt{cx+d}}{b} = 0$
$\implies\left(a+\sqrt{cx+d}\right)^4 - 2b^2x\left(a+\sqrt{cx+d}\right)^2 +b^4x^2 - b^4x - b^3\left(a+\sqrt{cx+d}\right) = 0$
$\implies a^4+4a^3\sqrt{cx+d}+6a^2(cx+d)+4a(cx+d)^{3/2}+(cx+d)^2-2a^2b^2x-4ab^2x\sqrt{cx+d}-2b^2x(cx+d)+b^4x^2-b^4x-ab^3-b^3\sqrt{cx+d}=0$
$\implies a^4+6a^2(cx+d)+(cx+d)^2-2a^2b^2x-2b^2x(cx+d)+b^4x^2-b^4x-ab^3=b^3\sqrt{cx+d}+4ab^2x\sqrt{cx+d}-4a^3\sqrt{cx+d}-4a(cx+d)^{3/2}$
$\implies a^4+6a^2(cx+d)+(cx+d)^2-2a^2b^2x-2b^2x(cx+d)+b^4x^2-b^4x-ab^3 = (b^3+4ab^2x-4a^3-4a(cx+d))\sqrt{cx+d}$
$\implies (a^4+6a^2(cx+d)+(cx+d)^2-2a^2b^2x-2b^2x(cx+d)+b^4x^2-b^4x-ab^3)^2 = (b^3+4ab^2x-4a^3-4a(cx+d))^2(cx+d)$
If you multiply this all out, and move everything to the left, you'll have a polynomial of the form $g(x)=M_4x^4 + M_3x^3 + M_2x^2 + M_1x + M_0 = 0$, where each of the $M_i$'s is a polynomial in $a$, $b$, $c$ and $d$. Since $g(x)=0$ for all $x>1$, we can say that $M_4=M_3=M_2=M_1=M_0=0$. If there is a solution to this system of five polynomials in four variables, then that's your answer.
In general, systems of polynomials can be difficult, but certain software packages allow for Groebner basis calculations, which can work, depending on the computer's capabilities and the user's patience.
One hopes there's an easier way to solve this, but the above is a sort of brute-force approach.
EDIT The solution posted here by Hagen von Eitzen using limits is considerably nicer.