I am trying to obtain a probability density function (PDF) from its moment generating function (MGF) and the problem is the nature of the MGF. My MGF is $M(t)=I_o^{\beta}(s)$, where $\beta$ is a positive integer: \begin{equation} M(t)=\left( \frac{1}{\pi}\int_0^{\pi} e^{t\cos \theta} d\theta \right)^{\beta} \end{equation}
I am stuck in finding the PDF from $M(t)$, as this is not a standard form for which a look-up table type approach can be used.
I tried replacing $t$ by $it$ and then taking the inverse Fourier Transform route, but it didn't yield anything. Even if exact solution cannot be found, some sort of approximation or bound will also be helpful.