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