Mathematics · Livro 4 · Bachelor Year 2

Matemática universitária — Graduação 2

Matemática universitária — Graduação 2 · Bachelor Year 2

9Integração

O volume do primeiro ano de graduação construiu a integral num segmento. Este capítulo a estende a intervalos arbitrários (integrais impróprias, com toda a caixa de ferramentas de comparação) e depois estuda integrais que dependem de um parâmetro — continuidade e derivação sob o sinal de integral — movidas pelo teorema da convergência dominada, o único resultado deste capítulo aceito por confiança. A função Γ\Gamma serve de exemplo condutor e de portal para metade das funções especiais da matemática.

9.1 Integrais num intervalo arbitrário

Definição 9.1

Seja ff contínua por partes em [a, b)\intco{a}{b} (b∈Rb \in \R ou +∞+\infty), com continuidade por partes. A integral converge quando lim⁡x→b−∫axf\lim_{x \to b^-} \int_a^x f existe; escreve-se então ∫abf\int_a^b f para o limite. (Analogamente em (a, b]\intoc{a}{b}, e em (a, b)\intoo{a}{b} separando num ponto interior — a escolha não importa, por Chasles.) A integral converge absolutamente quando ∫ab∣f∣\int_a^b \abs f converge; a convergência absoluta implica a convergência, pelo critério de Cauchy:

∣∫xyf∣≤∫xy∣f∣\Bigl| \int_x^{y} f \Bigr| \leq \int_x^{y} \abs f

e pela completude de R\R (a primitiva tem a propriedade de Cauchy). Em detalhe: sejam F(x)=∫axfF(x) = \int_a^x f e G(x)=∫ax∣f∣G(x) = \int_a^x \abs f. Se ∫b∣f∣\int^b\abs f converge, GG tem limite em b−b^-, de modo que, para todo ε>0\varepsilon > 0, existe c<bc < b com G(y)−G(x)≤εG(y) - G(x) \leq \varepsilon sempre que c≤x≤y<bc \leq x \leq y < b; a fórmula exibida transfere essa propriedade de Cauchy a FF. Para qualquer sequência xn→b−x_n \to b^-, os valores F(xn)F(x_n) formam então uma sequência de Cauchy de reais, convergente pela completude, e intercalar duas dessas sequências mostra que o limite é o mesmo para todas: FF tem limite em b−b^-.

Teorema 9.2 (Caixa de ferramentas de comparação positiva)

Para f,g≥0f, g \geq 0 contínua por partes em [a, b)\intco{a}{b}, com continuidade por partes:

  1. ∫abf\int_a^b f converge se e somente se a primitiva x↦∫axfx \mapsto \int_a^x f é limitada;
  2. f≤gf \leq g: a convergência de ∫g\int g força a de ∫f\int f; a divergência se transfere no outro sentido;
  3. f∼gf \sim g em bb: as duas integrais têm a mesma natureza;
  4. as escalas de referência: em +∞+\infty, ∫∞ ⁣dttα\int^{\infty} \frac{\dd t}{t^\alpha} converge se e somente se α>1\alpha > 1, e ∫∞ ⁣dtt(ln⁡t)β\int^\infty \frac{\dd t}{t(\ln t)^\beta} se e somente se β>1\beta > 1; numa extremidade finita bb, ∫b ⁣dt(b−t)α\int^b \frac{\dd t}{(b - t)^\alpha} converge se e somente se α<1\alpha < 1.

Demonstração. (1) A primitiva F(x)=∫axfF(x) = \int_a^x f é não decrescente (f≥0f \geq 0). Se ela é limitada, ℓ=sup⁡x<bF\ell = \sup_{x < b}F é finito e F(x)→ℓF(x) \to \ell: dado ε>0\varepsilon > 0, algum F(x0)>ℓ−εF(x_0) > \ell - \varepsilon, e a monotonicidade prende F(x)∈(ℓ−ε, ℓ]F(x) \in \intoc{\ell - \varepsilon}{\ell} para x0≤x<bx_0 \leq x < b. Se ela é ilimitada, F→+∞F \to +\infty: divergência.

(2) De f≤gf \leq g: ∫axf≤∫axg\int_a^x f \leq \int_a^x g para todo xx; se ∫bg\int^b g converge, o membro da direita é limitado, logo o da esquerda também, e (1) conclui. A contraposição transfere a divergência no outro sentido.

(3) f∼gf \sim g em bb fornece c<bc < b com

12 g(t)  ≤  f(t)  ≤  2 g(t)(c≤t<b):\tfrac12\,g(t) \;\leq\; f(t) \;\leq\; 2\,g(t) \qquad (c \leq t < b) :

por (2) aplicado nos dois sentidos em [c, b)\intco{c}{b}, as duas integrais têm a mesma natureza; a peça inicial [a, c]\intcc{a}{c} é uma integral própria e nada muda.

(4) Primitivas explícitas: para α≠1\alpha \neq 1 e β≠1\beta \neq 1,

∫cx ⁣dttα=x1−α−c1−α1−α,∫cx ⁣dtt(ln⁡t)β=(ln⁡x)1−β−(ln⁡c)1−β1−β,\int_c^x \frac{\dd t}{t^\alpha} = \frac{x^{1-\alpha} - c^{1-\alpha}}{1 - \alpha}, \qquad \int_c^x \frac{\dd t}{t(\ln t)^\beta} = \frac{(\ln x)^{1-\beta} - (\ln c)^{1-\beta}}{1 - \beta},

com logaritmos nos casos excluídos: limitadas quando x→+∞x \to +\infty exatamente se α>1\alpha > 1, resp. β>1\beta > 1. Numa extremidade finita, a substituição u=b−tu = b - t reduz à escala ∫0u−α  ⁣du\int_0 u^{-\alpha}\,\dd u, limitada se e somente se α<1\alpha < 1. Aplique (1) em cada caso. ∎

Exemplo 9.3 (Dois aquecimentos, levados até o fim)

(a) ∫01ln⁡t  ⁣dt\displaystyle\int_0^1 \ln t\,\dd t: o integrando explode em 0+0^+, mas ∣ln⁡t∣=o(t−1/2)\abs{\ln t} = o\bigl(t^{-1/2}\bigr) aí (os logaritmos perdem para as potências), e ∫0t−1/2\int_0 t^{-1/2} converge: convergência absoluta. O valor, por partes em [ε, 1]\intcc{\varepsilon}{1}:

∫ε1ln⁡t  ⁣dt=[tln⁡t−t]ε1=−1−εln⁡ε+ε→ε→0+−1.\int_\varepsilon^1 \ln t\,\dd t = \bigl[t\ln t - t\bigr]_\varepsilon^1 = -1 - \varepsilon\ln\varepsilon + \varepsilon \xrightarrow[\varepsilon\to0^+]{} -1 .

(b) ∫0∞ln⁡t1+t2  ⁣dt\displaystyle\int_0^\infty \frac{\ln t}{1 + t^2}\,\dd t: problema nas duas pontas, logo separe em 11. Perto de 00: ∣ln⁡t∣\abs{\ln t}, integrável como em (a); perto de ∞\infty: ln⁡t1+t2=o(t−3/2)\frac{\ln t}{1+t^2} = o(t^{-3/2}): absolutamente convergente. A substituição t=1ut = \frac1u leva (0, 1)\intoo{0}{1} em (1, ∞)\intoo{1}{\infty} e

∫01ln⁡t1+t2  ⁣dt=∫1∞−ln⁡u1+u−2⋅ ⁣duu2=−∫1∞ln⁡u1+u2  ⁣du:\int_0^1 \frac{\ln t}{1+t^2}\,\dd t = \int_1^{\infty} \frac{-\ln u}{1 + u^{-2}}\cdot \frac{\dd u}{u^2} = -\int_1^\infty \frac{\ln u}{1+u^2}\,\dd u :

as duas metades se cancelam, e a integral vale 00. Lição final: a simetria sob t↦1tt \mapsto \frac1t vale uma página de cálculo — o mesmo truque já movia o Exercício 9.3.

Exemplo 9.4 (Um valor, três integrais)

Estude I=∫0∞1−cos⁡tt2  ⁣dtI = \displaystyle\int_0^{\infty} \frac{1 - \cos t}{t^2}\,\dd t. Em 00: 1−cos⁡t∼t221 - \cos t \sim \frac{t^2}2, de modo que o integrando se estende continuamente pelo valor 12\frac12 — não há singularidade alguma. Em ∞\infty: 0≤1−cos⁡tt2≤2t20 \leq \frac{1 - \cos t}{t^2} \leq \frac{2}{t^2}: convergência absoluta (Teorema 9.2). Valor: integre por partes em [ε, M]\intcc{\varepsilon}{M} com u=1−cos⁡tu = 1 - \cos t, v′=t−2v' = t^{-2}:

∫εM1−cos⁡tt2  ⁣dt=[−1−cos⁡tt]εM+∫εMsin⁡tt  ⁣dt.\int_\varepsilon^M \frac{1 - \cos t}{t^2}\,\dd t = \Bigl[-\frac{1 - \cos t}{t}\Bigr]_\varepsilon^M + \int_\varepsilon^M \frac{\sin t}{t}\,\dd t .

O colchete se anula nas duas pontas (1−cos⁡εε∼ε2\frac{1 - \cos\varepsilon}{\varepsilon} \sim \frac\varepsilon2; numerador limitado em MM), e a integral tende ao valor de Dirichlet π2\frac\pi2 (Exercício 9.10): I=π2I = \frac\pi2. Lição final: com 1−cos⁡t=2sin⁡2t21 - \cos t = 2\sin^2\frac t2 e u=t2u = \frac t2,

I=∫0∞2sin⁡2u(2u)2  2  ⁣du=∫0∞(sin⁡uu) ⁣2 ⁣du:I = \int_0^\infty \frac{2\sin^2 u}{(2u)^2}\;2\,\dd u = \int_0^\infty \Bigl(\frac{\sin u}{u}\Bigr)^{\!2}\dd u :

os três clássicos ∫0∞sin⁡tt ⁣dt\int_0^\infty\frac{\sin t}{t}\dd t, ∫0∞(sin⁡tt)2 ⁣dt\int_0^\infty\bigl(\frac{\sin t}{t}\bigr)^2\dd t (Exercício 9.11) e II compartilham o valor π2\frac\pi2, passado adiante por partes e substituição — e só o primeiro é semiconvergente: a integração por partes trocou a convergência absoluta por um integrando mais simples.

Exemplo 9.5 (Uma integral semiconvergente)

∫1∞sin⁡tt  ⁣dt\displaystyle\int_1^{\infty} \frac{\sin t}{t}\,\dd t converge: integre por partes,

∫1xsin⁡tt ⁣dt=[−cos⁡tt]1x−∫1xcos⁡tt2 ⁣dt,\int_1^x \frac{\sin t}{t}\dd t = \Bigl[\frac{-\cos t}{t}\Bigr]_1^x - \int_1^x \frac{\cos t}{t^2}\dd t ,

em que o colchete tem limite e a última integral converge absolutamente (∣cos⁡t∣/t2≤t−2\abs{\cos t}/t^2 \leq t^{-2}). Mas não absolutamente: de ∣sin⁡t∣≥sin⁡2t\abs{\sin t} \geq \sin^2 t,

∫1x∣sin⁡t∣t  ⁣dt  ≥  ∫1xsin⁡2tt  ⁣dt=∫1x ⁣dt2t⏟= 12ln⁡x → ∞  −  ∫1xcos⁡2t2t  ⁣dt⏟convergente,\int_1^x \frac{\abs{\sin t}}{t}\,\dd t \;\geq\; \int_1^x \frac{\sin^2t}{t}\,\dd t = \underbrace{\int_1^x \frac{\dd t}{2t}}_{=\ \frac12\ln x \ \to\ \infty} \;-\; \underbrace{\int_1^x \frac{\cos 2t}{2t}\,\dd t}_{\text{convergente}} ,

em que a última integral converge pela mesma integração por partes acima (com sin⁡2t\sin 2t no colchete): uma peça divergente menos uma convergente diverge. Assim ∫1∞sin⁡tt ⁣dt\int_1^\infty\frac{\sin t}{t}\dd t converge sem convergir absolutamente — o análogo integral das séries alternadas, com a integração por partes fazendo o papel do critério das alternadas.

9.2 O teorema de convergência

Teorema 9.6 (Convergência dominada)

Sejam (fn)(f_n) contínuas por partes num intervalo II, com continuidade por partes, convergindo pontualmente para uma ff contínua por partes, com continuidade por partes, e suponha que exista uma função integrável fixa φ≥0\varphi \geq 0 (∫Iφ<∞\int_I \varphi < \infty) com

∣fn(t)∣≤φ(t)(t∈I, n∈N).\abs{f_n(t)} \leq \varphi(t) \qquad (t \in I,\ n \in \N).

Então todas as ∫Ifn\int_I f_n e ∫If\int_I f convergem absolutamente, e

∫Ifn→n→∞∫If.\int_I f_n \xrightarrow[n \to \infty]{} \int_I f .

Demonstração. Admitido neste nível. ∎

