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