Can anyone help me on deriving the closed form expression for the following definite integral of the modified Bessel function of the first kind, zero order, involving power and exponential functions
$$I={\displaystyle\int\limits_0^{a } {{x^{b}} \cdot \exp \left( { - \frac{{{x^2}}}{2}} \right){I_0}\left( {c x} \right)dx} }$$
where $a$, $b$, and $c$ are constant.
Many thanks for your help!