DSolve[Derivative[1][y][x] == y[x]^3/(E^(2*b*x)*(1 + y[x]/E^(b*x))),y[x],x]