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