Mathematics · Livro 3 · Bachelor Year 1

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

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

15Integração num Segmento

A integral do volume do ensino médio foi fundada em áreas tomadas intuitivamente. Este capítulo a constrói: primeiro para funções escada, em que a integral é uma soma finita, e depois para funções contínuas (e contínuas por partes) por aproximação uniforme — o lugar em que o teorema de Heine (Teorema 13.22) ganha o seu sustento. O teorema fundamental do cálculo conecta então a construção às primitivas, e as somas de Riemann a conectam às médias discretas.

Ao longo do capítulo, a<ba < b são reais.

15.1 Funções escada

Definição 15.1

φ ⁣:[a,b]R\varphi \colon \intcc{a}{b} \to \R é uma função escada quando existe uma subdivisão a=x0<x1<<xn=ba = x_0 < x_1 < \dots < x_n = b tal que φ\varphi é constante, igual a cic_i, em cada intervalo aberto (xi1,xi)\intoo{x_{i-1}}{x_i} (os valores nos nós não são restritos). A sua integral é

abφ=i=1nci(xixi1),\int_a^b \varphi = \sum_{i=1}^{n} c_i\,(x_i - x_{i-1}),

independente da subdivisão escolhida (refine duas subdivisões pela sua comum: cada lado fica inalterado por refinamento).

Proposição 15.2

Nas funções escada, a integral é linear, crescente (φψ    φψ\varphi \leq \psi \implies \int\varphi \leq \int\psi) e satisfaz a relação de Chasles ab=ac+cb\int_a^b = \int_a^c + \int_c^b para a<c<ba < c < b.

Demonstração. O motor é a invariância por refinamento, enunciada na definição: inserir um nó extra t(xi1,xi)t \in \intoo{x_{i-1}}{x_i} numa subdivisão substitui o termo ci(xixi1)c_i(x_i - x_{i-1}) por ci(txi1)+ci(xit)c_i(t - x_{i-1}) + c_i(x_i - t) — o mesmo número — de modo que a integral não muda por refinamento finito algum. Tome agora φ\varphi com subdivisão σ\sigma e ψ\psi com subdivisão σ\sigma': no refinamento comum σσ\sigma \cup \sigma', as duas são funções escada com os mesmos nós e, em cada peça, φ+λψ\varphi + \lambda\psi é constante igual a ci+λdic_i + \lambda d_i: a linearidade se reduz à linearidade de somas finitas. Crescimento: cidic_i \leq d_i em cada peça dá ciΔidiΔi\sum c_i \Delta_i \leq \sum d_i \Delta_i (comprimentos Δi0\Delta_i \geq 0). Chasles: insira o nó cc e separe a soma nele.

15.2 Integral de uma função contínua

Teorema 15.3 (Aproximação uniforme)

Seja ff contínua em [a,b]\intcc{a}{b}. Para todo ε>0\varepsilon > 0 existem funções escada φ,ψ\varphi, \psi com

φfψeψφε em[a,b].\varphi \leq f \leq \psi \qquad\text{e}\qquad \psi - \varphi \leq \varepsilon \text{ em} \intcc{a}{b}.

Demonstração. Pelo teorema de Heine (Teorema 13.22), ff é uniformemente contínua: tome δ\delta para ε\varepsilon e uma subdivisão de malha <δ< \delta (igualmente espaçada, digamos, com n>baδn > \frac{b - a}{\delta} peças). Em cada peça fechada [xi1,xi]\intcc{x_{i-1}}{x_i}, ff atinge um mínimo mim_i e um máximo MiM_i (Teorema 13.13), e MimiεM_i - m_i \leq \varepsilon (os dois pontos extremais estão a menos de δ\delta). Defina φ=mi\varphi = m_i e ψ=Mi\psi = M_i em (xi1,xi)\intoo{x_{i-1}}{x_i} (e φ=ψ=f\varphi = \psi = f nos nós).

Teorema 15.4 (Definição da integral)

Seja ff contínua em [a,b]\intcc{a}{b}. Os dois números

I(f)=sup{abφ:φ escada, φf},I+(f)=inf{abψ:ψ escada, ψf}I_-(f) = \sup\Bigl\{\int_a^b \varphi : \varphi \text{ escada},\ \varphi \leq f\Bigr\}, \qquad I_+(f) = \inf\Bigl\{\int_a^b \psi : \psi \text{ escada},\ \psi \geq f\Bigr\}

são iguais; o seu valor comum é a integral abf\int_a^b f (também escrita abf(t) ⁣dt\int_a^b f(t)\,\dd t). Ele coincide com a noção anterior nas funções escada e se estende às funções contínuas por partes, separando [a,b]\intcc{a}{b} nas descontinuidades (Chasles como definição ali).

Demonstração. Os dois conjuntos são não vazios (ff é limitada) e toda integral escada inferior é \leq toda superior (monotonicidade nas funções escada): assim, I(f)I+(f)I_-(f) \leq I_+(f). Pelo Teorema 15.3, para todo ε\varepsilon existe um par com ψφε(ba)\int\psi - \int\varphi \leq \varepsilon(b - a): o supremo e o ínfimo ficam comprimidos, I=I+I_- = I_+.

Exemplo 15.5 (Contínua por partes, sem drama)

A função parte inteira em [0,3]\intcc{0}{3} é uma função escada disfarçada: separando nos seus saltos,

03t ⁣dt=010+121+232=0+1+2=3,\int_0^3 \lfloor t \rfloor\,\dd t = \int_0^1 0 + \int_1^2 1 + \int_2^3 2 = 0 + 1 + 2 = 3 ,

e os valores nos pontos de salto 1,21, 2 são irrelevantes: mudar uma função em finitos pontos não muda integral alguma (as funções escada que a enquadram não são afetadas). É todo o conteúdo da extensão a “contínua por partes”: corte nas finitas descontinuidades, integre cada peça contínua e some — Chasles como definição.

Exemplo 15.6 (A definição calcula, uma vez)

Seja f(x)=xf(x) = x em [0,1]\intcc{0}{1} e corte em nn peças iguais. As melhores funções escada constantes nas peças são φ=k1n\varphi = \frac{k-1}{n} e ψ=kn\psi = \frac kn na kk-ésima peça, com

01φ=k=1nk1n1n=n12n,01ψ=k=1nkn1n=n+12n.\int_0^1 \varphi = \sum_{k=1}^{n} \frac{k-1}{n}\cdot\frac1n = \frac{n-1}{2n}, \qquad \int_0^1 \psi = \sum_{k=1}^{n} \frac{k}{n}\cdot\frac1n = \frac{n+1}{2n} .

Toda integral inferior é I(f)I+(f)\leq I_-(f) \leq I_+(f) \leq toda superior, de modo que n12nI(f)I+(f)n+12n\frac{n-1}{2n} \leq I_-(f) \leq I_+(f) \leq \frac{n+1}{2n} para todo nn: as duas se comprimem sobre 12\frac12, e 01x ⁣dx=12\int_0^1 x\,\dd x = \frac12 sai direto da definição. A ideia de fechamento: esta é a primeira e a última vez que integramos pela definição — o teorema fundamental, abaixo, substitui todos esses cálculos por uma consulta de primitiva, o que é todo o ponto econômico deste capítulo.

Teorema 15.7 (Propriedades)

Para f,gf, g contínuas (ou contínuas por partes) em [a,b]\intcc{a}{b} e λR\lambda \in \R:

  1. linearidade: (f+λg)=f+λg\int (f + \lambda g) = \int f + \lambda \int g;
  2. monotonicidade: fg    abfabgf \leq g \implies \int_a^b f \leq \int_a^b g; e abfabf(ba)supf\bigl|\int_a^b f\bigr| \leq \int_a^b \abs f \leq (b - a)\, \sup\abs f;
  3. Chasles: ab=ac+cb\int_a^b = \int_a^c + \int_c^b (com a convenção ba=ab\int_b^a = -\int_a^b, válida para qualquer ordem dos limites);
  4. positividade estrita: se ff é contínua, f0f \geq 0 e abf=0\int_a^b f = 0, então f=0f = 0 em toda parte em [a,b]\intcc{a}{b}.

Demonstração. (1)–(3) passam das funções escada ao limite pela definição por sup/inf. A linearidade merece os detalhes uma vez: dado ε>0\varepsilon > 0, enquadre φffψf\varphi_f \leq f \leq \psi_f e φggψg\varphi_g \leq g \leq \psi_g com folgas ε\leq \varepsilon (Teorema 15.3). Para λ0\lambda \geq 0, φf+λφgf+λgψf+λψg\varphi_f + \lambda\varphi_g \leq f + \lambda g \leq \psi_f + \lambda\psi_g é um enquadramento por funções escada de folga (1+λ)ε\leq (1 + \lambda)\varepsilon, e as suas integrais escada valem φf+λφg\int \varphi_f + \lambda\int\varphi_g, etc. (Proposição 15.2): fazer ε0\varepsilon \to 0 comprime (f+λg)\int(f + \lambda g) sobre f+λg\int f + \lambda\int g. Para λ<0\lambda < 0, multiplicar por λ\lambda inverte o enquadramento de gg — a função escada inferior de λg\lambda g é λψg\lambda \psi_g — e a mesma compressão roda com os papéis trocados. A estimativa ff\abs{\int f} \leq \int \abs f vem de fff-\abs f \leq f \leq \abs f e da monotonicidade.

(4) Contrapositiva: se f(x0)=m>0f(x_0) = m > 0, a continuidade fornece um subintervalo de comprimento η>0\eta > 0 no qual fm2f \geq \frac m2; a função escada que vale m2\frac m2 ali e 00 no resto é f\leq f, de modo que fmη2>0\int f \geq \frac{m\eta}{2} > 0.

Exemplo 15.8 (Chasles em ação: integrais com valores absolutos)

Para integrar um valor absoluto, corte onde o sinal muda.

02x1 ⁣dx=01(1x) ⁣dx+12(x1) ⁣dx=12+12=1,\int_0^2 \abs{x - 1}\,\dd x = \int_0^1 (1 - x)\,\dd x + \int_1^2 (x - 1)\,\dd x = \frac12 + \frac12 = 1 ,

e, cortando [0,2π]\intcc{0}{2\pi} em π\pi:

02πsint ⁣dt=0πsint ⁣dtπ2πsint ⁣dt=2+2=4,\int_0^{2\pi} \abs{\sin t}\,\dd t = \int_0^{\pi} \sin t\,\dd t - \int_{\pi}^{2\pi} \sin t\,\dd t = 2 + 2 = 4 ,

