Mathematics · Book 4 · Bachelor Year 2

Matemáticas universitarias — Grado 2

Matemáticas universitarias — Grado 2 · Bachelor Year 2

9Integración

El volumen del Año 1 construyó la integral en un segmento. este capitulo lo extiende a intervalos arbitrarios (integrales impropias, con el completo kit de herramientas de comparación), luego estudia las integrales dependiendo de un parámetrocontinuidad y la diferenciación bajo la integral signo — impulsado por el teorema de convergencia dominada, el que resultado de este capítulo tomado con confianza. La función Γ\Gamma sirve como ejemplo en ejecución, y como puerta de entrada a la mitad de la especial Funciones de las matemáticas.

9.1 Integrales en un intervalo arbitrario

Definición 9.1

Sea ff continuo por partes en [a,b)\intco{a}{b} (bRb \in \R o ++\infty). La integral converge cuando existe limxbaxf\lim_{x \to b^-} \int_a^x f; uno entonces escribe abf\int_a^b f para el límite. (De manera similar en (a,b]\intoc{a}{b}, y en (a,b)\intoo{a}{b} dividiendo en un punto interior — la elección no no importa, por Chasles.) La integral converges absolutamente cuando abf\int_a^b \abs f converge; la convergencia absoluta implica convergencia, por el criterio de Cauchy:

xyfxyf\Bigl| \int_x^{y} f \Bigr| \leq \int_x^{y} \abs f

y lo completo de R\R (la primitiva tiene la propiedad Cauchy). En detalle: sean F(x)=axfF(x) = \int_a^x f y G(x)=axfG(x) = \int_a^x \abs f. Si bf\int^b\abs fconverge, GGtiene un límite en bb^-, por lo que por cada ε>0\varepsilon > 0 hay c<bc < b con G(y)G(x)εG(y) - G(x) \leq \varepsilonsiempre que cxy<bc \leq x \leq y < b; la pantalla transfiere esta propiedad de Cauchy a FF. Para cualquier secuencia xnbx_n \to b^-, los valores F(xn)F(x_n) forman una secuencia de Cauchy de reales, convergentes por lo completo, y entrelazando dos de tales secuencias muestra que el límite es el mismo para todas ellas: FF tiene un límite en bb^-.

Teorema 9.2 (Kit de herramientas de comparación positiva)

Para f,g0f, g \geq 0 continuo por partes en [a,b)\intco{a}{b}:

  1. abf\int_a^b f converge si la primitiva xaxfx \mapsto \int_a^x f está acotada;
  2. fgf \leq g: la convergencia de g\int g fuerza a la de f\int f; la divergencia se transfiere en sentido contrario;
  3. fgf \sim g en bb: las dos integrales tienen la misma naturaleza;
  4. las escalas de referencia: en ++\infty,  ⁣dttα\int^{\infty} \frac{\dd t}{t^\alpha}converge iff α>1\alpha > 1, y  ⁣dtt(lnt)β\int^\infty \frac{\dd t}{t(\ln t)^\beta} y si β>1\beta > 1; en un punto final finito bb, b ⁣dt(bt)α\int^b \frac{\dd t}{(b - t)^\alpha}converge si α<1\alpha < 1.

Demostración. (1) La primitiva F(x)=axfF(x) = \int_a^x f no es decreciente (f0f \geq 0). Si es acotado, =supx<bF\ell = \sup_{x < b}F es finito y F(x)F(x) \to \ell: dado ε>0\varepsilon > 0, algunos F(x0)>εF(x_0) > \ell - \varepsilony trampas de monotonicidad F(x)(ε,]F(x) \in \intoc{\ell - \varepsilon}{\ell} para x0x<bx_0 \leq x < b. si es ilimitado, F+F \to +\infty: divergencia.

(2) De fgf \leq g: axfaxg\int_a^x f \leq \int_a^x g para todos los xx; si bg\int^b g converge, el lado derecho está acotado, por lo tanto también lo está la izquierda, y (1) concluye. La contraposición transfiere divergencia. al revés.

(3) fgf \sim g en bb proporciona a c<bc < b

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

por (2) aplicado en ambos sentidos en [c,b)\intco{c}{b}, las dos integrales tener la misma naturaleza; la pieza inicial [a,c]\intcc{a}{c} es una integral adecuada y no cambia nada.

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

cx ⁣dttα=x1αc1α1α,cx ⁣dtt(lnt)β=(lnx)1β(lnc)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},

con logaritmos en los casos excluidos: acotado como x+x \to +\inftyexactamente cuando α>1\alpha > 1, resp. β>1\beta > 1. en un punto final finito, la sustitución u=btu = b - t se reduce al escala 0uα ⁣du\int_0 u^{-\alpha}\,\dd u, acotada si α<1\alpha < 1. Aplicar (1) cada vez.

Ejemplo 9.3 (Dos calentamientos, trabajados hasta el final)

(a) 01lnt ⁣dt\displaystyle\int_0^1 \ln t\,\dd t: el integrando explota en 0+0^+, pero lnt=o(t1/2)\abs{\ln t} = o\bigl(t^{-1/2}\bigr) allí (los logaritmos pierden potencias), y 0t1/2\int_0 t^{-1/2} converge: convergencia absoluta. El valor, por partes [ε,1]\intcc{\varepsilon}{1}:

ε1lnt ⁣dt=[tlntt]ε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) 0lnt1+t2 ⁣dt\displaystyle\int_0^\infty \frac{\ln t}{1 + t^2}\,\dd t: problema en ambos extremos, por lo que se divide en 11. Cerca de 00: lnt\abs{\ln t} integrable como en (a); cerca de \infty: lnt1+t2=o(t3/2)\frac{\ln t}{1+t^2} = o(t^{-3/2}): absolutamente convergente. el la sustitución t=1ut = \frac1u asigna (0,1)\intoo{0}{1} a (1,)\intoo{1}{\infty} y

01lnt1+t2 ⁣dt=1lnu1+u2 ⁣duu2=1lnu1+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 :

las dos mitades se cancelan y la integral es 00. Cierre Insight: la simetría bajo t1tt \mapsto \frac1t vale una página de cálculo — el mismo truco ya está activado Ejercicio 9.3.

Ejemplo 9.4 (Un valor, tres integrales)

Estudio I=01costt2 ⁣dtI = \displaystyle\int_0^{\infty} \frac{1 - \cos t}{t^2}\,\dd t. En 00: 1costt221 - \cos t \sim \frac{t^2}2, entonces el integrando extiende continuamente por el valor 12\frac12 — no singularidad en absoluto. En \infty: 01costt22t20 \leq \frac{1 - \cos t}{t^2} \leq \frac{2}{t^2}: convergencia absoluta (Teorema 9.2). Valor: integrar por partes en [ε,M]\intcc{\varepsilon}{M} con u=1costu = 1 - \cos t, v=t2v' = t^{-2}:

εM1costt2 ⁣dt=[1costt]εM+εMsintt ⁣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 .

El soporte desaparece en ambos extremos (1cosεεε2\frac{1 - \cos\varepsilon}{\varepsilon} \sim \frac\varepsilon2; acotado numerador en MM), y la integral tiende al valor de Dirichlet π2\frac\pi2 (Ejercicio 9.10): I=π2I = \frac\pi2. Información final: con 1cost=2sin2t21 - \cos t = 2\sin^2\frac t2 y u=t2u = \frac t2,

I=02sin2u(2u)2  2 ⁣du=0(sinuu) ⁣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 :

los tres clásicos 0sintt ⁣dt\int_0^\infty\frac{\sin t}{t}\dd t, 0(sintt)2 ⁣dt\int_0^\infty\bigl(\frac{\sin t}{t}\bigr)^2\dd t (Ejercicio 9.11) y II comparten el valor π2\frac\pi2, circulado por piezas y sustitución — y Sólo el primero es semiconvergente: la integración por partes comercializadas. elimina la convergencia absoluta para un integrando más simple.

Ejemplo 9.5 (Una integral semiconvergente)

1sintt ⁣dt\displaystyle\int_1^{\infty} \frac{\sin t}{t}\,\dd t converge: integrar por partes,

1xsintt ⁣dt=[costt]1x1xcostt2 ⁣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 ,

donde el corchete tiene un límite y la última integral converge absolutamente (cost/t2t2\abs{\cos t}/t^2 \leq t^{-2}). Pero no absolutamente: de sintsin2t\abs{\sin t} \geq \sin^2 t,

1xsintt ⁣dt    1xsin2tt ⁣dt=1x ⁣dt2t= 12lnx      1xcos2t2t ⁣dtconvergent,\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{convergent}} ,

donde la última integral converge por la mismo integración por partes como arriba (con sin2t\sin 2t entre paréntesis): una pieza divergente menos una convergente diverge. entonces 1sintt ⁣dt\int_1^\infty\frac{\sin t}{t}\dd t converge sin convergente absolutamente — el análogo integral del series alternas, con la integración por partes desempeñando el papel de la prueba alterna.

9.2 El teorema de la convergencia

Teorema 9.6 (Convergencia dominada)

Sea (fn)(f_n) continuo por partes en un intervalo II, convergiendo puntualmente a un continuo ff por partes, y supongamos que hay un fijo integrable φ0\varphi \geq 0 (Iφ<\int_I \varphi < \infty) con

