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