ao passo que 02πsint ⁣dt=0\int_0^{2\pi} \sin t\,\dd t = 0: o cancelamento é real, e é por isso que o enunciado de positividade estrita (Teorema 15.7 (4)) carrega a hipótese f0f \geq 0 — sem ela, uma integral nula nada demonstra sobre ff. A ideia de fechamento: f\int \abs f mede área, e f\int f mede balanço com sinal; a desigualdade ff\abs{\int f} \leq \int\abs f é o registro exato do que o cancelamento pode destruir.

15.3 O teorema fundamental do cálculo

Teorema 15.9 (Teorema fundamental do cálculo)

Seja ff contínua num intervalo II e aIa \in I. A função

F(x)=axf(t) ⁣dtF(x) = \int_a^x f(t)\, \dd t

é de classe C1C^1 em II, com F=fF' = f: toda função contínua num intervalo tem primitivas. Consequentemente, para qualquer primitiva GG de ff:

abf(t) ⁣dt=G(b)G(a).\int_a^b f(t)\,\dd t = G(b) - G(a) .

Demonstração. Fixe x0Ix_0 \in I e ε>0\varepsilon > 0; a continuidade em x0x_0 fornece δ\delta com f(t)f(x0)ε\abs{f(t) - f(x_0)} \leq \varepsilon para tx0δ\abs{t - x_0} \leq \delta. Para 0<hδ0 < \abs{h} \leq \delta (e x0+hIx_0 + h \in I), Chasles dá

F(x0+h)F(x0)hf(x0)=1hx0x0+h(f(t)f(x0)) ⁣dt,\frac{F(x_0 + h) - F(x_0)}{h} - f(x_0) = \frac 1h \int_{x_0}^{x_0+h} \bigl(f(t) - f(x_0)\bigr)\dd t ,

cujo valor absoluto é 1hhε=ε\leq \frac{1}{\abs h}\cdot \abs h\, \varepsilon = \varepsilon (estimativa (2), válida para qualquer ordem dos limites). Assim, F(x0)=f(x0)F'(x_0) = f(x_0); F=fF' = f é contínua: FF é C1C^1. Se G=fG' = f também, então (GF)=0(G - F)' = 0 no intervalo, de modo que G=F+cG = F + c (Corolário 14.12), e G(b)G(a)=F(b)F(a)=abfG(b) - G(a) = F(b) - F(a) = \int_a^b f.

Exemplo 15.10 (Simetria antes do cálculo)

Num intervalo simétrico, a paridade faz o trabalho: se ff é ímpar, a substituição ttt \mapsto -t leva a0f\int_{-a}^{0} f a 0af-\int_0^a f, de modo que

aaf(t) ⁣dt=0;se f eˊ par,aaf=20af.\int_{-a}^{a} f(t)\,\dd t = 0 ; \qquad\text{se $f$ é par,}\quad \int_{-a}^{a} f = 2\int_0^a f .

Assim, 11t3cost1+t4 ⁣dt=0\int_{-1}^{1} \frac{t^3\cos t}{1 + t^4}\,\dd t = 0 sem primitiva à vista (o integrando é ímpar), e ππt2cost ⁣dt=20πt2cost ⁣dt\int_{-\pi}^{\pi} t^2\cos t\,\dd t = 2\int_0^\pi t^2\cos t\, \dd t. Verifique a simetria antes de recorrer a técnicas: a integral mais rápida é a que nunca é calculada.

Exemplo 15.11 (Reconhecendo uma derivada de imediato)

Calcule 0π/2 ⁣dx1+cosx\displaystyle\int_0^{\pi/2} \frac{\dd x}{1 + \cos x}. A identidade do arco metade 1+cosx=2cos2x21 + \cos x = 2\cos^2\frac x2 transforma o integrando em 12(1+tan2x2)\frac{1}{2}\bigl(1 + \tan^2\frac x2\bigr), que é exatamente a derivada de tanx2\tan\frac x2:

0π/2 ⁣dx1+cosx=[tanx2]0π/2=tanπ4=1.\int_0^{\pi/2} \frac{\dd x}{1 + \cos x} = \Bigl[\tan\frac x2\Bigr]_0^{\pi/2} = \tan\frac\pi4 = 1 .

Nenhuma maquinaria de substituição foi necessária — apenas o reflexo de ler um integrando como a derivada de alguém, com o teorema fundamental fazendo o resto. (A ferramenta sistemática por trás de tais integrais trigonométricas, a substituição t=tanx2t = \tan\frac x2, pertence ao instrumental padrão construído a partir do Teorema 15.15 (2).)

Exemplo 15.12 (Funções definidas por integrais)

O teorema fundamental fabrica funções. Seja

F(x)=0xet2 ⁣dt.F(x) = \int_0^x \eu^{-t^2}\,\dd t .

Nenhuma combinação de funções clássicas tem derivada et2\eu^{-t^2} (um teorema de Liouville, admitido); e, no entanto, FF existe, é C1C^1 com F(x)=ex2>0F'(x) = \eu^{-x^2} > 0, estritamente crescente, ímpar (substitua ttt \mapsto -t) e limitada: para x1x \geq 1,

F(x)F(1)=1xet2 ⁣dt1xet ⁣dte1,F(x) - F(1) = \int_1^x \eu^{-t^2}\dd t \leq \int_1^x \eu^{-t}\dd t \leq \eu^{-1} ,

de modo que FF(1)+e11+e1F \leq F(1) + \eu^{-1} \leq 1 + \eu^{-1}. (O limite exato, π2\frac{\sqrt\pi}{2}, é calculado com integrais duplas no volume do terceiro ano.) Regra da cadeia para limites móveis:  ⁣d ⁣dxxx2et2 ⁣dt=2xex4ex2\frac{\dd}{\dd x}\int_x^{x^2} \eu^{-t^2}\dd t = 2x\,\eu^{-x^4} - \eu^{-x^2}. A ideia de fechamento: a integração cria funções novas a partir de antigas, com todas as suas propriedades legíveis no integrando — a primitiva que você não consegue escrever é ainda uma função que você controla plenamente.

Exemplo 15.13 (Estimar sem avaliar)

As integrais Rn=01tn1+t ⁣dtR_n = \int_0^1 \frac{t^n}{1 + t}\,\dd t não têm forma fechada agradável e, no entanto, a monotonicidade as fixa com precisão: em [0,1]\intcc{0}{1}, 1211+t1\frac12 \leq \frac{1}{1+t} \leq 1, de modo que

12(n+1)=1201tn ⁣dt    Rn    01tn ⁣dt=1n+1:\frac{1}{2(n+1)} = \frac12\int_0^1 t^n\,\dd t \;\leq\; R_n \;\leq\; \int_0^1 t^n \,\dd t = \frac{1}{n+1} :

a ordem exata de decaimento (RnR_n \sim um múltiplo de 1n\frac1n, de fato Rn12nR_n \sim \frac{1}{2n}) com duas linhas e sem primitiva alguma. Os problemas de fim de semana deste capítulo e do próximo rodam exatamente sobre enquadramentos desses — o primeiro instinto do analista diante de uma integral deve ser estimá-la e só depois, se necessário, calculá-la.

Exemplo 15.14 (Valores médios)

A média de uma ff contínua em [a,b]\intcc{a}{b} é 1baabf\frac{1}{b-a}\int_a^b f. Para o arco do seno:

1π0πsint ⁣dt=1π[cost]0π=2π0.637:\frac{1}{\pi}\int_0^\pi \sin t\,\dd t = \frac{1}{\pi}\bigl[-\cos t\bigr]_0^\pi = \frac{2}{\pi} \approx 0.637 :

um arco positivo completo tem média não 12\frac12, mas 2π\frac2\pi — a curva passa mais tempo alta do que um triângulo passaria. Pelo Exercício 15.11 (teorema do valor médio para integrais, g=1g = 1), a média é um valor: sinc=2π\sin c = \frac2\pi para algum c(0,π)c \in \intoo{0}{\pi}. E, pelas somas de Riemann deste capítulo, a média é o limite de médias comuns de nn amostras — a ponte entre a média discreta de dados e a média contínua de um sinal, que é como a integral entra na física.

Teorema 15.15 (Integração por partes; substituição)

  1. Se u,vu, v são C1C^1 em [a,b]\intcc{a}{b}:

    abuv=[uv]ababuv.\int_a^b u'v = \bigl[uv\bigr]_a^b - \int_a^b uv' .
  2. Se φ\varphi é C1C^1 em [α,β]\intcc{\alpha}{\beta} e ff é contínua em φ([α,β])\varphi(\intcc{\alpha}{\beta}):

    αβf(φ(t))φ(t) ⁣dt=φ(α)φ(β)f(x) ⁣dx.\int_{\alpha}^{\beta} f\bigl(\varphi(t)\bigr)\,\varphi'(t)\, \dd t = \int_{\varphi(\alpha)}^{\varphi(\beta)} f(x)\, \dd x .

Demonstração. (1) (uv)=uv+uv(uv)' = u'v + uv'; integre em [a,b]\intcc{a}{b} e aplique o teorema fundamental à função C1C^1 uvuv.

(2) Seja FF uma primitiva de ff no intervalo imagem (Teorema 15.9). Então (Fφ)=(fφ)φ(F \circ \varphi)' = (f \circ \varphi)\,\varphi' (regra da cadeia), de modo que os dois lados valem F(φ(β))F(φ(α))F(\varphi(\beta)) - F(\varphi(\alpha)).

Exemplo 15.16

01tet ⁣dt=[tet]0101et ⁣dt=e(e1)=1\int_0^1 t\,\eu^t \dd t = \bigl[t\,\eu^t\bigr]_0^1 - \int_0^1 \eu^t\dd t = \eu - (\eu - 1) = 1. E, com a substituição x=sintx = \sin t (t[0,π2]t \in \intcc{0}{\frac\pi2}):

011x2 ⁣dx=0π/2costcost ⁣dt=0π/21+cos2t2 ⁣dt=π4,\int_0^1 \sqrt{1 - x^2}\, \dd x = \int_0^{\pi/2} \cos t \cdot \cos t \, \dd t = \int_0^{\pi/2} \frac{1 + \cos 2t}{2}\, \dd t = \frac\pi4 ,

— um quarto do disco unitário, como a geometria exige. (A linearização do Método 3.11 em ação.)

Observação 15.17 (Armadilhas frequentes no cálculo integral)

