Mathematics · Book 4 · Bachelor Year 2

Matemáticas universitarias — Grado 2

Matemáticas universitarias — Grado 2 · Bachelor Year 2

6Comparación de funciones

Análisis asintótico: el arte de reemplazar una cantidad complicada por uno simple más un error controlado — se inició en el Año 1 volumen con expansiones de Taylor. Este capítulo la convierte en una disciplina de propio: expansiones a lo largo de escalas generales, la serie–integral comparación con su fuerza asintótica total, la fórmula de Stirling (probada completamente), y el estudio sistemático de los conceptos implícitamente definidos. secuencias. Estas técnicas son el pan de cada día de los asintóticos. análisis, y cada capítulo posterior que estima algo — series, integrales, probabilidades — come de esta tabla.

6.1 Relaciones de comparación y escalas.

Definición 6.1

Cerca de un punto aa (aRa \in \R o ±\pm\infty), para funciones (o secuencias, con nn \to \infty): f=o(g)f = o(g), f=O(g)f = O(g), fgf \sim g como en el volumen del Año 1. Un comparación escala en aa es una familia de positivos funciones, comparables por pares, totalmente ordenadas por o()o(\cdot) — la escala estándar en ++\infty es

xα(lnx)β(α,βR),x^{\alpha} (\ln x)^{\beta} \qquad (\alpha, \beta \in \R),

ordenado lexicográficamente en (α,β)(\alpha, \beta), refinado cuando sea necesario por exponenciales eγx\eu^{\gamma x}.

Definición 6.2 (expansión asintótica)

ff admite el expansión asintótica

f=c1φ1+c2φ2++ckφk+o(φk)(φi+1=o(φi) in the scale)f = c_1 \varphi_1 + c_2\varphi_2 + \dots + c_k \varphi_k + o(\varphi_k) \qquad (\varphi_{i+1} = o(\varphi_i) \text{ in the scale})

cuando los restos sucesivos satisfacen las estimaciones mostradas. el los coeficientes son entonces únicos: c1=limf/φ1c_1 = \lim f/\varphi_1, y inductivamente ci+1=lim(fjicjφj)/φi+1c_{i+1} = \lim\,(f - \sum_{j \leq i} c_j\varphi_j)/\varphi_{i+1}.

Ejemplo 6.3

Las expansiones de Taylor son expansiones asintóticas a lo largo de la escala (xa)k(x - a)^ken aa. Pero la noción es estrictamente más amplia: en++\infty,

1xlnx=1x11lnxx=1x+lnxx2+o(lnxx2),\frac{1}{x - \ln x} = \frac1x \cdot \frac{1}{1 - \frac{\ln x}{x}} = \frac1x + \frac{\ln x}{x^2} + o\Bigl(\frac{\ln x}{x^2}\Bigr),

una expansión a lo largo de la escala mixta — no se aplica ningún teorema de Taylor, sólo la expansión geométrica y el cálculo de oo.

Ejemplo 6.4 (La báscula estándar está realmente ordenada)

El reclamo lexicográfico de Definición 6.1 necesita una línea de prueba por caso. Comparar xα(lnx)βx^{\alpha}(\ln x)^{\beta} y xα(lnx)βx^{\alpha'}(\ln x)^{\beta'} en ++\infty. Si α<α\alpha < \alpha': la relación es xαα(lnx)ββ0x^{\alpha - \alpha'}(\ln x)^{\beta - \beta'} \to 0, porque un La potencia negativa de xx aplasta cualquier potencia de lnx\ln x (establecer x=etx = \eu^t: e(αα)ttββ0\eu^{(\alpha - \alpha')t}\,t^{\beta - \beta'} \to 0 por el límite polinómico de latidos exponenciales del año 1 volumen). Si α=α\alpha = \alpha' y β<β\beta < \beta': la relación es (lnx)ββ0(\ln x)^{\beta - \beta'} \to 0 directamente. Entonces las parejas (α,β)(\alpha, \beta), ordenado lexicográficamente, ordena la escala por o()o(\cdot) — y la sustitución x=etx = \eu^t es la truco único para comparaciones de registros de potencia mixtos.

Ejemplo 6.5 (Clasificación de una colección de animales.)

Las escalas deben ser ordenado; Aquí está el taladro estándar. en ++\infty, comparar n10n^{10}, elnnn\eu^{\sqrt{\ln n}\,\cdot\,\sqrt n}, 2n2^n y nlnnn^{\ln n} tomando logaritmos:

10lnn    (lnn)2    nlnn    nln2,10\ln n \;\ll\; (\ln n)^2 \;\ll\; \sqrt{n\ln n} \;\ll\; n\ln 2 ,

donde anbna_n \ll b_n significa an=o(bn)a_n = o(b_n); la segunda entrada es ln(nlnn)\ln(n^{\ln n}). Los exponenciales preservan estas brechas estrictas (si lnunlnvn\ln u_n - \ln v_n \to -\infty luego un/vn0u_n/v_n \to 0), entonces

n10=o(nlnn),nlnn=o(enlnn),enlnn=o(2n).n^{10} = o\bigl(n^{\ln n}\bigr), \qquad n^{\ln n} = o\bigl(\eu^{\sqrt{n\ln n}}\bigr), \qquad \eu^{\sqrt{n\ln n}} = o(2^n) .

La moraleja, doble: comparar siempre mediante logaritmos (diferencias de registros, no proporciones de registros), y nunca concluir unvnu_n \sim v_n de lnunlnvn\ln u_n \sim \ln v_n — el par n10n^{10} y nlnnn^{\ln n} tiene una relación ln\ln que tiende a \infty, pero 2n2^n y 4n4^n tienen una relación ln\ln exactamente 22 y son tremendamente desigual.

6.2 Serie: comparación integral, asintóticamente

Teorema 6.6

Sea ff continuo, positivo, decreciente en [1,+)\intco{1}{+\infty}.

  1. Si 1f\int_1^{\infty} f converge, los restos satisfacen

    n+1f    k>nf(k)    nf.\int_{n+1}^{\infty} f \;\leq\; \sum_{k > n} f(k) \;\leq\; \int_{n}^{\infty} f .
  2. Si 1f\int_1^\infty f diverge, las sumas parciales satisfacen k=1nf(k)=1nf+C+o(1)\sum_{k=1}^{n} f(k) = \int_1^n f + C + o(1) para algunos constante CC: la diferencia knf(k)1nf\sum_{k \leq n} f(k) - \int_1^n f converge.

Demostración. El bracketing f(k+1)kk+1ff(k)f(k+1) \leq \int_k^{k+1} f \leq f(k) (disminución) fue el dispositivo del Año 1; sumando kn+1k \geq n+1, resp. knk \geq n, da (1). Para (2), configure uk=f(k)kk+1fu_k = f(k) - \int_k^{k+1} f: por el entre paréntesis, 0ukf(k)f(k+1)0 \leq u_k \leq f(k) - f(k+1), por lo que las sumas parciales de uk\sum u_k están delimitados por el f(1)f(n+1)f(1)f(1) - f(n+1) \leq f(1)telescópico: la serie converge. Además, la secuencia (nn+1f)n\bigl(\int_n^{n+1} f\bigr)_nno es creciente (ff disminuye) y no es negativa, por lo tanto convergente. escribiendo

k=1nf(k)1nf=k=1nuk+nn+1f,\sum_{k=1}^{n} f(k) - \int_1^n f = \sum_{k=1}^{n} u_k + \int_n^{n+1} f ,

el lado derecho converge como nn \to \infty: la diferencia converge a una constante CC, que es la declaración (2).

Ejemplo 6.7 (La expansión armónica)

Para f(t)=1tf(t) = \frac1t: Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1), recuperando Constante de Euler (volumen del año 1) con una prueba más limpia. Empujando uno ordenar más (Ejercicio 6.3):

Hn=lnn+γ+12n+o(1n).H_n = \ln n + \gamma + \frac{1}{2n} + o\Bigl(\frac1n\Bigr).

Los números hacen visible la ganancia en n=10n = 10: H10=2.928968H_{10} = 2.928968\dotsyln10=2.302585\ln 10 = 2.302585\dots, por lo que la estimación bruta de γ\gamma es H10ln10=0.626383H_{10} - \ln 10 = 0.626383, desviado por 0.0490.049; restando la corrección 120\frac1{20} se obtiene 0.5763830.576383, fuera de γ=0.577216\gamma = 0.577216 por solo 8.31048.3\cdot10^{-4} — que es a su vez el siguiente término 112100\frac{1}{12\cdot100} del expansión, como lo demuestra el problema del fin de semana (pregunta 8).

Ejemplo 6.8 (Un ln(n!)\ln(n!) tosco sin Stirling)

El dispositivo de fijación por sí solo ya localiza ln(n!)\ln(n!). desde ln\ln aumenta,

k1klnt ⁣dt    lnk    kk+1lnt ⁣dt,\int_{k-1}^{k}\ln t\,\dd t \;\leq\; \ln k \;\leq\; \int_{k}^{k+1}\ln t\,\dd t ,

y sumando k=2,,nk = 2, \dots, n (con 1nln=nlnnn+1\int_1^n\ln = n\ln n - n + 1):

nlnnn+1    ln(n!)    (n+1)ln(n+1)n.n\ln n - n + 1 \;\leq\; \ln(n!) \;\leq\; (n+1)\ln(n+1) - n .

Ambas vallas son nlnnn+O(lnn)n\ln n - n + O(\ln n): de ahí ln(n!)=nlnnn+O(lnn)\ln(n!) = n\ln n - n + O(\ln n), y en particular ln(n!)nlnn\ln(n!) \sim n\ln n. Lo que Stirling agrega son los siguientes dos peldaños: el 12lnn\frac12\ln n y la constante ln2π\ln\sqrt{2\pi} — que cuesta más fino telescópico de Teorema 6.13. sabiendo La precisión que adquiere cada herramienta es la mitad del arte asintótico.

Ejemplo 6.9 (Cancelación exige ampliaciones)

Calcule el límite de n2+nn\sqrt{n^2 + n} - n. Ambos términos son n\sim n y “nn\sim n - n” no tienen sentido: equivalentes no se puede restar. Expandir en su lugar:

n2+nn=n(1+1n1)=n(12n18n2+O(1n3))=1218n+O(1n2):\sqrt{n^2 + n} - n = n\Bigl(\sqrt{1 + \tfrac1n} - 1\Bigr) = n\Bigl(\frac{1}{2n} - \frac{1}{8n^2} + O\Bigl(\frac1{n^3}\Bigr)\Bigr) = \frac12 - \frac{1}{8n} + O\Bigl(\frac1{n^2}\Bigr) :

límite 12\frac12, con la velocidad de aproximación 18n\frac1{8n} como bonificación. El mecanismo merece un nombre: una diferencia de dos grandes cantidades equivalentes vive enteramente en su próximo términos, por lo que uno debe expandirse al primer orden en el que los dos lados difieren — y lleve el resto para certificar que nada más sobrevive en ese orden.

