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