Find analytical solution to first order inhomogeneous ODE with nonconstant coefficients I've given the following first order linear inhomogeneous ordinary differential equation (ODE) equation, unfortunately with nonconstant coefficients:
$$
\frac{A(t)+B(t)}{A(t)}\dot y + \frac{B(t)}{C(t)}y = \frac{1}{C(t)}x_0
$$
where $x_0$ is a constant and $A(t)$, $B(t)$ and $C(t)$ are continuous functions, given by
$$
A(t) = a_1\cdot \exp(a_2\cdot t^{a_3}) \\
B(t) = b_1\cdot \exp(b_2\cdot t^{b_3}) \\
C(t) = c_1\cdot \exp(c_2\cdot t^{c_3})
$$
where $a_i$, $b_i$ and $c_i$ are constants.
Now, I want to find an analytical solution for $y(t)$ while $t > 1$, if an IVP $y(t = 1) = y_0 = const.$ is given. And I hope that this is somehow possible.
I know that it is possible to separate the variables of the equation, so that $y(t) = y_\text{H} + y_\text{P}$. To find the solution of the homogeneous part should not be that problem, maybe "only" to solve the integrals. The bigger Problem for me is that I have no idea how I can find the particular solution. And unfortunately I didn't find any example of such equation-type on the web ...
Another Idea from a math-script was to use "integrating factors". Therefor, I rewrote the given ode to:
$$
\dot y + \underbrace{\frac{A(t)\cdot B(t)}{C(t)\cdot (A(t)+B(t))}}_{p(t)}\cdot y = \underbrace{\frac{A(t)}{C(t)\cdot (A(t) + B(t))}}_{g(t)}\cdot x_0
$$
what leads to the following Expression as solution for $y(t)$:
$$
y(t) = \frac{\int\exp(\int p(t) \text{d}t)\cdot g(t)\text{d}t + C_1}{\exp(\int p(t)\text{d}t)}
$$
But now I have the problem that I dont know how to solve the integrals of the solution for $y(t)$ that I can determine the constant $C_1$ for the solution of the IVP.
It would be fine if anyone of you has an idea how to deal with this problem. Or isn't it possible to solve this "thing" analytically?
Thank you in advance,
  pj
 A: With Maple 2017.1 we have:
eq := {(D(y))(t)+p(t)*y(t) = g(t)*x0, y(t0) = y0};
expand(subs(_z1 = s, dsolve(eq)))

$$y \left( t \right) ={\frac {{\it x0}\,\int_{{\it t0}}^{t}\!g \left( s
 \right) {{\rm e}^{\int_{{\it t0}}^{s}\!p \left( s \right) \,{\rm d}s}
}\,{\rm d}s}{{{\rm e}^{\int_{{\it t0}}^{t}\!p \left( s \right) 
\,{\rm d}s}}}}+{\frac {{\it y0}}{{{\rm e}^{\int_{{\it t0}}^{t}\!p
 \left( s \right) \,{\rm d}s}}}}
$$
y(t) = x0*(Int(g(s)*exp(Int(p(s), s = t0 .. s)), s = t0 .. t))/exp(Int(p(s), s = t0 .. t))+y0/exp(Int(p(s), s = t0 .. t))

With Mathematica 11.1 we have: 
 sol = First@DSolve[{y'[t] + p[t]*y[t] == g[t]*x0, y[t0] == y0}, y[t], t]
 Y[T] = y[t] /. sol /. K[1] -> s /. K[2] -> m

$$y(t)=e^{\int_1^t -p(s) \, ds-\int_1^{\text{t0}} -p(s) \, ds} \text{y0}+e^{\int_1^t -p(s) \,
   ds} \int_1^t e^{-\int_1^m -p(s) \, ds} \text{x0} g(m) \, dm-e^{\int_1^t -p(s) \, ds}
   \int_1^{\text{t0}} e^{-\int_1^m -p(s) \, ds} \text{x0} g(m) \, dm$$
  y[t] == E^(Integrate[-p[s], {s, 1, t}] - Integrate[-p[s], {s, 1,t0}])*y0 + 
  E^Integrate[-p[s], {s, 1, t}]*Integrate[(x0*g[m])/E^Integrate[-p[s], {s, 1, m}],{m, 1, t}] - E^Integrate[-p[s], {s, 1, t}]*
  Integrate[(x0*g[m])/E^Integrate[-p[s], {s, 1, m}], {m, 1, t0}]