fn(t)φ(t)(tI, nN).\abs{f_n(t)} \leq \varphi(t) \qquad (t \in I,\ n \in \N).

Entonces todos Ifn\int_I f_n y If\int_I f convergen absolutamente, y

IfnnIf.\int_I f_n \xrightarrow[n \to \infty]{} \int_I f .

Demostración. Admitido a este nivel.

Observación 9.7

La prueba honesta pertenece a la teoría de la integración de Lebesgue de Año 3; la declaración, sin embargo, se utiliza constantemente a partir de ahora. el La hipótesis dominación es el punto: puntual la convergencia por sí sola no es suficiente (fn=n1(0,1/n)f_n = n\,\mathbf{1}_{\intoo{0}{1/n}}, topes deslizantes: fn=1↛0=f\int f_n = 1 \not\to 0 = \int f). el teorema también es válido para un parámetro continuo (fλf_\lambda, λλ0\lambda \to \lambda_0), por la caracterización secuencial de límites.

Ejemplo 9.8 (Un límite gaussiano, por dominación)

Calcule limnIn\displaystyle\lim_{n\to\infty} I_n donde In=0(1+t2n) ⁣n ⁣dtI_n = \int_0^\infty \Bigl(1 + \frac{t^2}{n}\Bigr)^{\!-n}\dd t. Puntualmente, (1+t2/n)net2(1 + t^2/n)^n \to \eu^{t^2} (interés compuesto límite), por lo que los integrandos tienden a et2\eu^{-t^2}. Dominación: la La secuencia n(1+u/n)nn \mapsto (1 + u/n)^n no es decreciente para u0u \geq 0. (AM–GM en los factores n+1n + 1 1,1+un,,1+un1, 1 + \frac un, \dots, 1 + \frac unda (1+un+1)n+1(1+un)n(1 + \frac u{n+1})^{n+1} \geq (1 + \frac un)^n), entonces para n2n \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},

un dominador integrable (4t4\sim 4t^{-4} en el infinito). dominado convergencia:

Inn0et2 ⁣dt=π2I_n \xrightarrow[n\to\infty]{} \int_0^\infty \eu^{-t^2}\dd t = \frac{\sqrt\pi}{2}

(la integral gaussiana de Ejercicio 9.8). Cierre comprobar: la sustitución t=ntanθt = \sqrt n\tan\theta calcula InI_n exactamente, In=n0π/2cos2n2θ ⁣dθ=nW2n2I_n = \sqrt n\int_0^{\pi/2}\cos^{2n-2}\theta\,\dd\theta = \sqrt n\,W_{2n-2}, y las asintóticas de Wallis Wmπ/(2m)W_m \sim \sqrt{\pi/(2m)}(Lema 6.11) dan nuevamente nW2n2π2\sqrt n\,W_{2n-2} \to \frac{\sqrt\pi}2: los dos pilares de este El capítulo y el último están de acuerdo.

Ejemplo 9.9 (Convergencia dominada, parámetro continuo)

calcular

limx+0arctan(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 como xx \to \infty; y el dominación

arctan(xt)1+t2π/21+t2,integrable, independent of x,\Bigl|\frac{\arctan(xt)}{1+t^2}\Bigr| \leq \frac{\pi/2}{1+t^2}, \qquad\text{integrable, independent of } x,

es válido para todos los xx. Por la forma de parámetro continuo de Teorema 9.6 (secuencial caracterización: prueba a lo largo de cada xnx_n \to \infty),

0arctan(xt)1+t2 ⁣dtx+π20 ⁣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} .

Información final: el punto único t=0t = 0, donde el punto El límite es 00 en lugar de π2\frac\pi2, no cambia nada — el La función límite sólo entra por su integral, una de las misericordias silenciosas del teorema.

9.3 Integrales con un parámetro

Teorema 9.10 (Continuidad bajo el signo integral)

Sea f ⁣:A×IRf \colon A \times I \to \R (AA a espacio métrico, II un intervalo) con: tf(x,t)t \mapsto f(x, t) por partes continuo para cada xx; xf(x,t)x \mapsto f(x, t) continuo para cada tt; y un dominación f(x,t)φ(t)\abs{f(x,t)} \leq \varphi(t) (φ\varphi integrable en II, independiente de xx). entonces

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

está definido y continuo en AA.

Demostración. Definición: la dominación da convergencia absoluta. Continuidad en x0x_0: para cualquier secuencia xnx0x_n \to x_0, las funciones gn(t)=f(xn,t)g_n(t) = f(x_n, t)convergen puntualmente a f(x0,t)f(x_0, t)(continuidad en xx) bajo la dominación fija φ\varphi: la convergencia dominada da F(xn)F(x0)F(x_n) \to F(x_0); concluir por la caracterización secuencial de continuidad (Definición 4.5).

Teorema 9.11 (Diferenciación bajo el signo integral)

Deje que f ⁣:J×IRf \colon J \times I \to \R (JJ sea un intervalo de parámetros) con: tf(x,t)t \mapsto f(x,t) integrable en II para cada xx; xf(x,t)x \mapsto f(x,t)de la clase C1C^1 para cada tt, la derivada parcial fx\frac{\partial f}{\partial x} siendo continuo por partes en tt y dominado: fx(x,t)ψ(t)\bigl|\frac{\partial f}{\partial x}(x,t)\bigr| \leq \psi(t)con ψ\psiintegrable. Entonces F(x)=If(x,t) ⁣dtF(x) = \int_I f(x,t)\dd t es C1C^1 en JJ y

F(x)=Ifx(x,t) ⁣dt.F'(x) = \int_I \frac{\partial f}{\partial x}(x, t)\,\dd t .

Demostración. Repare xx y hn0h_n \to 0. Los cocientes de diferencias

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

tienen integrandos que convergen puntualmente a fx(x,t)\frac{\partial f}{\partial x}(x, t)y dominados por ψ(t)\psi(t): por el valor medio desigualdad aplicada en xx en tt fijo,

f(x+hn,t)f(x,t)hnsupξfx(ξ,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) .

La convergencia dominada da el límite Ifx(x,t) ⁣dt\int_I \frac{\partial f}{\partial x}(x,t)\dd tde los cocientes: FF es diferenciable con la derivada anunciada, que es continuo por Teorema 9.10 se aplicó a fx\frac{\partial f}{\partial x}.

Ejemplo 9.12 (Una integral de parámetro comparada con una fórmula)

Deja F(x)=0 ⁣dtt2+xF(x) = \displaystyle\int_0^\infty \frac{\dd t}{t^2 + x} para x>0x > 0. En cada [a,b](0,)\intcc{a}{b} \subset \intoo{0}{\infty}, el integrando está dominado por 1t2+a\frac{1}{t^2 + a}, integrable e independiente de xx: FF es continuo (Teorema 9.10). Aquí el teorema puede ser comparado con un valor explícito:

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

visiblemente continuo. Ahora diferencia bajo la integral: la xx-derivado 1(t2+x)2-\frac{1}{(t^2+x)^2} está dominado [a,b]\intcc ab por 1(t2+a)2\frac{1}{(t^2+a)^2}, integrables: Teorema 9.11 da

F(x)=0 ⁣dt(t2+x)2whileF(x)=π4x3/2,F'(x) = -\int_0^\infty \frac{\dd t}{(t^2 + x)^2} \qquad\text{while}\qquad F'(x) = -\frac{\pi}{4}\,x^{-3/2} ,

entonces tenemos calculado una nueva integral gratis: 0 ⁣dt(t2+x)2=π4x3/2\int_0^\infty\frac{\dd t}{(t^2+x)^2} = \frac{\pi}{4x^{3/2}}. Información final: diferenciar un integral de parámetros conocido es una fábrica de nuevas fórmulas — iterar da 0 ⁣dt(t2+1)n\int_0^\infty\frac{\dd t}{(t^2+1)^n} por cada nn, sin sustituciones trigonométricas.

Método 9.13 (Estudiando una integral impropia)

Dado abf\int_a^b f:

  1. Localice el problema: enumere los puntos finales (o el interior puntos) donde ff no está acotado o el intervalo es infinita y dividida para que cada pieza tenga exactamente una final problemático.
  2. Si ff tiene un signo constante cerca de ese extremo, encuentre un equivalente y comparar con las escalas de referencia de Teorema 9.2.
  3. Si ff oscila, pruebe primero f\abs f (valor absoluto). convergencia). Si f\int\abs f diverge, intente la integración por partes para cambiar la oscilación por decaimiento, como en Ejemplo 9.5; minorías como sintsin2t\abs{\sin t} \geq \sin^2t detectar genuino semiconvergencia.
  4. Para un valor, no sólo la naturaleza: partes, sustitución, o un parámetro (derivar una integral más simple, como en Ejemplo 9.12 y Ejemplo 9.21).
  5. Comprobaciones de cordura sobre cualquier valor calculado: signo y aproximado tamaño frente a un límite crudo (0et2 ⁣dt(0,1+1et)\int_0^\infty \eu^{-t^2}\dd t \in \intoo{0}{1 + \int_1^\infty \eu^{-t}}, por lo que π20.886\frac{\sqrt\pi}{2} \approx 0.886 es plausible); y consistencia dimensional bajo escala (tλtt \mapsto \lambda t debe reescalar ambos lados del de la misma manera — el detector más rápido de un factor perdido).

Observación 9.14 (Errores comunes)

