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