Una integral relacionada con la función digamma

Question

Estaba jugando con asintóticas y fórmulas integrales de la función digamma y la integral exponencial cuando me topé con esta:

$$\psi(x)=\int_0^\infty\left(\frac{e^{-t}}t-\frac{e^{-xt}}{1-e^{-t}}\right)~\mathrm dt,\quad\Re(x)>0$$

Trabajando un poco con ella, terminé con la siguiente integral:

$$I(a)=\int_1^\infty\frac{\ln(u-1)}{u^a}~\mathrm du,\quad\Re(a)>1$$

Y me preguntaba cómo evaluar esto. Al dejar que $u\mapsto u+1$,

$$I(a)=\int_0^\infty\frac{\ln(u)}{(u+1)^a}~\mathrm du$$

Al dejar que $u\mapsto e^u$,

$$I(a)=\int_{-\infty}^\infty\frac{ue^u}{(e^u+1)^a}~\mathrm du$$

Tal vez podamos aplicar IBP, pero eso parece ser complicado e involucra límites. Del mismo modo, los límites y el integrando no parecen ser muy propicios para una expansión en serie.

Preferiblemente, me gustaría resolver esta integral sin usar la función digamma.

Answer 1

$\newcommand{\bbx}[1]{\,\bbox[15px,border:1px groove navy]{\displaystyle{#1}}\,} \newcommand{\braces}[1]{\left\lbrace\,{#1}\,\right\rbrace} \newcommand{\bracks}[1]{\left\lbrack\,{#1}\,\right\rbrack} \newcommand{\dd}{\mathrm{d}} \newcommand{\ds}[1]{\displaystyle{#1}} \newcommand{\expo}[1]{\,\mathrm{e}^{#1}\,} \newcommand{\ic}{\mathrm{i}} \newcommand{\mc}[1]{\mathcal{#1}} \newcommand{\mrm}[1]{\mathrm{#1}} \newcommand{\pars}[1]{\left(\,{#1}\,\right)} \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}} \newcommand{\root}[2][]{\,\sqrt[#1]{\,{#2}\,}\,} \newcommand{\totald}[3][]{\frac{\mathrm{d}^{#1} #2}{\mathrm{d} #3^{#1}}} \newcommand{\verts}[1]{\left\vert\,{#1}\,\right\vert}$ \begin{align} &\left.\int_{0}^{\infty}\pars{{\expo{-t} \over t} - {\expo{-xt} \over 1 - \expo{-t}}}\,\dd t\,\right\vert_{\ \Re\pars{x}\ >\ 0} \,\,\,\stackrel{\exp\pars{-t}\ \mapsto\ t}{=} \int_{1}^{0}\bracks{{t \over -\ln\pars{t}} - {t^{x} \over 1 - t}}\, \pars{-\,{\dd t \over t}} \\[5mm] = &\ \int_{0}^{1}\bracks{{1 \over 1 - t}\int_{0}^{1}t^{y}\,\dd y - {t^{x - 1} \over 1 - t}}\,\dd t = \int_{0}^{1}\bracks{\int_{0}^{1}{1 - t^{x - 1} \over 1 - t}\,\dd t - \int_{0}^{1}{1 - t^{y} \over 1 - t}\,\dd t}\,\dd y \\[5mm] = &\ \int_{0}^{1}\bracks{\Psi\pars{x} - \Psi\pars{y + 1}}\,\dd y = \Psi\pars{x} - \ln\pars{\Gamma\pars{2} \over \Gamma\pars{1}} = \bbx{\Psi\pars{x}} \end{align}

Answer 2

No está claro lo que querías decir con "no usando $\psi(x)$".

Para $\Re(x) > 0$ $$\int_0^\infty\left(\frac{e^{-t}}t-\frac{e^{-xt}}{1-e^{-t}}\right)dt = \int_0^\infty e^{-t}\left(\frac{1}t-\frac{1}{1-e^{-t}}\right)dt+ \int_0^\infty\left(\frac{e^{-t}-e^{-xt}}{1-e^{-t}}\right)dt$$ $$ = C+\sum_{n=0}^\infty \int_0^\infty (e^{-(n+1)t}-e^{-(n+x) t})dt=C+\sum_{n=0}^\infty \left(\frac{1}{n+1}-\frac{1}{n+x}\right) = C+\gamma + \psi(x)$$ $C=- \gamma$ por ver a $x=1$

---

$$\int_1^\infty \ln(u-1)u^{-x}du = \frac{D}{x-1}+\int_1^\infty \ln(u-1)(u^{-x}-\frac{u^{-2}}{x-1})du \\= \frac{1}{x-1}(D+\int_1^\infty \frac{u^{1-x}-u^{-1}}{u-1} du)=\frac{1}{x-1}(D+\int_0^\infty \frac{e^{(1-x)t}-e^{-t}}{1-e^{-t}} dt )= \frac{\psi(x-1)+D}{x-1}$$

Answer 3

Deje $u\mapsto1/u$ en la definición de $I(a)$ para obtener

$$I(a)=\int_0^1\ln\left(\frac1u-1\right)u^{a-2}~\mathrm du$$

Ahora observe que:

$$\ln\left(\frac1u-1\right)=\ln(1-u)-\ln(u)$$

Así, por la expansión en series,

$$\begin{align}\int_0^1\ln(1-u)u^{a-2}~\mathrm du&=-\sum_{k=1}^\infty\frac1k\int_0^1u^{a+k-2}~\mathrm du\\&=-\sum_{k=1}^\infty\frac1{k(a+k-1)}\\&=-\frac{\gamma+\psi(a)}{a-1}\end{align}$$

Y por integración por partes,

$$\int_0^1\ln(u)u^{a-2}~\mathrm du=-\frac1{(a-1)^2}$$

Entonces la integral dada es

$$I(a)=\frac1{(a-1)^2}-\frac{\gamma+\psi(a)}{a-1}$$