(i) As substituições devem ser C1C^1 em todo o intervalo: a mudança x=1tx = \frac1t é ilegal atravessando 00; aplicada às cegas a 11 ⁣dx1+x2\int_{-1}^{1}\frac{\dd x}{1 + x^2}, ela “demonstra” que a integral é igual ao seu próprio oposto. Quando uma substituição tem singularidade, corte antes o intervalo (Chasles), substitua em cada peça e só então recombine. (ii) As primitivas logarítmicas precisam de valores absolutos:  ⁣dxx2=lnx2+C\int \frac{\dd x}{x - 2} = \ln\abs{x - 2} + C em cada lado de 22 separadamente — escrever ln(x2)\ln(x - 2) em (0,1)\intoo{0}{1} é escrever o logaritmo de um número negativo; e a constante CC pode diferir nos dois lados da singularidade. (iii) Uma integral nula não mata a função: 02πsin=0\int_0^{2\pi}\sin = 0; a positividade do integrando é exigida antes de concluir que f=0f = 0 (Exemplo 15.8). (iv) As somas de Riemann devem ser calibradas: em banf(a+kban)\frac{b-a}{n}\sum f\bigl(a + k\frac{b-a}{n}\bigr), o passo fora e os pontos dentro devem corresponder à mesma subdivisão — o erro frequente é uma soma k=1nf(kn)\sum_{k=1}^{n} f\bigl(\frac kn\bigr) sem o fator 1n\frac1n, que diverge em vez de convergir a 01f\int_0^1 f. Lista de verificação antes de invocar o Teorema 15.20: coloque 1n\frac1n em evidência, reescreva a parcela como ff de kn\frac kn, nomeie ff e verifique a sua continuidade.

Exemplo 15.18 (Adivinhe, derive, ajuste)

Quanto vale 1x(lnt)2 ⁣dt\int_1^x (\ln t)^2\,\dd t? Adivinhe uma primitiva da forma tP(lnt)t\,P(\ln t) com PP polinômio e derive:

(tP(lnt))=P(lnt)+P(lnt).\bigl(t\,P(\ln t)\bigr)' = P(\ln t) + P'(\ln t) .

Precisamos de P(u)+P(u)=u2P(u) + P'(u) = u^2: tome P(u)=u22u+2P(u) = u^2 - 2u + 2 (igualando coeficientes de cima para baixo a partir de u2u^2). Portanto,

1x(lnt)2 ⁣dt=[t((lnt)22lnt+2)]1x=x(lnx)22xlnx+2x2,\int_1^x (\ln t)^2\,\dd t = \bigl[t\bigl((\ln t)^2 - 2\ln t + 2\bigr)\bigr]_1^x = x(\ln x)^2 - 2x\ln x + 2x - 2 ,

resultado que de outro modo se alcançaria com duas integrações por partes. A ideia de fechamento: para integrandos da forma (polinômio em lnt\ln t) ou (polinômio vezes eλt\eu^{\lambda t}), a primitiva tem a mesma forma — derivar um palpite com forma converte a integração em álgebra linear sobre coeficientes, mais rápida e menos sujeita a erro que partes iteradas.

Exemplo 15.19 (A integral bumerangue)

Calcule I=0π/2excosx ⁣dxI = \int_0^{\pi/2} \eu^x \cos x\,\dd x. Integre por partes duas vezes, derivando o fator trigonométrico a cada vez:

I=[exsinx]0π/20π/2exsinx ⁣dx=eπ/2J,J=[excosx]0π/2+0π/2excosx ⁣dx=1+I.I = \bigl[\eu^x \sin x\bigr]_0^{\pi/2} - \int_0^{\pi/2} \eu^x \sin x\,\dd x = \eu^{\pi/2} - J, \qquad J = \bigl[-\eu^x\cos x\bigr]_0^{\pi/2} + \int_0^{\pi/2} \eu^x\cos x\,\dd x = 1 + I .

A integral voltou a si mesma: I=eπ/21II = \eu^{\pi/2} - 1 - I, donde

I=eπ/212.I = \frac{\eu^{\pi/2} - 1}{2} .

A ideia de fechamento: quando o integrando é um produto de duas funções que se reproduzem sob derivação (eax\eu^{ax}, cosbx\cos bx, sinbx\sin bx), duas integrações por partes produzem uma equação linear para a integral desconhecida — resolva-a em vez de integrar; equivalentemente, passe por e(a+ib)x\eu^{(a + \iu b)x} (Capítulo 3) e tome partes reais. Os dois caminhos dão a mesma resposta, e conferir que dão é um teste de sanidade gratuito.

15.4 Somas de Riemann

Teorema 15.20 (Somas de Riemann)

Seja ff contínua em [a,b]\intcc{a}{b}. Então

Sn=bank=0n1f(a+kban)nabf(t) ⁣dt,S_n = \frac{b - a}{n} \sum_{k=0}^{n-1} f\Bigl(a + k\,\frac{b-a}{n}\Bigr) \xrightarrow[n \to \infty]{} \int_a^b f(t)\, \dd t ,

e o mesmo com quaisquer pontos de avaliação dentro dos subintervalos.

Demonstração. SnS_n é a integral da função escada φn\varphi_n, igual a f(a+kban)f(a + k\frac{b-a}{n}) no kk-ésimo subintervalo. Dado ε>0\varepsilon > 0, a continuidade uniforme (Heine) fornece δ\delta; para n>baδn > \frac{b-a}{\delta}, todo ponto de um subintervalo está a menos de δ\delta do seu ponto de avaliação, de modo que fφnε\abs{f - \varphi_n} \leq \varepsilon em [a,b]\intcc{a}{b}, donde

abfSn=ab(fφn)(ba)ε.\Bigl| \int_a^b f - S_n \Bigr| = \Bigl| \int_a^b (f - \varphi_n) \Bigr| \leq (b-a)\,\varepsilon . \qedhere

Uma soma de Riemann à esquerda com n = 8 retângulos: quando a malha encolhe, a continuidade uniforme força a área da escada rumo a ∈t_ab f.
Uma soma de Riemann à esquerda com n=8n = 8 retângulos: quando a malha encolhe, a continuidade uniforme força a área da escada rumo a abf\int_a^b f.

Exemplo 15.21

k=1n1n+k=1nk=1n11+k/nn01 ⁣dx1+x=ln2\displaystyle\sum_{k=1}^{n} \frac{1}{n + k} = \frac 1n \sum_{k=1}^{n} \frac{1}{1 + k/n} \xrightarrow[n\to\infty]{} \int_0^1 \frac{\dd x}{1 + x} = \ln 2: um limite invisível a estimativas elementares, transparente como soma de Riemann.

Exemplo 15.22 (Uma segunda soma de Riemann, com calibração)

Encontre limnk=1nn(n+k)2\lim_{n\to\infty} \sum_{k=1}^{n} \dfrac{n}{(n+k)^2}. Calibre:

k=1nn(n+k)2=1nk=1nn2(n+k)2=1nk=1n1(1+kn)2,\sum_{k=1}^{n} \frac{n}{(n+k)^2} = \frac{1}{n}\sum_{k=1}^{n} \frac{n^2}{(n+k)^2} = \frac1n \sum_{k=1}^{n} \frac{1}{\bigl(1 + \frac kn\bigr)^2} ,

uma soma de Riemann da função contínua f(x)=1(1+x)2f(x) = \frac{1}{(1+x)^2} em [0,1]\intcc{0}{1}: o limite é

01 ⁣dx(1+x)2=[11+x]01=12.\int_0^1 \frac{\dd x}{(1 + x)^2} = \Bigl[-\frac{1}{1+x}\Bigr]_0^1 = \frac12 .

A ideia de fechamento: toda a arte está na linha do meio — force a parcela à forma f(kn)f(\frac kn) ao custo de extrair exatamente um fator 1n\frac1n; uma vez certa a forma, o teorema faz a análise e o teorema fundamental faz a aritmética.

Observação 15.23 (Onde a integral trabalha em seguida)

Cada construção do capítulo tem a sua sequência. As somas de Riemann voltam no Capítulo 17 como ponte entre séries e integrais (comparação de 1nα\sum \frac{1}{n^\alpha} com  ⁣dttα\int \frac{\dd t}{t^\alpha}); o resto integral é a forma mais fina da fórmula de Taylor no Capítulo 16; a definição baseada no sup\sup é o protótipo da integral de Lebesgue do volume do terceiro ano, em que as mesmas três propriedades (linearidade, monotonicidade, um teorema de convergência) são reconstruídas numa classe de funções muito maior. E o problema de fim de semana abaixo converte a integração por partes em aritmética: a irracionalidade de π2\pi^2.

15.5 Exercícios

Exercício 15.1

Calcule: 01 ⁣dxx24\displaystyle\int_0^1 \frac{\dd x}{x^2 - 4} (frações parciais, Capítulo 9); 1elnt ⁣dt\displaystyle\int_1^{\eu} \ln t \,\dd t; 0πtsint ⁣dt\displaystyle\int_0^{\pi} t \sin t\, \dd t.

Solução

Solução de Exercício 15.1.

1x24=1/4x21/4x+2\dfrac{1}{x^2-4} = \dfrac{1/4}{x - 2} - \dfrac{1/4}{x+2} (multiplicação-avaliação), de modo que

01 ⁣dxx24=14[lnx2lnx+2]01=14(ln13ln1)=ln34.\int_0^1 \frac{\dd x}{x^2 - 4} = \frac14\Bigl[\ln\abs{x-2} - \ln\abs{x+2}\Bigr]_0^1 = \frac14\Bigl(\ln\frac{1}{3} - \ln 1\Bigr) = -\frac{\ln 3}{4}.