Tres errores recurrentes. (i) Dependiente del parámetro dominadores: la dominación f(x,t)φ(t)\abs{f(x,t)} \leq \varphi(t) debe ser uniforme en xx en el conjunto considerado; normalmente aguanta segmentos [a,b]\intcc ab pero no globalmente — para 0ext ⁣dt\int_0^\infty\eu^{-xt}\dd t no hay dominador integrable válido para todo x>0x > 0, pero dominar en xa>0x \geq a > 0 es suficiente para trabajar en toda la media línea abierto, desde continuidad y los derivados son nociones locales. (ii) Comparación firmada integrandos: el kit de herramientas de comparación es para no negativo funciones; de fg\abs f \leq g con g\int g divergente se puede no concluir nada — 1sintt ⁣dt\int_1^\infty\frac{\sin t}t\,\dd t converge aunque todas las comparaciones con 1t\frac1t fallan. (iii) Olvidar la mitad del problema: en (0,)\intoo{0}{\infty} estudie siempre ambos extremos por separado; 0 ⁣dtt\int_0^\infty\frac{\dd t}{t} diverge en ambos, y un Una división de aspecto convergente puede cancelar silenciosamente dos infinitos. El reflejo seguro es la lista de verificación de Método 9.13.

Ejemplo 9.15 (Un caso límite de Bertrand, al dígito)

La escala  ⁣dtt(lnt)β\int^\infty\frac{\dd t}{t(\ln t)^\beta} de Teorema 9.2 se sienta exactamente en el borde de las escalas de potencia; sus casos límite merecen una revisión completa cálculo. Para β=2\beta = 2:

e ⁣dtt(lnt)2=[1lnt]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 ,

una integral convergente con un valor agradablemente exacto; mientras que para β=1\beta = 1,

ex ⁣dttlnt=[lnlnt]ex=lnlnx,\int_\eu^{x}\frac{\dd t}{t\ln t} = \bigl[\ln\ln t\bigr]_\eu^{x} = \ln\ln x \longrightarrow \infty ,

divergente — pero tan lentamente que llegando a lnlnx=10\ln\ln x = 10 requiere x=ee10109566x = \eu^{\eu^{10}} \approx 10^{9566}. Cierre insight: entre “todos los poderes t1εt^{-1-\varepsilon} converge” y “t1t^{-1} divergen” viven una escalera infinita de escalas logarítmicas, cada una refinando a la anterior; la sustitución u=lntu = \ln t colapsa cada renglón sobre el anterior, lo que Por eso los criterios de Bertrand se hacen eco de los de Riemann en un nivel arriba.

Observación 9.16 (Perspectivas dentro de este volumen)

Las herramientas de este capítulo están a punto de estar en todas partes. dominado La convergencia es el motor detrás de las identidades aproximadas de el siguiente capítulo (núcleos deslizantes, Bernstein y Fejér por igual); continuidad y diferenciación bajo el signo integral producir el cálculo de los coeficientes de Fourier en el sistema de Fourier capítulo, donde cada cn(f)c_n(f) es un integral de parámetros en disfrazarse. La función Γ\Gamma regresa dos veces: en el capítulo en integrales múltiples, donde una integral doble finalmente demuestra Fórmula Beta-Gamma de Euler en su totalidad y en la probabilidad capítulos, donde las integrales tipo Γ\Gamma normalizan el estándar densidades y calcular sus momentos. Y el semiconvergente sintt\int\frac{\sin t}{t} resurge como la constante de Gibbs del Capítulo de Fourier — la misma integral, midiendo el exceso de sumas parciales en un salto.

Definición 9.17 (La función Γ\Gamma)

Para x>0x > 0:

Γ(x)=0tx1et ⁣dt,\Gamma(x) = \int_0^{\infty} t^{x-1}\,\eu^{-t}\,\dd t ,

convergente en ambos extremos (tx1t^{x-1} integrable en 0+0^+ para x>0x > 0; decaimiento exponencial en \infty).

Teorema 9.18

Γ\Gamma es continuo en (0,+)\intoo{0}{+\infty}, satisface las ecuación funcional

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

y es de clase C1C^1 (de hecho CC^\infty) con Γ(x)=0tx1etlnt ⁣dt\Gamma'(x) = \int_0^\infty t^{x-1}\eu^{-t}\ln t\,\dd t.

Demostración. Ecuación funcional: integrar por partes en [ε,M]\intcc{\varepsilon}{M} y dejar ir los extremos: txet=[txet]+xtx1et\int t^{x}\eu^{-t} = [-t^x\eu^{-t}] + x\int t^{x-1}\eu^{-t}, los términos límite desaparecen — de hecho εxeε0\varepsilon^x\eu^{-\varepsilon} \to 0 como ε0+\varepsilon \to 0^+porque x>0x > 0y MxeM0M^x\eu^{-M} \to 0 como MM \to \infty porque la exponencial le gana a toda potencia; ambos truncados las integrales convergen a sus valores impropios por la convergencia establecido en Definición 9.17. Γ(1)=et=1\Gamma(1) = \int \eu^{-t} = 1; la inducción da el factorial.

Continuidad en [a,b](0,)\intcc{a}{b} \subset \intoo{0}{\infty}: dominar tx1ett^{x-1}\eu^{-t} por φ(t)=(ta1+tb1)et\varphi(t) = (t^{a-1} + t^{b-1})\eu^{-t}, integrable e independiente de x[a,b]x \in \intcc{a}{b}: Teorema 9.10 se aplica en cada uno de esos segmentos, por tanto, en toda la media línea. Diferenciabilidad: la derivada xx tx1etlntt^{x-1}\eu^{-t}\ln t está dominado sobre [a,b]\intcc{a}{b} por (ta1+tb1)etlnt(t^{a-1} + t^{b-1})\eu^{-t}\,\abs{\ln t}, aún integrable: Teorema 9.11; iterar lo da todo derivados (cada uno añade un poder de lnt\ln t, inofensivos).

Ejemplo 9.19 (Factoriales medio enteros)

La ecuación funcional y Γ(12)=π\Gamma\bigl(\frac12\bigr) = \sqrt\pi (una sustitución lejos de Ejercicio 9.8: establecer t=u2t = u^2 en la definición integral) genera todos los valores semienteros:

Γ(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} .

Desde Γ(n+1)=n!\Gamma(n+1) = n!, es justo decir “12!=π20.886\frac12! = \frac{\sqrt\pi}{2} \approx 0.886”: el factorial ha sido interpolado, y la curva de interpolación cae por debajo de 11 entre 0!=10! = 1 y 1!=11! = 1 (su mínimo 0.8856\approx 0.8856 en x1.4616x \approx 1.4616 coincide con la imagen de convexidad del Problema del fin de semana, Parte I). Perspectiva final: nada en el privilegios integrales 0tx1et ⁣dt\int_0^\infty t^{x-1}\eu^{-t}\dd t números enteros — la discreción del factorial fue un accidente de contando, y π\sqrt\pi es lo que vive entre 11 y 11.

Observación 9.20 (Adónde va Γ\Gamma desde aquí)

El problema del fin de semana de este capítulo construye todo el conjunto de Euler. cálculo alrededor de Γ\Gamma: el función beta, su integración por partes recursiones, las integrales de Wallis como valores Beta y las de Gauss. fórmula límite. El capítulo sobre integrales múltiples demuestra la afirmación de Euler. Beta: fórmula gamma para todos los argumentos mediante una integral doble; el Los capítulos de probabilidad se encuentran nuevamente con Γ\Gamma en la normalización de las densidades más habituales y en los momentos de tiempos de espera. El volumen del Año 3 reconstruye Γ\Gamma sobre los cimientos de Lebesgue, demuestra el teorema de unicidad de Bohr-Mollerup y extiende La fórmula de Stirling de números enteros a la media línea real por convergencia dominada.

Ejemplo 9.21 (Un cálculo clásico por diferenciación)

Para xRx \in \R, sea

F(x)=0et2cos(xt) ⁣dtF(x) = \int_0^{\infty} \eu^{-t^2}\cos(xt)\,\dd t

(absolutamente convergente, dominada por et2\eu^{-t^2}). Por Teorema 9.11 (dominación de la xx-derivada por tet2t\,\eu^{-t^2}, integrable):

