Resolvamos este problema mediante la variación de los parámetros.
Se nos da:
$$y''-4y'+3y=\dfrac{1}{1+e^{-x}}$$
La ecuación homogénea asociada es $y'' - 4 y' + 3 y = 0$ que tiene una solución complementaria de:
$$y_c = c_1 e^x + c_2 e^{3x}.$$
Para encontrar la función de Green, necesitamos el Wronskian de $(y_1(t), y_2(t))$ de la solución homogénea, por lo que dejamos que $y_1 = c_1 e^t$ y $y_2 = c_2 e^{3t}$ Así pues, tenemos:
$$\text{Wronskian} = W(t) = \text{Wronskian} (c_1 e^t, c_2 e^{3t}) = 2 c_1 c_2 e^{4 t}$$
Para formar la función de Green, tenemos:
$$G(x, t) = \dfrac{y_1(t)y_2(x) - y_1(x)y_2(t)}{W(t)} = \dfrac{c_1e^t c_2e^{3x} - c_1e^x c_2e^{3t}}{2 c_1 c_2 e^{4 t}} = \dfrac{1}{2}\left(e^{-3t}e^{3x} - e^x e^{-t}\right)$$
Para este problema, tenemos $f(t) = \dfrac{1}{1 + e^{-t}}$
Para encontrar la solución particular, establecemos y resolvemos:
$$y_p = \int G(x, t)f(t)dt = \int \dfrac{1}{2}\left(e^{-3t}e^{3x} - e^x e^{-t}\right)\left(\dfrac{1}{1 + e^{-t}}\right) dt$$
Estas dos integraciones dan como resultado:
$$\displaystyle y_p(x) = -\frac{e^x x}{2}+\frac{e^{2 x}}{2}-\frac{1}{2} e^{3 x} \ln(e^{-x}+1)+\frac{1}{2} e^x \ln(e^x+1) -\frac{1}{4}e^x$$
Por supuesto, para formar la solución final, tenemos:
$$\displaystyle y(x) = y_c + y_p = c_1 e^x+c_2 e^{3 x}-\frac{e^x x}{2}+\frac{e^{2 x}}{2}-\frac{1}{2} e^{3 x} \ln(e^{-x}+1)+\frac{1}{2} e^x \ln(e^x+1)-\frac{1}{4}e^x$$
Nota: Podemos combinar el $e^x$ en una única constante $\left(c_1 - \dfrac{1}{4}\right)e^x$ o simplemente $c_1 e^x$ pero lo he dejado de esta forma para que puedas trabajarlo y duplicar los pasos.