dsolve(diff(y(x),x) = y(x)*(y(x)^2+exp(-x^2)*y(x)+exp(-x^2)^2)/exp(-x^2)^2*x,y(x))