F(x)=0tet2sin(xt) ⁣dt=[et22sin(xt)]0x20et2cos(xt) ⁣dt=x2F(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 con u=tet2u' = t\eu^{-t^2}). La ecuación diferencial F=x2FF' = -\frac x2 Fse integra a F(x)=F(0)ex2/4F(x) = F(0)\,\eu^{-x^2/4}: la La integral de tipo gaussiano se reproduce a sí misma. La constante F(0)=0et2 ⁣dt=π2F(0) = \int_0^\infty \eu^{-t^2}\dd t = \frac{\sqrt\pi}{2} se calcula en Ejercicio 9.8 — y nuevamente, por doble integración, en Capítulo 20.

9.4 Ceremonias

Ejercicio 9.1

Naturaleza de: 01 ⁣dtt(1t)\displaystyle\int_0^1 \frac{\dd t}{\sqrt{t(1-t)}};   1lntt2 ⁣dt\;\displaystyle\int_1^\infty \frac{\ln t}{t^2}\dd t;   0 ⁣dt1+t2sin2t\;\displaystyle\int_0^\infty \frac{\dd t}{1 + t^2\sin^2 t} (compare with the divergent harmonic-type behavior near t=nπt = n\pi).

Solución

Solución de Ejercicio 9.1.

01 ⁣dtt(1t)\int_0^1 \frac{\dd t}{\sqrt{t(1-t)}}: cerca de 00, t1/2\sim t^{-1/2} (α=12<1\alpha = \frac12 < 1: converge); cerca de 11, (1t)1/2\sim (1-t)^{-1/2}: converge. Convergente (su valor es π\pi, por el sustitución t=sin2θt = \sin^2\theta).

1lntt2\int_1^\infty \frac{\ln t}{t^2}: lntt2=o(t3/2)\frac{\ln t}{t^2} = o(t^{-3/2}): convergente (valor 11 por partes).

0 ⁣dt1+t2sin2t\int_0^\infty \frac{\dd t}{1 + t^2\sin^2 t}: divergente. Cerca de t=nπt = n\pi, escriba t=nπ+ut = n\pi + u: sin2t=sin2uu2\sin^2 t = \sin^2 u \leq u^2, así sucesivamente u1n\abs u \leq \frac{1}{n}, 1+t2sin2t1+(nπ+1)2u2Cn2u2+11 + t^2\sin^2 t \leq 1 + (n\pi + 1)^2u^2 \leq C n^2 u^2 + 1; por lo tanto

nπ1/nnπ+1/n ⁣dt1+t2sin2t1/n1/n ⁣du1+Cn2u2=2arctanCC1n,\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} ,

un término de una serie de tipo armónico divergente: sumando nn, el primitivo es ilimitado.

Ejercicio 9.2

Calcular 0tneλt ⁣dt\displaystyle\int_0^\infty t^n \eu^{-\lambda t}\,\dd t (λ>0\lambda > 0) vía Γ\Gamma, y 01(lnt)n ⁣dt\displaystyle\int_0^1 (\ln t)^n \dd t mediante la sustitución t=eut = \eu^{-u}.

Solución

Solución de Ejercicio 9.2.

Sustituir u=λtu = \lambda t:

0tneλt ⁣dt=1λn+10uneu ⁣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}} .

Con t=eut = \eu^{-u} ( ⁣dt=eu ⁣du\dd t = -\eu^{-u}\dd u):

01(lnt)n ⁣dt=0(u)neu ⁣du=(1)nn!.\int_0^1 (\ln t)^n \dd t = \int_0^{\infty} (-u)^n \eu^{-u}\,\dd u = (-1)^n\, n! .

Ejercicio 9.3

Demuestre que 0 ⁣dt(1+t2)(1+tx)\displaystyle\int_0^{\infty} \frac{\dd t}{(1 + t^2)(1 + t^x)} está bien definido para cada xRx \in \R y independiente de xx. (Substitute t1tt \mapsto \frac1t and average the two expressions.) ¿Cuál es su valor?

Solución

Solución de Ejercicio 9.3.

Convergencia: el integrando es 11+t2\leq \frac{1}{1+t^2} cerca de \infty y acotado cerca de 00 (ambos factores acotados a continuación lejos de 00): absolutamente convergente, por cada xx. Sustituyendo t=1ut = \frac1u ( ⁣dt= ⁣duu2\dd t = -\frac{\dd u}{u^2}):