Ejemplo 6.10 (Una comparación divergente, funcionó.)

Para f(t)=1tlntf(t) = \frac{1}{t\ln t} en [2,+)\intco{2}{+\infty} (continuo, positivo, decreciente): 2xf=lnlnxlnln2\int_2^x f = \ln\ln x - \ln\ln 2 \to \infty, entonces por Teorema 6.6 (2),

k=2n1klnk=lnlnn+C+o(1)\sum_{k=2}^{n}\frac{1}{k\ln k} = \ln\ln n + C + o(1)

para alguna constante CC. Dos lecciones. Primero, la divergencia es real pero glacial: la suma parcial supera por primera vez 44 alrededor de nee4Cn \approx \eu^{\eu^{4 - C}}, astronómicamente grande. En segundo lugar, el forma lnlnn\ln\ln n fue entregado por una antiderivada, no adivinado: para términos monótonos, la integral es la canónica dispositivo sumador, y la constante CC — como γ\gamma de Euler — es la memoria de los términos iniciales.

6.3 La fórmula de Stirling.

Lema 6.11 (Integrales de Wallis, revisadas)

Deje Wn=0π/2sinnt ⁣dtW_n = \int_0^{\pi/2} \sin^n t\,\dd t. Luego nWnWn1=π2nW_nW_{n-1} = \frac\pi2 para n1n \geq 1, (Wn)(W_n)disminuye y Wnπ2nW_n \sim \sqrt{\dfrac{\pi}{2n}}.

Demostración. La integración por partes da nWn=(n1)Wn2nW_n = (n-1)W_{n-2} (n2n \geq 2), por lo que nWnWn1nW_nW_{n-1} es constante en nn, igual a 1W1W0=π21 \cdot W_1 W_0 = \frac\pi2. Disminución: sinn+1sinn\sin^{n+1} \leq \sin^n en [0,π2]\intcc{0}{\frac\pi2}. El apretón, en detalle: monotonicidad da Wn+1WnWn1W_{n+1} \leq W_n \leq W_{n-1}, y dividiendo por Wn1>0W_{n-1} > 0,

nn+1=Wn+1Wn1WnWn11,\frac{n}{n+1} = \frac{W_{n+1}}{W_{n-1}} \leq \frac{W_n}{W_{n-1}} \leq 1 ,

la identidad izquierda de la recurrencia en el índice n+1n + 1. ambos los límites tienden a 11: WnWn1W_n \sim W_{n-1}, de donde

nWn2nWnWn1=π2Wnπ2n.nW_n^2 \sim nW_nW_{n-1} = \frac\pi2 \qquad\Longrightarrow\qquad W_n \sim \sqrt{\frac{\pi}{2n}} .

Ejemplo 6.12 (Las primeras integrales de Wallis)

De W0=π2W_0 = \frac\pi2, W1=1W_1 = 1 y la recurrencia nWn=(n1)Wn2nW_n = (n-1)W_{n-2}:

W2=π4,W3=23,W4=3π16,W5=815,W6=5π32.W_2 = \frac\pi4, \qquad W_3 = \frac23, \qquad W_4 = \frac{3\pi}{16}, \qquad W_5 = \frac{8}{15}, \qquad W_6 = \frac{5\pi}{32}.

Los índices pares llevan un π\pi, los impares son racionales — los dos productos intercalados de las formas cerradas. Numéricamente W60.4909W_6 \approx 0.4909 contra lo asintótico π/120.5116\sqrt{\pi/12} \approx 0.5116: en n=6n = 6 el equivalente es ya dentro de 5%5\%, y la identidad del producto es exacta en cada nn: 6W6W5=65π32815=π26\,W_6W_5 = 6\cdot\frac{5\pi}{32}\cdot\frac8{15} = \frac\pi2. Las mesas pequeñas como esta son la forma más económica de atrapar un error de álgebra antes de que infecte un argumento asintótico.

Teorema 6.13 (Stirling)

n!    2πn(ne) ⁣n.n! \;\sim\; \sqrt{2\pi n}\, \Bigl(\frac{n}{\eu}\Bigr)^{\!n} .

Demostración. Step 1: n!Cn(n/e)nn! \sim C \sqrt n\, (n/\eu)^n for some constant C>0C > 0. Conjunto

dn=ln(n!)(n+12)lnn+n.d_n = \ln(n!) - \Bigl(n + \frac12\Bigr)\ln n + n .

entonces

dndn+1=(n+12)lnn+1n1=(n+12)(1n12n2+13n3+o(n3))1=112n2+o(1n2),d_n - d_{n+1} = \Bigl(n + \frac12\Bigr) \ln\frac{n+1}{n} - 1 = \Bigl(n + \frac12\Bigr)\Bigl(\frac1n - \frac{1}{2n^2} + \frac{1}{3n^3} + o\bigl(n^{-3}\bigr)\Bigr) - 1 = \frac{1}{12n^2} + o\Bigl(\frac{1}{n^2}\Bigr),

por la expansión de Taylor de ln(1+1n)\ln(1 + \frac1n). Por tanto, la serie (dndn+1)\sum (d_n - d_{n+1})converge absolutamente (comparación con n2\sum n^{-2}), por lo que (dn)(d_n)converge, digamos a dd; exponenciando, n!Cn(n/e)nn! \sim C\sqrt n\,(n/\eu)^ncon C=edC = \eu^{d}.

Step 2: C=2πC = \sqrt{2\pi} via Wallis. El formulario cerrado W2p=(2p)!4p(p!)2π2W_{2p} = \frac{(2p)!}{4^p (p!)^2}\cdot\frac\pi2 (de la recurrencia, Cálculo del año 1 rehecho en Lema 6.11 ajuste) se combina con el Paso 1:

W2pC2p(2p/e)2p4p(Cp(p/e)p)2π2=2pCpπ2=πC12p.W_{2p} \sim \frac{C\sqrt{2p}\,(2p/\eu)^{2p}} {4^p\,\bigl(C\sqrt p\,(p/\eu)^p\bigr)^2}\cdot\frac{\pi}{2} = \frac{\sqrt{2p}}{C\,p}\cdot\frac{\pi}{2} = \frac{\pi}{C}\cdot\frac{1}{\sqrt{2p}} .

Comparando con W2pπ4pW_{2p} \sim \sqrt{\frac{\pi}{4p}} (Lema 6.11): πC2p=π4p(1+o(1))\frac{\pi}{C\sqrt{2p}} = \sqrt{\frac{\pi}{4p}}\,(1 + o(1))fuerza a C=π4p2pπ=2πC = \pi \sqrt{\frac{4p}{2p\,\pi}} = \sqrt{2\pi}.

Ejemplo 6.14 (Coeficiente binomial central)

(2nn)=(2n)!(n!)24πn(2n/e)2n2πn(n/e)2n=4nπn:\binom{2n}{n} = \frac{(2n)!}{(n!)^2} \sim \frac{\sqrt{4\pi n}\,(2n/\eu)^{2n}}{2\pi n\,(n/\eu)^{2n}} = \frac{4^n}{\sqrt{\pi n}} :

la probabilidad de que un paseo aleatorio simétrico regrese a 00 en el momento 2n2n es 1πn\sim \frac{1}{\sqrt{\pi n}} — un anuncio de Capítulo 22.

Observación 6.15 (Perspectivas dentro de este volumen)

Cada capítulo cuantitativo que sigue habla de la esencia de este capítulo. idioma. Capítulo 7 clasifica series comparando términos a la escala nα(lnn)βn^{-\alpha}(\ln n)^{-\beta} — su El problema del fin de semana mapea esa frontera por completo. Capítulo 9 hace lo mismo para integrales impropias, con idéntica escala en la variable continuo. Capítulo 11 calcula radios de convergencia a partir de lim supan1/n\limsup\abs{a_n}^{1/n}, un equivalente de nn-ésimas raíces ejercicio donde Stirling es la clave estándar (n!nne\sqrt[n]{n!} \sim \frac n\eu, Ejercicio 6.4). y el capítulos de probabilidad en efectivo Stirling directamente: el local estimaciones de Capítulo 22 para coeficientes binomiales son Ejemplo 6.14 y Ejemplo 6.21 palabra por palabra. La asintótica no es un capítulo aquí; es el acento del volumen.

Método 6.16 (La lista de verificación de arranque)

Antes de confiar en una expansión iniciada, audite cuatro puntos. (1) La existencia primero: la raíz o secuencia debe estar fijada hacia abajo (monotonicidad, valores intermedios) antes de cualquier expansión — los símbolos sin referentes se expanden bellamente y significan nada. (2) Un pedido por pase: cada sustitución puede sólo se debe confiar en el orden de la estimación ingresada; extrayendo dos nuevos términos de una sola pasada es la fuente clásica de errores coeficientes. (3) Los restos siguen adelante: lleva el o()o(\cdot) a través de cada paso algebraico y dejar que la absorción (términos más pequeños absorbidos por restos más grandes) ocurren en el final, explícitamente. (4) Auditoría numérica: evaluar a la vez valor honesto de nn; un error de coeficiente sobrevive al algebraico La rederivación es sorprendentemente frecuente y casi nunca sobrevive. aritmética.

Observación 6.17 (Errores comunes)

(i) Equivalentes agregar mal: de unn+lnnu_n \sim n + \ln n y vnnv_n \sim -n se puede no concluir un+vnlnnu_n + v_n \sim \ln n; las cancelaciones exigen ampliaciones con restos explícitos, nunca desnudos equivalentes. (ii) Nunca exponenciar una equivalencia: n+1nn + 1 \sim n pero en+1≁en\eu^{n+1} \not\sim \eu^n; la caja fuerte dirección es tomar logaritmos de equivalentes que tienden a ++\infty (problema de fin de semana de este capítulo, pregunta 24). (iii) Un expansión asintótica está conectado a un escala: escritura f=1x+o(1x2)f = \frac1x + o\bigl(\frac1{x^2}\bigr) reclama más de f=1x+o(1x)f = \frac1x + o\bigl(\frac1x\bigr) y mezclar los dos invalida álgebra posterior. (iv) En bootstraps, sustituir el entero expansión actual, resto incluido — cayendo un pase medio o()o(\cdot) produce resultados plausibles pero incorrectos. coeficientes. (v) La serie: la comparación integral necesita monotonicidad: para términos oscilantes falla rotundamente (compárese sinkk\sum\frac{\sin k}k, Capítulo 7).

Ejemplo 6.18 (Stirling en números)

En n=10n = 10: la fórmula da 20π(10/e)103598696\sqrt{20\pi}\,(10/\eu)^{10} \approx 3\,598\,696 contra 10!=362880010! = 3\,628\,800: error relativo 8.31038.3\cdot10^{-3}, notable por una declaración "asintótica" en n=10n = 10. El error tiene una estructura — el refinamiento exacto n!=2πn(n/e)n(1+112n+O(n2))n! = \sqrt{2\pi n}\,(n/\eu)^n\bigl(1 + \frac1{12n} + O(n^{-2})\bigr)— cuya primera corrección 11208.3103\frac1{120} \approx 8.3\cdot10^{-3} explica casi exactamente la brecha observada. el La maquinaria Euler-Maclaurin del problema del fin de semana es precisamente la fuente sistemática de tales términos de corrección.

