The ODE has characteristic equation
[tex]r^2+6r+9=(r+3)^2=0[/tex]
with roots [tex]r=-3[/tex], and hence the characteristic solution
[tex]y_c=C_1e^{-3t}+C_2te^{-3t}[/tex]
For the particular solution, assume an ansatz of [tex]y_p=ae^{-t}[/tex], with derivatives
[tex]\dfrac{\mathrm dy_p}{\mathrm dt}=-ae^{-t}[/tex]
[tex]\dfrac{\mathrm d^2y_p}{\mathrm dt^2}=ae^{-t}[/tex]
Substituting these into the ODE gives
[tex]ae^{-t}-6ae^{-t}+9ae^{-t}=4ae^{-t}=4e^{-t}\implies a=1[/tex]
so that the particular solution is
[tex]\boxed{y(t)=C_1e^{-3t}+C_2te^{-3t}+e^{-t}}[/tex]