\[ y'(x)=x^2 \sqrt {x^2-4 y(x)+2 x+1}+\frac {x}{2}+\frac {1}{2} \] ✓ Mathematica : cpu = 0.246718 (sec), leaf count = 37
\[\left \{\left \{y(x)\to \frac {1}{36} \left (-4 x^6+9 x^2+24 c_1 x^3+18 x+9-36 c_1{}^2\right )\right \}\right \}\] ✓ Maple : cpu = 0.485 (sec), leaf count = 26
\[ \left \{ {\it \_C1}-{\frac {2\,{x}^{3}}{3}}-\sqrt {{x}^{2}+2\,x+1-4\,y \left ( x \right ) }=0 \right \} \]