Observación 6.19 (Dónde se utiliza este capítulo)

La comparación asintótica es la gramática de todo lo cuantitativo. downstream: las pruebas de convergencia y el panorama Bertrand de Capítulo 7, los criterios de integrabilidad de Capítulo 9, los cálculos del radio de convergencia de Capítulo 11, y los teoremas de límite de Capítulo 22 (donde Stirling dirige el de Moivre–estimaciones de Laplace). El volumen del Año 3 industrializa el una idea que probamos a mano aquí: extraiga el término principal, cota el resto — al método de Laplace y dominó la convergencia.

Ejemplo 6.20 (Una integral comparada consigo misma: 2x ⁣dtlnt\int_2^x \frac{\dd t}{\ln t})

La caja de herramientas de comparación también funciona con integrales. Sea F(x)=2x ⁣dtlntF(x) = \int_2^x\frac{\dd t}{\ln t} (el integrando es continuo en [2,)\intco2\infty). Integrar por partes:

F(x)=[tlnt]2x+2x ⁣dt(lnt)2=xlnx+O(2x ⁣dt(lnt)2)+O(1),F(x) = \Bigl[\frac{t}{\ln t}\Bigr]_2^x + \int_2^x\frac{\dd t}{(\ln t)^2} = \frac{x}{\ln x} + O\Bigl(\int_2^x\frac{\dd t}{(\ln t)^2}\Bigr) + O(1),

y la integral restante es o(xlnx)o\bigl(\frac{x}{\ln x}\bigr): dividirlo en x\sqrt x, delimitado por

2x ⁣dt(lnt)2xandxx ⁣dt(lnt)2x(lnx)2=4x(lnx)2.\int_2^{\sqrt x}\frac{\dd t}{(\ln t)^2} \leq \sqrt x \qquad\text{and}\qquad \int_{\sqrt x}^{x}\frac{\dd t}{(\ln t)^2} \leq \frac{x}{(\ln\sqrt x)^2} = \frac{4x}{(\ln x)^2} .

Por lo tanto F(x)xlnxF(x) \sim \frac{x}{\ln x}. Lectores que conocieron al mejor teorema de números en el problema de fin de semana de este capítulo reconocerá FF: es la integral logarítmica, el mejor estimador de π(x)\pi(x), y el cálculo muestra que concuerda con xlnx\frac{x}{\ln x} al primer pedido.

Ejemplo 6.21 (Stirling en un binomio torcido)

La misma rutina trifactorial que para Ejemplo 6.14 da, para (3nn)=(3n)!n!(2n)!\binom{3n}{n} = \frac{(3n)!}{n!\,(2n)!}:

(3nn)6πn(3n/e)3n2πn(n/e)n4πn(2n/e)2n=34πn(274) ⁣n.\binom{3n}{n} \sim \frac{\sqrt{6\pi n}\,(3n/\eu)^{3n}} {\sqrt{2\pi n}\,(n/\eu)^{n}\cdot\sqrt{4\pi n}\,(2n/\eu)^{2n}} = \sqrt{\frac{3}{4\pi n}}\, \Bigl(\frac{27}{4}\Bigr)^{\!n} .

La tasa exponencial 274=3322\frac{27}4 = \frac{3^3}{2^2} es e3nH(1/3)\eu^{3n\,H(1/3)} en la notación entrópica de información. teoría: los binomios asimétricos crecen estrictamente más lento que el central 4n4^n por dos pasos — aquí (27/4)1/31.89<2(27/4)^{1/3} \approx 1.89 < 2 por paso. Cada binomio asintótico en combinatoria y probabilidad. (Capítulo 22) es este cálculo con diferentes pesos.

6.4 Secuencias definidas implícitamente

Método 6.22

Para encontrar las asintóticas de las soluciones xnx_n de una ecuación F(x,n)=0F(x, n) = 0:

  1. Localizar: demostrar la existencia y unicidad de xnx_n en un intervalo definido (monotonicidad, valor intermedio teorema), y encontrar su comportamiento crudo (límite, orden de crecimiento).
  2. Oreja: sustituye la forma cruda xn=(main term)(1+εn)x_n = (\text{main term})(1 + \varepsilon_n) en la ecuación y resuelve el siguiente orden de εn\varepsilon_n; repetir, cada pase refinando un pedido.

Ejemplo 6.23

Para n1n \geq 1, la ecuación tanx=x\tan x = x tiene exactamente una solución xnx_n en (nππ2,nπ+π2)\intoo{n\pi - \frac\pi2}{n\pi + \frac\pi2} (la función tanxx\tan x - x aumenta de -\infty a ++\infty allí, es siendo el derivado tan2x0\tan^2 x \geq 0). Crudo: xn=nπ+π2ynx_n = n\pi + \frac\pi2 - y_ncon yn(0,π)y_n \in \intoo{0}{\pi}; ya que xnx_n \to \inftyytanxn=xn+\tan x_n = x_n \to +\infty, xnx_n se acerca al asíntota desde la izquierda: yn0y_n \to 0. Oreja: tanxn=cotyn=1tanyn1yn\tan x_n = \cot y_n = \frac{1}{\tan y_n} \sim \frac{1}{y_n}, y la ecuación cotyn=xnnπ\cot y_n = x_n \sim n\pi da yn1nπy_n \sim \frac{1}{n\pi}. Por lo tanto

xn=nπ+π21nπ+o(1n),x_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + o\Bigl(\frac1n\Bigr),

y el proceso continúa a cualquier orden (Ejercicio 6.6).

Ejemplo 6.24 (Una segunda ejecución del método)

Resuelva x+lnx=nx + \ln x = n asintóticamente. Localizar: xx+lnxx \mapsto x + \ln xaumenta de-\inftya++\infty en (0,+)\intoo{0}{+\infty}: una raíz única xnx_n y xnx_n \to \infty. Crudo: lnxn=o(xn)\ln x_n = o(x_n) da xnnx_n \sim n. Oreja: de xn=nlnxnx_n = n - \ln x_n y lnxn=lnn+o(1)\ln x_n = \ln n + o(1)(logaritmos de equivalentes, ambos lados \to \infty):

xn=nlnn+o(1);x_n = n - \ln n + o(1) ;

una pasada más, con lnxn=ln(nlnn+o(1))=lnnlnnn+o(lnnn)\ln x_n = \ln\bigl(n - \ln n + o(1)\bigr) = \ln n - \frac{\ln n}{n} + o\bigl(\frac{\ln n}n\bigr):

xn=nlnn+lnnn+o(lnnn).x_n = n - \ln n + \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr).

(Verifique en n=100n = 100: la raíz es x95.4415x \approx 95.4415; la La fórmula de tres términos da 1004.6052+0.0461=95.4409100 - 4.6052 + 0.0461 = 95.4409, el dos términos uno 95.394895.3948 — cada pase gana el pronóstico orden.) Mismo bucle, tercer paisaje: el método de Método 6.22 no le importa lo que el parece la ecuación, sólo que cada pasada aísla la dominante desconocido.

6.5 Ceremonias

Ejercicio 6.1

Expandir en ++\infty, dos términos más allá del principal:

x2+x+1,ln(x2+x)2lnx,x+sinxxlnx.\sqrt{x^2 + x + 1} , \qquad \ln(x^2 + x) - 2\ln x, \qquad \frac{x + \sin x}{x - \ln x} .
Solución

Solución de Ejercicio 6.1.

x2+x+1=x1+1x+1x2=x+12+381x+o(1x)\sqrt{x^2 + x + 1} = x\sqrt{1 + \tfrac1x + \tfrac{1}{x^2}} = x + \frac12 + \frac38\cdot\frac1x + o\bigl(\frac1x\bigr) (binomial expansión: 12u18u2\frac12 u - \frac18 u^2 con u=1x+1x2u = \frac1x + \frac{1}{x^2} da 12x+12x218x2=12x+38x2\frac{1}{2x} + \frac{1}{2x^2} - \frac{1}{8x^2} = \frac{1}{2x} + \frac{3}{8x^2}, luego multiplica por xx).

ln(x2+x)2lnx=ln(1+1x)=1x12x2+o(1x2)\ln(x^2 + x) - 2\ln x = \ln\bigl(1 + \tfrac1x\bigr) = \frac1x - \frac{1}{2x^2} + o\bigl(\frac{1}{x^2}\bigr).

Tercera función: expandir cada factor,

x+sinxxlnx=(1+sinxx)(1+lnxx+(lnx)2x2+O((lnx)3x3)).\frac{x + \sin x}{x - \ln x} = \Bigl(1 + \frac{\sin x}{x}\Bigr) \Bigl(1 + \frac{\ln x}{x} + \frac{(\ln x)^2}{x^2} + O\Bigl(\frac{(\ln x)^3}{x^3}\Bigr)\Bigr).

Ordene las aportaciones en la escala ++\infty: lnxx1xsinxx(lnx)2x2\frac{\ln x}{x} \gg \frac{1}{x} \geq \bigl|\frac{\sin x}{x}\bigr| \gg \frac{(\ln x)^2}{x^2}. Por lo tanto, los dos términos después del 11 inicial son lnxx\frac{\ln x}{x}, luego el término de oscilación acotada sinxx\frac{\sin x}{x}:

x+sinxxlnx=1+lnxx+sinxx+O((lnx)2x2).\frac{x + \sin x}{x - \ln x} = 1 + \frac{\ln x}{x} + \frac{\sin x}{x} + O\Bigl(\frac{(\ln x)^2}{x^2}\Bigr).

Ejercicio 6.2

Dé la naturaleza (convergencia/divergencia) y, cuando sea divergente, la asintóticas principales de knkα\sum_{k \leq n} k^\alpha para α>1\alpha > -1, α=1\alpha = -1, α<1\alpha < -1, vía Teorema 6.6.

Solución

Solución de Ejercicio 6.2.

f(t)=tαf(t) = t^\alpha (t1t \geq 1).

α>1\alpha > -1: divergencia, y por Teorema 6.6 (2), knkα=nα+1α+1+C+o(1)\sum_{k\leq n} k^\alpha = \frac{n^{\alpha+1}}{\alpha+1} + C + o(1)si α<0\alpha < 0(donde ffdisminuye); para α0\alpha \geq 0(ff en aumento) el El mismo corchetes con desigualdades invertidas da knkαnα+1α+1\sum_{k \leq n} k^\alpha \sim \frac{n^{\alpha + 1}}{\alpha + 1}.

α=1\alpha = -1: Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1) (Ejemplo 6.7).

α<1\alpha < -1: convergencia, con resto k>nkαnα+1(α+1)\sum_{k > n} k^\alpha \sim \frac{n^{\alpha+1}}{-(\alpha+1)} por el corchetes (1) (ambos los límites integrales son equivalentes a ese valor).