I(x)=01(1+1u2)(1+ux) ⁣duu2=0ux(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 .

Sumando las dos expresiones de I(x)I(x):

2I(x)=01+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, independiente de xx.

Ejercicio 9.4 ★★

(Integrales de Bertrand en un punto final finito) ¿Para qué (α,β)(\alpha, \beta)converge 01/2 ⁣dttαlntβ\displaystyle\int_0^{1/2} \frac{\dd t}{t^\alpha\,\abs{\ln t}^\beta}?

Solución

Solución de Ejercicio 9.4.

Cerca de 0+0^+, con u=lntu = \abs{\ln t} \to \infty. Si α<1\alpha < 1: convergencia independientemente de β\beta (compárese con tαt^{-\alpha'} para α<α<1\alpha < \alpha' < 1: se supera el factor logarítmico). Si α>1\alpha > 1: divergencia independientemente de β\beta (compárese con tαt^{-\alpha''}, 1<α<α1 < \alpha'' < \alpha). Si α=1\alpha = 1: sustituir t=eut = \eu^{-u}:

01/2 ⁣dttlntβ=ln2 ⁣duuβ,\int_0^{1/2} \frac{\dd t}{t\,\abs{\ln t}^\beta} = \int_{\ln 2}^{\infty} \frac{\dd u}{u^\beta},

convergente si β>1\beta > 1. Resumen: convergencia si α<1\alpha < 1, o (α=1\alpha = 1 y β>1\beta > 1) — el espejo del Bertrand serie.

Ejercicio 9.5 ★★

Deje F(x)=0ext1+t2 ⁣dtF(x) = \displaystyle\int_0^{\infty} \frac{\eu^{-xt}}{1 + t^2}\,\dd tpor x0x \geq 0. Demuestre que FF es continuo en [0,)\intco{0}{\infty}, C2C^2 en (0,)\intoo{0}{\infty}, satisface F+F=1xF'' + F = \frac1xallí, y ese F(x)0F(x) \to 0 como x+x \to +\infty.

Solución

Solución de Ejercicio 9.5.

Continuidad en [0,)\intco{0}{\infty}: dominación ext1+t211+t2\bigl|\frac{\eu^{-xt}}{1+t^2}\bigr| \leq \frac{1}{1+t^2}, integrable, uniforme en x0x \geq 0: Teorema 9.10.

C2C^2 en (0,)\intoo{0}{\infty}: en xa>0x \geq a > 0, los dos primeros xx-derivados text1+t2\frac{-t\,\eu^{-xt}}{1+t^2} y t2ext1+t2\frac{t^2\eu^{-xt}}{1+t^2} están dominados por teatt\,\eu^{-at} y eat\eu^{-at}: dos aplicaciones de Teorema 9.11. entonces

F(x)+F(x)=0t2+11+t2ext ⁣dt=0ext ⁣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 .

Límite: 0F(x)0ext ⁣dt=1x00 \leq F(x) \leq \int_0^\infty \eu^{-xt}\dd t = \frac1x \to 0.

Ejercicio 9.6 ★★

(Frullani) Sea ff continuo en [0,+)\intco{0}{+\infty} con un límite finito f()f(\infty) en ++\infty. Demuestre que para a,b>0a, b > 0:

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

(On [ε,M]\intcc{\varepsilon}{M}, substitute in each piece and regroup into aεbεaMbM\int_{a\varepsilon}^{b\varepsilon} - \int_{aM}^{bM} of f(u)u ⁣du\frac{f(u)}u\,\dd u; squeeze using the continuidad at 00 and the limit at \infty.) Computar 0ete2tt ⁣dt\int_0^\infty \frac{\eu^{-t} - \eu^{-2t}}{t}\dd t.

Solución

Solución de Ejercicio 9.6.

En [ε,M]\intcc{\varepsilon}{M}, sustituya u=atu = at y u=btu = bt en las dos mitades:

εMf(at)f(bt)t ⁣dt=aεaMf(u)u ⁣dubεbMf(u)u ⁣du=aεbεf(u)u ⁣duaMbMf(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 .

Primera pieza: f(u)=f(0)+o(1)f(u) = f(0) + o(1) cerca de 00, y aεbε ⁣duu=lnba\int_{a\varepsilon} ^{b\varepsilon} \frac{\dd u}{u} = \ln\frac ba: la pieza tiende a f(0)lnbaf(0)\ln\frac ba. Segunda pieza: f(u)f()f(u) \to f(\infty), igual Cálculo: tiende a f()lnbaf(\infty)\ln\frac ba. De ahí el inadecuado integral converge a (f(0)f())lnba\bigl(f(0) - f(\infty)\bigr)\ln\frac ba.

Con f(t)=etf(t) = \eu^{-t} (f(0)=1f(0) = 1, f()=0f(\infty) = 0), a=1a = 1, b=2b = 2:

0ete2tt ⁣dt=ln2.\int_0^\infty \frac{\eu^{-t} - \eu^{-2t}}{t}\dd t = \ln 2 .

Ejercicio 9.7 ★★

Justificar y calcular limn0n(1tn) ⁣ntx1 ⁣dt\lim_{n\to\infty} \displaystyle\int_0^n \Bigl(1 - \frac tn\Bigr)^{\!n} t^{x-1}\,\dd tpara x>0x > 0 (dominated convergence with φ(t)=ettx1\varphi(t) = \eu^{-t}t^{x-1}, using (1t/n)net(1 - t/n)^n \leq \eu^{-t}; the limit is Γ(x)\Gamma(x)).

Solución

Solución de Ejercicio 9.7.

Extienda el integrando en 00 más allá de t=nt = n: gn(t)=(1tn)ntx11tng_n(t) = (1 - \frac tn)^n t^{x-1}\mathbf{1}_{t \leq n}. Puntualmente, gn(t)ettx1g_n(t) \to \eu^{-t}t^{x-1} (el límite de interés compuesto, volumen del año 1). Dominación: ln(1u)u\ln(1 - u) \leq -u da (1tn)net(1 - \frac tn)^n \leq \eu^{-t} sobre [0,n]\intcc{0}{n}, por lo que gn(t)ettx1=φ(t)\abs{g_n(t)} \leq \eu^{-t}t^{x-1} = \varphi(t), integrable. Convergencia dominada:

0n(1tn)ntx1 ⁣dtn0ettx1 ⁣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) .

(Al calcular el lado izquierdo por partes repetidas se obtiene el producto de Euler formulario Γ(x)=limn!nxx(x+1)(x+n)\Gamma(x) = \lim \frac{n!\,n^x}{x(x+1)\cdots(x+n)}.)

Ejercicio 9.8 ★★★

(La integral gaussiana mediante un truco de parámetros) Para x0x \geq 0 establezca

G(x)=(0xet2 ⁣dt) ⁣2,H(x)=01ex2(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 .

Demuestre que G+H=0G' + H' = 0 (diferenciar HH bajo la integral y sustituya u=xtu = xt en la integral resultante), deduzca G(x)+H(x)=π4G(x) + H(x) = \frac\pi4 para todos los xx y concluya

0et2 ⁣dt=π2.\int_0^{\infty} \eu^{-t^2}\,\dd t = \frac{\sqrt\pi}{2} .
Solución

Solución de Ejercicio 9.8.

HH es diferenciable en xx (integrando C1C^1 en xx, derivada 2x(1+t2)ex2(1+t2)1+t2=2xex2ex2t2-2x(1+t^2)\cdot\frac{\eu^{-x^2(1+t^2)}}{1+t^2} = -2x\,\eu^{-x^2}\eu^{-x^2t^2}, continuo y acotado en compactos de xx, dominio sobre t[0,1]t \in \intcc{0}{1} trivial):

H(x)=2xex201ex2t2 ⁣dt=u=xt2ex20xeu2 ⁣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),

desde G(x)=2ex20xet2 ⁣dtG'(x) = 2\eu^{-x^2}\int_0^x \eu^{-t^2}\dd t (regla de cadena en el cuadrado, teorema fundamental del cálculo). Entonces G+HG + H es 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.

Como xx \to \infty: 0H(x)ex201 ⁣dt00 \leq H(x) \leq \eu^{-x^2}\int_0^1 \dd t \to 0, entonces G(x)π4G(x) \to \frac\pi4:

0et2 ⁣dt=π4=π2.\int_0^\infty \eu^{-t^2}\dd t = \sqrt{\frac\pi4} = \frac{\sqrt\pi}{2} .

(En consecuencia Γ(12)=20et2 ⁣dt=π\Gamma\bigl(\frac12\bigr) = 2\int_0^\infty \eu^{-t^2}\dd t = \sqrt\pi, por la sustitución t=ut = \sqrt u.)

Ejercicio 9.9 ★★★

Demuestre que Γ\Gamma es log-convexo: lnΓ\ln\Gamma es convexo en (0,)\intoo{0}{\infty}. (Cauchy–Schwarz for integrals applied to t(x+y)/21et=(tx1et)1/2(ty1et)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} gives Γ(x+y2)2Γ(x)Γ(y)\Gamma\bigl(\frac{x+y}{2}\bigr)^2 \leq \Gamma(x)\Gamma(y); combine with continuidad and Ejercicio 8.8.)

Solución

Solución de Ejercicio 9.9.

Cauchy–Schwarz (volumen del año 1, válido en [ε,M]\intcc{\varepsilon}{M} y pasado al límite) aplicado a la factorización tx+y21et=(tx1et)1/2(ty1et)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/2lnΓ(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 es el punto medio convexo; siendo continuo (Teorema 9.18), es convexo (Ejercicio 8.8). (La convexidad logarítmica fija Γ\Gamma hacia abajo únicamente entre las interpolaciones del factorial — el Teorema de Bohr-Mollerup, una perla del año 3.)

Ejercicio 9.10 ★★★

(Integral de Dirichlet) Configure F(x)=0sinttext ⁣dtF(x) = \displaystyle\int_0^{\infty} \frac{\sin t}{t}\,\eu^{-xt}\,\dd tpara x>0x > 0.

  1. Justifique F(x)=11+x2F'(x) = -\frac{1}{1 + x^2} (diferenciar en la integral; calcular 0extsint ⁣dt\int_0^\infty \eu^{-xt}\sin t\,\dd t mediante dos integraciones por partes).
  2. Demuestre que F(x)0F(x) \to 0 es x+x \to +\infty y deduzca F(x)=π2arctanxF(x) = \frac\pi2 - \arctan x.
  3. Admitiendo el continuidad de FF en 0+0^+ (un tipo Abel teorema), concluir el valor de la semiconvergente integrales:

    0sintt ⁣dt=π2.\int_0^{\infty} \frac{\sin t}{t}\,\dd t = \frac{\pi}{2}.
Solución

Solución de Ejercicio 9.10.

  1. En xa>0x \geq a > 0: la derivada xx del integrando es sintext-\sin t\,\eu^{-xt}, dominado por eat\eu^{-at}: Teorema 9.11 da F(x)=0extsint ⁣dtF'(x) = -\int_0^\infty \eu^{-xt}\sin t\,\dd t. Dos integraciones por partes (o el exponencial complejo):

    0extsint ⁣dt=0e(x+i)t ⁣dt=1xi=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)0ext ⁣dt=1x0\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: 0F(x)=(π2arctanx)0 - F(x) = -\bigl(\frac\pi2 - \arctan x\bigr), entonces F(x)=π2arctanxF(x) = \frac\pi2 - \arctan x.
  3. Dejemos x0+x \to 0^+ con el continuidad admitido: F(0+)=π2F(0^+) = \frac\pi2y F(0)=0sintt ⁣dtF(0) = \int_0^\infty \frac{\sin t}{t}\dd t (la integral de Dirichlet semiconvergente, Ejemplo 9.5): su valor es π2\frac\pi2.

Ejercicio 9.11 ★★

Justifique la convergencia de 0(sintt) ⁣2 ⁣dt\displaystyle\int_0^\infty \Bigl(\frac{\sin t}{t}\Bigr)^{\!2}\dd t, luego calculela por uno integración por partes y Ejercicio 9.10:

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

(El mismo valor que 0sintt ⁣dt\int_0^\infty \frac{\sin t}{t}\dd t — pero esta vez la convergencia es absoluta.)

Solución

Solución de Ejercicio 9.11.

Convergencia: cerca de 00 el integrando extiende a continuamente por el valor 11 (sintt\sin t \sim t); en el infinito es t2\leq t^{-2}: convergencia absoluta. En [ε,M]\intcc{\varepsilon}{M}, integrar por piezas con u=sin2tu = \sin^2 t, v=t2v' = t^{-2}:

εMsin2tt2 ⁣dt=[sin2tt]εM+εM2sintcostt ⁣dt=[sin2tt]εM+2ε2Msinuu ⁣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 en la última integral). El soporte tiende a 00 en ambos termina (sin2ε/εε\sin^2\varepsilon/\varepsilon \leq \varepsilon; sin2M/M1/M\sin^2 M/M \leq 1/M), y la última integral tiende a 0sinuu ⁣du=π2\int_0^\infty \frac{\sin u}{u}\dd u = \frac\pi2 (Ejercicio 9.10). Por lo tanto

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

Ejercicio 9.12 ★★★

(La cola gaussiana) Para x>0x > 0 establezca T(x)=xet2 ⁣dtT(x) = \displaystyle \int_x^\infty \eu^{-t^2}\dd t.

  1. Escribiendo et2=12t(2tet2)\eu^{-t^2} = \frac{1}{-2t}\cdot(-2t\,\eu^{-t^2}), integrar por partes dos veces para obtener

    T(x)=ex2(12x14x3)+34xet2t4 ⁣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. Cota el resto: 034xt4et2 ⁣dt38x5ex20 \leq \frac34\int_x^\infty t^{-4}\eu^{-t^2}\dd t \leq \frac{3}{8x^5}\,\eu^{-x^2}, y deducir el paréntesis

    ex2(12x14x3)T(x)ex22x,henceT(x)ex22x(x+).\eu^{-x^2}\Bigl(\frac{1}{2x} - \frac{1}{4x^3}\Bigr) \leq T(x) \leq \frac{\eu^{-x^2}}{2x}, \qquad\text{hence}\qquad T(x) \sim \frac{\eu^{-x^2}}{2x} \quad (x \to +\infty).
  3. ¿Por qué se puede obtener la serie alterna completa iterando? ¿Las piezas nunca convergen para xx fijo? (Compare the growth of the coefficients 13(2k1)1\cdot3\cdots(2k-1) with the powers (2x2)k(2x^2)^k.)
Solución

