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