Ejercicio 6.3 ★★

Prueba Hn=lnn+γ+12n+o(1n)H_n = \ln n + \gamma + \frac{1}{2n} + o\bigl(\frac1n\bigr). (Study vn=Hnlnnγv_n = H_n - \ln n - \gamma: show vnvn+1=12n2+O(n3)v_n - v_{n+1} = \frac{1}{2n^2} + O(n^{-3}) and sum the tail, comparing with kn12k212n\sum_{k \geq n} \frac{1}{2k^2} \sim \frac{1}{2n}Teorema 6.6 (1).)

Solución

Solución de Ejercicio 6.3.

Deje vn=Hnlnnγ0v_n = H_n - \ln n - \gamma \to 0. entonces

vnvn+1=lnn+1n1n+1=(1n12n2)(1n1n2)+O(1n3)=12n2+O(1n3),v_n - v_{n+1} = \ln\frac{n+1}{n} - \frac{1}{n+1} = \Bigl(\frac1n - \frac{1}{2n^2}\Bigr) - \Bigl(\frac1n - \frac{1}{n^2}\Bigr) + O\Bigl(\frac{1}{n^3}\Bigr) = \frac{1}{2n^2} + O\Bigl(\frac{1}{n^3}\Bigr),

utilizando 1n+1=1n1n2+O(n3)\frac{1}{n+1} = \frac1n - \frac{1}{n^2} + O(n^{-3}). desde vn0v_n \to 0, telescópicando la cola:

vn=kn(vkvk+1)=kn(12k2+O(k3))=12n+O(1n2),v_n = \sum_{k \geq n} (v_k - v_{k+1}) = \sum_{k\geq n} \Bigl(\frac{1}{2k^2} + O(k^{-3})\Bigr) = \frac{1}{2n} + O\Bigl(\frac{1}{n^2}\Bigr),

por Teorema 6.6 (1) aplicado a t2t^{-2} (resto 1n\sim \frac1n, reducido a la mitad) y a t3t^{-3}. Por lo tanto Hn=lnn+γ+12n+o(1n)H_n = \ln n + \gamma + \frac{1}{2n} + o(\frac1n).

Ejercicio 6.4 ★★

Usando Stirling, encuentre equivalentes de: (3n)!(n!)3\dfrac{(3n)!}{(n!)^3};   n!nn\;\dfrac{n!}{n^n};   n!n\;\sqrt[n]{n!} (como ne(1+o(1))\frac n\eu(1 + o(1)), precisa dos términos).

Solución

Solución de Ejercicio 6.4.

Stirling tres veces:

(3n)!(n!)36πn(3n/e)3n(2πn)3/2(n/e)3n=6  27n2πn12πn2πn  =327n2πn.\frac{(3n)!}{(n!)^3} \sim \frac{\sqrt{6\pi n}\,(3n/\eu)^{3n}} {(2\pi n)^{3/2}\,(n/\eu)^{3n}} = \frac{\sqrt{6}\; 27^{\,n}}{2\pi n} \cdot \frac{1}{\sqrt{2\pi n}}\cdot\sqrt{2\pi n}\; = \frac{\sqrt3\,27^n}{2\pi n} .

(Con atención: 6πn(2πn)3/2=6(2πn)2πnπn=32πn\frac{\sqrt{6\pi n}}{(2\pi n)^{3/2}} = \frac{\sqrt6}{(2\pi n)\sqrt{2\pi n}}\sqrt{\pi n} = \frac{\sqrt3}{2\pi n}.)

n!nn2πnen\dfrac{n!}{n^n} \sim \sqrt{2\pi n}\,\eu^{-n}.

n!n=exp(lnn!n)\sqrt[n]{n!} = \exp\bigl(\frac{\ln n!}{n}\bigr) con lnn!=nlnnn+12ln(2πn)+o(1)\ln n! = n\ln n - n + \frac12\ln(2\pi n) + o(1):

