I'm looking to solve the equation $$ x^2y'' + xy' - (x^2+p^2)y = 0 $$ in particular for $p = \frac12$. I've already determined the indicial equation to be $r^2 - p^2 = 0$, and so the roots are $r = \pm p$. I've also already determined a series solution (which seems to be a constant factor off, namely $a_0$, from the modified Bessel function $I_p$): $$ y_1 = a_0x^r\sum_{n=0}^\infty\left(\prod_{i=1}^n i(i+p)\right)^{-1}\left(\frac x2\right)^{2n} $$ Now, for $p = \frac12$, I noticed that $x^{-1/2}e^x$ was a solution (by brute-force). However, I can't for the life of me figure out why plugging $p = \frac12$ into the above series equation should give me something that looks remotely like $x^{-1/2}e^x$.
To go through my steps, I assumed $y = \sum_{n\ge 0}a_nx^{n+r}$, and substituted this into the DE to get $$ (r^2 - p^2)a_0x^r + ((r+1)^2 - p^2)a_1x^{r+1} + \sum_{n=0}^\infty ([(n+r)^2 - p^2]a_n - a_{n-2})x^{n+r} = 0 $$ Clearly $r = \pm p$ in order for $a_0\neq 0$. However, it's the $x^{r+1}$ coefficient which is bugging me. Namely, for $p = \frac12$, taking the root $r = -p = -\frac12$ gives $(r + 1)^2 - p^2 = 0$, which means $a_1\neq 0$ can be anything, meaning the odd powers of $x$ also have coefficients in the series solution. This doesn't effect the solution $y = x^{1/2}\sum(\cdots)$ since that still requires $a_1 = 0$, however the solution $y = x^{-1/2}\sum(\cdots)$ changes, which is why I think $x^{-1/2}e^x$ might still arise from the Frobenius solution.
That said, I simply can't figure out how to get to $x^{-1/2}e^x$. Even if you account for the coefficients of the odd powers the series doesn't look like $e^x$ at all as far as I can tell. Are you simply supposed to arrive at $x^{-1/2}e^x$ from another method?