\[ y'(x)=x^2 \sqrt {x^2+8 y(x)-2 x+1}+\sqrt {x^2+8 y(x)-2 x+1}+x^3 \sqrt {x^2+8 y(x)-2 x+1}-\frac {x}{4}+\frac {1}{4} \] ✓ Mathematica : cpu = 0.400379 (sec), leaf count = 77
\[\left \{\left \{y(x)\to \left (\frac {4}{3}-c_1\right ) x^4-\frac {4 c_1 x^3}{3}+\left (\frac {1}{4}-4 c_1\right ) x+2 c_1^2+\frac {x^8}{8}+\frac {x^7}{3}+\frac {2 x^6}{9}+x^5+\frac {15 x^2}{8}-\frac {1}{8}\right \}\right \}\]
✓ Maple : cpu = 0.29 (sec), leaf count = 32
\[ \left \{ {\it \_C1}+{x}^{4}+{\frac {4\,{x}^{3}}{3}}+4\,x-\sqrt {{x}^{2}-2\,x+1+8\,y \left ( x \right ) }=0 \right \} \]