n!n=exp(lnn1+ln(2πn)2n+o(lnnn))=ne(1+ln(2πn)2n+o(lnnn)).\sqrt[n]{n!} = \exp\Bigl(\ln n - 1 + \frac{\ln(2\pi n)}{2n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr) = \frac{n}{\eu}\Bigl(1 + \frac{\ln(2\pi n)}{2n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

Ejercicio 6.5 ★★

Para n2n \geq 2, demuestre que xn+x=1x^n + x = 1 tiene una solución única xn(0,1)x_n \in \intoo{0}{1}, que xn1x_n \to 1, y establezca

xn=1lnnn+o(lnnn).x_n = 1 - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr).

(From xnn=1xnx_n^n = 1 - x_n: take logarithms and bootstrap with xn=1εnx_n = 1 - \varepsilon_n.)

Solución

Solución de Ejercicio 6.5.

g(x)=xn+x1g(x) = x^n + x - 1 aumenta estrictamente en [0,1]\intcc{0}{1} desde 1-1 a 11: raíz única xnx_n. Desde xnn=1xn(0,1)x_n^n = 1 - x_n \in \intoo{0}{1}: si xnc<1x_n \leq c < 1a lo largo de una subsecuencia, entonces xnncn0x_n^n \leq c^n \to 0, entonces 1xn01 - x_n \to 0: contradicción con xncx_n \leq c. Por lo tanto xn1x_n \to 1.

Escriba xn=1εnx_n = 1 - \varepsilon_n, εn0+\varepsilon_n \to 0^+. el la ecuación dice (1εn)n=εn(1 - \varepsilon_n)^n = \varepsilon_n, es decir

nln(1εn)=lnεnnεn(1+o(1))=lnεn.n\ln(1 - \varepsilon_n) = \ln \varepsilon_n \quad\Longrightarrow\quad -n\varepsilon_n\bigl(1 + o(1)\bigr) = \ln\varepsilon_n .

Entonces nεn=lnεn(1+o(1))+n\varepsilon_n = -\ln\varepsilon_n\,(1 + o(1)) \to +\infty, y tomando logaritmos nuevamente: lnn+lnεn=ln(lnεn)+o(1)\ln n + \ln\varepsilon_n = \ln(-\ln\varepsilon_n) + o(1). Desde ln(lnεn)=o(ln(1/εn))\ln(-\ln \varepsilon_n) = o(\ln(1/\varepsilon_n)), esto da lnεnlnn\ln\varepsilon_n \sim -\ln n, de donde εn=lnεnn(1+o(1))lnnn\varepsilon_n = \frac{-\ln\varepsilon_n}{n}(1 + o(1)) \sim \frac{\ln n}{n}:

xn=1lnnn+o(lnnn).x_n = 1 - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr) .

Ejercicio 6.6 ★★

Empuje Ejemplo 6.23 un pedido más:

xn=nπ+π21nπ+12n2π+o(1n2).x_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + \frac{1}{2n^2\pi} + o\Bigl(\frac{1}{n^2}\Bigr).

(Write cotyn=xn\cot y_n = x_n exactly, expand coty=1yy3+o(y)\cot y = \frac1y - \frac y3 + o(y)and xn=nπ(1+12n)x_n = n\pi(1 + \frac{1}{2n} - \dots), and identify.)

Solución

Solución de Ejercicio 6.6.

Relación exacta: cotyn=xn=nπ+π2yn\cot y_n = x_n = n\pi + \frac\pi2 - y_n, con yn1nπy_n \sim \frac{1}{n\pi}(Ejemplo 6.23). Expandir coty=1yy3+O(y3)\cot y = \frac1y - \frac y3 + O(y^3):

1ynyn3+O(yn3)=nπ+π2yn1yn=nπ+π2+O(1n),\frac{1}{y_n} - \frac{y_n}{3} + O(y_n^3) = n\pi + \frac\pi2 - y_n \quad\Longrightarrow\quad \frac{1}{y_n} = n\pi + \frac\pi2 + O\Bigl(\frac1n\Bigr),

(los términos yn-y_n y yn3-\frac{y_n}{3} son O(1n)O(\frac1n)). Invertir:

yn=1nπ11+12n+O(n2)=1nπ(112n+O(1n2))=1nπ12n2π+O(1n3).y_n = \frac{1}{n\pi}\cdot\frac{1}{1 + \frac{1}{2n} + O(n^{-2})} = \frac{1}{n\pi}\Bigl(1 - \frac{1}{2n} + O\Bigl(\frac{1}{n^2}\Bigr)\Bigr) = \frac{1}{n\pi} - \frac{1}{2n^2\pi} + O\Bigl(\frac{1}{n^3}\Bigr).

Por lo tanto

xn=nπ+π2yn=nπ+π21nπ+12n2π+o(1n2).x_n = n\pi + \frac{\pi}{2} - y_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + \frac{1}{2n^2\pi} + o\Bigl(\frac{1}{n^2}\Bigr).

Ejercicio 6.7 ★★

Determinar limn1n!k=0nk!\lim_{n\to\infty} \dfrac{1}{n!}\sum_{k=0}^{n} k! (limita la suma de todos los términos excepto los dos últimos), y deducir el expansión asintótica knk!=n!(1+1n+O(n2))\sum_{k \leq n} k! = n!\bigl(1 + \frac1n + O(n^{-2})\bigr).

Solución

Solución de Ejercicio 6.7.

Separe los dos términos más grandes:

k=0nk!=n!+(n1)!+kn2k!,kn2k!(n1)(n2)!=(n1)!.\sum_{k=0}^{n} k! = n! + (n-1)! + \sum_{k \leq n-2} k! , \qquad \sum_{k\leq n-2} k! \leq (n-1)\,(n-2)! = (n-1)! .

Entonces 11n!k!1+2n1 \leq \frac{1}{n!}\sum k! \leq 1 + \frac{2}{n}: el límite es 11. Refinación: (n1)!n!=1n\frac{(n-1)!}{n!} = \frac1n y el crudo consolidado kn2k!(n1)!\sum_{k \leq n-2}k! \leq (n-1)! se puede afilar de la misma manera: kn2k!=(n2)!(1+O(1n))=O(n!n2)\sum_{k\leq n-2} k! = (n-2)!\,(1 + O(\frac1n)) = O\bigl(\frac{n!}{n^2}\bigr). Por lo tanto

k=0nk!=n!(1+1n+O(1n2)).\sum_{k=0}^{n} k! = n!\Bigl(1 + \frac1n + O\Bigl(\frac{1}{n^2}\Bigr)\Bigr).

Ejercicio 6.8 ★★★

Sean u0>0u_0 > 0 y un+1=un+1unu_{n+1} = u_n + \dfrac{1}{u_n}. Demuestre que unu_n \to \infty, luego que un2nu_n \sim \sqrt{2n}(study un2u_n^2: its increments are 2+un22 + u_n^{-2}; sum) y refine:

un=2n(1+lnn8n+o(lnnn)).u_n = \sqrt{2n}\Bigl(1 + \frac{\ln n}{8n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

(From un2=2n+k<nuk2+u02u_n^2 = 2n + \sum_{k<n} u_k^{-2} + u_0^2 and uk22ku_k^2 \sim 2k: the sum is 12lnn\sim \frac12\ln n by Teorema 6.6.)

Solución

Solución de Ejercicio 6.8.

(un)(u_n) aumenta; si está acotado, convergería a \ell con =+1\ell = \ell + \frac1\ell: absurdo. Entonces unu_n \to \infty.

Cuadrados: un+12=un2+2+un2u_{n+1}^2 = u_n^2 + 2 + u_n^{-2}, entonces

un2=u02+2n+k=0n11uk2.u_n^2 = u_0^2 + 2n + \sum_{k=0}^{n-1} \frac{1}{u_k^2} .

La suma es o(n)o(n) (los términos tienden a 00, Cesàro), por lo que un22nu_n^2 \sim 2n y un2nu_n \sim \sqrt{2n}.

Refinamiento: 1uk212k\frac{1}{u_k^2} \sim \frac{1}{2k}, así que en comparación (Teorema 6.6, o equivalentes de parcial sumas de series positivas) k<nuk212lnn\sum_{k<n} u_k^{-2} \sim \frac12 \ln n. Por lo tanto

un2=2n+lnn2(1+o(1))+O(1)un=2n1+lnn4n+o(lnnn)=2n(1+lnn8n+o(lnnn)).u_n^2 = 2n + \frac{\ln n}{2}\,(1 + o(1)) + O(1) \quad\Longrightarrow\quad u_n = \sqrt{2n}\sqrt{1 + \frac{\ln n}{4n} + o\Bigl(\frac{\ln n}{n}\Bigr)} = \sqrt{2n}\Bigl(1 + \frac{\ln n}{8n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

Ejercicio 6.9 ★★★

(Una suma de Riemann con un giro) Determinar el comportamiento asintótico de

Sn=k=1n1n+klnn.S_n = \sum_{k=1}^{n} \frac{1}{n + k\ln n} .

(Factor nn: Sn=1nk(1+klnnn)1S_n = \frac1n\sum_k \bigl(1 + \frac{k\ln n}{n}\bigr)^{-1}; recognize a Riemann-type sum with a slowly varying parameter t=lnnt = \ln n, compute 01 ⁣du1+tu=ln(1+t)t\int_0^1 \frac{\dd u}{1 + tu} = \frac{\ln(1+t)}{t}, and conclude SnlnlnnlnnS_n \sim \frac{\ln\ln n}{\ln n}.)

Solución

Solución de Ejercicio 6.9.

Factorice nn y configure t=lnnt = \ln n:

Sn=1nk=1n11+tkn.S_n = \frac1n \sum_{k=1}^{n} \frac{1}{1 + t\,\frac kn} .

Para tt fijo, la suma es una suma de Riemann de u11+tuu \mapsto \frac{1}{1 + tu} sobre [0,1]\intcc{0}{1}; la función es monótona en uu, por lo que la La suma de Riemann está entre paréntesis por la integral desplazada una malla:

01 ⁣du1+tu1nSn01 ⁣du1+tu+1n\int_0^1 \frac{\dd u}{1 + tu} - \frac1n \leq S_n \leq \int_0^1 \frac{\dd u}{1 + tu} + \frac1n

(comparación de las sumas de Riemann de una función monótona con su integral, válido para cada nn con su propio t=lnnt = \ln n). ahora 01 ⁣du1+tu=ln(1+t)t\int_0^1 \frac{\dd u}{1 + tu} = \frac{\ln(1 + t)}{t}, y 1n=o(lntt)\frac1n = o\bigl(\frac{\ln t}{t}\bigr): por lo tanto

Sn=ln(1+lnn)lnn+O(1n)    lnlnnlnn.S_n = \frac{\ln(1 + \ln n)}{\ln n} + O\Bigl(\frac 1n\Bigr) \;\sim\; \frac{\ln\ln n}{\ln n} .

Ejercicio 6.10

Demuestre la identidad (lnn)lnn=nlnlnn(\ln n)^{\ln n} = n^{\ln\ln n} y luego clasifique el siguiendo en orden creciente o()o(\cdot) al infinito, con pruebas: n2n^2, (lnn)lnn(\ln n)^{\ln n}, 2n2^n, n!n!, nnn^n.

Solución

Solución de Ejercicio 6.10.

Identidad: (lnn)lnn=elnnlnlnn=(elnn)lnlnn=nlnlnn(\ln n)^{\ln n} = \eu^{\ln n\,\ln\ln n} = \bigl(\eu^{\ln n}\bigr)^{\ln\ln n} = n^{\ln\ln n}. Clasificación: comparar logaritmos. ln(n2)=2lnn\ln(n^2) = 2\ln n; ln((lnn)lnn)=lnnlnlnn\ln\bigl((\ln n)^{\ln n}\bigr) = \ln n\ln\ln n; ln(2n)=nln2\ln(2^n) = n\ln2; ln(n!)=nlnnn+O(lnn)\ln(n!) = n\ln n - n + O(\ln n) (Stirling, o el más crudo entre corchetes lnn!nlnn\ln n! \sim n\ln n); ln(nn)=nlnn\ln(n^n) = n\ln n. desde 2lnn=o(lnnlnlnn)2\ln n = o(\ln n\ln\ln n), lnnlnlnn=o(n)\ln n\ln\ln n = o(n), nln2=o(nlnnn)n\ln 2 = o(n\ln n - n)y nlnnnnlnnn \ln n - n \sim n\ln npero n!/nn0n! / n^n \to 0(la diferencia de registros esn+O(lnn)-n + O(\ln n) \to -\infty):

n2=o((lnn)lnn),(lnn)lnn=o(2n),2n=o(n!),n!=o(nn).n^2 = o\bigl((\ln n)^{\ln n}\bigr),\quad (\ln n)^{\ln n} = o(2^n),\quad 2^n = o(n!),\quad n! = o(n^n).

(Para cada paso: la diferencia de logaritmos tiende a ++\infty, por lo que la relación tiende a 00.)

Ejercicio 6.11 ★★

(Cola de 1/k2\sum 1/k^2, dos términos) Usando el telescopio exacto k>n1k(k+1)=1n+1\sum_{k > n} \frac{1}{k(k+1)} = \frac{1}{n+1} y el descomposición 1k2=1k(k+1)+1k2(k+1)\frac1{k^2} = \frac{1}{k(k+1)} + \frac{1}{k^2(k+1)}, probar

k>n1k2=1n12n2+O(1n3).\sum_{k > n} \frac{1}{k^2} = \frac1n - \frac{1}{2n^2} + O\Bigl(\frac{1}{n^3}\Bigr).
Solución

Solución de Ejercicio 6.11.

Descomponer 1k2=1k(k+1)+1k2(k+1)\frac1{k^2} = \frac1{k(k+1)} + \frac1{k^2(k+1)} y suma para k>nk > n:

k>n1k2=1n+1+k>n1k2(k+1),\sum_{k>n}\frac1{k^2} = \frac1{n+1} + \sum_{k>n}\frac{1}{k^2(k+1)} ,

la primera suma se telescopa exactamente (1k(k+1)=1k1k+1\frac1{k(k+1)} = \frac1k - \frac1{k+1}). Por el segundo: 1k2(k+1)=1k3+O(1k4)\frac{1}{k^2(k+1)} = \frac1{k^3} + O\bigl(\frac1{k^4}\bigr)(desde 1k2(k+1)1k3=1k3(k+1)\frac{1}{k^2(k+1)} - \frac1{k^3} = \frac{-1}{k^3(k+1)}), y por el comparación integral k>n1k3=12n2+O(1n3)\sum_{k>n}\frac1{k^3} = \frac1{2n^2} + O\bigl(\frac1{n^3}\bigr), k>n1k4=O(1n3)\sum_{k>n}\frac1{k^4} = O\bigl(\frac1{n^3}\bigr). Por lo tanto

k>n1k2=1n+1+12n2+O(1n3)=1n1n2+12n2+O(1n3)=1n12n2+O(1n3),\sum_{k>n}\frac1{k^2} = \frac1{n+1} + \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr) = \frac1n - \frac1{n^2} + \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr) = \frac1n - \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr),

utilizando 1n+1=1n1n2+O(1n3)\frac1{n+1} = \frac1n - \frac1{n^2} + O\bigl(\frac1{n^3}\bigr).

Ejercicio 6.12 ★★★

Sean u0=12u_0 = \frac12 y un+1=un+eunu_{n+1} = u_n + \eu^{-u_n}. demostrar que unu_n \to \infty, luego — configurando vn=eunv_n = \eu^{u_n} y mostrando vn+1=vn+1+12vn+O(vn2)v_{n+1} = v_n + 1 + \frac{1}{2v_n} + O\bigl(v_n^{-2}\bigr) — establecer

un=lnn+lnn2n+O(1n).u_n = \ln n + \frac{\ln n}{2n} + O\Bigl(\frac1n\Bigr).
Solución

Solución de Ejercicio 6.12.

(un)(u_n) aumenta; si estuviera acotado convergería a un finito \ell con =+e\ell = \ell + \eu^{-\ell}: imposible. entonces unu_n \to \infty. Sea vn=eunv_n = \eu^{u_n} \to \infty: entonces

vn+1=eun+eun=vne1/vn=vn(1+1vn+12vn2+O(vn3))=vn+1+12vn+O(vn2).v_{n+1} = \eu^{u_n + \eu^{-u_n}} = v_n\,\eu^{1/v_n} = v_n\Bigl(1 + \frac1{v_n} + \frac1{2v_n^2} + O\bigl(v_n^{-3}\bigr)\Bigr) = v_n + 1 + \frac{1}{2v_n} + O\bigl(v_n^{-2}\bigr).

Sumar vk+1vk=1+O(1)v_{k+1} - v_k = 1 + O(1) primero da vn=n+O(n)v_n = n + O(n), por lo tanto vncnv_n \geq cn eventualmente; resumiendo con 12vk=O(1k)\frac1{2v_k} = O(\frac1k) da vn=n+O(lnn)v_n = n + O(\ln n). uno mas pasar: 12vk=12k(1+O(lnkk))\frac{1}{2v_k} = \frac{1}{2k}\bigl(1 + O\bigl(\tfrac{\ln k}k\bigr)\bigr), entonces

vn=n+k<n12k+O(1)=n+lnn2+O(1).v_n = n + \sum_{k<n}\frac1{2k} + O(1) = n + \frac{\ln n}2 + O(1).

Finalmente un=lnvn=lnn+ln(1+lnn2n+O(1n))=lnn+lnn2n+O(1n)u_n = \ln v_n = \ln n + \ln\Bigl(1 + \frac{\ln n}{2n} + O\bigl(\tfrac1n\bigr)\Bigr) = \ln n + \frac{\ln n}{2n} + O\bigl(\tfrac1n\bigr).

6.6 Problema: Bootstrapping, desde Euler-Maclaurin hasta el primos

Una cantidad implícita o acumulada rara vez entrega su asintóticos a la vez; uno los extrae en pasadas, cada pasada introducir la estimación anterior en la relación definitoria. Este fin de semana se entrena un problema que gira en torno a nuevas ecuaciones, lo demuestra el Fórmula de Euler-Maclaurin de primer orden (el trapezoide Actualización de la serie: comparación integral, con error riguroso. barras), invierte xlnx=nx\ln x = n y cobra el método más famoso. comprobar: del teorema de los números primos admitido, la asintótica ley pnnlnnp_n \sim n\ln n del nn-ésimo primo.

Problema 6.1

Problema de fin de semana — el Euler–Maclaurin corrección y las asintóticas del nn-ésimo primo

Parte I — The bootstrap loop on a fresh equation.

  1. Demostrar la afirmación de unicidad de Definición 6.2: si f=ikciφi+o(φk)=ikciφi+o(φk)f = \sum_{i\leq k} c_i\varphi_i + o(\varphi_k) = \sum_{i \leq k} c_i'\varphi_i + o(\varphi_k) está en la misma escala, entonces ci=cic_i = c_i' para todos los ii. Luego empuja el curso mixto ejemplo un peldaño más allá:

    1xlnx=1x+lnxx2+(lnx)2x3+o((lnx)2x3)(x+),\frac{1}{x - \ln x} = \frac1x + \frac{\ln x}{x^2} + \frac{(\ln x)^2}{x^3} + o\Bigl(\frac{(\ln x)^2}{x^3}\Bigr) \qquad (x \to +\infty),

    y explique por qué no aparece ningún término cx2\frac{c}{x^2}.

  2. Demuestre que para cada n1n \geq 1 la ecuación ex+x=n\eu^x + x = ntiene exactamente una solución real xnx_n, y que xn+x_n \to +\inftycon xnlnnx_n \sim \ln n.
  3. Arranca dos veces:

    xn=lnnlnnn(lnn)22n2+o((lnn)2n2).x_n = \ln n - \frac{\ln n}{n} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr).
  4. Verifique numéricamente en n=1000n = 1000: compare x10006.90083x_{1000} \approx 6.90083 con los de uno, dos y tres términos valores de la pregunta 3, a cinco decimales.

Parte II — Euler–Maclaurin, order one.

  1. Demuestre la identidad del núcleo trapezoidal: para gg de la clase C2C^2 en [0,1]\intcc{0}{1},

    01g(t) ⁣dt=g(0)+g(1)21201t(1t)g(t) ⁣dt\int_0^1 g(t)\,\dd t = \frac{g(0) + g(1)}{2} - \frac12\int_0^1 t(1 - t)\,g''(t)\,\dd t

    (integrate 12t(1t)g\frac12 t(1-t)g'' by parts twice).

  2. Sea ff C2C^2 en [1,+)\intco{1}{+\infty} con 1f<\int_1^\infty \abs{f''} < \infty. mostrar eso

    En=k=1nf(k)1nff(1)+f(n)2E_n = \sum_{k=1}^{n} f(k) - \int_1^n f - \frac{f(1) + f(n)}{2}

    converge a una constante EE, con la cola unida EEn18nf\abs{E - E_n} \leq \frac18\int_n^\infty\abs{f''}: el Fórmula de Euler-Maclaurin al primer pedido.

  3. Aplicar esto a f(t)=1tf(t) = \frac1t: probar

    Hn=lnn+γ+12n+εn,εn18n2,H_n = \ln n + \gamma + \frac{1}{2n} + \varepsilon_n, \qquad \abs{\varepsilon_n} \leq \frac{1}{8n^2},

    fortalecimiento Ejercicio 6.3 (identificar el constante con γ\gamma comparando con Ejemplo 6.7).

  4. Extrae el siguiente coeficiente: show εn=112n2+o(1n2)\varepsilon_n = -\frac{1}{12n^2} + o\bigl(\frac1{n^2}\bigr) (the increments of EnE_n are 1201t(1t)f(n+t) ⁣dt=112f(n)+o(f(n))\frac12\int_0^1t(1-t)f''(n+t)\dd t = \frac1{12}f''(n) + o(f''(n)); sum the tail with Teorema 6.6).
  5. Aplicar la pregunta 6 a f=lnf = \ln: volver a derivar en tres líneas la convergencia de dn=lnn!(n+12)lnn+nd_n = \ln n! - (n + \frac12)\ln n + n (Paso 1 de Teorema 6.13), con el error de bonificación tasa dn=d+O(1n)d_n = d + O\bigl(\frac1n\bigr).
  6. Aplicar la pregunta 6 a f(t)=1tf(t) = \frac{1}{\sqrt t}: mostrar

    k=1n1k=2n+c+12n+O(1n3/2)\sum_{k=1}^{n}\frac1{\sqrt k} = 2\sqrt n + c + \frac{1}{2\sqrt n} + O\Bigl(\frac{1}{n^{3/2}}\Bigr)

    para alguna constante cc y evalúe todos los términos en n=104n = 10^4(la constante es c1.4604c \approx -1.4604).

Parte III — Inversion: the equation xlnx=nx\ln x = n.

  1. Demuestre que xlnx=nx\ln x = n tiene exactamente una solución xn[1,+)x_n \in \intco{1}{+\infty} para n1n \geq 1, que xnx_n \to \inftyy que lnxnlnn\ln x_n \sim \ln n.
  2. Deduce la inversión de un término xnnlnnx_n \sim \dfrac{n}{\ln n}, luego arranca una vez más:

    lnxn=lnnlnlnn+o(1),xn=nlnn(1+lnlnnlnn+o(lnlnnlnn)).\ln x_n = \ln n - \ln\ln n + o(1), \qquad x_n = \frac{n}{\ln n}\Bigl(1 + \frac{\ln\ln n}{\ln n} + o\Bigl(\frac{\ln\ln n}{\ln n}\Bigr)\Bigr).
  3. Prueba en n=106n = 10^6: la raíz verdadera es x87848x \approx 87\,848; comparar con el de un término (72382\approx 72\,382) y valores de dos términos (86140\approx 86\,140), y explique los ganancia lenta (el parámetro de expansión es lnlnnlnn\frac{\ln\ln n}{\ln n}, solo 0.19\approx 0.19 en n=106n = 10^6).
  4. Ahora admitir teorema de los números primos: el número π(x)\pi(x) de números primos x\leq x satisface a π(x)xlnx\pi(x) \sim \frac{x}{\ln x} como xx \to \infty (probado honestamente en el volumen del año 3). Escribiendo pnp_n para el nn-ésimo primo, justifique π(pn)=n\pi(p_n) = n y ejecute la inversión de preguntas 11–12 para probar

    pnnlnn.p_n \sim n \ln n .
  5. Dividendos: (a) mostrar knpkn2lnn2\sum_{k \leq n} p_k \sim \frac{n^2\ln n}{2}(compare klnk\sum k\ln k with tlnt ⁣dt\int t\ln t\,\dd t); (b) calcule el valor aproximado probabilidad de que un entero uniformemente aleatorio con 100100 dígitos es primo (ln10100230.26\ln 10^{100} \approx 230.26: aproximadamente uno en 230230).

Parte IV — The method exported: xtanx=1x\tan x = 1.

  1. Demuestre que para cada n1n \geq 1 la ecuación tanx=1x\tan x = \frac1xtiene exactamente una solución xnx_n en (nπ,nπ+π2)\intoo{n\pi}{\,n\pi + \frac\pi2}, y que zn=xnnπ0+z_n = x_n - n\pi \to 0^+.
  2. Un término: zn1nπz_n \sim \dfrac{1}{n\pi}.
  3. Demuestre que la expansión de znz_n tiene no cn2\frac{c}{n^2} término: zn=1nπ+O(1n3)z_n = \frac1{n\pi} + O\bigl(\frac{1}{n^3}\bigr).
  4. Tres términos: usando arctanu=uu33+O(u5)\arctan u = u - \frac{u^3}3 + O(u^5)y 1xn=1nπzn(nπ)2+O(n3zn2)\frac1{x_n} = \frac{1}{n\pi} - \frac{z_n}{(n\pi)^2} + O(n^{-3}\cdot z_n^2), pruebe

    xn=nπ+1nπ43π3n3+o(1n3).x_n = n\pi + \frac{1}{n\pi} - \frac{4}{3\pi^3 n^3} + o\Bigl(\frac{1}{n^3}\Bigr).
  5. Verifique en n=3n = 3: raíz verdadera x39.5293344x_3 \approx 9.5293344; compare los valores de uno y tres términos, y contraste en una frase con el tanx=x\tan x = x del curso (Ejemplo 6.23): donde se ubica cada secuencia en su ventana y por qué.

Part V — A dynamical bootstrap, rules of the game, synthesis.

  1. Dejemos u0(0,π)u_0 \in \intoo{0}{\pi} y un+1=sinunu_{n+1} = \sin u_n. Muestre un0u_n \to 0 de forma decreciente y calcule el límite de 1un+121un2\dfrac{1}{u_{n+1}^2} - \dfrac{1}{u_n^2}(expand sin2\sin^{-2} via sinu=uu36+o(u3)\sin u = u - \frac{u^3}6 + o(u^3)).
  2. Deduce, vía Cesàro significa (volumen Año 1), el clásico

    un3n.u_n \sim \sqrt{\frac{3}{n}} .
  3. (Números certificados) Utilizando el límite riguroso de la pregunta 7, muestra que evaluar lnn+γ+12n\ln n + \gamma + \frac1{2n} en n=106n = 10^6 produce H106H_{10^6} con error como máximo 1.2510131.25\cdot10^{-13} — una suma millonaria calculada para trece dígitos por tres términos.
  4. (Reglas del juego) Probar o refutar, con pruebas o contraejemplos: (a) si unvn+u_n \sim v_n \to +\infty entonces lnunlnvn\ln u_n \sim \ln v_n; (b) si unvnu_n \sim v_n entonces eunevn\eu^{u_n} \sim \eu^{v_n}; (c) si fgf \sim g en ++\infty (f,gf, g diferenciable) y luego fgf' \sim g'.
  5. (Síntesis) En una frase cada uno: el bucle bootstrap de Método 6.22 como se usa en las Partes I, III, IV; ¿Qué aporta la corrección trapezoidal? Teorema 6.6; por qué inversión de xlnxx\ln x es exactamente el puente de π(x)\pi(x) a pnp_n; y cuál de las reglas de la pregunta 24 protegía cuál paso. Nombra las dos cumbres: la fórmula Euler-Maclaurin (primer orden), y la ley asintótica del nn-ésimo prima.
Solución

Solución de Problema 6.1.

1. Restando las dos expansiones: i(cici)φi=o(φk)\sum_i (c_i - c_i')\varphi_i = o(\varphi_k). Si algún coeficiente difiere, sea i0i_0 sea el primero: dividir por φi0\varphi_{i_0} y usar φj=o(φi0)\varphi_j = o(\varphi_{i_0}) para j>i0j > i_0 da ci0ci0=o(1)c_{i_0} - c_{i_0}' = o(1): cero, contradicción. Para la ampliación: con u=lnxx0u = \frac{\ln x}x \to 0,

1xlnx=1x11u=1x(1+u+u2+O(u3))=1x+lnxx2+(lnx)2x3+o((lnx)2x3).\frac{1}{x - \ln x} = \frac1x\cdot\frac{1}{1 - u} = \frac1x\bigl(1 + u + u^2 + O(u^3)\bigr) = \frac1x + \frac{\ln x}{x^2} + \frac{(\ln x)^2}{x^3} + o\Bigl(\frac{(\ln x)^2}{x^3}\Bigr).

No aparece ningún término cx2\frac c{x^2} porque la expansión es una serie geométrica en u=lnxxu = \frac{\ln x}{x}: cada término lleva como muchos poderes de lnx\ln x a partir de 1x\frac1x más allá del primero; el El peldaño de escala 1x2\frac1{x^2} (coeficiente de (lnx)0(\ln x)^0) es simplemente ausente, con coeficiente 00.

2. f(x)=ex+xf(x) = \eu^x + x es continuo, estrictamente creciente, con límites -\infty y ++\infty: una biyección RR\R \to \R, por lo que xn=f1(n)x_n = f^{-1}(n) existe y es único, y xn+x_n \to +\infty (f1f^{-1} aumenta a ++\infty). De exn=nxn\eu^{x_n} = n - x_n: xn=ln(nxn)lnnx_n = \ln(n - x_n) \leq \ln n, entonces xn/n0x_n/n \to 0 y xn=lnn+ln(1xn/n)=lnn+o(1)lnnx_n = \ln n + \ln(1 - x_n/n) = \ln n + o(1) \sim \ln n.

3. Escribe un=xn/nu_n = x_n/n. Segundo pase: un=lnn+o(1)nu_n = \frac{\ln n + o(1)}{n}, entonces

xn=lnn+ln(1un)=lnnun+O(un2)=lnnlnnn+o(lnnn).x_n = \ln n + \ln(1 - u_n) = \ln n - u_n + O(u_n^2) = \ln n - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}n\Bigr).

