E' semplice, ma poiché il secondo ordine é solo finto, faccio una poderosa abbreviazione di calcolo
Integrando rispetto a x abbiamo subito
y' + y = ln (ex + 5) + C1
e per l'omogenea associata a questa yo(x) = C e-x
Per la variazione delle costanti y = v(x) e-x
per cui omettendo alcune delle dipendenze da x per brevità
v' e-x + v e-x * (-1) + v e-x = ln (ex + 5) + C1
v' e-x = ln (ex + 5) + C1
v' = ex ln (ex + 5) + C1 ex
v(x) = S ln (ex + 5) d(ex) + C1 ex + C2 =
= [ t ln t - t ]_(t = ex + 5) + C1 ex + C2 =
= (ex + 5) ln (ex + 5) - (ex + 5) + C1 ex + C2
= (ex + 5) ln (ex + 5) + (C1 - 1) ex + C2 - 5 =
= (ex + 5) ln (ex + 5) + C1 ex + C2
y(x) = v(x) e-x = (1 + 5 e-x ln (ex + 5) + C1 + C2 e-x.