\[ y'(x)=x^2 \sqrt {x^2-4 y(x)+2 x+1}+\frac {x}{2}+\frac {1}{2} \] ✓ Mathematica : cpu = 0.265984 (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.294 (sec), leaf count = 26
\[c_{1}-\frac {2 x^{3}}{3}-\sqrt {x^{2}+2 x +1-4 y \relax (x )} = 0\]