Observação 9.7

A demonstração honesta pertence à teoria da integração de Lebesgue do terceiro ano; o enunciado, porém, é usado constantemente daqui em diante. A hipótese de dominação é o ponto todo: a convergência pontual sozinha não basta (fn=n 1(0, 1/n)f_n = n\,\mathbf{1}_{\intoo{0}{1/n}}, corcovas deslizantes: ∫fn=1↛0=∫f\int f_n = 1 \not\to 0 = \int f). O teorema também vale para um parâmetro contínuo (fλf_\lambda, λ→λ0\lambda \to \lambda_0), pela caracterização sequencial dos limites.

Exemplo 9.8 (Um limite gaussiano, por dominação)

Calcule lim⁡n→∞In\displaystyle\lim_{n\to\infty} I_n, em que In=∫0∞(1+t2n) ⁣−n ⁣dtI_n = \int_0^\infty \Bigl(1 + \frac{t^2}{n}\Bigr)^{\!-n}\dd t. Pontualmente, (1+t2/n)n→et2(1 + t^2/n)^n \to \eu^{t^2} (o limite dos juros compostos), de modo que os integrandos tendem a e−t2\eu^{-t^2}. Dominação: a sequência n↦(1+u/n)nn \mapsto (1 + u/n)^n é não decrescente para u≥0u \geq 0 (MA–MG nos n+1n + 1 fatores 1,1+un,…,1+un1, 1 + \frac un, \dots, 1 + \frac un dá (1+un+1)n+1≥(1+un)n(1 + \frac u{n+1})^{n+1} \geq (1 + \frac un)^n), logo, para n≥2n \geq 2:

(1+t2n) ⁣−n≤(1+t22) ⁣−2,\Bigl(1 + \frac{t^2}{n}\Bigr)^{\!-n} \leq \Bigl(1 + \frac{t^2}{2}\Bigr)^{\!-2},

um dominante integrável (∼4t−4\sim 4t^{-4} no infinito). Convergência dominada:

In→n→∞∫0∞e−t2 ⁣dt=π2I_n \xrightarrow[n\to\infty]{} \int_0^\infty \eu^{-t^2}\dd t = \frac{\sqrt\pi}{2}

(a integral gaussiana de Exercício 9.8). Verificação final: a substituição t=ntan⁡θt = \sqrt n\tan\theta calcula InI_n exatamente, In=n∫0π/2cos⁡2n−2θ  ⁣dθ=n W2n−2I_n = \sqrt n\int_0^{\pi/2}\cos^{2n-2}\theta\,\dd\theta = \sqrt n\,W_{2n-2}, e a assintótica de Wallis Wm∼π/(2m)W_m \sim \sqrt{\pi/(2m)} (Lema 6.11) dá n W2n−2→π2\sqrt n\,W_{2n-2} \to \frac{\sqrt\pi}2 de novo: os dois pilares deste capítulo e do anterior concordam.

Exemplo 9.9 (Convergência dominada, parâmetro contínuo)

Calcule

lim⁡x→+∞∫0∞arctan⁡(xt)1+t2  ⁣dt.\lim_{x\to+\infty}\int_0^\infty \frac{\arctan(xt)}{1+t^2}\,\dd t .

Para cada t>0t > 0, arctan⁡(xt)→π2\arctan(xt) \to \frac\pi2 quando x→∞x \to \infty; e a dominação

∣arctan⁡(xt)1+t2∣≤π/21+t2,integraˊvel, independente de x,\Bigl|\frac{\arctan(xt)}{1+t^2}\Bigr| \leq \frac{\pi/2}{1+t^2}, \qquad\text{integrável, independente de } x,

vale para todo xx. Pela forma com parâmetro contínuo do Teorema 9.6 (caracterização sequencial: teste ao longo de toda xn→∞x_n \to \infty),

∫0∞arctan⁡(xt)1+t2  ⁣dt→x→+∞π2∫0∞ ⁣dt1+t2=π24.\int_0^\infty\frac{\arctan(xt)}{1+t^2}\,\dd t \xrightarrow[x\to+\infty]{} \frac\pi2\int_0^\infty\frac{\dd t}{1+t^2} = \frac{\pi^2}{4} .

Lição final: o único ponto t=0t = 0, onde o limite pontual vale 00 em vez de π2\frac\pi2, nada muda — a função limite só entra por sua integral, uma das misericórdias silenciosas do teorema.

9.3 Integrais com parâmetro

Teorema 9.10 (Continuidade sob o sinal de integral)

Seja f ⁣:A×I→Rf \colon A \times I \to \R (AA um espaço métrico, II um intervalo) com: t↦f(x,t)t \mapsto f(x, t) contínua por partes para cada xx, com continuidade por partes; x↦f(x,t)x \mapsto f(x, t) contínua para cada tt; e uma dominação ∣f(x,t)∣≤φ(t)\abs{f(x,t)} \leq \varphi(t) (φ\varphi integrável em II, independente de xx). Então

F(x)=∫If(x,t)  ⁣dtF(x) = \int_I f(x, t)\,\dd t

está definida e é contínua em AA.

Demonstração. Boa definição: a dominação dá a convergência absoluta. Continuidade em x0x_0: para toda sequência xn→x0x_n \to x_0, as funções gn(t)=f(xn,t)g_n(t) = f(x_n, t) convergem pontualmente para f(x0,t)f(x_0, t) (continuidade em xx) sob a dominação fixa φ\varphi: a convergência dominada dá F(xn)→F(x0)F(x_n) \to F(x_0); conclua pela caracterização sequencial da continuidade (Definição 4.5). ∎

Teorema 9.11 (Derivação sob o sinal de integral)

Seja f ⁣:J×I→Rf \colon J \times I \to \R (JJ um intervalo de parâmetros) com: t↦f(x,t)t \mapsto f(x,t) integrável em II para cada xx; x↦f(x,t)x \mapsto f(x,t) de classe C1C^1 para cada tt, sendo a derivada parcial ∂f∂x\frac{\partial f}{\partial x} contínua por partes em tt, com continuidade por partes, e dominada: ∣∂f∂x(x,t)∣≤ψ(t)\bigl|\frac{\partial f}{\partial x}(x,t)\bigr| \leq \psi(t) com ψ\psi integrável. Então F(x)=∫If(x,t) ⁣dtF(x) = \int_I f(x,t)\dd t é C1C^1 em JJ e

F′(x)=∫I∂f∂x(x,t)  ⁣dt.F'(x) = \int_I \frac{\partial f}{\partial x}(x, t)\,\dd t .

Demonstração. Fixe xx e hn→0h_n \to 0. Os quocientes de diferenças

F(x+hn)−F(x)hn=∫If(x+hn,t)−f(x,t)hn  ⁣dt\frac{F(x + h_n) - F(x)}{h_n} = \int_I \frac{f(x + h_n, t) - f(x, t)}{h_n}\,\dd t

têm integrandos convergindo pontualmente para ∂f∂x(x,t)\frac{\partial f}{\partial x}(x, t) e dominados por ψ(t)\psi(t): pela desigualdade do valor médio aplicada em xx com tt fixo,

∣f(x+hn,t)−f(x,t)hn∣≤sup⁡ξ∣∂f∂x(ξ,t)∣≤ψ(t).\Bigl|\frac{f(x + h_n, t) - f(x,t)}{h_n}\Bigr| \leq \sup_{\xi} \Bigl|\frac{\partial f}{\partial x}(\xi, t)\Bigr| \leq \psi(t) .

A convergência dominada dá o limite ∫I∂f∂x(x,t) ⁣dt\int_I \frac{\partial f}{\partial x}(x,t)\dd t dos quocientes: FF é derivável com a derivada anunciada, que é contínua pelo Teorema 9.10 aplicado a ∂f∂x\frac{\partial f}{\partial x}. ∎

Exemplo 9.12 (Uma integral com parâmetro conferida contra uma fórmula)

Seja F(x)=∫0∞ ⁣dtt2+xF(x) = \displaystyle\int_0^\infty \frac{\dd t}{t^2 + x} para x>0x > 0. Em todo [a, b]⊂(0, ∞)\intcc{a}{b} \subset \intoo{0}{\infty}, o integrando é dominado por 1t2+a\frac{1}{t^2 + a}, integrável e independente de xx: FF é contínua (Teorema 9.10). Aqui o teorema pode ser conferido contra um valor explícito:

F(x)=[1xarctan⁡tx]0∞=π2x,F(x) = \Bigl[\frac{1}{\sqrt x}\arctan\frac{t}{\sqrt x}\Bigr]_0^\infty = \frac{\pi}{2\sqrt x} ,

visivelmente contínua. Derive agora sob a integral: a derivada em xx, −1(t2+x)2-\frac{1}{(t^2+x)^2}, é dominada em [a, b]\intcc ab por 1(t2+a)2\frac{1}{(t^2+a)^2}, integrável: o Teorema 9.11 dá

F′(x)=−∫0∞ ⁣dt(t2+x)2enquantoF′(x)=−π4 x−3/2,F'(x) = -\int_0^\infty \frac{\dd t}{(t^2 + x)^2} \qquad\text{enquanto}\qquad F'(x) = -\frac{\pi}{4}\,x^{-3/2} ,

de modo que calculamos uma integral nova de graça: ∫0∞ ⁣dt(t2+x)2=π4x3/2\int_0^\infty\frac{\dd t}{(t^2+x)^2} = \frac{\pi}{4x^{3/2}}. Lição final: derivar uma integral com parâmetro conhecida é uma fábrica de fórmulas novas — iterando obtém-se ∫0∞ ⁣dt(t2+1)n\int_0^\infty\frac{\dd t}{(t^2+1)^n} para todo nn, sem nenhuma substituição trigonométrica.

Método 9.13 (Estudar uma integral imprópria)

Dada ∫abf\int_a^b f:

  1. Localize o problema: liste as extremidades (ou pontos interiores) em que ff é ilimitada ou o intervalo é infinito, e separe de modo que cada peça tenha exatamente uma ponta problemática.
  2. Se ff tem sinal constante perto dessa ponta, ache um equivalente e compare com as escalas de referência do Teorema 9.2.
  3. Se ff oscila, teste ∣f∣\abs f primeiro (convergência absoluta). Se ∫∣f∣\int\abs f diverge, tente a integração por partes para trocar a oscilação por decaimento, como no Exemplo 9.5; minorações do tipo ∣sin⁡t∣≥sin⁡2t\abs{\sin t} \geq \sin^2t detectam a semiconvergência genuína.
  4. Para um valor, e não apenas a natureza: partes, substituição ou um parâmetro (derive uma integral mais simples, como no Exemplo 9.12 e no Exemplo 9.21).
  5. Verificações de bom senso em qualquer valor calculado: sinal e ordem de grandeza contra uma cota bruta (∫0∞e−t2 ⁣dt∈(0, 1+∫1∞e−t)\int_0^\infty \eu^{-t^2}\dd t \in \intoo{0}{1 + \int_1^\infty \eu^{-t}}, logo π2≈0.886\frac{\sqrt\pi}{2} \approx 0.886 é plausível); e coerência dimensional por reescalamento (t↦λtt \mapsto \lambda t deve reescalar os dois lados do mesmo jeito — o detector mais rápido de um fator perdido).

Observação 9.14 (Armadilhas comuns)

Três erros recorrentes. (i) Dominantes que dependem do parâmetro: a dominação ∣f(x,t)∣≤φ(t)\abs{f(x,t)} \leq \varphi(t) tem de ser uniforme em xx no conjunto considerado; ela costuma valer em segmentos [a, b]\intcc ab mas não globalmente — para ∫0∞e−xt ⁣dt\int_0^\infty\eu^{-xt}\dd t não há dominante integrável válido para todo x>0x > 0, e no entanto dominar em x≥a>0x \geq a > 0 basta para trabalhar em toda a semirreta aberta, pois continuidade e derivadas são noções locais. (ii) Comparar integrandos com sinal: a caixa de ferramentas de comparação é para funções não negativas; de ∣f∣≤g\abs f \leq g com ∫g\int g divergente nada se pode concluir — ∫1∞sin⁡tt  ⁣dt\int_1^\infty\frac{\sin t}t\,\dd t converge embora toda comparação com 1t\frac1t falhe. (iii) Esquecer metade do problema: em (0, ∞)\intoo{0}{\infty} estude sempre as duas pontas separadamente; ∫0∞ ⁣dtt\int_0^\infty\frac{\dd t}{t} diverge em ambas, e uma separação de aparência convergente pode cancelar silenciosamente dois infinitos. O reflexo seguro é a lista de verificação do Método 9.13.

Exemplo 9.15 (Um caso de fronteira de Bertrand, até o algarismo)

A escala ∫∞ ⁣dtt(ln⁡t)β\int^\infty\frac{\dd t}{t(\ln t)^\beta} do Teorema 9.2 fica exatamente na borda das escalas de potências; seus casos de fronteira merecem um cálculo completo. Para β=2\beta = 2:

