  In this paper we deal with Mellin convolution of generalized Gamma densities
which leads to integrals of modified Bessel functions of the second kind. Such
convolutions allow us to explicitly write the solutions of the time-fractional
diffusion equations involving the adjoint operators of a square Bessel process
and a Bessel process.
