\[ 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 = 0.422984 (sec), leaf count = 69
\[\left \{\left \{y(x)\to \frac {1}{144} \left (72 c_1 x^4+96 c_1 x^3+288 c_1 x-144 c_1^2-9 x^8-24 x^7-16 x^6-72 x^5-96 x^4-108 x^2+72 x+36\right )\right \}\right \}\]
✓ Maple : cpu = 1.476 (sec), leaf count = 34
\[ \left \{ {\it \_C1}-{\frac {{x}^{4}}{2}}-{\frac {2\,{x}^{3}}{3}}-2\,x-\sqrt {{x}^{2}+2\,x+1-4\,y \left ( x \right ) }=0 \right \} \]