\[ y'(x)=\frac {x^6 \sqrt {4 x^2 y(x)+1}+x^5 \sqrt {4 x^2 y(x)+1}+x^3 \sqrt {4 x^2 y(x)+1}+\frac {1}{2}}{x^3} \] ✓ Mathematica : cpu = 0.451552 (sec), leaf count = 81
\[\left \{\left \{y(x)\to -\frac {2 c_1 x^5}{5}-\frac {1}{4} \left (2 c_1-1\right ) x^4-c_1 x^2+c_1^2+\frac {x^{10}}{25}+\frac {x^9}{10}+\frac {x^8}{16}+\frac {x^7}{5}+\frac {x^6}{4}-\frac {1}{4 x^2}\right \}\right \}\]
✓ Maple : cpu = 0.312 (sec), leaf count = 34
\[ \left \{ {\it \_C1}-{\frac {1}{x}\sqrt {4\,{x}^{2}y \left ( x \right ) +1}}+{x}^{2}+{\frac {{x}^{4}}{2}}+{\frac {2\,{x}^{5}}{5}}=0 \right \} \]