∫e∞ ⁣dtt(ln⁡t)2=[−1ln⁡t]e∞=0−(−1)=1,\int_\eu^{\infty}\frac{\dd t}{t(\ln t)^2} = \Bigl[-\frac{1}{\ln t}\Bigr]_\eu^{\infty} = 0 - (-1) = 1 ,

uma integral convergente de valor agradavelmente exato; ao passo que, para β=1\beta = 1,

∫ex ⁣dttln⁡t=[ln⁡ln⁡t]ex=ln⁡ln⁡x⟶∞,\int_\eu^{x}\frac{\dd t}{t\ln t} = \bigl[\ln\ln t\bigr]_\eu^{x} = \ln\ln x \longrightarrow \infty ,

divergente — mas tão devagar que atingir ln⁡ln⁡x=10\ln\ln x = 10 exige x=ee10≈109566x = \eu^{\eu^{10}} \approx 10^{9566}. Lição final: entre “toda potência t−1−εt^{-1-\varepsilon} converge” e “t−1t^{-1} diverge” mora uma escada infinita de escalas logarítmicas, cada uma refinando a anterior; a substituição u=ln⁡tu = \ln t colapsa cada degrau sobre o anterior, e é por isso que os critérios de Bertrand ecoam os de Riemann um nível acima.

Observação 9.16 (Perspectivas dentro deste volume)

As ferramentas deste capítulo estão prestes a aparecer em toda parte. A convergência dominada é o motor por trás das identidades aproximadas do capítulo seguinte (núcleos deslizantes, tanto Bernstein quanto Fejér); a continuidade e a derivação sob o sinal de integral produzem o cálculo dos coeficientes de Fourier no capítulo de Fourier, onde todo cn(f)c_n(f) é uma integral com parâmetro disfarçada. A função Γ\Gamma volta duas vezes: no capítulo sobre integrais múltiplas, em que uma integral dupla finalmente demonstra a fórmula Beta–Gama de Euler por inteiro, e nos capítulos de probabilidade, em que integrais do tipo Γ\Gamma normalizam as densidades padrão e calculam seus momentos. E a semiconvergente ∫sin⁡tt\int\frac{\sin t}{t} ressurge como a constante de Gibbs do capítulo de Fourier — a mesma integral, medindo a ultrapassagem das somas parciais num salto.

Definição 9.17 (A função Γ\Gamma)

Para x>0x > 0:

Γ(x)=∫0∞tx−1 e−t  ⁣dt,\Gamma(x) = \int_0^{\infty} t^{x-1}\,\eu^{-t}\,\dd t ,

convergente nas duas pontas (tx−1t^{x-1} integrável em 0+0^+ para x>0x > 0; decaimento exponencial em ∞\infty).

Teorema 9.18

Γ\Gamma é contínua em (0, +∞)\intoo{0}{+\infty}, satisfaz a equação funcional

Γ(x+1)=x Γ(x),Γ(1)=1,logoΓ(n+1)=n!,\Gamma(x + 1) = x\,\Gamma(x), \qquad \Gamma(1) = 1, \qquad\text{logo}\qquad \Gamma(n + 1) = n! ,

e é de classe C1C^1 (na verdade C∞C^\infty) com Γ′(x)=∫0∞tx−1e−tln⁡t  ⁣dt\Gamma'(x) = \int_0^\infty t^{x-1}\eu^{-t}\ln t\,\dd t.

Demonstração. Equação funcional: integre por partes em [ε, M]\intcc{\varepsilon}{M} e faça as pontas irem ao limite: ∫txe−t=[−txe−t]+x∫tx−1e−t\int t^{x}\eu^{-t} = [-t^x\eu^{-t}] + x\int t^{x-1}\eu^{-t}, anulando-se os termos de bordo — de fato, εxe−ε→0\varepsilon^x\eu^{-\varepsilon} \to 0 quando ε→0+\varepsilon \to 0^+, pois x>0x > 0, e Mxe−M→0M^x\eu^{-M} \to 0 quando M→∞M \to \infty, pois a exponencial vence toda potência; as duas integrais truncadas convergem para seus valores impróprios pela convergência estabelecida na Definição 9.17. Γ(1)=∫e−t=1\Gamma(1) = \int \eu^{-t} = 1; a indução dá o fatorial.

Continuidade em [a, b]⊂(0, ∞)\intcc{a}{b} \subset \intoo{0}{\infty}: domine tx−1e−tt^{x-1}\eu^{-t} por φ(t)=(ta−1+tb−1)e−t\varphi(t) = (t^{a-1} + t^{b-1})\eu^{-t}, integrável e independente de x∈[a, b]x \in \intcc{a}{b}: o Teorema 9.10 se aplica em todo segmento desse tipo, logo em toda a semirreta. Derivabilidade: a derivada em xx, tx−1e−tln⁡tt^{x-1}\eu^{-t}\ln t, é dominada em [a, b]\intcc{a}{b} por (ta−1+tb−1)e−t ∣ln⁡t∣(t^{a-1} + t^{b-1})\eu^{-t}\,\abs{\ln t}, ainda integrável: Teorema 9.11; iterando obtêm-se todas as derivadas (cada uma acrescenta uma potência de ln⁡t\ln t, inofensiva). ∎

Exemplo 9.19 (Fatoriais semi-inteiros)

A equação funcional e Γ(12)=π\Gamma\bigl(\frac12\bigr) = \sqrt\pi (a uma substituição de Exercício 9.8: ponha t=u2t = u^2 na integral que define) geram todos os valores semi-inteiros:

Γ(32)=12 Γ(12)=π2,Γ(52)=32⋅π2=3π4,Γ(72)=15π8.\Gamma\Bigl(\frac32\Bigr) = \frac12\,\Gamma\Bigl(\frac12\Bigr) = \frac{\sqrt\pi}{2}, \qquad \Gamma\Bigl(\frac52\Bigr) = \frac32\cdot\frac{\sqrt\pi}{2} = \frac{3\sqrt\pi}{4}, \qquad \Gamma\Bigl(\frac72\Bigr) = \frac{15\sqrt\pi}{8} .

Como Γ(n+1)=n!\Gamma(n+1) = n!, é justo dizer “12!=π2≈0.886\frac12! = \frac{\sqrt\pi}{2} \approx 0.886”: o fatorial foi interpolado, e a curva interpoladora mergulha abaixo de 11 entre 0!=10! = 1 e 1!=11! = 1 (seu mínimo ≈0.8856\approx 0.8856 em x≈1.4616x \approx 1.4616 combina com o retrato de convexidade da Parte I do problema de fim de semana). Lição final: nada na integral ∫0∞tx−1e−t ⁣dt\int_0^\infty t^{x-1}\eu^{-t}\dd t privilegia os inteiros — a discretude do fatorial era um acidente da contagem, e π\sqrt\pi é o que vive entre 11 e 11.

Observação 9.20 (Para onde vai Γ\Gamma daqui)

O problema de fim de semana deste capítulo constrói todo o cálculo de Euler em torno de Γ\Gamma: a função Beta, suas recursões por integração por partes, as integrais de Wallis como valores de Beta e a fórmula-limite de Gauss. O capítulo sobre integrais múltiplas demonstra a fórmula Beta–Gama de Euler para todos os argumentos por uma integral dupla; os capítulos de probabilidade reencontram Γ\Gamma na normalização das densidades mais comuns e nos momentos dos tempos de espera. O volume do terceiro ano de graduação reconstrói Γ\Gamma sobre fundamentos de Lebesgue, demonstra o teorema de unicidade de Bohr–Mollerup e estende a fórmula de Stirling dos inteiros à semirreta real por convergência dominada.

Exemplo 9.21 (Um cálculo clássico por derivação)

Para x∈Rx \in \R, seja F(x)=∫0∞e−t2cos⁡(xt)  ⁣dtF(x) = \int_0^{\infty} \eu^{-t^2}\cos(xt)\,\dd t (absolutamente convergente, dominada por e−t2\eu^{-t^2}). Pelo Teorema 9.11 (dominação da derivada em xx por t e−t2t\,\eu^{-t^2}, integrável):

F′(x)=−∫0∞t e−t2sin⁡(xt)  ⁣dt=[e−t22sin⁡(xt)]0∞−x2∫0∞e−t2cos⁡(xt)  ⁣dt=−x2 F(x),F'(x) = -\int_0^\infty t\,\eu^{-t^2}\sin(xt)\,\dd t = \Bigl[\frac{\eu^{-t^2}}{2}\sin(xt)\Bigr]_0^\infty - \frac x2\int_0^\infty \eu^{-t^2}\cos(xt)\,\dd t = -\frac x2\,F(x),