Solución de Ejercicio 9.12.

  1. Piezas con u=12tu = \frac{-1}{2t}, v=2tet2v' = -2t\,\eu^{-t^2} (entonces v=et2v = \eu^{-t^2}):

    T(x)=[et22t]xxet22t2 ⁣dt=ex22xxet22t2 ⁣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 .

    Mismo dispositivo en la nueva integral (u=14t3u = \frac{-1}{4t^3}, v=2tet2v' = -2t\,\eu^{-t^2}):

    xet22t2 ⁣dt=ex24x334xet2t4 ⁣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 ,

    De ahí la identidad anunciada.

  2. Una integración más por partes acota el resto:

    xet2t4 ⁣dt=ex22x552xet2t6 ⁣dtex22x5,\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},

    entonces 034xt4et2 ⁣dt38x5ex20 \leq \frac34\int_x^\infty t^{-4}\eu^{-t^2}\dd t \leq \frac{3}{8x^5}\eu^{-x^2}. Dejando caer el (positivo) resto en la identidad de la pregunta 1 da el menor atado; eliminando el segundo término (negativo) del primero piezas da T(x)ex22xT(x) \leq \frac{\eu^{-x^2}}{2x}. Dividiendo el bracketing por ex22x\frac{\eu^{-x^2}}{2x}: la relación es comprimido entre 112x21 - \frac{1}{2x^2} y 11, por lo que T(x)ex22xT(x) \sim \frac{\eu^{-x^2}}{2x}.

  3. La iteración de las piezas produce la serie formal.

    T(x)ex22x(112x2+13(2x2)2135(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),

    cuyo kk-ésimo coeficiente 13(2k1)=(2k)!2kk!1\cdot3\cdots(2k-1) = \frac{(2k)!}{2^k k!} crece más rápido que cualquier geométrico secuencia: para xx fijos los términos 13(2k1)(2x2)k\frac{1\cdot3\cdots(2k-1)}{(2x^2)^k} tiende al infinito (su relación es 2k+12x2\frac{2k+1}{2x^2} \to \infty), por lo que el La serie diverge para cada xx. es un Expansión asintótico: truncada en cualquier punto fijo orden, el error es del orden del primero omitido término como xx \to \infty — pero nunca convergente serie. (Esta estimación de cola es la norma gaussiana límite final de los capítulos de probabilidad.)

9.5 Problema: integrales de Euler — Beta, Gamma y Gauss fórmula límite

Problema 9.1

La función Γ\Gamma de Definición 9.17 es una la mitad del cálculo de integrales de Euler; la otra mitad es la función beta

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

Este problema desarrolla el par (Γ,B)(\Gamma, B) con las herramientas de sólo en este capítulo — integración por partes, sustitución, convergencia dominada — y culmina en Euler Fórmula Beta-Gamma B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \frac{\Gamma(x)\Gamma(y)} {\Gamma(x+y)} en los semienteros y en El límite de Gauss fórmula para Γ\Gamma. En el camino reaparecen las integrales de Wallis de Lema 6.11 como valores Beta, y la fórmula de duplicación de Legendre desaparece.

Parte I — Fine structure of Γ\Gamma.

  1. Recuerde por qué Γ(x)=0tx1et ⁣dt\Gamma(x) = \int_0^\infty t^{x-1}\eu^{-t}\dd tconverge exactamente para x>0x > 0 y muestre

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

    (functional equation plus continuidad of Γ\Gamma at 11).

  2. Demostrar Γ(12)=π\Gamma\bigl(\tfrac12\bigr) = \sqrt\pi (substitute t=u2t = u^2 and invoke Ejercicio 9.8), y deducir Reu2/2 ⁣du=2π\int_\R \eu^{-u^2/2}\dd u = \sqrt{2\pi}.
  3. Mostrar por inducción, para nNn \in \N:

    Γ(n+12)=(2n)!4nn!π.\Gamma\Bigl(n + \frac12\Bigr) = \frac{(2n)!}{4^n\,n!}\,\sqrt\pi .
  4. Justifica Γ(x)=0tx1et(lnt)2 ⁣dt>0\Gamma''(x) = \int_0^\infty t^{x-1}\eu^{-t}(\ln t)^2\dd t > 0 y deduce que Γ\Gamma es estrictamente convexo, alcanza un mínimo único en algunos x0(1,2)x_0 \in \intoo{1}{2} (Γ(1)=Γ(2)=1\Gamma(1) = \Gamma(2) = 1 and Rolle), disminuye en (0,x0)\intoo{0}{x_0} y aumenta en (x0,)\intoo{x_0}{\infty}.
  5. Demuestre que Γ\Gamma supera todos los poderes: por cada kNk \in \N, xk=o(Γ(x))x^k = o\bigl(\Gamma(x)\bigr)como x+x \to +\infty (squeeze xx between integers and use Γ(n+1)=n!\Gamma(n+1) = n! with the monotonicity of question 4).

Parte II — The función beta, by parts.

  1. Muestre que B(x,y)B(x,y) converge exactamente para x>0x > 0 y y>0y > 0, y que B(x,y)=B(y,x)B(x,y) = B(y,x).
  2. Calcule B(x,1)=1xB(x, 1) = \frac1x y pruebe por integración mediante piezas, para x,y>0x, y > 0:

    B(x,y+1)=yxB(x+1,y).B(x, y+1) = \frac{y}{x}\,B(x+1, y) .
  3. De la división tx1(1t)y1=tx(1t)y1+tx1(1t)yt^{x-1}(1-t)^{y-1} = t^{x}(1-t)^{y-1} + t^{x-1}(1-t)^{y} deduzca B(x,y)=B(x+1,y)+B(x,y+1)B(x,y) = B(x+1,y) + B(x,y+1) y combine con la pregunta 7 para obtener el relaciones de ascendencia

    B(x,y+1)=yx+yB(x,y),B(x+1,y)=xx+yB(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. Deducir, para números enteros m,n1m, n \geq 1:

    B(m,n)=(m1)!(n1)!(m+n1)!=1(m+n1)(m+n2m1).B(m, n) = \frac{(m-1)!\,(n-1)!}{(m+n-1)!} = \frac{1}{(m+n-1)\binom{m+n-2}{m-1}} .
  5. Demuestre Fórmula de Euler para un argumento entero: para cada x>0x > 0 y nNn \in \N^*,

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

    (induction on nn: both sides equal 1x\frac1x at n=1n = 1 and obey the same descent relation).

Parte III — Wallis integrals as Beta values.

  1. Sustituya t=sin2θt = \sin^2\theta para obtener la ecuación trigonométrica. forma

    B(x,y)=20π/2sin2x1θcos2y1θ ⁣dθ.B(x, y) = 2\int_0^{\pi/2} \sin^{2x-1}\theta\,\cos^{2y-1}\theta\,\dd\theta .
  2. Deducir Wn=12B(n+12,12)W_n = \frac12\,B\bigl(\frac{n+1}2, \frac12\bigr) para la integral de Wallis Wn=0π/2sinnθ ⁣dθW_n = \int_0^{\pi/2}\sin^n \theta\,\dd\theta, y recuperar la recurrencia Wn=n1nWn2W_n = \frac{n-1}{n}W_{n-2} de Lema 6.11 de las relaciones de descendencia de la pregunta 8 únicamente.
  3. Calcular B(12,12)=2W0=πB\bigl(\frac12, \frac12\bigr) = 2W_0 = \pi y compruébelo con Γ(12)2/Γ(1)\Gamma\bigl(\frac12\bigr)^2/\Gamma(1): La fórmula de Euler se mantiene en (12,12)\bigl(\frac12, \frac12\bigr).
  4. Derive la forma cerrada W2n=π2(2n)!4n(n!)2W_{2n} = \frac\pi2\, \frac{(2n)!}{4^n(n!)^2} de la recurrencia y 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)} .

    Concluimos, por inducción con las relaciones de descendencia, que La fórmula de Euler B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} se cumple siempre que 2x2x y 2y2y son números enteros positivos.

  5. Sustituya u=t1tu = \frac{t}{1-t} para obtener el tercero. forma clásica

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

    y verifique el caso x=y=12x = y = \frac12 directamente (u=v2u = v^2 lo reduce a 02 ⁣dv1+v2\int_0^\infty\frac{2\,\dd v}{1+v^2}).

Parte IV — Gauss’s limit formula.

  1. Para x>0x > 0 y nNn \in \N^*, probar por nn sucesivos integraciones por partes:

    0n(1tn) ⁣ntx1 ⁣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. Concluye con Ejercicio 9.7 (dominado convergencia) Fórmula del límite de Gauss:

    Γ(x)=limnn!  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, muestra que para x>0x > 0:

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

    donde γ\gamma es la constante de Euler (Ejemplo 6.7); justificar el convergencia de la serie (the general term is x22k2\sim \frac{x^2}{2k^2}).

  4. Utilice la fórmula de Gauss en x=12x = \frac12 y la central asintóticas binomiales (2nn)4nπn\binom{2n}{n} \sim \frac{4^n}{\sqrt{\pi n}} (Ejemplo 6.14) para volver a calcular Γ(12)=π\Gamma\bigl(\frac12\bigr) = \sqrt\pi: Stirling constante y la integral gaussiana son el mismo número en dos disfraces.
  5. Comprobar que la fórmula de Gauss reprueba la funcionalidad ecuación: de la identidad exacta

    n!nx+1(x+1)(x+n+1)=n!nxx(x+1)(x+n)nxx+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},

    Concluya Γ(x+1)=xΓ(x)\Gamma(x+1) = x\,\Gamma(x) nuevamente. (Gauss la fórmula determina Γ\Gamma directamente; el volumen del año 3 demuestra el teorema de Bohr-Mollerup, más agudo: el ecuación funcional más log-convexidad ya fijada Γ\Gamma abajo.)

