Solution to a higher order ordinary differential equation Let $q$ be a strictly positive integer and let $\beta \neq q$ be a real number. Consider a following Ordinary Differential Equation(ODE):
\begin{equation}
\frac{d^q r_t}{d t^q} + \frac{1}{t^\beta} r_t = 0
\end{equation}
Using Mathematica I have checked that the fundamental solutions to that ODE reads:
\begin{equation}
r_t^{(\eta)} = t^\eta F_{0,q-1}\left[\begin{array}{c} \{\} \\ \left\{ 1+ \frac{\xi}{q-\beta} \right\}_{\xi=1}^\eta  \cup  \left\{ 1- \frac{\xi}{q-\beta} \right\}_{\xi=1}^{q-1-\eta} \end{array}; -\frac{t^{q-\beta}}{(q-\beta)^q}\right]
\end{equation}
where $\eta=0,\cdots,q-1$. Here $F_{p,q}$ is the hypergeometric function. The quantity on the bottom of the square brackets is understood as a concatenation of two sequences the first one of length $\eta$ and the second one of length $q-1-\eta$.
In the special case $q=2$ the fundamental solutions read:
\begin{equation}
\left(
\sqrt{t} J_{\frac{1}{-2+\beta}}\left[\frac{t^{1-\frac{\beta}{2}}}{\frac{\beta}{2}-1}\right]
,
\sqrt{t} J_{\frac{-1}{-2+\beta}}\left[\frac{t^{1-\frac{\beta}{2}}}{\frac{\beta}{2}-1}\right]
\right)
\end{equation}
where $J_\beta$ is the Bessel function. The solutions match those given in Solution to a second order ordinary differential equation .
Is there a way to derive the solutions using some different method?
 A: Note, that integrating our ODE $q$ times and  solving the resulting integral equation by successive iterations we end up with the following solution in terms of infinite series:
\begin{equation}
r_t = \sum\limits_{\eta=0}^{q-1} r^{(\eta)}_0 \sum\limits_{p=0}^\infty \frac{1}{((q-1)!)^p} I^{(t,t_ 0)}_{\vec{q}_\eta,p}
\end{equation}
where $\vec{q}_\eta := (q-1,q-1,\cdots,q-1,\eta)$ where the number $q-1$ is repeated $p$ times. The initial conditions of the ODE are 
\begin{equation}
\frac{d^\eta  r_t }{ d t^\eta}(t= t_0) = r^{(\eta)}_0
\end{equation}
for $\eta=0,\cdots,q-1$ and the quantities $I^{(t,t_ 0)}_{\vec{q}_\eta,p}$ are multivariable integrals over a simplex, defined in Yet another multivariable integral over a simplex . 
For simplicity we only consider the last ($\eta=q-1$ ) term only. Using the closed form solution for the multivariable integral we easily arrive at the following special solution:
\begin{eqnarray}
&&\sum\limits_{l=0}^{q-1} \binom{q-1}{l} (-1)^l \\&& 
F_{0,q-1}\left[\begin{array}{c} \left\{\right\} \\ \left\{ 1 - \frac{l-\eta}{q-\beta}\right\}_{\stackrel{\eta=0}{\eta\neq l}}^{q-1}\end{array}; - \frac{t^{q-\beta}}{(q-\beta)^q}\right]
F_{0,q-1}\left[\begin{array}{c} \left\{\right\} \\ \left\{ 1 + \frac{l+\eta}{q-\beta}\right\}_{\stackrel{\eta=0}{\eta\neq l}}^{q-1}\end{array}; (-1)^q \frac{t_0^{q-\beta}}{(q-\beta)^q}\right] t^{q-1-l} t_0^l
\end{eqnarray}
We will further simplify the result soon.
