\[ y'(x)=x^2 \sqrt {x^2-4 y(x)+2 x+1}+\sqrt {x^2-4 y(x)+2 x+1}+x^3 \sqrt {x^2-4 y(x)+2 x+1}+\frac {x}{2}+\frac {1}{2} \] ✓ Mathematica : cpu = 1.44182 (sec), leaf count = 103
\[\text {Solve}\left [\frac {x^4}{4}+\frac {x^3}{3}+\frac {1}{2} \sqrt {x^2-4 y(x)+2 x+1}+\frac {1}{2} \log \left (\sqrt {x^2-4 y(x)+2 x+1}+x+1\right )-\frac {1}{2} \tanh ^{-1}\left (\frac {2 x+2}{2 \sqrt {x^2-4 y(x)+2 x+1}}\right )-\frac {1}{4} \log (y(x))+x=c_1,y(x)\right ]\] ✓ Maple : cpu = 0.384 (sec), leaf count = 34
\[\left \{-\frac {x^{4}}{2}-\frac {2 x^{3}}{3}+c_{1}-2 x -\sqrt {x^{2}+2 x -4 y \left (x \right )+1} = 0\right \}\]