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