\[ y'(x)=\frac {x^5 \sqrt {4 x^2 y(x)+1}+\frac {1}{2}}{x^3} \] ✓ Mathematica : cpu = 0.190646 (sec), leaf count = 31
\[\left \{\left \{y(x)\to \frac {1}{16} \left (-8 c_1 x^4+16 c_1^2+x^8-\frac {4}{x^2}\right )\right \}\right \}\]
✓ Maple : cpu = 0.248 (sec), leaf count = 26
\[ \left \{ {\it \_C1}-{\frac {1}{x}\sqrt {4\,{x}^{2}y \left ( x \right ) +1}}+{\frac {{x}^{4}}{2}}=0 \right \} \]