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