Part V — Dividends.

  1. Para a>0a > 0 muestre 0eta ⁣dt=Γ(1+1a)\int_0^\infty \eu^{-t^a}\dd t = \Gamma\bigl(1 + \frac1a\bigr) y calcule el límite como a+a \to +\infty por convergencia dominada (pointwise limit 1t<1\mathbf 1_{t < 1}; dominate by 11 on (0,1]\intoc{0}{1} and by et2\eu^{-t^2} beyond, for a2a \geq 2). Verifique la respuesta con el continuidad de Γ\Gamma.
  2. Para espectáculo n1n \geq 1

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

    y recuperar los valores 22 (n=1n = 1) y π2\frac\pi2 (n=2n = 2). (Para n=4n = 4 este es el lemniscata constante, que no tiene forma cerrada elemental; su La historia pertenece a la teoría de las integrales elípticas.)

  3. (Momentos) Para x>0x > 0 y kNk \in \N, mostrar

    1Γ(x)0tktx1et ⁣dt=Γ(x+k)Γ(x)=x(x+1)(x+k1),\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),

    el factorial en ascenso; comprobar que x=1x = 1 da k!k!. (En En los capítulos de probabilidad, este es el kk-ésimo momento de un densidad de tiempo de espera estándar.)

  4. Probar la identidad Beta, válida para todo x>0x > 0:

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

    (substitute t=1+s2t = \frac{1+s}2, exploit the symmetry in ss, then set s=vs = \sqrt v). Deducir, por 2xN2x \in \N^*, La fórmula de duplicación de Legendre

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

    y verifíquelo directamente en x=nx = n a través de la pregunta 3. (Para general xx se sigue de la misma identidad una vez La fórmula de Euler es conocida por todos los argumentos — el prueba de doble integral en el capítulo sobre múltiples integrales.)

  5. Síntesis. En una oración cada uno: (i) cuando la integración por partes llevó la totalidad de la Parte II; (ii) donde la convergencia dominada entró en las Partes IV y V; (iii) cuál Las entradas asintóticas se importaron de la comparación. capítulo; (iv) lo que ahora se demuestra de la fórmula de Euler B(x,y)=Γ(x)Γ(y)/Γ(x+y)B(x,y) = \Gamma(x)\Gamma(y)/\Gamma(x+y), y que Queda por resolver la integral doble.
Solución

Solución de Problema 9.1.

1. En 0+0^+ el integrando es tx1\sim t^{x-1}: el escala de punto final finito converge iff 1x<11 - x < 1, es decir, x>0x > 0 (y para x0x \leq 0, tx1t1t^{x-1} \geq t^{-1} diverge); en ++\infty, tx1et=o(t2)t^{x-1}\eu^{-t} = o(t^{-2}) converge para cada xx. Luego Γ(x)=Γ(x+1)x\Gamma(x) = \frac{\Gamma(x+1)}{x} y Γ(x+1)Γ(1)=1\Gamma(x+1) \to \Gamma(1) = 1 como x0+x \to 0^+ (continuidad, Teorema 9.18): Γ(x)1x\Gamma(x) \sim \frac1x.

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

Γ(12)=0t1/2et ⁣dt=0eu2u2u ⁣du=20eu2 ⁣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

por Ejercicio 9.8. Con u=v/2u = v/\sqrt2:

Rev2/2 ⁣dv=220eu2 ⁣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. Verdadero para n=0n = 0 (ambos lados π\sqrt\pi). si Γ(n+12)=(2n)!4nn!π\Gamma(n + \frac12) = \frac{(2n)!}{4^n n!}\sqrt\pi, el ecuación funcional da

Γ(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 ,

el último paso porque (2n+2)!(2n)!=(2n+2)(2n+1)\frac{(2n+2)!}{(2n)!} = (2n+2)(2n+1) y 2n+12=(2n+2)(2n+1)4(n+1)\frac{2n+1}{2} = \frac{(2n+2)(2n+1)}{4(n+1)}.

4. Teorema 9.18 da Γ(x)=0tx1et(lnt)2 ⁣dt\Gamma''(x) = \int_0^\infty t^{x-1}\eu^{-t}(\ln t)^2\dd t (dos aplicaciones de la regla de Leibniz, dominaciones como en el demostración del teorema); el integrando es 0\geq 0 y no es idéntico cero, entonces Γ>0\Gamma'' > 0: Γ\Gamma es estrictamente convexo y Γ\Gamma' es estrictamente creciente. Desde Γ(1)=Γ(2)=1\Gamma(1) = \Gamma(2) = 1, Rolle proporciona x0(1,2)x_0 \in \intoo12 con Γ(x0)=0\Gamma'(x_0) = 0; La estricta monotonicidad de Γ\Gamma'hace que x0x_0 sea único. cero, con Γ<0\Gamma' < 0 antes y Γ>0\Gamma' > 0 después: Γ\Gamma disminuye en (0,x0)\intoo0{x_0}, aumenta en (x0,)\intoo{x_0}\infty y x0x_0 es el mínimo único.

5. Dejemos que kNk \in \N y x3x \geq 3; elige el numero entero nn con n+1x<n+2n + 1 \leq x < n + 2 (es decir, n1n \geq 1). por el monotonía de la pregunta 4 (válida desde x0<2x_0 < 2 en adelante): Γ(x)Γ(n+1)=n!\Gamma(x) \geq \Gamma(n + 1) = n!, mientras que xk(n+2)kx^k \leq (n+2)^k. desde n!(n+2)k\frac{n!}{(n+2)^k} \to \infty (factoriales vencen a las potencias, Año 1 volumen), Γ(x)xkn!(n+2)k\frac{\Gamma(x)}{x^k} \geq \frac{n!}{(n+2)^k} \to \inftycomo xx \to \infty: xk=o(Γ(x))x^k = o(\Gamma(x)).

6. Cerca de 00 el integrando es tx1\sim t^{x-1} (iff convergente x>0x > 0), cerca de 11 es (1t)y1\sim (1-t)^{y-1} (iff y>0y > 0); ambas comparaciones son entre funciones positivas, entonces B(x,y)B(x,y) converge exactamente para x,y>0x, y > 0. la sustitucion t1tt \mapsto 1 - t intercambia los dos factores: B(x,y)=B(y,x)B(x,y) = B(y,x).

7. B(x,1)=01tx1 ⁣dt=1xB(x,1) = \int_0^1 t^{x-1}\dd t = \frac1x. Piezas en [ε,1ε]\intcc\varepsilon{1-\varepsilon} con u=(1t)yu = (1-t)^y, v=txxv = \frac{t^x}{x}:

tx1(1t)y ⁣dt=[tx(1t)yx]+yxtx(1t)y1 ⁣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 ;

el soporte desaparece en ambos extremos como ε0\varepsilon \to 0 (x>0x > 0 en 00, y>0y > 0 en 11), dejando B(x,y+1)=yxB(x+1,y)B(x, y+1) = \frac yx\,B(x+1, y).

8. Desde t+(1t)=1t + (1-t) = 1:

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

entonces B(x,y)=B(x+1,y)+B(x,y+1)B(x,y) = B(x+1,y) + B(x,y+1). La pregunta 7 dice B(x+1,y)=xyB(x,y+1)B(x+1,y) = \frac xy B(x,y+1); sustituyendo,

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

es decir B(x,y+1)=yx+yB(x,y)B(x,y+1) = \frac{y}{x+y}B(x,y); la relación gemela sigue por la simetría de la pregunta 6.

9. Inducción en nn en mm fijo: B(m,1)=1m=(m1)!0!m!B(m,1) = \frac1m = \frac{(m-1)!\,0!}{m!}, y si la fórmula se mantiene en nn,

B(m,n+1)=nm+nB(m,n)=nm+n(m1)!(n1)!(m+n1)!=(m1)!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)!} .

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

10. Ambos lados de B(x,n)=Γ(x)Γ(n)Γ(x+n)B(x,n) = \frac{\Gamma(x)\Gamma(n)}{\Gamma(x+n)} son iguales a 1x\frac1xyn=1n = 1(Γ(1)=1\Gamma(1) = 1, Γ(x+1)=xΓ(x)\Gamma(x+1) = x\Gamma(x)). si estan de acuerdo en nn, luego por la relación de descenso y el funcional ecuación:

B(x,n+1)=nx+nB(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)} :

las dos secuencias obedecen a la misma recursividad de la misma semilla, por lo tanto, acepta para todos los nNn \in \N^* y todos los x>0x > 0.

11. Con t=sin2θt = \sin^2\theta (θ(0,π/2)\theta \in \intoo0{\pi/2},  ⁣dt=2sinθcosθ ⁣dθ\dd t = 2\sin\theta\cos\theta\,\dd\theta), tx1=sin2x2θt^{x-1} = \sin^{2x-2}\theta y (1t)y1=cos2y2θ(1-t)^{y-1} = \cos^{2y-2}\theta:

B(x,y)=0π/2sin2x2θcos2y2θ2sinθcosθ ⁣dθ=20π/2sin2x1θcos2y1θ ⁣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 (eliminando el factor coseno) y 2x1=n2x - 1 = n: B(n+12,12)=2WnB\bigl(\frac{n+1}2, \frac12\bigr) = 2W_n, es decir Wn=12B(n+12,12)W_n = \frac12 B\bigl(\frac{n+1}2,\frac12\bigr). el descenso relación en la primera variable da

WnWn2=B(n12+1,12)B(n12,12)=n12n12+12=n1n:\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} :

la recurrencia de Wallis, esta vez sin integración por partes en senos — La Parte II hizo el trabajo de una vez por todas.

