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