I am trying to solve this integral:

\begin{equation} \frac{1}{c^{b}}\int_{0}^{\infty} x^{n}\, e^{-a x}\, \gamma(b,c(-d+x)) \ \mathrm{d}x \end{equation}

where, $n>0$ is an integer, and $a$, $b$, $c$, $-d$ are all positive. I tried using $u = -d+x$ and then binomial expansion for the resulting $(u+d)^{n}$ term to get the following form:

\begin{equation} \frac{1}{c^{b}}\sum_{j=0}^{n} \ (\mathrm{constant}) \int_{-d}^{\infty} u^{j}\, e^{-a u}\, \gamma(b,cu) \ \mathrm{d}u \end{equation}

Which gets me to the same form of the integral in Gradshteyn and Ryzhik's Table (see entry 6.455):

\begin{equation} \int_{0}^{\infty} x^{\mu-1}\, e^{-\beta x}\, \gamma(\nu,\, \alpha x) \ \mathrm{d}x = \frac{\alpha^{\nu}\, \Gamma(\mu+\nu)}{\nu(\alpha+\beta)^{\mu+\nu}}\, {_{2}}F_{1}\left(1,\, \mu+\nu;\, \nu+1;\frac{\alpha}{\alpha+\beta}\right) \end{equation}

Which of course has different limits of integration. Any ideas on how I might solve this? I am looking for suggestions on how to approach the problem. Any methods of integration and/or helpful identities would also be great.



\begin{equation} \frac{1}{a_{2}^{b}}\int_{0}^{\infty} x^{n}\, e^{-a_{1}x} \, \gamma(b, a_{2}(-p+x)) \ \mathrm{d}x \end{equation}

Use this fact: $\frac{\partial}{\partial x}\, \gamma(a,f(x)) = e^{-f(x)}\, f(x)^{a-1} \, f'(x) \ \mathrm{d}x$ and do integration by parts.

\begin{equation*} \begin{aligned} u &= \gamma(b, a_{2}(-p+x)) &&v = -\frac{\Gamma(n+1,a_{1}x)}{a_{1}^{n+1}}\\ du &= a_{2}^{b}\, e^{-a_{2}(-p+x)}\, (-p+x)^{b-1} \ \mathrm{d}x &&dv = x^{n}\, e^{-a_{1}x} \ \mathrm{d}x \end{aligned} \end{equation*}

Which gives:

\begin{equation} -\frac{\gamma(b, a_{2}(-p+x))\,\Gamma(n+1,a_{1}x)}{a_{1}^{n+1}a_{2}^{b}}\Bigg|_{0}^{\infty} %%%%%%%%%%%%%%%%%%%% +\frac{e^{a_{2}p}}{a_{1}^{n+1}}\int_{0}^{\infty} (-p+x)^{b-1}\, e^{-a_{2}x}\, \Gamma(n+1,a_{1}x) \ \mathrm{d}x \end{equation}

For the first term we need these limits: $\lim_{x\to\infty}\Gamma(a,b x)=0$, and $\lim_{x\to 0^{+}}\Gamma(a,b x)=\Gamma(a)$. Also, use binomial theorem to expand the $(-p+x)^{b-1}$ inside the integral.

\begin{equation} \frac{\gamma(b, -a_{2}p)\,\Gamma(n+1)}{a_{1}^{n+1}a_{2}^{b}} +\frac{e^{a_{2}p}}{a_{1}^{n+1}} \sum_{j=0}^{b-1} \binom{b-1}{j} (-p)^{b-1-j} \int_{0}^{\infty} x^{j}\, e^{-a_{2}x}\, \Gamma(n+1,a_{1}x) \ \mathrm{d}x \end{equation}

Now we can see the integral is of the form in the question:

\begin{equation*} \int_{0}^{\infty} x^{j}\, e^{-a_{2} x}\, \Gamma(n+1,\, a_{1} x) \ \mathrm{d}x = \frac{a_{1}^{n+1}\, \Gamma(j+n+2)}{(j+1)(a_{1}+a_{2})^{j+n+2}}\, {_{2}}F_{1}\left(1,\, j+n+2;\, j+2;\frac{a_{2}}{a_{1}+a_{2}}\right) \end{equation*}

So the final answer is:

\begin{equation} \frac{\gamma(b, -a_{2}p)\,\Gamma(n+1)}{a_{1}^{n+1}a_{2}^{b}} +e^{a_{2}p} \sum_{j=0}^{b-1} \binom{b-1}{j} \frac{(-p)^{b-1-j}\, \Gamma(j+n+2)}{(j+1)(a_{1}+a_{2})^{j+n+2}}\, {_{2}}F_{1}\left(1,\, j+n+2;\, j+2;\frac{a_{2}}{a_{1}+a_{2}}\right) \end{equation}

Not sure if the summation can be solved to simplify the solution but my hunch is that it can be with a bit more work.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.