13. B(12,12)=2W0=2π2=πB\bigl(\frac12,\frac12\bigr) = 2W_0 = 2\cdot\frac\pi2 = \pi, mientras Γ(12)2/Γ(1)=(π)2=π\Gamma\bigl(\frac12\bigr)^2/\Gamma(1) = (\sqrt\pi)^2 = \pi: La fórmula de Euler se mantiene en (12,12)\bigl(\frac12,\frac12\bigr).

14. Iterando W2n=2n12nW2n2W_{2n} = \frac{2n-1}{2n}W_{2n-2} desde W0=π2W_0 = \frac\pi2:

W2n=π2k=1n2k12k=π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},

desde (2k1)=(2n)!2nn!\prod(2k-1) = \frac{(2n)!}{2^n n!} y 2k=2nn!\prod 2k = 2^n n!. Por lo tanto, utilizando la pregunta 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)} .

Ahora arregle x12Nx \in \frac12\N^*. La fórmula de Euler se cumple en (x,12)(x, \frac12): para el número entero xx esta es la pregunta 10 (con simetría), para x=n+12x = n + \frac12 es la pantalla de arriba. ambos lados de La fórmula de Euler obedece a la recursividad descendente yy+1y \mapsto y + 1. (pregunta 8 a la izquierda, la ecuación funcional a la derecha, como en la pregunta 10): la inducción propaga la fórmula de y=12y = \frac12y y=1y = 1a cada y12Ny \in \frac12\N^*. Euler por lo tanto, la fórmula es válida siempre que 2x,2yN2x, 2y \in \N^*.

15. Con u=t1tu = \frac{t}{1-t}, es decir t=u1+ut = \frac{u}{1+u}, 1t=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)x1(11+u)y1 ⁣du(1+u)2=0ux1(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 .

En x=y=12x = y = \frac12, con u=v2u = v^2:

0u1/21+u ⁣du=02 ⁣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. Una integración por piezas, para 1kn1 \leq k \leq n y s>0s > 0 (u=(1t/n)ku = (1 - t/n)^k, v=ts/sv = t^s/s; los términos límite desaparecer):

0n(1tn) ⁣kts1 ⁣dt=kns0n(1tn) ⁣k1ts ⁣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 .

A partir de k=nk = n, s=xs = x e iterando nn veces:

0n(1tn) ⁣ntx1 ⁣dt=n(n1)1nnx(x+1)(x+n1)0ntx+n1 ⁣dt=n!nnnx+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 es n!nxx(x+1)(x+n)\dfrac{n!\,n^x}{x(x+1)\cdots(x+n)}.

17. Por Ejercicio 9.7 el lado izquierdo tiende a Γ(x)\Gamma(x) (convergencia dominada con dominador tx1ett^{x-1}\eu^{-t}); el lado derecho es el cociente de Gauss:

Γ(x)=limnn!nxx(x+1)(x+n).\Gamma(x) = \lim_{n\to\infty} \frac{n!\,n^x}{x(x+1)\cdots(x+n)} .

18. Tomando logaritmos en el cociente de la pregunta 16 Gn(x)G_n(x) y dividiendo ln(x+k)=lnk+ln(1+x/k)\ln(x+k) = \ln k + \ln(1 + x/k) para k1k \geq 1:

lnGn(x)=xlnnlnxk=1nln(1+xk)=lnx+x(lnnHn)+k=1n(xkln(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 u0u \geq 0, uu22ln(1+u)uu - \frac{u^2}2 \leq \ln(1+u) \leq u, entonces el término general se encuentra en [0,x2/(2k2)]\intcc{0}{x^2/(2k^2)}: la serie converge (comparación con k2\sum k^{-2}). Desde lnnHnγ\ln n - H_n \to -\gamma(Ejemplo 6.7) y lnGn(x)lnΓ(x)\ln G_n(x) \to \ln\Gamma(x)(pregunta 17 y continuidad de ln\ln):

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

19. En x=12x = \frac12, el denominador es k=0n(k+12)=(2n+1)!22n+1n!\prod_{k=0}^n\bigl(k+\frac12\bigr) = \frac{(2n+1)!}{2^{2n+1}n!} (multiplica las mitades), entonces

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}} .

Con (2nn)4nπn\binom{2n}{n} \sim \frac{4^n}{\sqrt{\pi n}} (Ejemplo 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) .

El coeficiente binomial central π\sqrt\pi (que proviene de Wallis, por lo tanto de la constante de Stirling) y el gaussiano Los π\sqrt\pi de la integral son el mismo número.

20. La identidad es álgebra directa: multiplicar 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} y absorba xx en el producto, nn en nxn^x. Dejando nn \to \infty: el lado izquierdo tiende a Γ(x+1)\Gamma(x+1)(Gauss en x+1x+1), el lado derecho a Γ(x)x1\Gamma(x)\cdot x\cdot 1 desde nx+n+11\frac{n}{x+n+1} \to 1: Γ(x+1)=xΓ(x)\Gamma(x+1) = x\Gamma(x) — recuperado sin una sola integración por partes.

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

0eta ⁣dt=1a0u1a1eu ⁣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) .

Como a+a \to +\infty (a lo largo de cualquier secuencia): eta1\eu^{-t^a} \to 1 para 0<t<10 < t < 1, e1\to \eu^{-1} en t=1t = 1, 0\to 0 para t>1t > 1; para a2a \geq 2 domina por 1t1+et21t>1\mathbf 1_{t \leq 1} + \eu^{-t^2}\mathbf 1_{t > 1}(tat2t^a \geq t^2 para t1t \geq 1), integrables. Convergencia dominada: la integral tiende a 011 ⁣dt=1\int_0^1 1\,\dd t = 1 — como debe, desde Γ(1+1a)Γ(1)=1\Gamma(1 + \frac1a) \to \Gamma(1) = 1 por continuidad.

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

01 ⁣dt1tn=1n01u1n1(1u)1/2 ⁣du=1nB(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, coincidente 01 ⁣dt1t=2\int_0^1\frac{\dd t}{\sqrt{1-t}} = 2. n=2n = 2: 12B(12,12)=π2=arcsin1\frac12 B\bigl(\frac12,\frac12\bigr) = \frac\pi2 = \arcsin 1. Para n=4n = 4 el valor 14B(14,12)\frac14 B\bigl(\frac14,\frac12\bigr) es la constante de lemniscata: no hay forma cerrada elemental.

23. Iterando la ecuación funcional:

1Γ(x)0tx+k1et ⁣dt=Γ(x+k)Γ(x)=(x+k1)(x+k2)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 ,

el factorial ascendente con factores kk. En x=1x = 1: Γ(1+k)/Γ(1)=k!\Gamma(1+k)/\Gamma(1) = k!, los momentos de et\eu^{-t} de Ejercicio 9.2.

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

B(x,x)=11(1s24)x1 ⁣ds2=41x01(1s2)x1 ⁣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

(el integrando es par). Entonces s=vs = \sqrt v ( ⁣ds= ⁣dv2v\dd s = \frac{\dd v}{2\sqrt v}):

B(x,x)=41x201v1/2(1v)x1 ⁣dv=212xB(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 2xN2x \in \N^* cada argumento a la vista radica en 12N\frac12\N^*, por lo que la fórmula de Euler (pregunta 14) se aplica a ambos lados:

Γ(x)2Γ(2x)=212xΓ(12)Γ(x)Γ(x+12)Γ(x)Γ(x+12)=212xπ  Γ(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) .

Verificación directa en x=nx = n: el lado izquierdo es (n1)!(2n)!π4nn!=(2n)!π4nn(n-1)!\cdot \frac{(2n)!\sqrt\pi}{4^n n!} = \frac{(2n)!\sqrt\pi}{4^n n}, el lado derecho 24nπ(2n1)!=(2n)!π4nn2\cdot4^{-n}\sqrt\pi\,(2n-1)! = \frac{(2n)!\sqrt\pi}{4^n n}: igual.

25. (i) Integración por piezas producidas B(x,y+1)=yxB(x+1,y)B(x,y+1) = \frac yx B(x+1,y), identidad única a partir de la cual cada relación de descendencia, los valores enteros y semienteros, y la La recurrencia de Wallis todo fluye. (ii) La convergencia dominada se volvió las integrales elementales 0n(1t/n)ntx1\int_0^n(1-t/n)^n t^{x-1} en Γ(x)\Gamma(x) (fórmula de Gauss, pregunta 17) y calculó el límite aa \to \infty en la pregunta 21. (iii) De la comparación En el capítulo importamos la constante de Euler (lnnHnγ\ln n - H_n \to -\gamma, pregunta 18) y las asintóticas binomiales centrales. (pregunta 19) — es decir, la fórmula de Stirling disfrazada. (v) La fórmula de Euler B(x,y)=Γ(x)Γ(y)/Γ(x+y)B(x,y) = \Gamma(x)\Gamma(y)/\Gamma(x+y) es ahora probado para yNy \in \N^* con x>0x > 0 arbitrario (pregunta 10) y para todos los pares de semienteros (pregunta 14); el general el caso x,y>0x, y > 0 espera el cálculo de la integral doble del capítulo sobre integrales múltiples, que factoriza Γ(x)Γ(y)\Gamma(x)\Gamma(y) sobre un cuarto de plano.