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