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