\[ y'(x)=x^2 \sqrt {x^2+8 y(x)-2 x+1}-\frac {x}{4}+\frac {1}{4} \] ✓ Mathematica : cpu = 0.259776 (sec), leaf count = 37
DSolve[Derivative[1][y][x] == 1/4 - x/4 + x^2*Sqrt[1 - 2*x + x^2 + 8*y[x]],y[x],x]
\[\left \{\left \{y(x)\to \frac {1}{72} \left (16 x^6-96 c_1 x^3-9 x^2+18 x-9+144 c_1{}^2\right )\right \}\right \}\] ✓ Maple : cpu = 0.287 (sec), leaf count = 26
dsolve(diff(y(x),x) = -1/4*x+1/4+x^2*(x^2-2*x+1+8*y(x))^(1/2),y(x))
\[c_{1}+\frac {4 x^{3}}{3}-\sqrt {x^{2}-2 x +1+8 y \left (x \right )} = 0\]