Tercer pase: ahora un=lnnnlnnn2+o(lnnn2)u_n = \frac{\ln n}{n} - \frac{\ln n}{n^2} + o\bigl(\frac{\ln n}{n^2}\bigr)y ln(1un)=unun22+O(un3)\ln(1 - u_n) = -u_n - \frac{u_n^2}2 + O(u_n^3):

xn=lnnlnnn+lnnn2(lnn)22n2+o((lnn)2n2)=lnnlnnn(lnn)22n2+o((lnn)2n2),x_n = \ln n - \frac{\ln n}n + \frac{\ln n}{n^2} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr) = \ln n - \frac{\ln n}{n} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr),

el término lnnn2\frac{\ln n}{n^2} es absorbido por o((lnn)2n2)o\bigl(\frac{(\ln n)^2}{n^2}\bigr).

4. En n=1000n = 1000: ln10006.90776\ln 1000 \approx 6.90776 (error 71037\cdot10^{-3}); dos términos: 6.900856.90085 (error 21052\cdot10^{-5}); tres términos: 6.900826.90082 (error a continuación 10510^{-5}), contra x10006.90083x_{1000} \approx 6.90083. Cada pase compra aproximadamente el factor previsto lnnn\frac{\ln n}{n}.

5. Dos integraciones por partes, empezando por la derecha: con  ⁣d ⁣dt[12t(1t)]=12t\frac{\dd}{\dd t}\bigl[\tfrac12t(1-t)\bigr] = \tfrac12 - tyt(1t)t(1-t) desapareciendo en ambos extremos,

