\[ 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.606804 (sec), leaf count = 107
\[\text {Solve}\left [\frac {x^4}{4}+\frac {x^3}{3}-\frac {1}{4} \sqrt {x^2+8 y(x)-2 x+1}+\frac {1}{4} \log \left (-\sqrt {x^2+8 y(x)-2 x+1}-x+1\right )-\frac {1}{4} \tanh ^{-1}\left (\frac {2 x-2}{2 \sqrt {x^2+8 y(x)-2 x+1}}\right )-\frac {1}{8} \log (y(x))+x=c_1,y(x)\right ]\] ✓ Maple : cpu = 0.213 (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 \} \]