It is known that, for $y^2=x^4-2x^3+3x^2+4x+1$, $$\int\frac{x\,\mathrm{d}x}{y}=\frac13\mathrm{arctanh}\frac{x^2-x+1}{y}+\frac13\mathrm{arctanh}\frac{x^2-x-1}{y}+\frac16\mathrm{arctanh}\frac{x^2-x-1/2}{y}+C.$$
To get this result, I tried to apply Risch's algorithm as described in Manuel Bronstein's Symbolic Integration Tutorial, as follows. The integrand is $xy/D$ where $y^2=D=x^4-2x^3+3x^2+4x+1$ is square-free, so the indefinite integral has only a logarithm part, and we then get to section 2.5. There a procedure is described, which in step 3 computes the polynomial $$R(t)=\mathrm{Resultant}_x(\mathrm{Resultant}_y(xy-t\frac{\mathrm{d}D}{\mathrm{d}x},y^2-D),D).$$ From this formula, we have $R(t)=186624t^8$ for the given integrand. This has no nonzero roots, so by step 6 the integral is not elementary. This contradicts the antiderivative result given above. What is wrong?