1201t(1t)g(t) ⁣dt=01(12t)g(t) ⁣dt=[(12t)g]0101g=g(0)+g(1)201g.\frac12\int_0^1 t(1-t)g''(t)\dd t = -\int_0^1\Bigl(\frac12 - t\Bigr)g'(t)\dd t = -\Bigl[\Bigl(\frac12 - t\Bigr)g\Bigr]_0^1 - \int_0^1 g = \frac{g(0) + g(1)}2 - \int_0^1 g .

Reorganizada, esta es la identidad declarada.

6. Calcule el incremento, luego aplique la pregunta 5 a g(t)=f(n+t)g(t) = f(n + t):

En+1En=f(n+1)nn+1 ⁣ff(n+1)f(n)2=f(n)+f(n+1)2nn+1 ⁣f=1201t(1t)f(n+t) ⁣dt.\begin{align*} E_{n+1} - E_n &= f(n{+}1) - \int_n^{n+1}\!f - \frac{f(n{+}1) - f(n)}2 \\ &= \frac{f(n) + f(n{+}1)}2 - \int_n^{n+1}\!f = \frac12\int_0^1 t(1-t)f''(n+t)\dd t . \end{align*}

Desde 0t(1t)140 \leq t(1-t) \leq \frac14: En+1En18nn+1f\abs{E_{n+1} - E_n} \leq \frac18\int_n^{n+1}\abs{f''}, cuya suma sobre nn converge por hipótesis: (En)(E_n) converge (incrementos absolutamente sumables) a algunos EE, con

EEnknEk+1Ek18nf.\abs{E - E_n} \leq \sum_{k\geq n}\abs{E_{k+1} - E_k} \leq \frac18\int_n^\infty\abs{f''} .

7. f(t)=1tf(t) = \frac1t: f(t)=2t3f''(t) = \frac2{t^3}, 1f=1<\int_1^\infty\abs{f''} = 1 < \infty. Pregunta 6:

Hn=lnn+1+1n2+E+(EnE)=lnn+(E+12)+12n+εn,H_n = \ln n + \frac{1 + \frac1n}{2} + E + (E_n - E) = \ln n + \Bigl(E + \frac12\Bigr) + \frac1{2n} + \varepsilon_n,

con εn=EnE18n2 ⁣dtt3=18n2\abs{\varepsilon_n} = \abs{E_n - E} \leq \frac18\int_n^\infty\frac{2\dd t}{t^3} = \frac1{8n^2}. Comparando con Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1) (Ejemplo 6.7) identifica E+12=γE + \frac12 = \gamma.

8. De la fórmula de incremento de la pregunta 6,

εn=EnE=kn1201t(1t)2 ⁣dt(k+t)3=kn(1k301t(1t) ⁣dt+O(1k4)),\varepsilon_n = E_n - E = -\sum_{k\geq n}\frac12\int_0^1 t(1-t)\,\frac{2\,\dd t}{(k+t)^3} = -\sum_{k \geq n}\Bigl(\frac1{k^3}\int_0^1t(1-t)\dd t + O\Bigl(\frac1{k^4}\Bigr)\Bigr),

usando 1(k+t)3=1k3+O(1k4)\frac{1}{(k+t)^3} = \frac1{k^3} + O\bigl(\frac1{k^4}\bigr)uniformemente para t[0,1]t \in \intcc01. con 01t(1t)=16\int_0^1 t(1-t) = \frac16 y kn1k312n2\sum_{k\geq n}\frac1{k^3} \sim \frac{1}{2n^2} (Teorema 6.6):

εn=1612n2+o(1n2)=112n2+o(1n2).\varepsilon_n = -\frac16\cdot\frac{1}{2n^2} + o\Bigl(\frac1{n^2}\Bigr) = -\frac{1}{12n^2} + o\Bigl(\frac{1}{n^2}\Bigr).

9. f=lnf = \ln: f(t)=1t2f''(t) = -\frac1{t^2}, absolutamente integrables. La pregunta 6 da

lnn!=1nlnt ⁣dt+lnn2+E+O(18n ⁣dtt2)=(n+12)lnnn+1+E+O(1n),\ln n! = \int_1^n\ln t\,\dd t + \frac{\ln n}2 + E + O\Bigl( \frac1{8}\int_n^\infty\frac{\dd t}{t^2}\Bigr) = \Bigl(n + \frac12\Bigr)\ln n - n + 1 + E + O\Bigl(\frac1n\Bigr),

entonces dn=1+E+O(1n)d_n = 1 + E + O\bigl(\frac1n\bigr): convergencia de (dn)(d_n) — Paso 1 de Teorema 6.13 — más el tasa O(1/n)O(1/n). (El valor del límite de Stirling da E=ln2π1E = \ln\sqrt{2\pi} - 1).

10. f(t)=t1/2f(t) = t^{-1/2}: f(t)=34t5/2f''(t) = \frac34 t^{-5/2}, absolutamente integrable. Pregunta 6:

k=1n1k=2n2+1+1n2+E+O(n3/2)=2n+c+12n+O(n3/2),\sum_{k=1}^n \frac1{\sqrt k} = 2\sqrt n - 2 + \frac{1 + \frac1{\sqrt n}}2 + E + O\bigl(n^{-3/2}\bigr) = 2\sqrt n + c + \frac{1}{2\sqrt n} + O\bigl(n^{-3/2}\bigr),

con c=E32c = E - \frac32. En n=104n = 10^4: 2n=2002\sqrt n = 200, c1.46035c \approx -1.46035, 12n=0.005\frac1{2\sqrt n} = 0.005: previsto 198.54465198.54465 y, de hecho, k104k1/2=198.544645\sum_{k\leq10^4}k^{-1/2} = 198.544645\dots — tres términos, siete dígitos.

11. ttlntt \mapsto t\ln t es continuo y estrictamente aumentando en [1,)\intco1\infty (derivada lnt+11\ln t + 1 \geq 1), de 00 a ++\infty: existe un xnx_n único y un xnx_n \to \infty(de lo contrario, xnlnxnx_n\ln x_n permanecería limitado). tomando logaritmos en xnlnxn=nx_n\ln x_n = n: lnxn+lnlnxn=lnn\ln x_n + \ln\ln x_n = \ln n; desde lnlnxn=o(lnxn)\ln\ln x_n = o(\ln x_n), al dividir por lnxn\ln x_n se obtiene lnnlnxn1\frac{\ln n}{\ln x_n} \to 1: lnxnlnn\ln x_n \sim \ln n.

12. De xn=nlnxnx_n = \frac{n}{\ln x_n} y lnxnlnn\ln x_n \sim \ln n: xnnlnnx_n \sim \frac{n}{\ln n}. Siguiente paso: lnlnxn=ln(lnn(1+o(1)))=lnlnn+o(1)\ln\ln x_n = \ln\bigl(\ln n\,(1 + o(1))\bigr) = \ln\ln n + o(1), entonces lnxn=lnnlnlnn+o(1)\ln x_n = \ln n - \ln\ln n + o(1) y

xn=nlnnlnlnn+o(1)=nlnn11lnlnn+o(1)lnn=nlnn(1+lnlnnlnn+o(lnlnnlnn)).x_n = \frac{n}{\ln n - \ln\ln n + o(1)} = \frac{n}{\ln n}\cdot\frac{1}{1 - \frac{\ln\ln n + o(1)}{\ln n}} = \frac{n}{\ln n}\Bigl(1 + \frac{\ln\ln n}{\ln n} + o\Bigl(\frac{\ln\ln n}{\ln n}\Bigr)\Bigr).

