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