[tex]
\mu'(t) + 9 \mu(t) = 15[/tex]
is an inhomogeneous 1st order linear ODE. The integrating factor is:
[tex]
A(t) = \exp \left(\int{9 \, dt} \right) = e^{9 t}[/tex]
Then, we have:
[tex]
\frac{d}{dt} \left( e^{9 t} \, \mu(t) \right) = 15 \, e^{9 t}[/tex]
Integrate once:
[tex]
e^{9 t} \, \mu(t) = \frac{15}{9} e^{9 t} + C_1[/tex]
where [itex]C_1[/itex] is an arbitrary integration constant. We have a general solution:
[tex]
\mu(t) = \frac{5}{3} + C_1 \, e^{-9 t}[/tex]
To find [itex]C_1[/itex], we need an inital conditon. Look at the integral equation:
[tex]
\mu_x(t) = 2 + \int_{0}^{t}{\left(15 - 9 \, \mu_x(s) \right) \, ds}[/tex]
Substitute [itex]t = 0[/itex]. The integral is zero because the upper and lower bound coincide! We have:
[tex]
\mu_x(0) =2[/tex]
From here, we have:
[tex]
2 = \frac{5}{3} + C_1 \Rightarrow C_1 = \frac{1}{3}[/tex]
Thus, the mean is:
[tex]
\mu_x(t) = E[X_t] = \frac{5 + e^{-9 t}}{3}[/tex]