DSolve[Derivative[1][y][x] == (1/2 + x^4*Sqrt[1 + 4*x^2*y[x]])/x^3,y[x],x]