Por partes (u=1u' = 1, v=lntv = \ln t): 1elnt ⁣dt=[tlnt]1e1e ⁣dt=e(e1)=1\int_1^{\eu} \ln t\,\dd t = \bigl[t\ln t\bigr]_1^{\eu} - \int_1^{\eu} \dd t = \eu - (\eu - 1) = 1.

Por partes (u=sintu' = \sin t, v=tv = t): 0πtsint ⁣dt=[tcost]0π+0πcost ⁣dt=π+0=π\int_0^\pi t\sin t\,\dd t = \bigl[-t\cos t\bigr]_0^\pi + \int_0^\pi \cos t\,\dd t = \pi + 0 = \pi.

Exercício 15.2

Calcule 01t(t2+1)2 ⁣dt\displaystyle\int_0^{1} \frac{t}{(t^2+1)^2}\,\dd t (substituição) e 0π/2cos3t ⁣dt\displaystyle\int_0^{\pi/2} \cos^3 t\, \dd t (escreva cos3=cos(1sin2)\cos^3 = \cos(1 - \sin^2)).

Solução

Solução de Exercício 15.2.

Substituição u=t2+1u = t^2 + 1,  ⁣du=2t ⁣dt\dd u = 2t\,\dd t:

01t ⁣dt(t2+1)2=1212 ⁣duu2=12[1u]12=14.\int_0^1 \frac{t\,\dd t}{(t^2+1)^2} = \frac12 \int_1^2 \frac{\dd u}{u^2} = \frac12\Bigl[-\frac1u\Bigr]_1^2 = \frac14 .

Com u=sintu = \sin t: 0π/2cos3t ⁣dt=0π/2(1sin2t)cost ⁣dt=[sintsin3t3]0π/2=113=23\int_0^{\pi/2} \cos^3 t\,\dd t = \int_0^{\pi/2} (1 - \sin^2 t)\cos t\,\dd t = \bigl[\sin t - \frac{\sin^3 t}{3}\bigr]_0^{\pi/2} = 1 - \frac13 = \frac23.

Exercício 15.3

Encontre os limites, como somas de Riemann:

un=k=1nnn2+k2,vn=1n(2n)!n!nnn (take logarithms).u_n = \sum_{k=1}^{n} \frac{n}{n^2 + k^2}, \qquad v_n = \frac{1}{n}\sqrt[n]{\frac{(2n)!}{n!\,n^n}} \ \emph{(take logarithms)}.
Solução

Solução de Exercício 15.3.

un=1nk=1n11+(k/n)2u_n = \dfrac1n \sum_{k=1}^{n} \dfrac{1}{1 + (k/n)^2}: uma soma de Riemann de x11+x2x \mapsto \frac{1}{1+x^2} em [0,1]\intcc{0}{1}, de modo que un01 ⁣dx1+x2=arctan1=π4u_n \to \int_0^1 \frac{\dd x}{1+x^2} = \arctan 1 = \dfrac\pi4.

lnvn=1nk=1nlnn+kn=1nk=1nln(1+kn)01ln(1+x) ⁣dx=[(1+x)ln(1+x)x]01=2ln21\ln v_n = \dfrac1n \sum_{k=1}^{n} \ln\dfrac{n+k}{n} = \dfrac 1n \sum_{k=1}^n \ln\Bigl(1 + \dfrac kn\Bigr) \to \int_0^1 \ln(1+x)\,\dd x = \bigl[(1+x)\ln(1+x) - x\bigr]_0^1 = 2\ln 2 - 1. Portanto, vne2ln21=4ev_n \to \eu^{2\ln 2 - 1} = \dfrac 4\eu. (Verificação da identificação: (2n)!n!nn=k=1nn+kn\frac{(2n)!}{n!\,n^n} = \prod_{k=1}^{n} \frac{n+k}{n}.)

Exercício 15.4

Seja ff contínua em [0,1]\intcc{0}{1}. Calcule limn01xnf(x) ⁣dx\lim_{n\to\infty} \int_0^1 x^n f(x)\,\dd x. (Corte [0,1]\intcc{0}{1} em 1δ1 - \delta.)

Solução

Solução de Exercício 15.4.

O limite é 00. Sejam M=supfM = \sup \abs f e ε(0,1)\varepsilon \in \intoo{0}{1}. Corte em 1ε1 - \varepsilon:

01xnf01εxnf+1ε1xnfM(1ε)n+Mε.\Bigl| \int_0^1 x^n f \Bigr| \leq \int_0^{1 - \varepsilon} x^n \abs f + \int_{1-\varepsilon}^1 x^n \abs f \leq M\,(1-\varepsilon)^n + M\varepsilon .

Como (1ε)n0(1 - \varepsilon)^n \to 0 (Exercício 11.3), o limite superior do lado esquerdo é Mε\leq M\varepsilon para todo ε\varepsilon: a integral tende a 00.

Exercício 15.5 ★★

(Cauchy–Schwarz) Para f,gf, g contínuas em [a,b]\intcc{a}{b}, demonstre que

(abfg) ⁣2abf2abg2,\Bigl(\int_a^b fg\Bigr)^{\!2} \leq \int_a^b f^2 \cdot \int_a^b g^2 ,

expandindo ab(f+λg)20\int_a^b (f + \lambda g)^2 \geq 0 como um polinômio do segundo grau em λ\lambda. Quando há igualdade?

Solução

Solução de Exercício 15.5.

Q(λ)=ab(f+λg)2=f2+2λfg+λ2g20Q(\lambda) = \int_a^b (f + \lambda g)^2 = \int f^2 + 2\lambda \int fg + \lambda^2 \int g^2 \geq 0 para todo λ\lambda. Se g2=0\int g^2 = 0, então g=0g = 0 (positividade estrita, Teorema 15.7 (4)), e a desigualdade é 000 \leq 0. Caso contrário, QQ é um genuíno polinômio do segundo grau, sempre 0\geq 0: o seu discriminante é 0\leq 0, isto é, (fg)2f2g2\bigl(\int fg\bigr)^2 \leq \int f^2 \int g^2.

Igualdade se, e somente se, o discriminante se anula, se, e somente se, Q(λ0)=0Q(\lambda_0) = 0 para algum λ0\lambda_0, isto é, (f+λ0g)2=0\int (f + \lambda_0 g)^2 = 0, ou seja (de novo por positividade estrita), f=λ0gf = -\lambda_0 g: a igualdade vale exatamente quando ff e gg são proporcionais.

Exercício 15.6 ★★

Seja ff contínua em R\R e TT-periódica. Demonstre que aa+Tf\int_a^{a+T} f não depende de aa, e que 1x0xf(t) ⁣dt1T0Tf\frac1x \int_0^x f(t)\,\dd t \to \frac 1T \int_0^T f quando x+x \to +\infty.

Solução

Solução de Exercício 15.6.

Seja Φ(a)=aa+Tf\Phi(a) = \int_a^{a+T} f. Pelo teorema fundamental (Teorema 15.9), Φ\Phi é derivável com Φ(a)=f(a+T)f(a)=0\Phi'(a) = f(a + T) - f(a) = 0: constante.

Para x>0x > 0, escreva x=nT+rx = nT + r, 0r<T0 \leq r < T (n=x/Tn = \lfloor x/T \rfloor). Chasles:

0xf=n0Tf+nTnT+rf,nTnT+rfTsup[0,T]f=C.\int_0^x f = n \int_0^T f + \int_{nT}^{nT + r} f, \qquad \Bigl| \int_{nT}^{nT+r} f \Bigr| \leq T \sup_{\intcc{0}{T}} \abs f = C .

Então 1x0xf=nTx1T0Tf+O(1x)\frac 1x \int_0^x f = \frac{nT}{x}\cdot\frac 1T \int_0^T f + O\bigl(\frac 1x\bigr), e nTx1\frac{nT}{x} \to 1: o limite é 1T0Tf\frac1T \int_0^T f.

Exercício 15.7 ★★

Para ff contínua em [0,1]\intcc{0}{1} com 01f=12\int_0^1 f = \frac12, demonstre que ff tem um ponto fixo em [0,1]\intcc{0}{1}. (Integre f(x)xf(x) - x e use a positividade estrita, Teorema 15.7 (4), pela sua contrapositiva combinada com o teorema do valor intermediário.)

Solução

Solução de Exercício 15.7.

Seja g(x)=f(x)xg(x) = f(x) - x: contínua, com

01g=01f12=0.\int_0^1 g = \int_0^1 f - \frac12 = 0 .

Se gg nunca se anulasse, o teorema do valor intermediário forçaria um sinal constante (uma função contínua num intervalo que assume os dois sinais se anula); digamos g>0g > 0. Então, pela positividade estrita (Teorema 15.7 (4) aplicada a g>0g > 0, dando g>0\int g > 0): contradição com g=0\int g = 0. Logo, g(c)=0g(c) = 0 para algum cc: f(c)=cf(c) = c.

Exercício 15.8 ★★★

(Integrais de Wallis) Seja Wn=0π/2sinnt ⁣dtW_n = \int_0^{\pi/2} \sin^n t\,\dd t.

  1. Demonstre a recorrência nWn=(n1)Wn2n W_n = (n-1) W_{n-2} (n2n \geq 2) por partes, e calcule W0,W1W_0, W_1, e depois W2pW_{2p} e W2p+1W_{2p+1} em forma fechada.
  2. Demonstre que (Wn)(W_n) é decrescente com Wn+1Wn1\frac{W_{n+1}}{W_n} \to 1; demonstre que a quantidade (n+1)Wn+1Wn(n+1)\,W_{n+1} W_n é constante, igual a π2\frac\pi2; e deduza a equivalência Wnπ2nW_n \sim \sqrt{\dfrac{\pi}{2n}}.
Solução

Solução de Exercício 15.8.

  1. Por partes, com u=sintu' = \sin t, v=sinn1tv = \sin^{n-1} t:

    Wn=[costsinn1t]0π/2+(n1)0π/2cos2tsinn2t ⁣dt=(n1)(Wn2Wn),W_n = \bigl[-\cos t\sin^{n-1}t\bigr]_0^{\pi/2} + (n-1)\int_0^{\pi/2} \cos^2 t\,\sin^{n-2} t\,\dd t = (n-1)(W_{n-2} - W_n),

    de modo que nWn=(n1)Wn2nW_n = (n-1)W_{n-2}. De W0=π2W_0 = \frac\pi2, W1=1W_1 = 1:

    W2p=(2p1)(2p3)1(2p)(2p2)2π2=(2p)!4p(p!)2π2,W2p+1=(2p)(2p2)2(2p+1)(2p1)3=4p(p!)2(2p+1)!.W_{2p} = \frac{(2p-1)(2p-3)\cdots 1}{(2p)(2p-2)\cdots 2}\, \frac{\pi}{2} = \frac{(2p)!}{4^p (p!)^2}\,\frac\pi2, \qquad W_{2p+1} = \frac{(2p)(2p-2)\cdots 2}{(2p+1)(2p-1)\cdots 3} = \frac{4^p (p!)^2}{(2p+1)!} .
  2. Em (0,π2)\intoo{0}{\frac\pi2}, 0<sint<10 < \sin t < 1, de modo que sinn+1<sinn\sin^{n+1} < \sin^n e (Wn)(W_n) é (estritamente) decrescente e positiva. Comprimindo com a recorrência:

    nn+1=Wn+1Wn1Wn+1Wn1    Wn+1Wn1.\frac{n}{n+1} = \frac{W_{n+1}}{W_{n-1}} \leq \frac{W_{n+1}}{W_n} \leq 1 \quad\implies\quad \frac{W_{n+1}}{W_n} \to 1 .

    Invariante: an=(n+1)Wn+1Wna_n = (n+1)W_{n+1}W_n satisfaz an=an1a_n = a_{n-1} pela recorrência (n+1)Wn+1=nWn1(n+1)W_{n+1} = nW_{n-1}, de modo que an=a0=1W1W0=π2a_n = a_0 = 1 \cdot W_1 W_0 = \frac\pi2. Então

    nWn2(n+1)Wn+1Wn=π2    Wnπ2n.n W_n^2 \sim (n+1) W_{n+1} W_n = \frac\pi2 \quad\implies\quad W_n \sim \sqrt{\frac{\pi}{2n}} .

Exercício 15.9 ★★★

(Niven: π\pi é irracional) Suponha π=ab\pi = \frac ab com a,bNa, b \in \N^*, e ponha, para um nn a ser escolhido,

P(x)=xn(abx)nn!,In=0πP(x)sinx ⁣dx.P(x) = \frac{x^n (a - bx)^n}{n!}, \qquad I_n = \int_0^{\pi} P(x)\sin x\, \dd x .
  1. Demonstre que 0<Inπ(πa)nn!0 < I_n \leq \pi\,\frac{(\pi a)^n}{n!}, que é <1< 1 para nn grande.
  2. Demonstre que PP e todas as suas derivadas assumem valores inteiros em 00 e em π=ab\pi = \frac ab. (Expansão binomial: os coeficientes de PP vezes k!k! são inteiros para knk \geq n; e P(πx)=P(x)P(\pi - x) = P(x).)
  3. Ponha Q=PP+P(4)Q = P - P'' + P^{(4)} - \dots (uma soma finita). Verifique que (QsinxQcosx)=Psinx\bigl(Q'\sin x - Q\cos x\bigr)' = P \sin x, e deduza que In=Q(π)+Q(0)I_n = Q(\pi) + Q(0) é um inteiro.
  4. Conclua.
Solução

Solução de Exercício 15.9.

  1. Em (0,π)\intoo{0}{\pi}: x>0x > 0, abx=b(abx)=b(πx)>0a - bx = b(\frac ab - x) = b(\pi - x) > 0 e sinx>0\sin x > 0, de modo que o integrando é >0> 0 e In>0I_n > 0 (positividade estrita). Estimativa: em [0,π]\intcc{0}{\pi}, xπx \leq \pi e abxaa - bx \leq a, de modo que Pπnann!P \leq \frac{\pi^n a^n}{n!} e Inπ(πa)nn!I_n \leq \pi\,\frac{(\pi a)^n}{n!}, que tende a 00 (o fatorial vence o termo geométrico: é o termo geral da série exponencial convergente, cf. Exemplo 11.12); em particular, In<1I_n < 1 para nn grande.
  2. Expanda xn(abx)n=j=0n(nj)anj(b)jxn+jx^n(a - bx)^n = \sum_{j=0}^{n} \binom nj a^{n-j} (-b)^j x^{n+j}: assim, P=1n!jcjxn+jP = \frac{1}{n!}\sum_j c_j x^{n+j} com cjc_j inteiros. Então P(k)(0)=0P^{(k)}(0) = 0 para k<nk < n (valoração) e, para nk2nn \leq k \leq 2n, P(k)(0)=k!n!cknP^{(k)}(0) = \frac{k!}{n!} c_{k-n}, um inteiro, pois n!k!n! \mid k!. Além disso, P(πx)=P(x)P(\pi - x) = P(x) (substitua: πx\pi - x troca os fatores, usando ab(πx)=bxa - b(\pi - x) = bx), de modo que P(k)(π)=±P(k)(0)P^{(k)}(\pi) = \pm P^{(k)}(0): também inteiros.
  3. Com Q=PP+P(4)Q = P - P'' + P^{(4)} - \dots (finita: PP tem grau 2n2n): Q+Q=PQ + Q'' = P, e

    (QsinxQcosx)=(Q+Q)sinx=Psinx.\bigl(Q'\sin x - Q\cos x\bigr)' = (Q + Q'')\sin x = P\sin x .

    Portanto, In=[QsinxQcosx]0π=Q(π)+Q(0)I_n = \bigl[Q'\sin x - Q\cos x\bigr]_0^{\pi} = Q(\pi) + Q(0), uma soma de valores P(2k)P^{(2k)} em 00 e em π\pi: um inteiro, por (2).

  4. Para nn grande, InI_n é um inteiro com 0<In<10 < I_n < 1: impossível. A suposição π=ab\pi = \frac ab falha: π\pi é irracional.

Exercício 15.10 ★★★

Seja ff de classe C1C^1 em [a,b]\intcc{a}{b}. Demonstre o limite do tipo Riemann–Lebesgue

abf(t)sin(λt) ⁣dtλ+0\int_a^b f(t)\sin(\lambda t)\,\dd t \xrightarrow[\lambda \to +\infty]{} 0

integrando por partes. Depois demonstre-o de novo para ff apenas contínua, por aproximação uniforme com funções escada (Teorema 15.3).

Solução

Solução de Exercício 15.10.

Caso C1C^1: por partes,

abf(t)sinλt ⁣dt=[f(t)cosλtλ]ab+1λabf(t)cosλt ⁣dt,\int_a^b f(t)\sin\lambda t\,\dd t = \Bigl[-f(t)\frac{\cos\lambda t}{\lambda}\Bigr]_a^b + \frac{1}{\lambda}\int_a^b f'(t)\cos\lambda t\,\dd t,

limitado em valor absoluto por 2supf+(ba)supfλ0\frac{2\sup\abs f + (b - a)\sup\abs{f'}} {\lambda} \to 0.

Caso contínuo: seja ε>0\varepsilon > 0 e tome uma função escada φ\varphi com fφε\abs{f - \varphi} \leq \varepsilon (o Teorema 15.3 fornece φfψ\varphi \leq f \leq \psi com folga ε\leq\varepsilon; tome φ\varphi). Então

abfsinλtabfφ+abφsinλt(ba)ε+icixi1xisinλt ⁣dt,\Bigl|\int_a^b f\sin\lambda t\Bigr| \leq \int_a^b \abs{f - \varphi} + \Bigl|\int_a^b \varphi \sin\lambda t\Bigr| \leq (b-a)\varepsilon + \sum_i \abs{c_i}\,\Bigl|\int_{x_{i-1}}^{x_i} \sin\lambda t\,\dd t\Bigr| ,

e cada sinλt=cosλxi1cosλxiλ2λ\bigl|\int \sin \lambda t\bigr| = \bigl|\frac{\cos\lambda x_{i-1} - \cos\lambda x_i}{\lambda}\bigr| \leq \frac{2}{\lambda}: o segundo termo tende a 00. Portanto, o limite superior é (ba)ε\leq (b-a)\varepsilon para todo ε\varepsilon: o limite é 00.

Exercício 15.11 ★★

(Teorema do valor médio para integrais) Sejam f,gf, g contínuas em [a,b]\intcc{a}{b} com g0g \geq 0. Demonstre que existe c[a,b]c \in \intcc{a}{b} com

abf(t)g(t) ⁣dt=f(c)abg(t) ⁣dt,\int_a^b f(t)\,g(t)\,\dd t = f(c)\int_a^b g(t)\,\dd t ,

e mostre, com um exemplo, que a hipótese g0g \geq 0 não pode ser abandonada.

Solução

Solução de Exercício 15.11.

Sejam m=minfm = \min f e M=maxfM = \max f, atingidos pelo teorema de Weierstrass. Como g0g \geq 0: mgfgMgm\,g \leq fg \leq M\,g, de modo que, por monotonicidade,

mabg    abfg    Mabg.m \int_a^b g \;\leq\; \int_a^b fg \;\leq\; M \int_a^b g .

Se abg=0\int_a^b g = 0: a positividade estrita (Teorema 15.7 (4)) força g0g \equiv 0, os dois lados se anulam, e qualquer cc serve. Caso contrário, t=fggt = \frac{\int fg}{\int g} está em [m,M]=f([a,b])\intcc{m}{M} = f(\intcc{a}{b}) (Teoremas 13.13 e 13.10), de modo que t=f(c)t = f(c) para algum cc.

O sinal importa: em [1,1]\intcc{-1}{1}, com f(t)=g(t)=tf(t) = g(t) = t: fg=11t2=23\int fg = \int_{-1}^1 t^2 = \frac23, ao passo que f(c)11t ⁣dt=0f(c)\int_{-1}^1 t\,\dd t = 0 para todo cc.

Exercício 15.12 ★★★

(Os momentos forçam zeros) Seja ff contínua em [a,b]\intcc{a}{b} com

abf(t)tk ⁣dt=0parak=0,1,,n.\int_a^b f(t)\,t^k\,\dd t = 0 \qquad \text{para} k = 0, 1, \dots, n .

Demonstre que ff se anula em n+1n + 1 pontos distintos de (a,b)\intoo{a}{b}. (Se ff muda de sinal apenas em z1<<zmz_1 < \dots < z_m com mnm \leq n, integre ff contra P(t)=(tz1)(tzm)P(t) = (t - z_1) \cdots (t - z_m) e use a positividade estrita.)

Solução

Solução de Exercício 15.12.

Se f0f \equiv 0, a afirmação é vazia (todo ponto é zero). Suponha, então, f≢0f \not\equiv 0 e que ela tenha no máximo nn zeros distintos em (a,b)\intoo{a}{b}; sejam z1<<zmz_1 < \dots < z_m (mnm \leq n) os zeros em que ff muda de sinal (possivelmente nenhum). Ponha P(t)=i=1m(tzi)P(t) = \prod_{i=1}^{m}(t - z_i) (produto vazio =1= 1), de grau mnm \leq n. Em cada subintervalo cortado pelos ziz_i, tanto ff quanto PP têm sinal constante, e ambos trocam de sinal ao cruzar algum ziz_i: o produto fPfP tem um único sinal constante em todo (a,b)\intoo{a}{b}. Sendo contínuo, não identicamente nulo e de sinal constante, ele tem abfP>0\bigl|\int_a^b fP\bigr| > 0 (positividade estrita aplicada a fP\abs{fP}). Mas fP\int f P é uma combinação linear dos momentos ftk\int f\,t^k, knk \leq n, todos nulos: contradição. Portanto, ff tem ao menos n+1n + 1 zeros distintos em (a,b)\intoo{a}{b}.

Observação 15.24 (Perspectivas dentro deste volume)

Três capítulos adiante se apoiam diretamente neste. O Capítulo 16 carrega o resto integral — a mais fina das três fórmulas de Taylor é uma integração por partes iterada nn vezes. O Capítulo 17 converte o enquadramento de somas por integrais no teste decisivo para nα\sum n^{-\alpha}, e o seu problema de fim de semana refina esse enquadramento na constante de Euler. O Capítulo 24 torna a integral geométrica: o comprimento de um arco parametrizado é x(t)2+y(t)2 ⁣dt\int \sqrt{x'(t)^2 + y'(t)^2}\,\dd t, integral de uma função contínua num segmento — precisamente o objeto aqui construído, sem necessidade de teoria imprópria alguma. O fato mais reutilizado será o mais humilde: f(ba)supf\bigl|\int f\bigr| \leq (b - a)\sup\abs{f}, a desigualdade que converte toda estimativa pontual numa estimativa integral.

15.6 Problema: A máquina integral de irracionalidade

Problema 15.1

Problema de fim de semana — e\eu e π2\pi^2 são irracionais, e\eu com seis casas decimais, e 227>π\frac{22}{7} > \pi com demonstração

Um único mecanismo alimenta todo este problema: uma expressão que deve ser um inteiro positivo e que, no entanto, é comprovadamente menor que 11, não pode existir. O Exercício 15.9 (Niven) o rodou uma vez para demonstrar que πQ\pi \notin \Q; aqui nós o industrializamos. A máquina precisa de três partes: uma entrada de integralidade (valores nas extremidades de polinômios bem escolhidos), uma entrada de pequenez (um fator 1n!\frac{1}{n!} esmagando a integral) e uma ponte (a integração por partes) que as conecta. Demonstramos que e\eu é irracional e o calculamos com erro certificado, demonstramos o teorema mais fino de Legendre de que π2\pi^2 é irracional, e terminamos com a integral mais encantadora da análise: 01x4(1x)41+x2 ⁣dx=227π\int_0^1 \frac{x^4(1 - x)^4}{1 + x^2}\dd x = \frac{22}{7} - \pi, que enquadra π\pi à mão.

Parte I — Combustível.

  1. Demonstre que cnn!0\dfrac{c^{\,n}}{n!} \to 0 para todo c>0c > 0 fixado (a partir de n2cn \geq 2c, cada passo divide ao menos ao meio o termo).
  2. (Integrais beta) Demonstre, por indução em mm com integração por partes:

    01xk(1x)m ⁣dx=k!m!(k+m+1)!(k,mN).\int_0^1 x^{\,k}\,(1 - x)^{\,m}\,\dd x = \frac{k!\,m!}{(k + m + 1)!} \qquad (k, m \in \N).
  3. Deduza que 01(x(1x))n ⁣dx=1(2n+1)(2nn)\displaystyle\int_0^1 \bigl(x(1-x)\bigr)^n \dd x = \frac{1}{(2n+1)\binom{2n}{n}} e — comparando com a estimativa x(1x)14x(1 - x) \leq \frac14 — a estimativa (2nn)4n2n+1\binom{2n}{n} \geq \dfrac{4^n}{2n+1}, que corresponde a (2nn)1/n4\binom{2n}{n}^{1/n} \to 4 do Problema 11.1.
  4. Demonstre o lema da pequenez usado duas vezes abaixo: para toda gg contínua e >0> 0 em (0,1)\intoo{0}{1},

    0<01(x(1x))ng(x) ⁣dxsup[0,1]g4n.0 < \int_0^1 \bigl(x(1-x)\bigr)^n g(x)\,\dd x \leq \frac{\sup_{\intcc{0}{1}}\abs g}{4^{\,n}} .

Parte II — e\eu: irracionalidade e depois seis casas decimais. Ponha An=01xnex ⁣dxA_n = \displaystyle\int_0^1 x^n \eu^x \dd x.

  1. Calcule A0A_0 e A1A_1, demonstre a recorrência An=enAn1A_n = \eu - n A_{n-1} e as estimativas 0<Anen+10 < A_n \leq \dfrac{\eu}{n+1}.
  2. Mostre por indução que An=αn+βneA_n = \alpha_n + \beta_n \eu com αn,βnZ\alpha_n, \beta_n \in \Z.
  3. Deduza que e\eu é irracional (se e=pq\eu = \frac pq, então qAnq A_n é um inteiro preso em (0,1)\intoo{0}{1} para nn grande). Compare com a demonstração do Exercício 11.9: mesma conclusão, combustível diferente.
  4. Demonstre, por indução e integração por partes, a fórmula exata do resto

    e=k=0n1k!+Rn,Rn=1n!01(1t)net ⁣dt,1(n+1)!Rne(n+1)!.\eu = \sum_{k=0}^{n} \frac{1}{k!} + R_n, \qquad R_n = \frac{1}{n!}\int_0^1 (1 - t)^n\,\eu^{\,t}\,\dd t, \qquad \frac{1}{(n+1)!} \leq R_n \leq \frac{\eu}{(n+1)!} .
  5. Tome n=9n = 9: estime R9R_9 usando e<2.75\eu < 2.75 (a partir de b2=2.75b_2 = 2.75 no Exemplo 11.12), avalie a soma e conclua o enquadramento certificado 2.7182818e2.71828232.7182818 \leq \eu \leq 2.7182823 — seis casas decimais, e2.718282\eu \approx 2.718282, com demonstração.

Parte III — O teorema de Legendre: π2\pi^2 é irracional. Seja f(x)=xn(1x)nn!f(x) = \dfrac{x^n (1 - x)^n}{n!}, e suponha que π2=ab\pi^2 = \frac ab com a,bNa, b \in \N^*.

  1. Mostre que f(1x)=f(x)f(1 - x) = f(x) e 0<f14nn!0 < f \leq \dfrac{1}{4^n\,n!} em (0,1)\intoo{0}{1}.
  2. Mostre que f(k)(0)f^{(k)}(0) e f(k)(1)f^{(k)}(1) são inteiros para todo k0k \geq 0 (expanda xn(1x)nx^n(1-x)^n com coeficientes inteiros; k!n!Z\frac{k!}{n!} \in \Z para knk \geq n; depois use a simetria).
  3. Defina

    G=bnk=0n(1)kπ2n2kf(2k).G = b^{\,n} \sum_{k=0}^{n} (-1)^k\, \pi^{2n - 2k} f^{(2k)} .

    Mostre que G(0)G(0) e G(1)G(1) são inteiros (cada bnπ2n2k=ankbkb^n \pi^{2n-2k} = a^{\,n-k}\,b^{\,k}).

  4. Verifique o telescopamento G+π2G=bnπ2n+2f=π2anfG'' + \pi^2 G = b^n \pi^{2n+2} f = \pi^2 a^n f e depois

     ⁣d ⁣dx(G(x)sinπxπG(x)cosπx)=π2anf(x)sinπx.\frac{\dd}{\dd x}\Bigl(G'(x)\sin \pi x - \pi\,G(x)\cos \pi x\Bigr) = \pi^2 a^n f(x)\sin \pi x .
  5. Integre em [0,1]\intcc{0}{1} e conclua que

    πan01f(x)sin(πx) ⁣dx=G(0)+G(1)Z,\pi\,a^n \int_0^1 f(x)\sin(\pi x)\,\dd x = G(0) + G(1) \in \Z ,

    é um inteiro positivo limitado por πan4nn!\dfrac{\pi a^n}{4^n\,n!}.

  6. Conclua, com a questão 1, que π2\pi^2 é irracional (Legendre, 1794) e que isso reforça o Exercício 15.9: por que a irracionalidade de π2\pi^2 implica a de π\pi, e não o contrário?

Parte IV — Entendendo a máquina.

  1. Localize as duas forças opostas (integralidade dos dados nas extremidades; pequenez analítica da integral) e a ponte, nas Partes II e III. Depois explique por que o fator 1n!\frac{1}{n!} em ff é o ponto crucial: se ele for removido, a integralidade sobrevive, mas que desigualdade morre, e para que frações ab\frac ab alegadas a demonstração passa então a falhar?
  2. Efetividade: suponha que alguém alegue que π2=ab\pi^2 = \frac ab com a10a \leq 10. Mostre que a contradição já aparece em n=7n = 7: calcule π(10/4)7/7!0.38<1\pi\,(10/4)^7/7! \approx 0.38 < 1. A máquina não apenas refuta; ela refuta num estágio fixo e computável.
  3. Verifique a ponte incondicionalmente: demonstre, por duas integrações por partes, que

    01x(1x)sin(πx) ⁣dx=4π3,\int_0^1 x(1 - x)\sin(\pi x)\,\dd x = \frac{4}{\pi^3},

    e reconcilie com a questão 14 em n=1n = 1 (mantenha π2\pi^2 simbólico: a identidade telescopada diz π301f1sinπx=(f1(0)+f1(1))=4\pi^3 \int_0^1 f_1 \sin \pi x = -(f_1''(0) + f_1''(1)) = 4).

  4. O que torna ex\eu^x e sinπx\sin \pi x elegíveis como núcleos da máquina? Identifique a propriedade (cada uma satisfaz uma equação diferencial linear com coeficientes constantes, de modo que integrações por partes repetidas retornam ao início) e nomeie a fronteira: a mesma máquina, refinada por Hermite e Lindemann, demonstra que e\eu e π\pi são transcendentes — além deste volume.

Parte V — 227\frac{22}{7} contra π\pi, e a moral.

  1. Estabeleça a divisão polinomial

    x4(1x)41+x2=x64x5+5x44x2+441+x2,\frac{x^4(1-x)^4}{1 + x^2} = x^6 - 4x^5 + 5x^4 - 4x^2 + 4 - \frac{4}{1 + x^2},

    e deduza a célebre identidade

    01x4(1x)41+x2 ⁣dx=227π.\int_0^1 \frac{x^4 (1-x)^4}{1 + x^2}\,\dd x = \frac{22}{7} - \pi .
  2. O integrando é positivo: conclua que π<227\pi < \frac{22}{7}. Depois, limitando 11+x2\frac{1}{1+x^2} entre 12\frac12 e 11 e usando 01(x(1x))4=1630\int_0^1 (x(1-x))^4 = \frac{1}{630} (questão 3), demonstre que

    2271630    π    22711260,\frac{22}{7} - \frac{1}{630} \;\leq\; \pi \;\leq\; \frac{22}{7} - \frac{1}{1260} ,

    isto é, 3.14126π3.142073.14126 \leq \pi \leq 3.14207: duas casas decimais corretas, à mão.

  3. Generalize: dividindo x4m(1x)4mx^{4m}(1-x)^{4m} por 1+x21 + x^2, mostre que o resto é a constante (4)m(-4)^m (trabalhe módulo x2+1x^2 + 1: (1x)22x(1-x)^2 \equiv -2x), deduza racionais rmr_m com

    πrm415m,\abs{\pi - r_m} \leq 4^{\,1 - 5m} ,

    e verifique que m=1m = 1 reproduz as questões 20–21.

  4. Confronte esses racionais com a teoria de aproximação do Problema 14.1: calcule π2271.26103\abs{\pi - \frac{22}{7}} \approx 1.26\cdot10^{-3} contra a garantia de Dirichlet 149\frac{1}{49}, e cite π3551132.7107\abs{\pi - \frac{355}{113}} \approx 2.7\cdot10^{-7} contra 111327.8105\frac{1}{113^2} \approx 7.8\cdot10^{-5}: existem aproximações racionais excepcionalmente boas para π\pi — coerente, pois não se sabe que π\pi seja mal aproximável.
  5. (A armadilha do inteiro, abstraída) Demonstre o lema que unifica tudo: se xRx \in \R e existem inteiros an,bna_n, b_n com 0<an+bnx00 < \abs{a_n + b_n x} \to 0, então xx é irracional. Liste as suas instâncias neste problema, no Exercício 15.9, no Exercício 11.9 e no Problema 14.1.
  6. Síntese, uma frase para cada: (i) as três partes da máquina e onde cada uma vive no instrumental deste capítulo; (ii) o que a integral contribui que o teorema do valor médio do Problema 14.1 não podia; (iii) o inventário dos resultados extraídos (duas irracionalidades, uma constante com seis casas decimais, um enquadramento de π\pi, uma estimativa binomial); (iv) a fronteira (Hermite, Lindemann; e a mesma armadilha, rodada em ζ(2)\zeta(2) e ζ(3)\zeta(3), na aritmética do século XX).
Solução

Solução de Problema 15.1.

1. Seja N=2cN = \lceil 2c \rceil. Para nNn \geq N: cn+1/(n+1)!cn/n!=cn+112\frac{c^{n+1}/(n+1)!}{c^n/n!} = \frac{c}{n+1} \leq \frac12, de modo que 0<cnn!cNN!2(nN)00 < \frac{c^n}{n!} \leq \frac{c^N}{N!}\,2^{-(n - N)} \to 0: confronto.

2. Fixe kk; indução em mm. Para m=0m = 0: 01xk=1k+1=k!0!(k+1)!\int_0^1 x^k = \frac{1}{k+1} = \frac{k!\,0!}{(k+1)!}. Passo, por partes (u=(1x)mu = (1-x)^m, v=xkv' = x^k):

01xk(1x)m ⁣dx=[xk+1k+1(1x)m]01+mk+101xk+1(1x)m1 ⁣dx=mk+1(k+1)!(m1)!(k+m+1)!,\int_0^1 x^k (1-x)^m \dd x = \Bigl[\frac{x^{k+1}}{k+1}(1-x)^m\Bigr]_0^1 + \frac{m}{k+1}\int_0^1 x^{k+1}(1-x)^{m-1}\dd x = \frac{m}{k+1}\cdot\frac{(k+1)!\,(m-1)!}{(k+m+1)!} ,

que é k!m!(k+m+1)!\frac{k!\,m!}{(k+m+1)!}.

3. k=m=nk = m = n: 01(x(1x))n=(n!)2(2n+1)!=1(2n+1)(2nn)\int_0^1 (x(1-x))^n = \frac{(n!)^2}{(2n+1)!} = \frac{1}{(2n+1)\binom{2n}{n}}. Como x(1x)14x(1-x) \leq \frac14 em [0,1]\intcc{0}{1}, a integral é 4n\leq 4^{-n}, donde (2nn)4n2n+1\binom{2n}{n} \geq \frac{4^n}{2n+1} — coerente com (2nn)1/n4\binom{2n}{n}^{1/n} \to 4 (Problema 11.1).

4. O integrando é contínuo, 0\geq 0 e positivo em (0,1)\intoo{0}{1}, logo não identicamente nulo: a sua integral é >0> 0 (Teorema 15.7 (4)). Estimativa superior: (x(1x))n4n(x(1-x))^n \leq 4^{-n} e gsupgg \leq \sup\abs g, e depois monotonicidade.

5. A0=e1A_0 = \eu - 1; A1=[xex]0101ex=e(e1)=1A_1 = [x\eu^x]_0^1 - \int_0^1 \eu^x = \eu - (\eu - 1) = 1. Por partes: An=[xnex]01n01xn1ex=enAn1A_n = [x^n \eu^x]_0^1 - n\int_0^1 x^{n-1}\eu^x = \eu - n A_{n-1}. Estimativas: o integrando é positivo, de modo que An>0A_n > 0; e exe\eu^x \leq \euAne01xn=en+1A_n \leq \eu\int_0^1 x^n = \frac{\eu}{n+1}.

6. A0=1+1eA_0 = -1 + 1\cdot\eu. Se An1=αn1+βn1eA_{n-1} = \alpha_{n-1} + \beta_{n-1}\eu com entradas inteiras, então

An=enαn1nβn1e=(nαn1)αn+(1nβn1)βne,A_n = \eu - n\alpha_{n-1} - n\beta_{n-1}\eu = \underbrace{(-n\,\alpha_{n-1})}_{\alpha_n} + \underbrace{(1 - n\,\beta_{n-1})}_{\beta_n}\,\eu ,

ambos inteiros.

7. Se e=pq\eu = \frac pq: qAn=qαn+pβnZq A_n = q\alpha_n + p\beta_n \in \Z, e 0<qAnqen+1<10 < qA_n \leq \frac{q\eu}{n+1} < 1 para nn grande: um inteiro estritamente entre 00 e 11 — impossível. Logo, eQ\eu \notin \Q. No Exercício 11.9, o inteiro preso era q!pqq!aqq!\,\frac pq - q!\,a_q; aqui é qAnqA_n: mesma armadilha, combustível integral.

8. n=0n = 0: R0=01et=e1R_0 = \int_0^1 \eu^t = \eu - 1, de modo que e=1+R0\eu = 1 + R_0. Por partes (u=etu = \eu^t, v=(1t)n+1n+1v = -\frac{(1-t)^{n+1}}{n+1}):

Rn=1n!([(1t)n+1n+1et]01+1n+101(1t)n+1et)=1(n+1)!+Rn+1,R_n = \frac{1}{n!}\Bigl(\Bigl[-\frac{(1-t)^{n+1}}{n+1} \eu^t\Bigr]_0^1 + \frac{1}{n+1}\int_0^1 (1-t)^{n+1}\eu^t\Bigr) = \frac{1}{(n+1)!} + R_{n+1} ,

de modo que a fórmula se propaga de nn a n+1n + 1. Estimativas: 1ete1 \leq \eu^t \leq \eu em [0,1]\intcc{0}{1} e 01(1t)n=1n+1\int_0^1 (1-t)^n = \frac{1}{n+1} dão 1(n+1)!Rne(n+1)!\frac{1}{(n+1)!} \leq R_n \leq \frac{\eu}{(n+1)!}.

9. k=091k!=986410362880=2.71828152\sum_{k=0}^{9} \frac{1}{k!} = \frac{986410}{362880} = 2.71828152\dots, e

110!=2.76107R92.7510!=7.58107,\frac{1}{10!} = 2.76\cdot10^{-7} \leq R_9 \leq \frac{2.75}{10!} = 7.58\cdot10^{-7} ,

de modo que 2.7182818e2.71828232.7182818 \leq \eu \leq 2.7182823: com demonstração, e=2.718282\eu = 2.718282 com seis casas decimais (valor verdadeiro 2.71828182.7182818\dots).

10. f(1x)=(1x)nxnn!=f(x)f(1 - x) = \frac{(1-x)^n x^n}{n!} = f(x). Em (0,1)\intoo{0}{1}: 0<x(1x)140 < x(1-x) \leq \frac14, de modo que 0<f14nn!0 < f \leq \frac{1}{4^n\,n!}.

11. xn(1x)n=j=0n(1)j(nj)xn+jx^n(1-x)^n = \sum_{j=0}^{n} (-1)^j\binom nj\,x^{n+j}, de modo que f=1n!jcjxn+jf = \frac{1}{n!}\sum_j c_j\,x^{n+j} com cjZc_j \in \Z. Portanto, f(k)(0)=0f^{(k)}(0) = 0 para k<nk < n ou k>2nk > 2n, e, para nk2nn \leq k \leq 2n: f(k)(0)=k!n!cknf^{(k)}(0) = \frac{k!}{n!}\, c_{k-n}, um inteiro, pois n!k!n! \mid k!. A simetria dá f(k)(1)=(1)kf(k)(0)Zf^{(k)}(1) = (-1)^k f^{(k)}(0) \in \Z.

12. bnπ2n2k=bn(ab)nk=ankbkZb^n \pi^{2n-2k} = b^n\bigl(\frac ab\bigr)^{\,n-k} = a^{\,n-k}\,b^{\,k} \in \Z, de modo que G(0)=k(1)kankbkf(2k)(0)G(0) = \sum_k (-1)^k a^{n-k} b^k f^{(2k)}(0) e, do mesmo modo, G(1)G(1) são inteiros pela questão 11.

13. Em π2G+G\pi^2 G + G'', o termo kk de π2G\pi^2 G carrega π2n2k+2f(2k)\pi^{2n-2k+2}f^{(2k)}, e o termo j=k1j = k - 1 de GG'' carrega (1)k1π2n2k+2f(2k)(-1)^{k-1}\pi^{2n-2k+2}f^{(2k)}: tudo se cancela, exceto k=0k = 0 na primeira soma e j=nj = n na segunda, isto é,

G+π2G=bn(π2n+2f+(1)nf(2n+2))=bnπ2n+2f=π2anfG'' + \pi^2 G = b^n\bigl(\pi^{2n+2} f + (-1)^n f^{(2n+2)}\bigr) = b^n \pi^{2n+2} f = \pi^2 a^n f

(ff tem grau 2n2n, de modo que f(2n+2)=0f^{(2n+2)} = 0; e bnπ2n=anb^n\pi^{2n} = a^n). Então

(GsinπxπGcosπx)=(G+π2G)sinπx=π2anfsinπx.\bigl(G'\sin\pi x - \pi G\cos\pi x\bigr)' = (G'' + \pi^2 G)\sin \pi x = \pi^2 a^n f\sin\pi x .

14. Integrando em [0,1]\intcc{0}{1}:

π2an01fsin(πx) ⁣dx=[GsinπxπGcosπx]01=π(G(1)+G(0)),\pi^2 a^n \int_0^1 f\sin(\pi x)\,\dd x = \bigl[G'\sin\pi x - \pi G\cos\pi x\bigr]_0^1 = \pi\bigl(G(1) + G(0)\bigr) ,

de modo que πan01fsinπx=G(0)+G(1)Z\pi a^n \int_0^1 f\sin\pi x = G(0) + G(1) \in \Z. Em (0,1)\intoo{0}{1}, f>0f > 0 e sinπx>0\sin\pi x > 0: o lado esquerdo é positivo, logo G(0)+G(1)1G(0) + G(1) \geq 1; e sin1\sin \leq 1, com a questão 10, o limita por πan4nn!\frac{\pi a^n}{4^n\,n!}.

15. Pela questão 1 (com c=a4c = \frac a4), πan4nn!0\frac{\pi a^n}{4^n n!} \to 0: para nn grande, ele é <1< 1, contradizendo G(0)+G(1)1G(0) + G(1) \geq 1. Assim, fração ab\frac ab alguma é igual a π2\pi^2: o teorema de Legendre. Se π\pi fosse racional, π2\pi^2 também seria: logo, πQ\pi \notin \Q — e a implicação só corre nesse sentido (2\sqrt2 é irracional com quadrado racional), e é por isso que π2Q\pi^2 \notin \Q é estritamente mais forte que o Exercício 15.9.

16. Integralidade: questões 11–12 (derivadas nas extremidades); pequenez: questões 10 e 1; ponte: questões 13–14 (a dupla integração por partes telescopada). Sem 1n!\frac{1}{n!}, os dados nas extremidades continuam inteiros (com ainda mais facilidade), mas a estimativa se torna πan4n\frac{\pi a^n}{4^n}, que tende a 00 apenas quando a<4a < 4 — e todo candidato tem a=bπ2>9a = b\,\pi^2 > 9. O fatorial é exatamente o que ultrapassa o crescimento geométrico ana^n: sem fatorial, sem teorema.

17. Para a10a \leq 10, o inteiro G(0)+G(1)G(0) + G(1) é positivo e no máximo π(10/4)n/n!\pi\,(10/4)^n/n!. Em n=7n = 7: 2.57=610.352.5^7 = 610.35\dots, de modo que a estimativa é π×610.3550400.38<1\frac{\pi \times 610.35}{5040} \approx 0.38 < 1 (em n=6n = 6 ela ainda vale 1,071{,}07): a contradição aparece no sétimo estágio, explicitamente.

18. Duas integrações por partes:

01x(1x)sinπx ⁣dx=1π01(12x)cosπx ⁣dx=2π201sinπx ⁣dx=2π22π=4π3\int_0^1 x(1-x)\sin\pi x\,\dd x = \frac1\pi\int_0^1 (1 - 2x)\cos\pi x\,\dd x = \frac{2}{\pi^2}\int_0^1 \sin\pi x\,\dd x = \frac{2}{\pi^2}\cdot\frac{2}{\pi} = \frac{4}{\pi^3}

(os termos de fronteira se anulam: x(1x)x(1-x) em 0,10, 1, e sinπx\sin\pi x em 0,10, 1). Simbolicamente, o telescópio n=1n = 1 (sem suposição sobre π\pi) diz π301f1sinπx=(f1(0)+f1(1))\pi^3\int_0^1 f_1\sin\pi x = -(f_1''(0) + f_1''(1)) com f1=x(1x)f_1 = x(1 - x), f1=2f_1'' = -2: lado direito 44 — os dois cálculos concordam.

19. ex\eu^x resolve y=yy' = y e sinπx\sin\pi x resolve y=π2yy'' = -\pi^2 y: equações lineares com coeficientes constantes, de modo que a integração por partes faz o núcleo circular de volta a si mesmo e mantém todos os dados de fronteira dentro de Z+Ze\Z + \Z\eu (resp. polinômios inteiros em π2\pi^2). Essa propriedade de fechamento é o que a máquina precisa. Refinado com núcleos adaptados a vários pontos ao mesmo tempo, o mesmo mecanismo produz o teorema de Hermite (e\eu transcendente, 1873) e o de Lindemann (π\pi transcendente, 1882) — além deste volume.

20. Divisão polinomial (ou multiplique de volta e confira):

x4(1x)4=(x64x5+5x44x2+4)(1+x2)4.x^4(1-x)^4 = (x^6 - 4x^5 + 5x^4 - 4x^2 + 4)(1 + x^2) - 4 .

Integrando a identidade exibida, dividida por 1+x21 + x^2:

01x4(1x)41+x2 ⁣dx=(1746+143+4)4arctan1=227π,\int_0^1 \frac{x^4(1-x)^4}{1+x^2}\dd x = \Bigl(\frac17 - \frac46 + 1 - \frac43 + 4\Bigr) - 4\arctan 1 = \frac{22}{7} - \pi ,

usando 1723+143+4=3+17\frac17 - \frac23 + 1 - \frac43 + 4 = 3 + \frac17 e arctan1=π4\arctan 1 = \frac\pi4.

21. O integrando é contínuo e positivo em (0,1)\intoo{0}{1}: a integral é >0> 0, de modo que π<227\pi < \frac{22}{7}. Além disso, 1211+x21\frac12 \leq \frac{1}{1+x^2} \leq 1 em [0,1]\intcc{0}{1} e 01(x(1x))4=(4!)29!=1630\int_0^1 (x(1-x))^4 = \frac{(4!)^2}{9!} = \frac{1}{630} (questão 3):

11260227π16303.14126<2271630π22711260<3.14207.\frac{1}{1260} \leq \frac{22}{7} - \pi \leq \frac{1}{630} \quad\Longrightarrow\quad 3.14126 < \frac{22}{7} - \frac{1}{630} \leq \pi \leq \frac{22}{7} - \frac{1}{1260} < 3.14207 .

22. Módulo x2+1x^2 + 1: x21x^2 \equiv -1, de modo que x4m=(x2)2m1x^{4m} = (x^2)^{2m} \equiv 1 e (1x)2=12x+x22x(1 - x)^2 = 1 - 2x + x^2 \equiv -2x, donde (1x)4m(2x)2m=4m(x2)m(4)m(1-x)^{4m} \equiv (-2x)^{2m} = 4^m (x^2)^m \equiv (-4)^m: o resto é a constante (4)m(-4)^m, e o quociente QmQ_m tem coeficientes inteiros (divisão por um polinômio inteiro mônico). Dividindo a identidade por 1+x21 + x^2 e integrando:

Jm:=01(x(1x))4m1+x2 ⁣dx=sm+(4)mπ4,sm=01QmQ.J_m := \int_0^1 \frac{(x(1-x))^{4m}}{1+x^2}\dd x = s_m + (-4)^m\,\frac{\pi}{4}, \qquad s_m = \int_0^1 Q_m \in \Q .

Resolvendo para π\pi: com rm=(1)m+141msmQr_m = (-1)^{m+1}\,4^{\,1-m} s_m \in \Q,   πrm=41mJm41m44m=415m\;\abs{\pi - r_m} = 4^{\,1-m} J_m \leq 4^{\,1-m}\cdot4^{-4m} = 4^{\,1-5m}. Para m=1m = 1: s1=227s_1 = \frac{22}{7}, r1=227r_1 = \frac{22}{7}, estimativa 44=12564^{-4} = \frac{1}{256} — de novo as questões 20–21.

23. π227=227π1.26103\bigl|\pi - \frac{22}{7}\bigr| = \frac{22}{7} - \pi \approx 1.26\cdot10^{-3}, dezesseis vezes melhor que o padrão de ordem 22, 1722.0102\frac{1}{7^2} \approx 2.0\cdot10^{-2}, garantido por Dirichlet (Problema 14.1, questão 4); e π3551132.7107\bigl|\pi - \frac{355}{113}\bigr| \approx 2.7\cdot10^{-7} supera 111327.8105\frac{1}{113^2} \approx 7.8\cdot10^{-5} por um fator 300\approx 300. Nenhuma contradição com o que foi demonstrado: as desigualdades de Liouville limitam inferiormente os erros de aproximação apenas para números algébricos, e nenhuma estimativa dessas para π\pi está disponível neste nível — π\pi está livre para ser aproximado espetacularmente bem.

24. Lema: suponha x=pqx = \frac pq e 0<an+bnx00 < \abs{a_n + b_n x} \to 0. Então an+bnx=qan+pbnq\abs{a_n + b_n x} = \frac{\abs{q a_n + p b_n}}{q}, com qan+pbnq a_n + p b_n um inteiro não nulo (não nulo, porque o valor absoluto é >0> 0): assim, an+bnx1q\abs{a_n + b_n x} \geq \frac1q para todo nn, contradizendo a convergência a 00. Instâncias: questão 7 (x=ex = \eu, an=αna_n = \alpha_n, bn=βnb_n = \beta_n); Exercício 11.9 (x=ex = \eu de novo, com aq=q!kq1k!a_q = -q!\sum_{k \leq q}\frac{1}{k!}, bq=q!b_q = q!); e o Problema 14.1, questão 1, é a sua forma geométrica. Na Parte III e no Exercício 15.9, a armadilha roda dentro da contradição: supor a racionalidade converte uma expressão num inteiro, que a análise em seguida comprime em (0,1)\intoo{0}{1} — o mesmo princípio, transposto.

25. (i) A integralidade vive no cálculo nas extremidades de polinômios (questões 6, 11–12), a pequenez nas estimativas que a monotonicidade do sup\sup dá (questões 4 e 10), e a ponte na integração por partes (questões 8, 13–14) — os três são teoremas deste capítulo. (ii) A integral fornece o que o teorema do valor médio não podia: uma identidade exata entre o objeto analítico e os dados aritméticos (igualdade, e não apenas uma desigualdade com um cc desconhecido), e é por isso que a máquina alcança π2\pi^2, ao passo que o Problema 14.1 alcançou apenas expoentes de aproximação. (iii) Extraídos: eQ\eu \notin \Q, π2Q\pi^2 \notin \Q (logo, πQ\pi \notin \Q), e=2.718282\eu = 2.718282 certificado, 2271630π22711260\frac{22}{7} - \frac{1}{630} \leq \pi \leq \frac{22}{7} - \frac{1}{1260}, e (2nn)4n2n+1\binom{2n}{n} \geq \frac{4^n}{2n+1}. (iv) Fronteira: Hermite e Lindemann empurram a mesma máquina até a transcendência; e Apéry (1979) rodou a armadilha do inteiro em ζ(3)\zeta(3) — a máquina ainda está produzindo matemática do século XX.