Integral involving a confluent hypergeometric function I have the following integral involving a confluent hypergeometric function:
$$\int_{0}^{\infty}x^3e^{-ax^2}{}_1F_1(1+n,1,bx^2)dx$$
where $a>b>0$ are real constants, and $n\geq 0$ is an integer.
Wolfram Mathematica returns the following solution: $\frac{a^{n-1}(a+bn)}{(a-b)^{n+2}}$.  However, I can't figure out how it arrived at it (I always try to check the solutions "on paper" that Mathematica gives me -- or at least using Gradshteyn and Ryzhik).  Can anyone help?
 A: You can probably resolve it by recognizing that $_{1}F_{1}(1+n,1,bx^{2})$
is the Confluent HyperGeometric function $M(1+n,1,bx^{2})$ 
substitute $y=x^{2}$.   You have to keep the b for the next step.
Render the first term negative by the first  Kummer transform 
$M(1+n,1,by)=e^{y}M(-n,1,-by)$ and set $-b=k$.  This makes the series representation have a finite ($n$) number of terms if you want; it might not be necessary .
Now we can use http://dlmf.nist.gov/13.10.E3 for a direct integration/series or we can recognize that $M(-n,1,-by)$ is the 0 th Laguerre Polynomial like so:
$M(-n,1,-by)= \frac{n!}{(1+n)_{n}} L^{0}_{n}(-by)$
Either way you have a finite series to sum.
I leave summing the series as an exercise for the reader (no I haven't done it!).
A: Let's start with the hypergeometric function. We have:
\begin{eqnarray}
F_{2,1}[1+n,1;b x^2] &=& 
\sum\limits_{m=0}^\infty  \frac{(1+n)^{(m)}}{m!} \cdot \frac{(b x^2)^m}{m!} \\
&=& \sum\limits_{m=0}^\infty \frac{(m+1)^{(n)}}{n!} \cdot \frac{(b x^2)^m}{m!} \\
&=&\left. \frac{1}{n!} \frac{d^n}{d t^n}  \left( t^n \cdot e^{b x^2 \cdot t} \right)\right|_{t=1} \\
&=& \frac{1}{n!} \sum\limits_{p=0}^n \binom{n}{p} n_{(p)} (b x^2) ^{n-p} \cdot e^{b x^2}
\end{eqnarray}
Therefore the integral in question reads:
\begin{eqnarray}
&&\int\limits_0^\infty x^3 e^{-a x^2} F_{2,1}[1+n,1;b x^2] dx =\\ 
&& \frac{1}{2} \sum\limits_{p=0}^n \frac{n^{(p)}}{p! (n-p)!} b^{n-p} \frac{(n-p+1)!}{(a-b)^{n-p+2}} \\
&&\frac{1}{2} \frac{1}{(a-b)^{n+2}} \sum\limits_{p=0}^n \binom{n}{p} \frac{(n-p+1)^{(p)}}{(n-p+2)^{(p-1)}} \cdot b^{n-p} (a-b)^p \\
&& \frac{1}{2} \frac{1}{(a-b)^{n+2}} \sum\limits_{p=0}^n \binom{n}{p} (n-p+1) \cdot b^{n-p} (a-b)^p \\
&& \frac{1}{2} \frac{1}{(a-b)^{n+2}} a^{n-1} (a+b n)
\end{eqnarray}
