Let $\displaystyle Z = X + Y$ where $\displaystyle X \sim N\left(\mu, \sigma^2 \right)$ and $\displaystyle Y \sim \Gamma\left(k, \theta \right)$ using this parametrization of the Gamma distribution. Also assume $\displaystyle X$ and $\displaystyle Y$ are independent. Then what is the distribution (pdf) of $\displaystyle Z$?

----------

This question doesn't seem as straightforward as it sounds. For example, I have tried using the convolution formulas here, but can't seem to find a closed form expression for the integral. I have also tried multiplying the moment generating functions (mgfs) of $\displaystyle X$ and $\displaystyle Y$, but it does not seem to match up to any known mgfs.

Does anyone have any ideas on how to find the distribution for $\displaystyle Z$?