(partes com u′=te−t2u' = t\eu^{-t^2}). A equação diferencial F′=−x2FF' = -\frac x2 F integra-se em F(x)=F(0) e−x2/4F(x) = F(0)\,\eu^{-x^2/4}: a integral do tipo gaussiano se reproduz. A constante F(0)=∫0∞e−t2 ⁣dt=π2F(0) = \int_0^\infty \eu^{-t^2}\dd t = \frac{\sqrt\pi}{2} é calculada no Exercício 9.8 — e de novo, por integração dupla, no Capítulo 20.

9.4 Exercícios

Exercício 9.1 ★

Natureza de: ∫01 ⁣dtt(1−t)\displaystyle\int_0^1 \frac{\dd t}{\sqrt{t(1-t)}};   ∫1∞ln⁡tt2 ⁣dt\;\displaystyle\int_1^\infty \frac{\ln t}{t^2}\dd t;   ∫0∞ ⁣dt1+t2sin⁡2t\;\displaystyle\int_0^\infty \frac{\dd t}{1 + t^2\sin^2 t} (compare com o comportamento divergente do tipo harmônico perto de t=nπt = n\pi).

Solução

Solução de Exercício 9.1.

∫01 ⁣dtt(1−t)\int_0^1 \frac{\dd t}{\sqrt{t(1-t)}}: perto de 00, ∼t−1/2\sim t^{-1/2} (α=12<1\alpha = \frac12 < 1: converge); perto de 11, ∼(1−t)−1/2\sim (1-t)^{-1/2}: converge. Convergente (seu valor é π\pi, pela substituição t=sin⁡2θt = \sin^2\theta).

∫1∞ln⁡tt2\int_1^\infty \frac{\ln t}{t^2}: ln⁡tt2=o(t−3/2)\frac{\ln t}{t^2} = o(t^{-3/2}): convergente (valor 11 por partes).

∫0∞ ⁣dt1+t2sin⁡2t\int_0^\infty \frac{\dd t}{1 + t^2\sin^2 t}: divergente. Perto de t=nπt = n\pi, escreva t=nπ+ut = n\pi + u: sin⁡2t=sin⁡2u≤u2\sin^2 t = \sin^2 u \leq u^2, logo, em ∣u∣≤1n\abs u \leq \frac{1}{n}, 1+t2sin⁡2t≤1+(nπ+1)2u2≤Cn2u2+11 + t^2\sin^2 t \leq 1 + (n\pi + 1)^2u^2 \leq C n^2 u^2 + 1; portanto

∫nπ−1/nnπ+1/n ⁣dt1+t2sin⁡2t≥∫−1/n1/n ⁣du1+Cn2u2=2arctan⁡CC⋅1n,\int_{n\pi - 1/n}^{n\pi + 1/n} \frac{\dd t}{1 + t^2\sin^2 t} \geq \int_{-1/n}^{1/n} \frac{\dd u}{1 + Cn^2u^2} = \frac{2\arctan\sqrt C}{\sqrt C}\cdot\frac{1}{n} ,

um termo de uma série divergente do tipo harmônico: somando sobre nn, a primitiva é ilimitada.

Exercício 9.2 ★

Calcule ∫0∞tne−λt  ⁣dt\displaystyle\int_0^\infty t^n \eu^{-\lambda t}\,\dd t (λ>0\lambda > 0) via Γ\Gamma, e ∫01(ln⁡t)n ⁣dt\displaystyle\int_0^1 (\ln t)^n \dd t pela substituição t=e−ut = \eu^{-u}.

Solução

Solução de Exercício 9.2.

Substitua u=λtu = \lambda t:

∫0∞tne−λt ⁣dt=1λn+1∫0∞une−u ⁣du=Γ(n+1)λn+1=n!λn+1.\int_0^\infty t^n \eu^{-\lambda t}\dd t = \frac{1}{\lambda^{n+1}}\int_0^\infty u^n\eu^{-u}\dd u = \frac{\Gamma(n+1)}{\lambda^{n+1}} = \frac{n!}{\lambda^{n+1}} .

Com t=e−ut = \eu^{-u} ( ⁣dt=−e−u ⁣du\dd t = -\eu^{-u}\dd u):

∫01(ln⁡t)n ⁣dt=∫0∞(−u)ne−u  ⁣du=(−1)n n!.\int_0^1 (\ln t)^n \dd t = \int_0^{\infty} (-u)^n \eu^{-u}\,\dd u = (-1)^n\, n! .

Exercício 9.3 ★

Prove que ∫0∞ ⁣dt(1+t2)(1+tx)\displaystyle\int_0^{\infty} \frac{\dd t}{(1 + t^2)(1 + t^x)} está bem definida para todo x∈Rx \in \R e é independente de xx. (Substitua t↦1tt \mapsto \frac1t e faça a média das duas expressões.) Qual é seu valor?

Solução

Solução de Exercício 9.3.

Convergência: o integrando é ≤11+t2\leq \frac{1}{1+t^2} perto de ∞\infty e limitado perto de 00 (os dois fatores são minorados longe de 00): absolutamente convergente, para todo xx. Substituindo t=1ut = \frac1u ( ⁣dt=− ⁣duu2\dd t = -\frac{\dd u}{u^2}):

I(x)=∫0∞1(1+1u2)(1+u−x)⋅ ⁣duu2=∫0∞ux(1+u2)(1+ux)  ⁣du.I(x) = \int_0^\infty \frac{1}{\bigl(1 + \frac1{u^2}\bigr)\bigl(1 + u^{-x}\bigr)}\cdot\frac{\dd u}{u^2} = \int_0^\infty \frac{u^x}{(1 + u^2)(1 + u^x)}\,\dd u .

Somando as duas expressões de I(x)I(x):

2I(x)=∫0∞1+tx(1+t2)(1+tx) ⁣dt=∫0∞ ⁣dt1+t2=π2:2I(x) = \int_0^\infty \frac{1 + t^x}{(1+t^2)(1+t^x)}\dd t = \int_0^\infty \frac{\dd t}{1 + t^2} = \frac{\pi}{2} :

I(x)=π4I(x) = \frac\pi4, independente de xx.

Exercício 9.4 ★★

(Integrais de Bertrand numa extremidade finita) Para quais (α,β)(\alpha, \beta) a integral ∫01/2 ⁣dttα ∣ln⁡t∣β\displaystyle\int_0^{1/2} \frac{\dd t}{t^\alpha\,\abs{\ln t}^\beta} converge?

Solução

Solução de Exercício 9.4.

Perto de 0+0^+, com u=∣ln⁡t∣→∞u = \abs{\ln t} \to \infty. Se α<1\alpha < 1: convergência qualquer que seja β\beta (compare com t−α′t^{-\alpha'} para α<α′<1\alpha < \alpha' < 1: o fator logarítmico é vencido). Se α>1\alpha > 1: divergência qualquer que seja β\beta (compare com t−α′′t^{-\alpha''}, 1<α′′<α1 < \alpha'' < \alpha). Se α=1\alpha = 1: substitua t=e−ut = \eu^{-u}:

∫01/2 ⁣dtt ∣ln⁡t∣β=∫ln⁡2∞ ⁣duuβ,\int_0^{1/2} \frac{\dd t}{t\,\abs{\ln t}^\beta} = \int_{\ln 2}^{\infty} \frac{\dd u}{u^\beta},

convergente se e somente se β>1\beta > 1. Resumo: convergência se e somente se α<1\alpha < 1, ou (α=1\alpha = 1 e β>1\beta > 1) — o espelho das séries de Bertrand.

Exercício 9.5 ★★

Seja F(x)=∫0∞e−xt1+t2  ⁣dtF(x) = \displaystyle\int_0^{\infty} \frac{\eu^{-xt}}{1 + t^2}\,\dd t para x≥0x \geq 0. Prove que FF é contínua em [0, ∞)\intco{0}{\infty}, de classe C2C^2 em (0, ∞)\intoo{0}{\infty}, satisfaz F′′+F=1xF'' + F = \frac1x aí, e que F(x)→0F(x) \to 0 quando x→+∞x \to +\infty.

Solução

Solução de Exercício 9.5.

Continuidade em [0, ∞)\intco{0}{\infty}: dominação ∣e−xt1+t2∣≤11+t2\bigl|\frac{\eu^{-xt}}{1+t^2}\bigr| \leq \frac{1}{1+t^2}, integrável, uniforme em x≥0x \geq 0: Teorema 9.10.

C2C^2 em (0, ∞)\intoo{0}{\infty}: em x≥a>0x \geq a > 0, as duas primeiras derivadas em xx, −t e−xt1+t2\frac{-t\,\eu^{-xt}}{1+t^2} e t2e−xt1+t2\frac{t^2\eu^{-xt}}{1+t^2}, são dominadas por t e−att\,\eu^{-at} e e−at\eu^{-at}: duas aplicações do Teorema 9.11. Então

F′′(x)+F(x)=∫0∞t2+11+t2 e−xt ⁣dt=∫0∞e−xt ⁣dt=1x.F''(x) + F(x) = \int_0^\infty \frac{t^2 + 1}{1 + t^2}\,\eu^{-xt}\dd t = \int_0^\infty \eu^{-xt}\dd t = \frac1x .

Limite: 0≤F(x)≤∫0∞e−xt ⁣dt=1x→00 \leq F(x) \leq \int_0^\infty \eu^{-xt}\dd t = \frac1x \to 0.

Exercício 9.6 ★★

(Frullani) Seja ff contínua em [0, +∞)\intco{0}{+\infty} com limite finito f(∞)f(\infty) em +∞+\infty. Prove que, para a,b>0a, b > 0:

∫0∞f(at)−f(bt)t  ⁣dt=(f(0)−f(∞)) ln⁡ba.\int_0^{\infty} \frac{f(at) - f(bt)}{t}\,\dd t = \bigl(f(0) - f(\infty)\bigr)\,\ln\frac ba .

(Em [ε, M]\intcc{\varepsilon}{M}, substitua em cada peça e reagrupe em ∫aεbε−∫aMbM\int_{a\varepsilon}^{b\varepsilon} - \int_{aM}^{bM} de f(u)u  ⁣du\frac{f(u)}u\,\dd u; confronte usando a continuidade em 00 e o limite em ∞\infty.) Calcule ∫0∞e−t−e−2tt ⁣dt\int_0^\infty \frac{\eu^{-t} - \eu^{-2t}}{t}\dd t.

Solução

Solução de Exercício 9.6.

Em [ε, M]\intcc{\varepsilon}{M}, substitua u=atu = at e u=btu = bt nas duas metades:

∫εMf(at)−f(bt)t ⁣dt=∫aεaMf(u)u ⁣du−∫bεbMf(u)u ⁣du=∫aεbεf(u)u ⁣du−∫aMbMf(u)u ⁣du.\int_\varepsilon^M \frac{f(at) - f(bt)}{t}\dd t = \int_{a\varepsilon}^{aM}\frac{f(u)}{u}\dd u - \int_{b\varepsilon}^{bM}\frac{f(u)}{u}\dd u = \int_{a\varepsilon}^{b\varepsilon} \frac{f(u)}{u}\dd u - \int_{aM}^{bM} \frac{f(u)}{u}\dd u .

Primeira peça: f(u)=f(0)+o(1)f(u) = f(0) + o(1) perto de 00, e ∫aεbε ⁣duu=ln⁡ba\int_{a\varepsilon} ^{b\varepsilon} \frac{\dd u}{u} = \ln\frac ba: a peça tende a f(0)ln⁡baf(0)\ln\frac ba. Segunda peça: f(u)→f(∞)f(u) \to f(\infty), mesmo cálculo: tende a f(∞)ln⁡baf(\infty)\ln\frac ba. Logo a integral imprópria converge para (f(0)−f(∞))ln⁡ba\bigl(f(0) - f(\infty)\bigr)\ln\frac ba.

Com f(t)=e−tf(t) = \eu^{-t} (f(0)=1f(0) = 1, f(∞)=0f(\infty) = 0), a=1a = 1, b=2b = 2:

∫0∞e−t−e−2tt ⁣dt=ln⁡2.\int_0^\infty \frac{\eu^{-t} - \eu^{-2t}}{t}\dd t = \ln 2 .

Exercício 9.7 ★★

Justifique e calcule lim⁡n→∞∫0n(1−tn) ⁣ntx−1  ⁣dt\lim_{n\to\infty} \displaystyle\int_0^n \Bigl(1 - \frac tn\Bigr)^{\!n} t^{x-1}\,\dd t para x>0x > 0 (convergência dominada com φ(t)=e−ttx−1\varphi(t) = \eu^{-t}t^{x-1}, usando (1−t/n)n≤e−t(1 - t/n)^n \leq \eu^{-t}; o limite é Γ(x)\Gamma(x)).

Solução

Solução de Exercício 9.7.

Estenda o integrando por 00 além de t=nt = n: gn(t)=(1−tn)ntx−11t≤ng_n(t) = (1 - \frac tn)^n t^{x-1}\mathbf{1}_{t \leq n}. Pontualmente, gn(t)→e−ttx−1g_n(t) \to \eu^{-t}t^{x-1} (o limite dos juros compostos, volume do primeiro ano de graduação). Dominação: ln⁡(1−u)≤−u\ln(1 - u) \leq -u dá (1−tn)n≤e−t(1 - \frac tn)^n \leq \eu^{-t} em [0, n]\intcc{0}{n}, logo ∣gn(t)∣≤e−ttx−1=φ(t)\abs{g_n(t)} \leq \eu^{-t}t^{x-1} = \varphi(t), integrável. Convergência dominada:

∫0n(1−tn)ntx−1 ⁣dt→n→∞∫0∞e−ttx−1 ⁣dt=Γ(x).\int_0^n \Bigl(1 - \frac tn\Bigr)^n t^{x-1}\dd t \xrightarrow[n\to\infty]{} \int_0^\infty \eu^{-t}t^{x-1}\dd t = \Gamma(x) .

(Calcular o membro da esquerda por partes repetidas dá a forma de produto de Euler Γ(x)=lim⁡n! nxx(x+1)⋯(x+n)\Gamma(x) = \lim \frac{n!\,n^x}{x(x+1)\cdots(x+n)}.)

Exercício 9.8 ★★★

(A integral gaussiana por um truque de parâmetro) Para x≥0x \geq 0 ponha

G(x)=(∫0xe−t2 ⁣dt) ⁣2,H(x)=∫01e−x2(1+t2)1+t2  ⁣dt.G(x) = \Bigl(\int_0^x \eu^{-t^2}\dd t\Bigr)^{\!2}, \qquad H(x) = \int_0^1 \frac{\eu^{-x^2(1+t^2)}}{1 + t^2}\,\dd t .

Prove que G′+H′=0G' + H' = 0 (derive HH sob a integral e substitua u=xtu = xt na integral resultante), deduza G(x)+H(x)=π4G(x) + H(x) = \frac\pi4 para todo xx e conclua

∫0∞e−t2  ⁣dt=π2.\int_0^{\infty} \eu^{-t^2}\,\dd t = \frac{\sqrt\pi}{2} .
Solução

Solução de Exercício 9.8.

HH é derivável em xx (integrando C1C^1 em xx, derivada −2x(1+t2)⋅e−x2(1+t2)1+t2=−2x e−x2e−x2t2-2x(1+t^2)\cdot\frac{\eu^{-x^2(1+t^2)}}{1+t^2} = -2x\,\eu^{-x^2}\eu^{-x^2t^2}, contínua e limitada nos compactos de xx, dominação sobre t∈[0, 1]t \in \intcc{0}{1} trivial):

H′(x)=−2x e−x2∫01e−x2t2  ⁣dt=u=xt−2 e−x2∫0xe−u2  ⁣du=−G′(x),H'(x) = -2x\,\eu^{-x^2}\int_0^1 \eu^{-x^2t^2}\,\dd t \overset{u = xt}{=} -2\,\eu^{-x^2}\int_0^x \eu^{-u^2}\,\dd u = -G'(x),

pois G′(x)=2e−x2∫0xe−t2 ⁣dtG'(x) = 2\eu^{-x^2}\int_0^x \eu^{-t^2}\dd t (regra da cadeia no quadrado, teorema fundamental do cálculo). Logo G+HG + H é constante, igual a G(0)+H(0)=0+∫01 ⁣dt1+t2=π4G(0) + H(0) = 0 + \int_0^1 \frac{\dd t}{1+t^2} = \frac\pi4.

Quando x→∞x \to \infty: 0≤H(x)≤e−x2∫01 ⁣dt→00 \leq H(x) \leq \eu^{-x^2}\int_0^1 \dd t \to 0, logo G(x)→π4G(x) \to \frac\pi4:

∫0∞e−t2 ⁣dt=π4=π2.\int_0^\infty \eu^{-t^2}\dd t = \sqrt{\frac\pi4} = \frac{\sqrt\pi}{2} .

(Por consequência, Γ(12)=2∫0∞e−t2 ⁣dt=π\Gamma\bigl(\frac12\bigr) = 2\int_0^\infty \eu^{-t^2}\dd t = \sqrt\pi, pela substituição t=ut = \sqrt u.)

Exercício 9.9 ★★★

Prove que Γ\Gamma é log-convexa: ln⁡Γ\ln\Gamma é convexa em (0, ∞)\intoo{0}{\infty}. (A desigualdade de Cauchy–Schwarz para integrais aplicada a t(x+y)/2−1e−t=(tx−1e−t)1/2(ty−1e−t)1/2t^{(x+y)/2 - 1}\eu^{-t} = \bigl(t^{x-1}\eu^{-t}\bigr)^{1/2} \bigl(t^{y-1}\eu^{-t}\bigr)^{1/2} dá Γ(x+y2)2≤Γ(x)Γ(y)\Gamma\bigl(\frac{x+y}{2}\bigr)^2 \leq \Gamma(x)\Gamma(y); combine com a continuidade e Exercício 8.8.)

Solução

Solução de Exercício 9.9.

Cauchy–Schwarz (volume do primeiro ano de graduação, válida em [ε, M]\intcc{\varepsilon}{M} e passada ao limite) aplicada à fatoração tx+y2−1e−t=(tx−1e−t)1/2(ty−1e−t)1/2t^{\frac{x+y}{2}-1}\eu^{-t} = \bigl(t^{x-1}\eu^{-t}\bigr)^{1/2} \bigl(t^{y-1}\eu^{-t}\bigr)^{1/2}:

Γ(x+y2)≤Γ(x)1/2 Γ(y)1/2⟹ln⁡Γ(x+y2)≤ln⁡Γ(x)+ln⁡Γ(y)2:\Gamma\Bigl(\frac{x+y}{2}\Bigr) \leq \Gamma(x)^{1/2}\,\Gamma(y)^{1/2} \quad\Longrightarrow\quad \ln\Gamma\Bigl(\frac{x+y}{2}\Bigr) \leq \frac{\ln\Gamma(x) + \ln\Gamma(y)}{2} :

ln⁡Γ\ln\Gamma é convexa no ponto médio; sendo contínua (Teorema 9.18), ela é convexa (Exercício 8.8). (A log-convexidade fixa Γ\Gamma de modo único entre as interpolações do fatorial — o teorema de Bohr–Mollerup, uma pérola do terceiro ano.)

Exercício 9.10 ★★★

(Integral de Dirichlet) Ponha F(x)=∫0∞sin⁡tt e−xt  ⁣dtF(x) = \displaystyle\int_0^{\infty} \frac{\sin t}{t}\,\eu^{-xt}\,\dd t para x>0x > 0.

  1. Justifique F′(x)=−11+x2F'(x) = -\frac{1}{1 + x^2} (derive sob a integral; calcule ∫0∞e−xtsin⁡t  ⁣dt\int_0^\infty \eu^{-xt}\sin t\,\dd t por duas integrações por partes).
  2. Prove que F(x)→0F(x) \to 0 quando x→+∞x \to +\infty e deduza F(x)=π2−arctan⁡xF(x) = \frac\pi2 - \arctan x.
  3. Admitindo a continuidade de FF em 0+0^+ (um teorema do tipo Abel), conclua o valor da integral semiconvergente:

    ∫0∞sin⁡tt  ⁣dt=π2.\int_0^{\infty} \frac{\sin t}{t}\,\dd t = \frac{\pi}{2}.
Solução

Solução de Exercício 9.10.

  1. Em x≥a>0x \geq a > 0: a derivada em xx do integrando é −sin⁡t e−xt-\sin t\,\eu^{-xt}, dominada por e−at\eu^{-at}: o Teorema 9.11 dá F′(x)=−∫0∞e−xtsin⁡t  ⁣dtF'(x) = -\int_0^\infty \eu^{-xt}\sin t\,\dd t. Duas integrações por partes (ou a exponencial complexa):

    ∫0∞e−xtsin⁡t  ⁣dt=ℑ∫0∞e(−x+i)t ⁣dt=ℑ1x−i=11+x2.\int_0^\infty \eu^{-xt}\sin t\,\dd t = \Im \int_0^\infty \eu^{(-x+\iu)t}\dd t = \Im\frac{1}{x - \iu} = \frac{1}{1 + x^2} .
  2. ∣F(x)∣≤∫0∞e−xt ⁣dt=1x→0\abs{F(x)} \leq \int_0^\infty \eu^{-xt}\dd t = \frac1x \to 0. Integrando F′=−11+x2F' = -\frac{1}{1+x^2} de xx a ∞\infty: 0−F(x)=−(π2−arctan⁡x)0 - F(x) = -\bigl(\frac\pi2 - \arctan x\bigr), logo F(x)=π2−arctan⁡xF(x) = \frac\pi2 - \arctan x.
  3. Fazendo x→0+x \to 0^+ com a continuidade admitida: F(0+)=π2F(0^+) = \frac\pi2, e F(0)=∫0∞sin⁡tt ⁣dtF(0) = \int_0^\infty \frac{\sin t}{t}\dd t (a integral semiconvergente de Dirichlet, Exemplo 9.5): seu valor é π2\frac\pi2.

Exercício 9.11 ★★

Justifique a convergência de ∫0∞(sin⁡tt) ⁣2 ⁣dt\displaystyle\int_0^\infty \Bigl(\frac{\sin t}{t}\Bigr)^{\!2}\dd t e depois calcule-a por uma integração por partes e Exercício 9.10:

∫0∞(sin⁡tt) ⁣2 ⁣dt=π2.\int_0^\infty \Bigl(\frac{\sin t}{t}\Bigr)^{\!2}\dd t = \frac{\pi}{2} .

(O mesmo valor do ∫0∞sin⁡tt ⁣dt\int_0^\infty \frac{\sin t}{t}\dd t — mas desta vez a convergência é absoluta.)

Solução

Solução de Exercício 9.11.

Convergência: perto de 00 o integrando se estende continuamente pelo valor 11 (sin⁡t∼t\sin t \sim t); no infinito ele é ≤t−2\leq t^{-2}: convergência absoluta. Em [ε, M]\intcc{\varepsilon}{M}, integre por partes com u=sin⁡2tu = \sin^2 t, v′=t−2v' = t^{-2}:

∫εMsin⁡2tt2 ⁣dt=[−sin⁡2tt]εM+∫εM2sin⁡tcos⁡tt ⁣dt=[−sin⁡2tt]εM+∫2ε2Msin⁡uu ⁣du\int_\varepsilon^M \frac{\sin^2 t}{t^2}\dd t = \Bigl[-\frac{\sin^2 t}{t}\Bigr]_\varepsilon^M + \int_\varepsilon^M \frac{2\sin t\cos t}{t}\dd t = \Bigl[-\frac{\sin^2 t}{t}\Bigr]_\varepsilon^M + \int_{2\varepsilon}^{2M} \frac{\sin u}{u}\dd u

(u=2tu = 2t na última integral). O colchete tende a 00 nas duas pontas (sin⁡2ε/ε≤ε\sin^2\varepsilon/\varepsilon \leq \varepsilon; sin⁡2M/M≤1/M\sin^2 M/M \leq 1/M), e a última integral tende a ∫0∞sin⁡uu ⁣du=π2\int_0^\infty \frac{\sin u}{u}\dd u = \frac\pi2 (Exercício 9.10). Logo

∫0∞(sin⁡tt) ⁣2 ⁣dt=π2.\int_0^\infty \Bigl(\frac{\sin t}{t}\Bigr)^{\!2}\dd t = \frac\pi2 .

Exercício 9.12 ★★★

(A cauda gaussiana) Para x>0x > 0 ponha T(x)=∫x∞e−t2 ⁣dtT(x) = \displaystyle \int_x^\infty \eu^{-t^2}\dd t.

  1. Escrevendo e−t2=1−2t⋅(−2t e−t2)\eu^{-t^2} = \frac{1}{-2t}\cdot(-2t\,\eu^{-t^2}), integre por partes duas vezes para obter

    T(x)=e−x2(12x−14x3)+34∫x∞e−t2t4  ⁣dt.T(x) = \eu^{-x^2}\Bigl(\frac{1}{2x} - \frac{1}{4x^3}\Bigr) + \frac34\int_x^\infty \frac{\eu^{-t^2}}{t^4}\,\dd t .
  2. Majore o resto: 0≤34∫x∞t−4e−t2 ⁣dt≤38x5 e−x20 \leq \frac34\int_x^\infty t^{-4}\eu^{-t^2}\dd t \leq \frac{3}{8x^5}\,\eu^{-x^2}, e deduza o enquadramento

    e−x2(12x−14x3)≤T(x)≤e−x22x,logoT(x)∼e−x22x(x→+∞).\eu^{-x^2}\Bigl(\frac{1}{2x} - \frac{1}{4x^3}\Bigr) \leq T(x) \leq \frac{\eu^{-x^2}}{2x}, \qquad\text{logo}\qquad T(x) \sim \frac{\eu^{-x^2}}{2x} \quad (x \to +\infty).
  3. Por que a série alternada completa obtida iterando as partes nunca pode convergir para xx fixo? (Compare o crescimento dos coeficientes 1⋅3⋯(2k−1)1\cdot3\cdots(2k-1) com as potências (2x2)k(2x^2)^k.)
Solução

Solução de Exercício 9.12.

  1. Partes com u=−12tu = \frac{-1}{2t}, v′=−2t e−t2v' = -2t\,\eu^{-t^2} (de modo que v=e−t2v = \eu^{-t^2}):

    T(x)=[−e−t22t]x∞−∫x∞e−t22t2 ⁣dt=e−x22x−∫x∞e−t22t2 ⁣dt.T(x) = \Bigl[\frac{-\eu^{-t^2}}{2t}\Bigr]_x^\infty - \int_x^\infty \frac{\eu^{-t^2}}{2t^2}\dd t = \frac{\eu^{-x^2}}{2x} - \int_x^\infty \frac{\eu^{-t^2}}{2t^2}\dd t .

    Mesmo dispositivo na nova integral (u=−14t3u = \frac{-1}{4t^3}, v′=−2t e−t2v' = -2t\,\eu^{-t^2}):

    ∫x∞e−t22t2 ⁣dt=e−x24x3−34∫x∞e−t2t4 ⁣dt,\int_x^\infty \frac{\eu^{-t^2}}{2t^2}\dd t = \frac{\eu^{-x^2}}{4x^3} - \frac34\int_x^\infty \frac{\eu^{-t^2}}{t^4}\dd t ,

    donde a identidade anunciada.

  2. Mais uma integração por partes majora o resto:

    ∫x∞e−t2t4 ⁣dt=e−x22x5−52∫x∞e−t2t6 ⁣dt≤e−x22x5,\int_x^\infty \frac{\eu^{-t^2}}{t^4}\dd t = \frac{\eu^{-x^2}}{2x^5} - \frac52\int_x^\infty\frac{\eu^{-t^2}}{t^6}\dd t \leq \frac{\eu^{-x^2}}{2x^5},

    logo 0≤34∫x∞t−4e−t2 ⁣dt≤38x5e−x20 \leq \frac34\int_x^\infty t^{-4}\eu^{-t^2}\dd t \leq \frac{3}{8x^5}\eu^{-x^2}. Descartar o resto (positivo) na identidade da questão 1 dá a cota inferior; descartar o segundo termo (negativo) da primeira integração por partes dá T(x)≤e−x22xT(x) \leq \frac{\eu^{-x^2}}{2x}. Dividindo o enquadramento por e−x22x\frac{\eu^{-x^2}}{2x}: a razão fica espremida entre 1−12x21 - \frac{1}{2x^2} e 11, logo T(x)∼e−x22xT(x) \sim \frac{\eu^{-x^2}}{2x}.

  3. Iterar as partes produz a série formal

    T(x)≈e−x22x(1−12x2+1⋅3(2x2)2−1⋅3⋅5(2x2)3+⋯ ),T(x) \approx \frac{\eu^{-x^2}}{2x}\Bigl(1 - \frac{1}{2x^2} + \frac{1\cdot3}{(2x^2)^2} - \frac{1\cdot3\cdot5}{(2x^2)^3} + \cdots\Bigr),

    cujo kk-ésimo coeficiente 1⋅3⋯(2k−1)=(2k)!2kk!1\cdot3\cdots(2k-1) = \frac{(2k)!}{2^k k!} cresce mais depressa do que qualquer sequência geométrica: para xx fixo os termos 1⋅3⋯(2k−1)(2x2)k\frac{1\cdot3\cdots(2k-1)}{(2x^2)^k} tendem a infinito (sua razão é 2k+12x2→∞\frac{2k+1}{2x^2} \to \infty), de modo que a série diverge para todo xx. Trata-se de um desenvolvimento assintótico: truncado em qualquer ordem fixa, o erro é da ordem do primeiro termo omitido quando x→∞x \to \infty — mas nunca de uma série convergente. (Essa estimativa de cauda é a cota padrão da cauda gaussiana dos capítulos de probabilidade.)

9.5 Problema: as integrais de Euler — Beta, Gama e a fórmula-limite de Gauss

Problema 9.1

A função Γ\Gamma da Definição 9.17 é uma metade do cálculo de integrais de Euler; a outra metade é a função Beta

B(x,y)=∫01tx−1(1−t)y−1  ⁣dt.B(x, y) = \int_0^1 t^{x-1}(1 - t)^{y-1}\,\dd t .

Este problema desenvolve o par (Γ,B)(\Gamma, B) apenas com as ferramentas deste capítulo — integração por partes, substituição, convergência dominada — e culmina na fórmula Beta–Gama de Euler B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \frac{\Gamma(x)\Gamma(y)} {\Gamma(x+y)} nos semi-inteiros e na fórmula-limite de Gauss para Γ\Gamma. Pelo caminho, as integrais de Wallis do Lema 6.11 reaparecem como valores de Beta, e a fórmula de duplicação de Legendre cai no colo.

Parte I — Estrutura fina de Γ\Gamma.

  1. Relembre por que Γ(x)=∫0∞tx−1e−t ⁣dt\Gamma(x) = \int_0^\infty t^{x-1}\eu^{-t}\dd t converge exatamente para x>0x > 0 e mostre que

    Γ(x)∼1x(x→0+)\Gamma(x) \sim \frac1x \qquad (x \to 0^+)

    (equação funcional mais continuidade de Γ\Gamma em 11).

  2. Prove Γ(12)=π\Gamma\bigl(\tfrac12\bigr) = \sqrt\pi (substitua t=u2t = u^2 e invoque o Exercício 9.8) e deduza ∫Re−u2/2 ⁣du=2π\int_\R \eu^{-u^2/2}\dd u = \sqrt{2\pi}.
  3. Mostre por indução, para n∈Nn \in \N:

    Γ(n+12)=(2n)!4n n! π.\Gamma\Bigl(n + \frac12\Bigr) = \frac{(2n)!}{4^n\,n!}\,\sqrt\pi .
  4. Justifique Γ′′(x)=∫0∞tx−1e−t(ln⁡t)2 ⁣dt>0\Gamma''(x) = \int_0^\infty t^{x-1}\eu^{-t}(\ln t)^2\dd t > 0 e deduza que Γ\Gamma é estritamente convexa, atinge um mínimo único em algum x0∈(1, 2)x_0 \in \intoo{1}{2} (Γ(1)=Γ(2)=1\Gamma(1) = \Gamma(2) = 1 e Rolle), decresce em (0, x0)\intoo{0}{x_0} e cresce em (x0, ∞)\intoo{x_0}{\infty}.
  5. Mostre que Γ\Gamma vence toda potência: para cada k∈Nk \in \N, xk=o(Γ(x))x^k = o\bigl(\Gamma(x)\bigr) quando x→+∞x \to +\infty (confronte xx entre inteiros e use Γ(n+1)=n!\Gamma(n+1) = n! com a monotonicidade da questão 4).

Parte II — A função Beta, por partes.

  1. Mostre que B(x,y)B(x,y) converge exatamente para x>0x > 0 e y>0y > 0, e que B(x,y)=B(y,x)B(x,y) = B(y,x).
  2. Calcule B(x,1)=1xB(x, 1) = \frac1x e prove, por integração por partes, para x,y>0x, y > 0:

    B(x,y+1)=yx B(x+1,y).B(x, y+1) = \frac{y}{x}\,B(x+1, y) .
  3. Da decomposição tx−1(1−t)y−1=tx(1−t)y−1+tx−1(1−t)yt^{x-1}(1-t)^{y-1} = t^{x}(1-t)^{y-1} + t^{x-1}(1-t)^{y} deduza B(x,y)=B(x+1,y)+B(x,y+1)B(x,y) = B(x+1,y) + B(x,y+1) e combine com a questão 7 nas relações de descida

    B(x,y+1)=yx+y B(x,y),B(x+1,y)=xx+y B(x,y).B(x, y+1) = \frac{y}{x+y}\,B(x,y), \qquad B(x+1, y) = \frac{x}{x+y}\,B(x,y) .
  4. Deduza, para m,n≥1m, n \geq 1 inteiros:

    B(m,n)=(m−1)! (n−1)!(m+n−1)!=1(m+n−1)(m+n−2m−1).B(m, n) = \frac{(m-1)!\,(n-1)!}{(m+n-1)!} = \frac{1}{(m+n-1)\binom{m+n-2}{m-1}} .
  5. Demonstre a fórmula de Euler com um argumento inteiro: para todo x>0x > 0 e n∈N∗n \in \N^*,

    B(x,n)=Γ(x) Γ(n)Γ(x+n)B(x, n) = \frac{\Gamma(x)\,\Gamma(n)}{\Gamma(x + n)}

    (indução sobre nn: os dois lados valem 1x\frac1x em n=1n = 1 e obedecem à mesma relação de descida).

Parte III — As integrais de Wallis como valores de Beta.

  1. Substitua t=sin⁡2θt = \sin^2\theta para obter a forma trigonométrica

    B(x,y)=2∫0π/2sin⁡2x−1θ cos⁡2y−1θ  ⁣dθ.B(x, y) = 2\int_0^{\pi/2} \sin^{2x-1}\theta\,\cos^{2y-1}\theta\,\dd\theta .
  2. Deduza Wn=12 B(n+12,12)W_n = \frac12\,B\bigl(\frac{n+1}2, \frac12\bigr) para a integral de Wallis Wn=∫0π/2sin⁡nθ  ⁣dθW_n = \int_0^{\pi/2}\sin^n \theta\,\dd\theta e recupere a recorrência Wn=n−1nWn−2W_n = \frac{n-1}{n}W_{n-2} do Lema 6.11 só a partir das relações de descida da questão 8.
  3. Calcule B(12,12)=2W0=πB\bigl(\frac12, \frac12\bigr) = 2W_0 = \pi e confira contra Γ(12)2/Γ(1)\Gamma\bigl(\frac12\bigr)^2/\Gamma(1): a fórmula de Euler vale em (12,12)\bigl(\frac12, \frac12\bigr).
  4. Deduza a forma fechada W2n=π2 (2n)!4n(n!)2W_{2n} = \frac\pi2\, \frac{(2n)!}{4^n(n!)^2} da recorrência e verifique

    B(n+12,12)=Γ(n+12)Γ(12)Γ(n+1).B\Bigl(n + \frac12, \frac12\Bigr) = \frac{\Gamma\bigl(n + \frac12\bigr)\Gamma\bigl( \frac12\bigr)}{\Gamma(n+1)} .

    Conclua, por indução com as relações de descida, que a fórmula de Euler B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} vale sempre que 2x2x e 2y2y são inteiros positivos.

  5. Substitua u=t1−tu = \frac{t}{1-t} para obter a terceira forma clássica

    B(x,y)=∫0∞ux−1(1+u)x+y  ⁣du,B(x,y) = \int_0^\infty \frac{u^{x-1}}{(1+u)^{x+y}}\,\dd u ,

    e confira o caso x=y=12x = y = \frac12 diretamente (u=v2u = v^2 reduz a ∫0∞2  ⁣dv1+v2\int_0^\infty\frac{2\,\dd v}{1+v^2}).

Parte IV — A fórmula-limite de Gauss.

  1. Para x>0x > 0 e n∈N∗n \in \N^*, prove por nn integrações por partes sucessivas:

    ∫0n(1−tn) ⁣ntx−1  ⁣dt=n!  nxx(x+1)⋯(x+n).\int_0^n \Bigl(1 - \frac tn\Bigr)^{\!n} t^{x-1}\,\dd t = \frac{n!\;n^x}{x(x+1)\cdots(x+n)} .
  2. Conclua com Exercício 9.7 (convergência dominada) a fórmula-limite de Gauss:

    Γ(x)=lim⁡n→∞n!  nxx(x+1)⋯(x+n)(x>0).\Gamma(x) = \lim_{n\to\infty} \frac{n!\;n^x}{x(x+1)\cdots(x+n)} \qquad (x > 0).
  3. Tomando logaritmos, mostre que, para x>0x > 0:

    ln⁡Γ(x)=−ln⁡x−γx+∑k=1∞(xk−ln⁡(1+xk)),\ln\Gamma(x) = -\ln x - \gamma x + \sum_{k=1}^{\infty}\Bigl(\frac xk - \ln\Bigl(1 + \frac xk\Bigr)\Bigr),

    em que γ\gamma é a constante de Euler (Exemplo 6.7); justifique a convergência da série (o termo geral é ∼x22k2\sim \frac{x^2}{2k^2}).

  4. Use a fórmula de Gauss em x=12x = \frac12 e a assintótica do binomial central (2nn)∼4nπn\binom{2n}{n} \sim \frac{4^n}{\sqrt{\pi n}} (Exemplo 6.14) para recalcular Γ(12)=π\Gamma\bigl(\frac12\bigr) = \sqrt\pi: a constante de Stirling e a integral gaussiana são o mesmo número em dois disfarces.
  5. Verifique que a fórmula de Gauss redemonstra a equação funcional: da identidade exata

    n! nx+1(x+1)⋯(x+n+1)=n! nxx(x+1)⋯(x+n)⋅n xx+n+1,\frac{n!\,n^{x+1}}{(x+1)\cdots(x+n+1)} = \frac{n!\,n^{x}}{x(x+1)\cdots(x+n)}\cdot \frac{n\,x}{x+n+1},

    conclua Γ(x+1)=x Γ(x)\Gamma(x+1) = x\,\Gamma(x) mais uma vez. (A fórmula de Gauss determina Γ\Gamma de imediato; o volume do terceiro ano de graduação demonstra o teorema mais fino de Bohr–Mollerup: a equação funcional mais a log-convexidade já fixam Γ\Gamma.)

Parte V — Dividendos.

  1. Para a>0a > 0 mostre que ∫0∞e−ta ⁣dt=Γ(1+1a)\int_0^\infty \eu^{-t^a}\dd t = \Gamma\bigl(1 + \frac1a\bigr) e calcule o limite quando a→+∞a \to +\infty por convergência dominada (limite pontual 1t<1\mathbf 1_{t < 1}; domine por 11 em (0, 1]\intoc{0}{1} e por e−t2\eu^{-t^2} além, para a≥2a \geq 2). Confira a resposta contra a continuidade de Γ\Gamma.
  2. Para n≥1n \geq 1 mostre que

    ∫01 ⁣dt1−tn=1n B(1n,12),\int_0^1 \frac{\dd t}{\sqrt{1 - t^n}} = \frac1n\,B\Bigl(\frac1n, \frac12\Bigr),

    e recupere os valores 22 (n=1n = 1) e π2\frac\pi2 (n=2n = 2). (Para n=4n = 4 essa é a constante da lemniscata, que não tem forma fechada elementar; sua história pertence à teoria das integrais elípticas.)

  3. (Momentos) Para x>0x > 0 e k∈Nk \in \N, mostre que

    1Γ(x)∫0∞tk tx−1e−t  ⁣dt=Γ(x+k)Γ(x)=x(x+1)⋯(x+k−1),\frac{1}{\Gamma(x)}\int_0^\infty t^{k}\,t^{x-1}\eu^{-t}\,\dd t = \frac{\Gamma(x+k)}{\Gamma(x)} = x(x+1)\cdots(x+k-1),

    o fatorial ascendente; confira que x=1x = 1 dá k!k!. (Nos capítulos de probabilidade esse é o kk-ésimo momento de uma densidade padrão de tempo de espera.)

  4. Demonstre a identidade Beta, válida para todos x>0x > 0:

    B(x,x)=21−2x B(x,12)B(x, x) = 2^{1-2x}\,B\Bigl(x, \frac12\Bigr)

    (substitua t=1+s2t = \frac{1+s}2, explore a simetria em ss e depois ponha s=vs = \sqrt v). Deduza, para 2x∈N∗2x \in \N^*, a fórmula de duplicação de Legendre

    Γ(x) Γ(x+12)=21−2x π  Γ(2x),\Gamma(x)\,\Gamma\Bigl(x + \frac12\Bigr) = 2^{1-2x}\,\sqrt\pi\;\Gamma(2x),

    e verifique-a diretamente em x=nx = n pela questão 3. (Para xx geral ela decorre da mesma identidade uma vez conhecida a fórmula de Euler para todos os argumentos — a demonstração por integral dupla no capítulo sobre integrais múltiplas.)

  5. Síntese. Uma frase para cada: (i) onde a integração por partes carregou toda a Parte II; (ii) onde a convergência dominada entrou nas Partes IV e V; (iii) quais insumos assintóticos foram importados do capítulo de comparação; (iv) o que está agora demonstrado da fórmula de Euler B(x,y)=Γ(x)Γ(y)/Γ(x+y)B(x,y) = \Gamma(x)\Gamma(y)/\Gamma(x+y), e o que resta para a integral dupla resolver.
Solução

Solução de Problema 9.1.

1. Em 0+0^+ o integrando é ∼tx−1\sim t^{x-1}: a escala de extremidade finita converge se e somente se 1−x<11 - x < 1, isto é, x>0x > 0 (e, para x≤0x \leq 0, tx−1≥t−1t^{x-1} \geq t^{-1} diverge); em +∞+\infty, tx−1e−t=o(t−2)t^{x-1}\eu^{-t} = o(t^{-2}) converge para todo xx. Então Γ(x)=Γ(x+1)x\Gamma(x) = \frac{\Gamma(x+1)}{x} e Γ(x+1)→Γ(1)=1\Gamma(x+1) \to \Gamma(1) = 1 quando x→0+x \to 0^+ (continuidade, Teorema 9.18): Γ(x)∼1x\Gamma(x) \sim \frac1x.

2. Com t=u2t = u^2,  ⁣dt=2u  ⁣du\dd t = 2u\,\dd u:

Γ(12)=∫0∞t−1/2e−t ⁣dt=∫0∞e−u2u 2u  ⁣du=2∫0∞e−u2 ⁣du=π\Gamma\Bigl(\frac12\Bigr) = \int_0^\infty t^{-1/2}\eu^{-t}\dd t = \int_0^\infty \frac{\eu^{-u^2}}{u}\,2u\,\dd u = 2\int_0^\infty \eu^{-u^2}\dd u = \sqrt\pi

pelo Exercício 9.8. Com u=v/2u = v/\sqrt2:

∫Re−v2/2 ⁣dv=22∫0∞e−u2 ⁣du=2 π=2π.\int_\R \eu^{-v^2/2}\dd v = 2\sqrt2\int_0^\infty \eu^{-u^2}\dd u = \sqrt2\,\sqrt\pi = \sqrt{2\pi} .

3. Verdadeiro para n=0n = 0 (os dois lados valem π\sqrt\pi). Se Γ(n+12)=(2n)!4nn!π\Gamma(n + \frac12) = \frac{(2n)!}{4^n n!}\sqrt\pi, a equação funcional dá

Γ(n+1+12)=(n+12)Γ(n+12)=2n+12⋅(2n)!4nn!π=(2n+2)!4n+1(n+1)!π,\Gamma\Bigl(n + 1 + \frac12\Bigr) = \Bigl(n + \frac12\Bigr)\Gamma\Bigl(n + \frac12\Bigr) = \frac{2n+1}{2}\cdot\frac{(2n)!}{4^n n!}\sqrt\pi = \frac{(2n+2)!}{4^{n+1}(n+1)!}\sqrt\pi ,

o último passo porque (2n+2)!(2n)!=(2n+2)(2n+1)\frac{(2n+2)!}{(2n)!} = (2n+2)(2n+1) e 2n+12=(2n+2)(2n+1)4(n+1)\frac{2n+1}{2} = \frac{(2n+2)(2n+1)}{4(n+1)}.

4. O Teorema 9.18 dá Γ′′(x)=∫0∞tx−1e−t(ln⁡t)2 ⁣dt\Gamma''(x) = \int_0^\infty t^{x-1}\eu^{-t}(\ln t)^2\dd t (duas aplicações da regra de Leibniz, com dominações como na demonstração do teorema); o integrando é ≥0\geq 0 e não identicamente nulo, logo Γ′′>0\Gamma'' > 0: Γ\Gamma é estritamente convexa e Γ′\Gamma' é estritamente crescente. Como Γ(1)=Γ(2)=1\Gamma(1) = \Gamma(2) = 1, Rolle fornece x0∈(1, 2)x_0 \in \intoo12 com Γ′(x0)=0\Gamma'(x_0) = 0; a monotonicidade estrita de Γ′\Gamma' faz de x0x_0 seu único zero, com Γ′<0\Gamma' < 0 antes e Γ′>0\Gamma' > 0 depois: Γ\Gamma decresce em (0, x0)\intoo0{x_0}, cresce em (x0, ∞)\intoo{x_0}\infty, e x0x_0 é o mínimo único.

5. Sejam k∈Nk \in \N e x≥3x \geq 3; escolha o inteiro nn com n+1≤x<n+2n + 1 \leq x < n + 2 (de modo que n≥1n \geq 1). Pela monotonicidade da questão 4 (válida a partir de x0<2x_0 < 2): Γ(x)≥Γ(n+1)=n!\Gamma(x) \geq \Gamma(n + 1) = n!, enquanto xk≤(n+2)kx^k \leq (n+2)^k. Como n!(n+2)k→∞\frac{n!}{(n+2)^k} \to \infty (os fatoriais vencem as potências, volume do primeiro ano de graduação), Γ(x)xk≥n!(n+2)k→∞\frac{\Gamma(x)}{x^k} \geq \frac{n!}{(n+2)^k} \to \infty quando x→∞x \to \infty: xk=o(Γ(x))x^k = o(\Gamma(x)).

6. Perto de 00 o integrando é ∼tx−1\sim t^{x-1} (convergente se e somente se x>0x > 0), perto de 11 ele é ∼(1−t)y−1\sim (1-t)^{y-1} (se e somente se y>0y > 0); as duas comparações são entre funções positivas, logo B(x,y)B(x,y) converge exatamente para x,y>0x, y > 0. A substituição t↦1−tt \mapsto 1 - t troca os dois fatores: B(x,y)=B(y,x)B(x,y) = B(y,x).

7. B(x,1)=∫01tx−1 ⁣dt=1xB(x,1) = \int_0^1 t^{x-1}\dd t = \frac1x. Partes em [ε, 1−ε]\intcc\varepsilon{1-\varepsilon} com u=(1−t)yu = (1-t)^y, v=txxv = \frac{t^x}{x}:

∫tx−1(1−t)y ⁣dt=[tx(1−t)yx]+yx∫tx(1−t)y−1 ⁣dt;\int t^{x-1}(1-t)^{y}\dd t = \Bigl[\frac{t^x(1-t)^y}{x}\Bigr] + \frac{y}{x}\int t^{x}(1-t)^{y-1}\dd t ;

o colchete se anula nas duas pontas quando ε→0\varepsilon \to 0 (x>0x > 0 em 00, y>0y > 0 em 11), restando B(x,y+1)=yx B(x+1,y)B(x, y+1) = \frac yx\,B(x+1, y).

8. Como t+(1−t)=1t + (1-t) = 1:

tx−1(1−t)y−1=tx(1−t)y−1+tx−1(1−t)y,t^{x-1}(1-t)^{y-1} = t^{x}(1-t)^{y-1} + t^{x-1}(1-t)^{y},

logo B(x,y)=B(x+1,y)+B(x,y+1)B(x,y) = B(x+1,y) + B(x,y+1). A questão 7 se lê B(x+1,y)=xyB(x,y+1)B(x+1,y) = \frac xy B(x,y+1); substituindo,

B(x,y)=(xy+1)B(x,y+1)=x+yy B(x,y+1),B(x,y) = \Bigl(\frac xy + 1\Bigr)B(x,y+1) = \frac{x+y}{y}\,B(x,y+1),

isto é, B(x,y+1)=yx+yB(x,y)B(x,y+1) = \frac{y}{x+y}B(x,y); a relação gêmea segue da simetria da questão 6.

9. Indução sobre nn com mm fixo: B(m,1)=1m=(m−1)! 0!m!B(m,1) = \frac1m = \frac{(m-1)!\,0!}{m!} e, se a fórmula vale em nn,

B(m,n+1)=nm+n B(m,n)=nm+n⋅(m−1)!(n−1)!(m+n−1)!=(m−1)! n!(m+n)!.B(m, n+1) = \frac{n}{m+n}\,B(m,n) = \frac{n}{m+n}\cdot\frac{(m-1)!(n-1)!}{(m+n-1)!} = \frac{(m-1)!\,n!}{(m+n)!} .

Reescrevendo: B(m,n)=(m−1)!(n−1)!(m+n−1)!=[(m+n−1)(m+n−2m−1)]−1B(m,n) = \frac{(m-1)!(n-1)!}{(m+n-1)!} = \bigl[(m+n-1)\binom{m+n-2}{m-1}\bigr]^{-1}.

10. Os dois lados de B(x,n)=Γ(x)Γ(n)Γ(x+n)B(x,n) = \frac{\Gamma(x)\Gamma(n)}{\Gamma(x+n)} valem 1x\frac1x em n=1n = 1 (Γ(1)=1\Gamma(1) = 1, Γ(x+1)=xΓ(x)\Gamma(x+1) = x\Gamma(x)). Se eles coincidem em nn, então, pela relação de descida e pela equação funcional:

B(x,n+1)=nx+n B(x,n),Γ(x)Γ(n+1)Γ(x+n+1)=nx+n⋅Γ(x)Γ(n)Γ(x+n):B(x, n+1) = \frac{n}{x+n}\,B(x,n), \qquad \frac{\Gamma(x)\Gamma(n+1)}{\Gamma(x+n+1)} = \frac{n}{x+n}\cdot \frac{\Gamma(x)\Gamma(n)}{\Gamma(x+n)} :

as duas sequências obedecem à mesma recursão a partir da mesma semente, logo coincidem para todo n∈N∗n \in \N^* e todo x>0x > 0.

11. Com t=sin⁡2θt = \sin^2\theta (θ∈(0, π/2)\theta \in \intoo0{\pi/2},  ⁣dt=2sin⁡θcos⁡θ  ⁣dθ\dd t = 2\sin\theta\cos\theta\,\dd\theta), tx−1=sin⁡2x−2θt^{x-1} = \sin^{2x-2}\theta e (1−t)y−1=cos⁡2y−2θ(1-t)^{y-1} = \cos^{2y-2}\theta:

B(x,y)=∫0π/2sin⁡2x−2θ cos⁡2y−2θ⋅2sin⁡θcos⁡θ  ⁣dθ=2∫0π/2sin⁡2x−1θ cos⁡2y−1θ  ⁣dθ.B(x,y) = \int_0^{\pi/2}\sin^{2x-2}\theta\,\cos^{2y-2}\theta \cdot 2\sin\theta\cos\theta\,\dd\theta = 2\int_0^{\pi/2}\sin^{2x-1}\theta\,\cos^{2y-1}\theta\, \dd\theta .

12. Tome y=12y = \frac12 (matando o fator cosseno) e 2x−1=n2x - 1 = n: B(n+12,12)=2WnB\bigl(\frac{n+1}2, \frac12\bigr) = 2W_n, isto é, Wn=12B(n+12,12)W_n = \frac12 B\bigl(\frac{n+1}2,\frac12\bigr). A relação de descida na primeira variável dá

WnWn−2=B(n−12+1,12)B(n−12,12)=n−12n−12+12=n−1n:\frac{W_n}{W_{n-2}} = \frac{B\bigl(\frac{n-1}2 + 1, \frac12\bigr)} {B\bigl(\frac{n-1}2, \frac12\bigr)} = \frac{\frac{n-1}2}{\frac{n-1}2 + \frac12} = \frac{n-1}{n} :

a recorrência de Wallis, desta vez sem nenhuma integração por partes sobre senos — a Parte II fez o trabalho de uma vez por todas.

13. B(12,12)=2W0=2⋅π2=πB\bigl(\frac12,\frac12\bigr) = 2W_0 = 2\cdot\frac\pi2 = \pi, enquanto Γ(12)2/Γ(1)=(π)2=π\Gamma\bigl(\frac12\bigr)^2/\Gamma(1) = (\sqrt\pi)^2 = \pi: a fórmula de Euler vale em (12,12)\bigl(\frac12,\frac12\bigr).

14. Iterando W2n=2n−12nW2n−2W_{2n} = \frac{2n-1}{2n}W_{2n-2} a partir de W0=π2W_0 = \frac\pi2:

W2n=π2∏k=1n2k−12k=π2⋅(2n)!4n(n!)2,W_{2n} = \frac\pi2\prod_{k=1}^{n}\frac{2k-1}{2k} = \frac\pi2\cdot\frac{(2n)!}{4^n(n!)^2},

pois ∏(2k−1)=(2n)!2nn!\prod(2k-1) = \frac{(2n)!}{2^n n!} e ∏2k=2nn!\prod 2k = 2^n n!. Logo, usando a questão 3:

B(n+12,12)=2W2n=π (2n)!4n(n!)2=(2n)!π4nn!⋅πn!=Γ(n+12)Γ(12)Γ(n+1).B\Bigl(n+\frac12, \frac12\Bigr) = 2W_{2n} = \pi\,\frac{(2n)!}{4^n(n!)^2} = \frac{(2n)!\sqrt\pi}{4^n n!}\cdot\frac{\sqrt\pi}{n!} = \frac{\Gamma\bigl(n+\frac12\bigr)\Gamma\bigl(\frac12\bigr)} {\Gamma(n+1)} .

Fixe agora x∈12N∗x \in \frac12\N^*. A fórmula de Euler vale em (x,12)(x, \frac12): para xx inteiro isso é a questão 10 (com a simetria), para x=n+12x = n + \frac12 é a fórmula acima. Os dois lados da fórmula de Euler obedecem à recursão de descida y↦y+1y \mapsto y + 1 (questão 8 à esquerda, equação funcional à direita, como na questão 10): a indução propaga a fórmula de y=12y = \frac12 e y=1y = 1 a todo y∈12N∗y \in \frac12\N^*. A fórmula de Euler vale, portanto, sempre que 2x,2y∈N∗2x, 2y \in \N^*.

15. Com u=t1−tu = \frac{t}{1-t}, isto é, t=u1+ut = \frac{u}{1+u}, 1−t=11+u1 - t = \frac{1}{1+u},  ⁣dt= ⁣du(1+u)2\dd t = \frac{\dd u}{(1+u)^2}:

B(x,y)=∫0∞(u1+u)x−1(11+u)y−1 ⁣du(1+u)2=∫0∞ux−1(1+u)x+y  ⁣du.B(x,y) = \int_0^\infty \Bigl(\frac{u}{1+u}\Bigr)^{x-1} \Bigl(\frac{1}{1+u}\Bigr)^{y-1} \frac{\dd u}{(1+u)^2} = \int_0^\infty \frac{u^{x-1}}{(1+u)^{x+y}}\,\dd u .

Em x=y=12x = y = \frac12, com u=v2u = v^2:

∫0∞u−1/21+u ⁣du=∫0∞2  ⁣dv1+v2=π=B(12,12).✓\int_0^\infty \frac{u^{-1/2}}{1+u}\dd u = \int_0^\infty \frac{2\,\dd v}{1+v^2} = \pi = B\Bigl(\frac12,\frac12\Bigr) . \checkmark

16. Uma integração por partes, para 1≤k≤n1 \leq k \leq n e s>0s > 0 (u=(1−t/n)ku = (1 - t/n)^k, v=ts/sv = t^s/s; os termos de bordo se anulam):

∫0n(1−tn) ⁣kts−1 ⁣dt=kns∫0n(1−tn) ⁣k−1ts ⁣dt.\int_0^n \Bigl(1-\frac tn\Bigr)^{\!k} t^{s-1}\dd t = \frac{k}{ns}\int_0^n \Bigl(1-\frac tn\Bigr)^{\!k-1} t^{s}\dd t .

Partindo de k=nk = n, s=xs = x e iterando nn vezes:

∫0n(1−tn) ⁣ntx−1 ⁣dt=n(n−1)⋯1nn x(x+1)⋯(x+n−1)∫0ntx+n−1 ⁣dt=n!nn⋅nx+nx(x+1)⋯(x+n),\int_0^n \Bigl(1-\frac tn\Bigr)^{\!n} t^{x-1}\dd t = \frac{n(n-1)\cdots1}{n^n\,x(x+1)\cdots(x+n-1)} \int_0^n t^{x+n-1}\dd t = \frac{n!}{n^n}\cdot \frac{n^{x+n}}{x(x+1)\cdots(x+n)} ,

que é n! nxx(x+1)⋯(x+n)\dfrac{n!\,n^x}{x(x+1)\cdots(x+n)}.

17. Pelo Exercício 9.7 o membro da esquerda tende a Γ(x)\Gamma(x) (convergência dominada com dominante tx−1e−tt^{x-1}\eu^{-t}); o membro da direita é o quociente de Gauss:

Γ(x)=lim⁡n→∞n! nxx(x+1)⋯(x+n).\Gamma(x) = \lim_{n\to\infty} \frac{n!\,n^x}{x(x+1)\cdots(x+n)} .

18. Tomando logaritmos no quociente da questão 16, Gn(x)G_n(x), e separando ln⁡(x+k)=ln⁡k+ln⁡(1+x/k)\ln(x+k) = \ln k + \ln(1 + x/k) para k≥1k \geq 1:

ln⁡Gn(x)=xln⁡n−ln⁡x−∑k=1nln⁡(1+xk)=−ln⁡x+x(ln⁡n−Hn)+∑k=1n(xk−ln⁡(1+xk)).\ln G_n(x) = x\ln n - \ln x - \sum_{k=1}^n \ln\Bigl(1+\frac xk\Bigr) = -\ln x + x(\ln n - H_n) + \sum_{k=1}^n\Bigl(\frac xk - \ln\Bigl(1+\frac xk\Bigr)\Bigr).

Para u≥0u \geq 0, u−u22≤ln⁡(1+u)≤uu - \frac{u^2}2 \leq \ln(1+u) \leq u, de modo que o termo geral está em [0, x2/(2k2)]\intcc{0}{x^2/(2k^2)}: a série converge (comparação com ∑k−2\sum k^{-2}). Como ln⁡n−Hn→−γ\ln n - H_n \to -\gamma (Exemplo 6.7) e ln⁡Gn(x)→ln⁡Γ(x)\ln G_n(x) \to \ln\Gamma(x) (questão 17 e continuidade de ln⁡\ln):

ln⁡Γ(x)=−ln⁡x−γx+∑k=1∞(xk−ln⁡(1+xk)).\ln\Gamma(x) = -\ln x - \gamma x + \sum_{k=1}^\infty\Bigl(\frac xk - \ln\Bigl(1+\frac xk\Bigr)\Bigr) .

19. Em x=12x = \frac12, o denominador é ∏k=0n(k+12)=(2n+1)!22n+1n!\prod_{k=0}^n\bigl(k+\frac12\bigr) = \frac{(2n+1)!}{2^{2n+1}n!} (desenvolvendo as metades), logo

Gn(12)=n! n  22n+1n!(2n+1)!=2n  4n(2n+1)(2nn).G_n\Bigl(\frac12\Bigr) = \frac{n!\,\sqrt n\;2^{2n+1}n!}{(2n+1)!} = \frac{2\sqrt n\;4^n}{(2n+1)\binom{2n}{n}} .

Com (2nn)∼4nπn\binom{2n}{n} \sim \frac{4^n}{\sqrt{\pi n}} (Exemplo 6.14):

Gn(12)∼2n πn2n+1⟶π=Γ(12).G_n\Bigl(\frac12\Bigr) \sim \frac{2\sqrt n\,\sqrt{\pi n}}{2n+1} \longrightarrow \sqrt\pi = \Gamma\Bigl(\frac12\Bigr) .

O π\sqrt\pi do coeficiente binomial central (que veio de Wallis, logo da constante de Stirling) e o π\sqrt\pi da integral gaussiana são o mesmo número.

20. A identidade é álgebra direta: multiplique n! nxx(x+1)⋯(x+n)\frac{n!\,n^x}{x(x+1)\cdots(x+n)} por nxx+n+1\frac{nx}{x+n+1} e absorva xx no produto e nn em nxn^x. Fazendo n→∞n \to \infty: o membro da esquerda tende a Γ(x+1)\Gamma(x+1) (Gauss em x+1x+1), o da direita a Γ(x)⋅x⋅1\Gamma(x)\cdot x\cdot 1, pois nx+n+1→1\frac{n}{x+n+1} \to 1: Γ(x+1)=xΓ(x)\Gamma(x+1) = x\Gamma(x) — recuperado sem uma única integração por partes.

21. Com u=tau = t^a, t=u1/at = u^{1/a},  ⁣dt=1au1/a−1 ⁣du\dd t = \frac1a u^{1/a - 1}\dd u:

∫0∞e−ta ⁣dt=1a∫0∞u1a−1e−u ⁣du=1a Γ(1a)=Γ(1+1a).\int_0^\infty \eu^{-t^a}\dd t = \frac1a\int_0^\infty u^{\frac1a - 1}\eu^{-u}\dd u = \frac1a\,\Gamma\Bigl(\frac1a\Bigr) = \Gamma\Bigl(1 + \frac1a\Bigr) .

Quando a→+∞a \to +\infty (ao longo de qualquer sequência): e−ta→1\eu^{-t^a} \to 1 para 0<t<10 < t < 1, →e−1\to \eu^{-1} em t=1t = 1, →0\to 0 para t>1t > 1; para a≥2a \geq 2 domine por 1t≤1+e−t21t>1\mathbf 1_{t \leq 1} + \eu^{-t^2}\mathbf 1_{t > 1} (ta≥t2t^a \geq t^2 para t≥1t \geq 1), integrável. Convergência dominada: a integral tende a ∫011  ⁣dt=1\int_0^1 1\,\dd t = 1 — como tinha de ser, pois Γ(1+1a)→Γ(1)=1\Gamma(1 + \frac1a) \to \Gamma(1) = 1 pela continuidade.

22. Com u=tnu = t^n,  ⁣dt=1nu1/n−1 ⁣du\dd t = \frac1n u^{1/n - 1}\dd u:

∫01 ⁣dt1−tn=1n∫01u1n−1(1−u)−1/2 ⁣du=1n B(1n,12).\int_0^1 \frac{\dd t}{\sqrt{1-t^n}} = \frac1n\int_0^1 u^{\frac1n-1}(1-u)^{-1/2}\dd u = \frac1n\,B\Bigl(\frac1n, \frac12\Bigr) .

n=1n = 1: B(1,12)=B(12,1)=2B\bigl(1,\frac12\bigr) = B\bigl(\frac12,1\bigr) = 2, coincidindo com ∫01 ⁣dt1−t=2\int_0^1\frac{\dd t}{\sqrt{1-t}} = 2. n=2n = 2: 12B(12,12)=π2=arcsin⁡1\frac12 B\bigl(\frac12,\frac12\bigr) = \frac\pi2 = \arcsin 1. Para n=4n = 4 o valor 14B(14,12)\frac14 B\bigl(\frac14,\frac12\bigr) é a constante da lemniscata: nenhuma forma fechada elementar.

23. Iterando a equação funcional:

1Γ(x)∫0∞tx+k−1e−t ⁣dt=Γ(x+k)Γ(x)=(x+k−1)(x+k−2)⋯x,\frac{1}{\Gamma(x)}\int_0^\infty t^{x+k-1}\eu^{-t}\dd t = \frac{\Gamma(x+k)}{\Gamma(x)} = (x+k-1)(x+k-2)\cdots x ,

o fatorial ascendente com kk fatores. Em x=1x = 1: Γ(1+k)/Γ(1)=k!\Gamma(1+k)/\Gamma(1) = k!, os momentos de e−t\eu^{-t} vindos do Exercício 9.2.

24. Substitua t=1+s2t = \frac{1+s}2 (s∈(−1, 1)s \in \intoo{-1}1,  ⁣dt= ⁣ds2\dd t = \frac{\dd s}2, t(1−t)=1−s24t(1-t) = \frac{1-s^2}4):

B(x,x)=∫−11(1−s24)x−1 ⁣ds2=41−x∫01(1−s2)x−1 ⁣dsB(x,x) = \int_{-1}^{1}\Bigl(\frac{1-s^2}{4}\Bigr)^{x-1} \frac{\dd s}{2} = 4^{1-x}\int_0^1 (1-s^2)^{x-1}\dd s

(o integrando é par). Depois s=vs = \sqrt v ( ⁣ds= ⁣dv2v\dd s = \frac{\dd v}{2\sqrt v}):

B(x,x)=41−x2∫01v−1/2(1−v)x−1 ⁣dv=21−2x B(12,x).B(x,x) = \frac{4^{1-x}}{2}\int_0^1 v^{-1/2}(1-v)^{x-1}\dd v = 2^{1-2x}\,B\Bigl(\frac12, x\Bigr) .

Para 2x∈N∗2x \in \N^* todos os argumentos à vista estão em 12N∗\frac12\N^*, de modo que a fórmula de Euler (questão 14) se aplica aos dois lados:

Γ(x)2Γ(2x)=21−2x Γ(12)Γ(x)Γ(x+12)⟺Γ(x) Γ(x+12)=21−2xπ  Γ(2x).\frac{\Gamma(x)^2}{\Gamma(2x)} = 2^{1-2x}\, \frac{\Gamma\bigl(\frac12\bigr)\Gamma(x)} {\Gamma\bigl(x+\frac12\bigr)} \quad\Longleftrightarrow\quad \Gamma(x)\,\Gamma\Bigl(x+\frac12\Bigr) = 2^{1-2x}\sqrt\pi\;\Gamma(2x) .

Verificação direta em x=nx = n: o membro da esquerda é (n−1)!⋅(2n)!π4nn!=(2n)!π4nn(n-1)!\cdot \frac{(2n)!\sqrt\pi}{4^n n!} = \frac{(2n)!\sqrt\pi}{4^n n}, o da direita 2⋅4−nπ (2n−1)!=(2n)!π4nn2\cdot4^{-n}\sqrt\pi\,(2n-1)! = \frac{(2n)!\sqrt\pi}{4^n n}: iguais.

25. (i) A integração por partes produziu B(x,y+1)=yxB(x+1,y)B(x,y+1) = \frac yx B(x+1,y), a única identidade da qual decorrem toda relação de descida, os valores inteiros e semi-inteiros e a recorrência de Wallis. (ii) A convergência dominada transformou as integrais elementares ∫0n(1−t/n)ntx−1\int_0^n(1-t/n)^n t^{x-1} em Γ(x)\Gamma(x) (fórmula de Gauss, questão 17) e calculou o limite a→∞a \to \infty na questão 21. (iii) Do capítulo de comparação importamos a constante de Euler (ln⁡n−Hn→−γ\ln n - H_n \to -\gamma, questão 18) e a assintótica do binomial central (questão 19) — isto é, a fórmula de Stirling disfarçada. (iv) A fórmula de Euler B(x,y)=Γ(x)Γ(y)/Γ(x+y)B(x,y) = \Gamma(x)\Gamma(y)/\Gamma(x+y) está agora demonstrada para y∈N∗y \in \N^* com x>0x > 0 arbitrário (questão 10) e para todos os pares semi-inteiros (questão 14); o caso geral x,y>0x, y > 0 aguarda o cálculo por integral dupla do capítulo sobre integrais múltiplas, que fatoriza Γ(x)Γ(y)\Gamma(x)\Gamma(y) sobre um quarto de plano.

Termos definidos neste capítulo

Ver todos os 395 termos do glossário