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