\[ y'(x)=x \sqrt {x^2+8 y(x)-2 x+1}-\frac {x}{4}+\frac {1}{4} \] ✓ Mathematica : cpu = 0.268263 (sec), leaf count = 37
\[\left \{\left \{y(x)\to \frac {1}{8} \left (4 x^4-x^2-16 c_1 x^2+2 x-1+16 c_1{}^2\right )\right \}\right \}\] ✓ Maple : cpu = 0.334 (sec), leaf count = 27
\[c_{1}+2 x^{2}+\frac {1}{4}-\sqrt {x^{2}-2 x +1+8 y \relax (x )} = 0\]