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