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