integration - Solving defined integral

Is there an analytical solution to the following integral:

$$ I = \iint\limits_{\mathcal{D}} \exp\left(-kx\right) \mathrm{d}x \mathrm{d}y $$

Where:

$$ \mathcal{D}(x,y) \equiv x^2 + y^2 \leq R^2 $$

And:

$$ R, k \in \mathbb{R}^+_0$$

This integral arose in a simple problem (detector response to a cylindrical cell using Beer-Lambert law) but I am struggling to solve it. I did a polar coordinate change, but then I got:

$$ I = \iint\limits_{\mathcal{D}} \rho \exp \left( -k\rho \cos(\theta) \right) \mathrm{d}\rho \mathrm{d}\theta $$

Which seems not easier to solve.

I looked for this integral form in Handbook of Integrals but I cannot found it. I tried to solve it with a symbolic solver without success. I tried to apply the Green's Theorem but the mixed exponential/trigonometric term reappeared.

Edit:

As pointed out by Santosh Linkha and JJacquelin (featuring Mathematica), this integral can be solved using First Kind modified Bessel function. Readers which are not used to Bessel functions - as I was, might be intersted to know that resolution of this integral requires the following identities:

$$ I_n(x) = \frac{1}{\pi}\int\limits_0^\pi \cos(n\theta)\exp(x\cos(\theta))\mathrm{d}\theta $$

And:

$$ \frac{\mathrm{d}}{\mathrm{d}x} \left( x^\nu I_\nu(x) \right) = x^\nu I_{\nu-1}(x) $$

Answer

I cannot do better than Mathematica !