Considere la siguiente integral: $$\mathcal{A}=\int_0^\infty\exp\left(-\sqrt2\,x^2\right)\,\operatorname{erfi}(x)\,\log(x)\,x^3\,dx,\tag1$$ donde $\operatorname{erfi}(x)$ denota la imaginaria de la función de error: $$\begin{align}\operatorname{erfi}(x)&=\frac2{\sqrt\pi}\int_0^x\exp \left(z^2\right)dz\tag2\\&={_1F_1}\left(\frac12;\frac32;x^2\right)\cdot\frac{2\,x}{\sqrt\pi}.\end{align}\tag3$$
Para evaluar esto, he intentado reemplazar el poder $3$ con un parámetro de $\alpha$ y eliminar el logarítmica del factor, y tengo esto: $$\int_0^\infty\exp\left(-\sqrt2\,x^2\right)\,\operatorname{erfi}(x)\,x^\alpha\,dx=\frac{{_2F_1}\left(\frac12,\frac\alpha2+1;\frac32;\frac1{\sqrt{2}}\right)\,\Gamma\left(\frac\alpha2+1\right)}{2^{\alpha/4}\,\sqrt{2\pi}}.\tag4$$ Ahora tomando un derivado con respecto al parámetro de $\alpha$ y la sustitución de $\alpha=3$ me dan el resultado deseado: $$\mathcal{A}=\frac{\sqrt{7+5\sqrt2}}{48\sqrt2}\left(\left(3-\sqrt2\right)\cdot(16-15\log2-6\gamma)+9\sqrt{10-7\sqrt2}\cdot\mathcal{D}\right),\tag5$$ donde $$\mathcal{D}=\left.\partial_\beta\Bigg({_2F_1}\left(\tfrac12,\beta;\tfrac32;\tfrac1{\sqrt2}\right)\Bigg)\right|_{\beta=5/2}\tag6$$ El problema es que no sé cómo evaluar el $\mathcal{D}$ plazo — cómo tomar un derivado de la función hipergeométrica con respecto al parámetro de $\beta$. Podría usted por favor me ayude con esta parte o sugerir otro enfoque que podría dar un formato cerrado resultado de $\mathcal{A}$ (si es posible)?