13. En n=106n = 10^6: nlnn72382\frac{n}{\ln n} \approx 72\,382 (desactivado por 18%18\%), dos términos dan 86140\approx 86\,140 (desactivado por 1.9%1.9\%), contra el verdadero x87848x \approx 87\,848. La ganancia por pasa solo el factor lnlnnlnn2.6313.80.19\frac{\ln\ln n}{\ln n} \approx \frac{2.63}{13.8} \approx 0.19: las escalas logarítmicas convergen con una lentitud enloquecedora: una realidad dondequiera que se encuentren los números primos. involucrados.

14. Hay exactamente nn números primos pn\leq p_n (es decir, p1,,pnp_1, \dots, p_n): π(pn)=n\pi(p_n) = n. El teorema de los números primos (admitido; volumen del año 3) da n=π(pn)pnlnpnn = \pi(p_n) \sim \frac{p_n}{\ln p_n}, es decir pnnlnpnp_n \sim n\ln p_n: este es el ecuación xlnxnx\ln x \approx n leída al revés. Tomando logaritmos: lnpn=lnn+lnlnpn+o(1)\ln p_n = \ln n + \ln\ln p_n + o(1) y lnlnpn=o(lnpn)\ln\ln p_n = o(\ln p_n)fuerzan a lnpnlnn\ln p_n \sim \ln n como en la pregunta 11. Sustituyendo de nuevo:

pnnlnpn=nlnnlnpnlnnnlnn.p_n \sim n\ln p_n = n\,\ln n\,\frac{\ln p_n}{\ln n} \sim n\ln n .

15. (a) Reparar ε>0\varepsilon > 0; para kk grande, (1ε)klnkpk(1+ε)klnk(1 - \varepsilon)k\ln k \leq p_k \leq (1 + \varepsilon)k\ln k. Por comparación con el creciente tlntt\ln t (horquillado tipo Teorema 6.6), knklnk=1ntlnt ⁣dt+O(nlnn)=n2lnn2n24+O(nlnn)n2lnn2\sum_{k\leq n}k\ln k = \int_1^n t\ln t\,\dd t + O(n\ln n) = \frac{n^2\ln n}2 - \frac{n^2}4 + O(n\ln n) \sim \frac{n^2\ln n}2. Por lo tanto knpk=n2lnn2(1+O(ε)+o(1))\sum_{k\leq n}p_k = \frac{n^2\ln n}{2}(1 + O(\varepsilon) + o(1))para cada ε\varepsilon: knpkn2lnn2\sum_{k\leq n}p_k \sim \frac{n^2\ln n}2. (b) Por la prima teorema de los números, entre los números enteros hasta 1010010^{100} a proporción 1ln10100=1230.26\sim \frac{1}{\ln 10^{100}} = \frac1{230.26\dots} son primos: un entero uniformemente aleatorio de dígitos 100100 es primo con probabilidad aproximadamente 1230\frac1{230}.

16. En (nπ,nπ+π2)\intoo{n\pi}{n\pi + \frac\pi2}, g(x)=tanx1xg(x) = \tan x - \frac1xes continuo y aumenta estrictamente (g=1+tan2x+1x2>0g' = 1 + \tan^2x + \frac1{x^2} > 0), con g1nπ<0g \to -\frac1{n\pi} < 0 en el extremo izquierdo y g+g \to +\infty en el derecho: exactamente uno raíz xnx_n. Desde tanzn=tanxn=1xn0\tan z_n = \tan x_n = \frac1{x_n} \to 0 con zn(0,π2)z_n \in \intoo{0}{\frac\pi2}: zn0+z_n \to 0^+.

17. tanznzn\tan z_n \sim z_n y 1xn1nπ\frac1{x_n} \sim \frac1{n\pi}: zn1nπz_n \sim \frac1{n\pi}.

18. zn=arctan1xnz_n = \arctan\frac1{x_n} y arctanu=u+O(u3)\arctan u = u + O(u^3). Con zn=O(1n)z_n = O(\frac1n):

1xn=1nπ11+znnπ=1nπznn2π2+O(1n4)=1nπ+O(1n3),\frac1{x_n} = \frac{1}{n\pi}\cdot\frac1{1 + \frac{z_n}{n\pi}} = \frac1{n\pi} - \frac{z_n}{n^2\pi^2} + O\Bigl(\frac1{n^4}\Bigr) = \frac1{n\pi} + O\Bigl(\frac1{n^3}\Bigr),

entonces zn=1nπ+O(1n3)z_n = \frac1{n\pi} + O\bigl(\frac1{n^3}\bigr): el El renglón cn2\frac{c}{n^2} lleva el coeficiente 00, porque el La primera corrección a 1xn\frac1{x_n} es en sí misma de tamaño. znn2=O(n3)\frac{z_n}{n^2} = O(n^{-3}).

19. Inserte zn=1nπ+O(n3)z_n = \frac1{n\pi} + O(n^{-3}) en el visualización anterior:

1xn=1nπ1n3π3+O(1n5),\frac{1}{x_n} = \frac{1}{n\pi} - \frac{1}{n^3\pi^3} + O\Bigl(\frac1{n^5}\Bigr),

entonces zn=arctan1xn=1xn13(1xn)3+O(1n5)=1nπ1n3π313n3π3+O(1n5)z_n = \arctan\frac1{x_n} = \frac1{x_n} - \frac{1}{3}\Bigl(\frac1{x_n}\Bigr)^3 + O\Bigl(\frac1{n^5}\Bigr) = \frac1{n\pi} - \frac{1}{n^3\pi^3} - \frac{1}{3n^3\pi^3} + O\Bigl(\frac1{n^5}\Bigr):

xn=nπ+1nπ43π3n3+O(1n5).x_n = n\pi + \frac{1}{n\pi} - \frac{4}{3\pi^3n^3} + O\Bigl(\frac1{n^5}\Bigr).

20. En n=3n = 3: un término 9.530889.53088, tres términos 9.529299.52929, raíz verdadera 9.529339.52933: errores 1.51031.5\cdot10^{-3} y 51055\cdot10^{-5}. Contraste: para tanx=x\tan x = x la raíz debe hacer tan\tan enorme, por lo que abraza el extremo bien nπ+π2n\pi + \frac\pi2 de la ventana, a distancia 1nπ\sim\frac1{n\pi} antes de la asíntota; para xtanx=1x\tan x = 1 la raíz debe hacer tan\tan pequeño, por lo que se encuentra justo después del extremo izquierda nπn\pi, a la distancia 1nπ\sim\frac1{n\pi} después del cero. Mismo método, Geografía espejo.

21. sinu<u\sin u < u en los mapas (0,π)\intoo0\pi y sin\sin (0,π)\intoo0\pi en (0,1](0,π)\intoc01 \subseteq \intoo0\pi: después del uno paso u1(0,1]u_1 \in \intoc{0}{1}, luego (un)(u_n) disminuye y es acotado abajo por 00: converge a un punto fijo de sin\sin, es decir, a 00. Expansión: sinu=u(1u26+o(u2))\sin u = u(1 - \frac{u^2}6 + o(u^2)), entonces

1un+121un2=1un2((1un26+o(un2))21)=1un2(un23+o(un2))13.\frac{1}{u_{n+1}^2} - \frac1{u_n^2} = \frac{1}{u_n^2}\Bigl(\bigl(1 - \tfrac{u_n^2}6 + o(u_n^2)\bigr)^{-2} - 1\Bigr) = \frac{1}{u_n^2}\Bigl(\frac{u_n^2}{3} + o(u_n^2)\Bigr) \longrightarrow \frac13 .

22. Por Cesàro (volumen Año 1), la media de los incrementos converge al mismo límite:

1n1un2=1n(1u02+k=0n1(1uk+121uk2))13,\frac{1}{n}\cdot\frac{1}{u_n^2} = \frac1n\Bigl(\frac1{u_0^2} + \sum_{k=0}^{n-1} \Bigl(\frac1{u_{k+1}^2} - \frac1{u_k^2}\Bigr)\Bigr) \longrightarrow \frac13 ,

entonces un23nu_n^2 \sim \frac3n y, siendo todos los términos positivos, un3/nu_n \sim \sqrt{3/n}.

23. Por la pregunta 7, Hnlnnγ12n18n2\abs{H_n - \ln n - \gamma - \frac1{2n}} \leq \frac1{8n^2}. En n=106n = 10^6 este límite es 181012=1.251013\frac{1}{8\cdot10^{12}} = 1.25\cdot10^{-13}: tres computados términos entregan la suma armónica de un millón de términos a trece dígitos, con un certificado de error totalmente riguroso — el conjunto punto de una fórmula asintótica con resto explícito.

24. (a) Verdadero: lnunlnvn=lnunvn0\ln u_n - \ln v_n = \ln\frac{u_n}{v_n} \to 0 mientras que lnvn+\ln v_n \to +\infty, por lo que la relación de logaritmos tiende a 11. (b) Falso: un=n+1vn=nu_n = n + 1 \sim v_n = n, pero eun/evn=e1\eu^{u_n}/\eu^{v_n} = \eu \neq 1. La equivalencia tolera errores aditivos o(1)o(1) en el exponente, no O(1)O(1). (c) Falso: f(x)=x+sin(x2)g(x)=xf(x) = x + \sin(x^2) \sim g(x) = x en ++\infty, pero f(x)=1+2xcos(x2)f'(x) = 1 + 2x\cos(x^2)oscila ilimitadamente mientras que g=1g' = 1: Las derivadas de funciones equivalentes no necesitan ser comparables en todos.

25. Se ejecutó el bucle de Método 6.22. idénticamente tres veces: localizar la raíz, extraer un crudo plazo, retroalimentarlo para el siguiente pedido — en ex+x=n\eu^x + x = n (Parte I), en xlnx=nx\ln x = n (Parte III), en xtanx=1x\tan x = 1 (Parte IV). La corrección trapezoidal actualiza la serie: integral comparación de "la diferencia converge" a una definición explícita Término f(1)+f(n)2\frac{f(1) + f(n)}2 con un resto certificado O(nf)O(\int_n^\infty \abs{f''}) — constantes y barras de error en lugar de mera convergencia. El puente hacia los números primos es pura inversión: el El teorema de los números primos dice π(x)lnxx\pi(x)\ln x \approx x, por lo que pnp_n, definido por π(pn)=n\pi(p_n) = n, resuelve una ecuación xlnx=nx\ln x = n — y hereda sus asintóticas. Regla a) de la pregunta 24 legitimó cada pasaje de unvnu_n \sim v_n a lnunlnvn\ln u_n \sim \ln v_n (preguntas 11, 14); la falsedad de (b) es la razón por la que nunca exponenciar equivalencias. Cumbres: Euler-Maclaurin fórmula de primer orden (pregunta 6), y la ley asintótica pnnlnnp_n \sim n\ln n del nn-ésimo primo (pregunta 14).