\[ 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.365323 (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 = 0.398 (sec), leaf count = 36
\[ \left \{ -{\frac {1}{2}\sqrt {{x}^{2}+2\,x+1-4\,y \left ( x \right ) }}-{\frac {{x}^{4}}{4}}-{\frac {{x}^{3}}{3}}-x-{\it \_C1}=0 \right \} \]