Matemáticas universitarias — Grado 2 · Bachelor Year 2
6Comparación de funciones
El análisis asintótico —el arte de sustituir una cantidad complicada por otra sencilla más un error controlado— se inició en el volumen del primer año con los desarrollos de Taylor. Este capítulo lo convierte en disciplina propia: desarrollos sobre escalas generales, la comparación serie–integral con toda su fuerza asintótica, la fórmula de Stirling (demostrada por completo) y el estudio sistemático de las sucesiones definidas implícitamente. Estas técnicas son el pan de cada día del análisis asintótico, y todo capítulo posterior que estime algo —series, integrales, probabilidades— come de esta mesa.
6.1 Relaciones de comparación y escalas
Definición 6.1
Cerca de un punto a (a∈R o ±∞), para funciones (o para sucesiones, con n→∞): f=o(g), f=O(g), f∼g como en el volumen del primer año. Una escala de comparación en a es una familia de funciones positivas, comparables dos a dos y totalmente ordenada por o(⋅); la escala estándar en +∞ es
xα(lnx)β(α,β∈R),
ordenada lexicográficamente en (α,β) y refinada cuando hace falta con exponenciales eγx.
Definición 6.2(Desarrollo asintótico)
f admite el desarrollo asintótico
f=c1φ1+c2φ2+⋯+ckφk+o(φk)(φi+1=o(φi) en la escala)
cuando los restos sucesivos cumplen las estimaciones indicadas. Entonces los coeficientes son únicos: c1=limf/φ1 y, por inducción, ci+1=lim(f−∑j≤icjφj)/φi+1.
Ejemplo 6.3
Los desarrollos de Taylor son desarrollos asintóticos sobre la escala (x−a)k en a. Pero la noción es estrictamente más amplia: en +∞,
x−lnx1=x1⋅1−xlnx1=x1+x2lnx+o(x2lnx),
un desarrollo sobre la escala mixta: no se aplica ningún teorema de Taylor, solo el desarrollo geométrico y el cálculo con o.
Ejemplo 6.4(La escala estándar está realmente ordenada)
La afirmación lexicográfica de la Definición 6.1 necesita una línea de demostración por caso. Comparemos xα(lnx)β y xα′(lnx)β′ en +∞. Si α<α′: el cociente es xα−α′(lnx)β−β′→0, porque una potencia negativa de x aplasta cualquier potencia de lnx (póngase x=et: e(α−α′)ttβ−β′→0 por el límite exponencial-gana-a-polinomio del volumen del primer año). Si α=α′ y β<β′: el cociente es (lnx)β−β′→0 directamente. Así pues, los pares (α,β), ordenados lexicográficamente, ordenan la escala por o(⋅); y la sustitución x=et es el truco universal para las comparaciones mixtas de potencias y logaritmos.
Ejemplo 6.5(Ordenar una fauna variada)
Las escalas han de estar ordenadas; he aquí el ejercicio estándar. En +∞, comparemos n10, elnn⋅n, 2n y nlnn tomando logaritmos:
10lnn≪(lnn)2≪nlnn≪nln2,
donde an≪bn significa an=o(bn); la segunda entrada es ln(nlnn). Las exponenciales conservan esos saltos estrictos (si lnun−lnvn→−∞, entonces un/vn→0), de modo que
n10=o(nlnn),nlnn=o(enlnn),enlnn=o(2n).
La moraleja, por partida doble: compárese siempre a través de logaritmos (diferencias de logaritmos, no cocientes de logaritmos), y no se concluya nunca un∼vn a partir de lnun∼lnvn; el par n10 y nlnn tiene cociente de logaritmos que tiende a ∞, pero 2n y 4n tienen cociente de logaritmos exactamente 2 y son tremendamente no equivalentes.
6.2 Comparación serie–integral, en clave asintótica
Teorema 6.6
Sea fcontinua, positiva y decreciente sobre [1,+∞).
Si ∫1∞f converge, los restos cumplen
∫n+1∞f≤k>n∑f(k)≤∫n∞f.
Si ∫1∞f diverge, las sumas parciales cumplen ∑k=1nf(k)=∫1nf+C+o(1) para cierta constante C: la diferencia ∑k≤nf(k)−∫1nfconverge.
Demostración. El encuadre f(k+1)≤∫kk+1f≤f(k) (por ser decreciente) fue el recurso del primer año; sumando sobre k≥n+1, o bien sobre k≥n, se obtiene (1). Para (2), pongamos uk=f(k)−∫kk+1f: por el encuadre, 0≤uk≤f(k)−f(k+1), de modo que las sumas parciales de ∑uk están acotadas por la suma telescópica f(1)−f(n+1)≤f(1): la serie converge. Además, la sucesión (∫nn+1f)n es no creciente (f decrece) y no negativa, luego convergente. Escribiendo
k=1∑nf(k)−∫1nf=k=1∑nuk+∫nn+1f,
el miembro derecho converge cuando n→∞: la diferencia converge a una constante C, que es el enunciado (2). ∎
Ejemplo 6.7(El desarrollo armónico)
Para f(t)=t1: Hn=lnn+γ+o(1), recuperando la constante de Euler (volumen del primer año) con una demostración más limpia. Un orden más allá (Ejercicio 6.3):
Hn=lnn+γ+2n1+o(n1).
Los números hacen visible la ganancia en n=10: H10=2.928968… y ln10=2.302585…, de modo que la estimación bruta de γ es H10−ln10=0.626383, con un error de 0.049; restando la corrección 201 se obtiene 0.576383, que dista de γ=0.577216 solo 8.3⋅10−4, cantidad que es a su vez el término siguiente 12⋅1001 del desarrollo, como demuestra el problema de fin de semana (pregunta 8).
Ejemplo 6.8(Un ln(n!) tosco, sin Stirling)
El recurso del encuadre basta ya para localizar ln(n!). Como ln es creciente,
∫k−1klntdt≤lnk≤∫kk+1lntdt,
y sumando sobre k=2,…,n (con ∫1nln=nlnn−n+1):
nlnn−n+1≤ln(n!)≤(n+1)ln(n+1)−n.
Ambas vallas valen nlnn−n+O(lnn), de donde ln(n!)=nlnn−n+O(lnn) y, en particular, ln(n!)∼nlnn. Lo que Stirling añade son los dos peldaños siguientes —el 21lnn y la constante ln2π—, que cuestan la suma telescópica más fina del Teorema 6.13. Saber qué precisión compra cada herramienta es la mitad del oficio asintótico.
Ejemplo 6.9(Las cancelaciones exigen desarrollos)
Calculemos el límite de n2+n−n. Ambos términos son ∼n, y “∼n−n” no significa nada: los equivalentes no se pueden restar. Desarrollemos en su lugar:
el límite es 21, con la velocidad de aproximación 8n1 de regalo. El mecanismo merece nombre: una diferencia de dos cantidades grandes equivalentes vive por completo en sus términos siguientes, así que hay que desarrollar hasta el primer orden en el que ambos miembros difieren, arrastrando el resto para certificar que a ese orden no sobrevive nada más.
Para f(t)=tlnt1 sobre [2,+∞) (continua, positiva y decreciente): ∫2xf=lnlnx−lnln2→∞, luego por el Teorema 6.6 (2),
k=2∑nklnk1=lnlnn+C+o(1)
para cierta constante C. Dos lecciones. Primera: la divergencia es real pero glacial; la suma parcial supera 4 por primera vez hacia n≈ee4−C, astronómicamente grande. Segunda: la formalnlnn la ha entregado una primitiva, no una conjetura; para términos monótonos, la integral es el dispositivo canónico de suma, y la constante C —como la γ de Euler— es la memoria de los primeros términos.
6.3 La fórmula de Stirling
Lema 6.11(Integrales de Wallis, revisitadas)
Sea Wn=∫0π/2sinntdt. Entonces nWnWn−1=2π para n≥1, la sucesión (Wn) decrece y Wn∼2nπ.
Demostración. La integración por partes da nWn=(n−1)Wn−2 (n≥2), luego nWnWn−1 es constante en n e igual a 1⋅W1W0=2π. Decrecimiento: sinn+1≤sinn sobre [0,2π]. El encaje, en detalle: la monotonía da Wn+1≤Wn≤Wn−1 y, dividiendo por Wn−1>0,
n+1n=Wn−1Wn+1≤Wn−1Wn≤1,
donde la identidad de la izquierda sale de la recurrencia en el índice n+1. Ambas cotas tienden a 1: Wn∼Wn−1, de donde
nWn2∼nWnWn−1=2π⟹Wn∼2nπ.
∎
Ejemplo 6.12(Las primeras integrales de Wallis)
A partir de W0=2π, W1=1 y la recurrencia nWn=(n−1)Wn−2:
W2=4π,W3=32,W4=163π,W5=158,W6=325π.
Los índices pares llevan un π y los impares son racionales: los dos productos entrelazados de las formas cerradas. Numéricamente, W6≈0.4909 frente al valor asintótico π/12≈0.5116: ya en n=6 el equivalente está a menos del 5%, y la identidad del producto es exacta para todo n: 6W6W5=6⋅325π⋅158=2π. Tablas pequeñas como esta son la manera más barata de cazar un desliz algebraico antes de que infecte un argumento asintótico.
Teorema 6.13(Stirling)
n!∼2πn(en)n.
Demostración.Paso 1: n!∼Cn(n/e)n para cierta constante C>0. Pongamos
por el desarrollo de Taylor de ln(1+n1). La serie ∑(dn−dn+1) converge, pues, absolutamente (comparación con ∑n−2), de modo que (dn) converge, digamos a d; exponenciando, n!∼Cn(n/e)n con C=ed.
Paso 2: C=2π mediante Wallis. La forma cerrada W2p=4p(p!)2(2p)!⋅2π (procedente de la recurrencia, cálculo del primer año rehecho en el marco del Lema 6.11) se combina con el paso 1:
la probabilidad de que un paseo aleatorio simétrico vuelva a 0 en el instante 2n es ∼πn1, un anuncio del Capítulo 22.
Observación 6.15(Perspectivas dentro de este volumen)
Todos los capítulos cuantitativos que vienen hablan la lengua de este. El Capítulo 7 clasifica las series comparando sus términos con la escala n−α(lnn)−β, y su problema de fin de semana cartografía por completo esa frontera. El Capítulo 9 hace lo mismo con las integrales impropias, con la escala idéntica en la variable continua. El Capítulo 11 calcula radios de convergencia a partir de limsup∣an∣1/n, un ejercicio de equivalentes de raíces n-ésimas donde Stirling es la clave habitual (nn!∼en, Ejercicio 6.4). Y los capítulos de probabilidad cobran Stirling directamente: las estimaciones locales del Capítulo 22 para coeficientes binomiales son el Ejemplo 6.14 y el Ejemplo 6.21 palabra por palabra. La asintótica no es aquí un capítulo: es el acento del volumen.
Método 6.16(Lista de comprobación del arranque sucesivo)
Antes de fiarse de un desarrollo obtenido por arranques sucesivos, audítense cuatro puntos. (1) Primero la existencia: la raíz o la sucesión ha de quedar fijada (monotonía, valores intermedios) antes de todo desarrollo; los símbolos sin referente se desarrollan maravillosamente y no significan nada. (2) Un orden por pasada: cada sustitución solo es fiable hasta el orden de la estimación introducida; extraer dos términos nuevos de una sola pasada es la fuente clásica de coeficientes erróneos. (3) Los restos viajan con nosotros: arrástrese el o(⋅) por todos los pasos algebraicos y déjese la absorción (los términos menores tragados por los restos mayores) para el final, de forma explícita. (4) Auditoría numérica: evalúese en un valor honesto de n; un error de coeficiente sobrevive a una nueva deducción algebraica con sorprendente frecuencia, y casi nunca sobrevive a la aritmética.
Observación 6.17(Errores frecuentes)
(i) Los equivalentes se suman mal: de un∼n+lnn y vn∼−nno puede concluirse un+vn∼lnn; las cancelaciones exigen desarrollos con restos explícitos y nunca equivalentes a secas. (ii) No hay que exponenciar nunca una equivalencia: n+1∼n pero en+1∼en; la dirección segura es tomar logaritmos de equivalentes que tienden a +∞ (problema de fin de semana de este capítulo, pregunta 24). (iii) Un desarrollo asintótico va unido a una escala: escribir f=x1+o(x21) afirma más que f=x1+o(x1), y mezclar ambas cosas invalida el álgebra posterior. (iv) En los arranques sucesivos hay que sustituir el desarrollo actual entero, resto incluido; descartar un o(⋅) a mitad de pasada produce coeficientes plausibles pero falsos. (v) La comparación serie–integral necesita la monotonía: para términos oscilantes falla sin más (compárese ∑ksink, Capítulo 7).
Ejemplo 6.18(Stirling en números)
Para n=10: la fórmula da 20π(10/e)10≈3598696 frente a 10!=3628800, con error relativo 8.3⋅10−3, notable para un enunciado “asintótico” en n=10. El error tiene estructura —el refinamiento exacto n!=2πn(n/e)n(1+12n1+O(n−2))— y su primera corrección 1201≈8.3⋅10−3 explica casi exactamente la discrepancia observada. La maquinaria de Euler–Maclaurin del problema de fin de semana es precisamente la fuente sistemática de esos términos correctores.
Observación 6.19(Dónde se usa este capítulo)
La comparación asintótica es la gramática de todo lo cuantitativo que viene después: los criterios de convergencia y el panorama de Bertrand del Capítulo 7, los criterios de integrabilidad del Capítulo 9, los cálculos de radio de convergencia del Capítulo 11 y los teoremas límite del Capítulo 22 (donde Stirling mueve las estimaciones de de Moivre–Laplace). El volumen del tercer año industrializa la única idea que aquí demostramos a mano —extraer el término principal y acotar el resto— convirtiéndola en el método de Laplace y en la convergencia dominada.
Ejemplo 6.20(Una integral comparada consigo misma: ∫2xlntdt)
La caja de herramientas de comparación también funciona con integrales. Sea F(x)=∫2xlntdt (el integrando es continuo sobre [2,∞)). Integrando por partes:
Por tanto F(x)∼lnxx. Quien haya visto el teorema de los números primos en el problema de fin de semana de este capítulo reconocerá a F: es la integral logarítmica, el mejor estimador de π(x), y el cálculo muestra que coincide con lnxx en primer orden.
Ejemplo 6.21(Stirling sobre un binomial desequilibrado)
La misma rutina de tres factoriales que en el Ejemplo 6.14 da, para (n3n)=n!(2n)!(3n)!:
La tasa exponencial 427=2233 es e3nH(1/3) en la notación de entropía de la teoría de la información: los binomiales desequilibrados crecen estrictamente más despacio que el central 4n por cada dos pasos —aquí, (27/4)1/3≈1.89<2 por paso—. Toda asintótica de binomiales en combinatoria y en probabilidad (Capítulo 22) es este mismo cálculo con otros pesos.
6.4 Sucesiones definidas implícitamente
Método 6.22
Para hallar la asintótica de las soluciones xn de una ecuación F(x,n)=0:
Localizar: demuéstrese la existencia y la unicidad de xn en un intervalo determinado (monotonía, teorema del valor intermedio) y hállese su comportamiento tosco (límite, orden de crecimiento).
Arrancar: sustitúyase la forma tosca xn=(teˊrmino principal)(1+εn) en la ecuación y despéjese el orden siguiente de εn; repítase, refinando un orden por pasada.
Ejemplo 6.23
Para n≥1, la ecuación tanx=x tiene exactamente una solución xn en (nπ−2π,nπ+2π) (la función tanx−x crece allí de −∞ a +∞, pues su derivada es tan2x≥0). Tosco:xn=nπ+2π−yn con yn∈(0,π); como xn→∞ y tanxn=xn→+∞, xn se acerca a la asíntota por la izquierda: yn→0. Arranque:tanxn=cotyn=tanyn1∼yn1, y la ecuación cotyn=xn∼nπ da yn∼nπ1. Por tanto
xn=nπ+2π−nπ1+o(n1),
y el proceso continúa hasta cualquier orden (Ejercicio 6.6).
Ejemplo 6.24(Una segunda vuelta del método)
Resolvamos asintóticamente x+lnx=n. Localizar:x↦x+lnx crece de −∞ a +∞ sobre (0,+∞): hay una única raíz xn, y xn→∞. Tosco: de lnxn=o(xn) resulta xn∼n. Arranque: de xn=n−lnxn y lnxn=lnn+o(1) (logaritmos de equivalentes, con ambos miembros →∞):
xn=n−lnn+o(1);
una pasada más, con lnxn=ln(n−lnn+o(1))=lnn−nlnn+o(nlnn):
xn=n−lnn+nlnn+o(nlnn).
(Comprobación en n=100: la raíz es x≈95.4415; la fórmula de tres términos da 100−4.6052+0.0461=95.4409 y la de dos términos 95.3948; cada pasada gana el orden previsto.) Mismo bucle, tercer paisaje: al método del Método 6.22 le da igual qué aspecto tenga la ecuación, con tal de que cada pasada aísle la incógnita dominante.
6.5 Ejercicios
Ejercicio 6.1★
Desarrolla en +∞, dos términos más allá del principal:
x2+x+1,ln(x2+x)−2lnx,x−lnxx+sinx.
Solución
Solución de Ejercicio 6.1.
x2+x+1=x1+x1+x21=x+21+83⋅x1+o(x1) (desarrollo binomial: 21u−81u2 con u=x1+x21 da 2x1+2x21−8x21=2x1+8x23, y después se multiplica por x).
Ordenemos las contribuciones en la escala de +∞: xlnx≫x1≥xsinx≫x2(lnx)2. Los dos términos que siguen al 1 principal son, por tanto, xlnx y después el término de oscilación acotada xsinx:
x−lnxx+sinx=1+xlnx+xsinx+O(x2(lnx)2).
Ejercicio 6.2★
Da el carácter (convergencia o divergencia) y, en caso divergente, la asintótica principal de ∑k≤nkα para α>−1, α=−1 y α<−1, mediante el Teorema 6.6.
Solución
Solución de Ejercicio 6.2.
f(t)=tα (t≥1).
α>−1: divergencia y, por el Teorema 6.6 (2), ∑k≤nkα=α+1nα+1+C+o(1) si α<0 (donde f decrece); para α≥0 (f creciente), el mismo encuadre con las desigualdades invertidas da ∑k≤nkα∼α+1nα+1.
α<−1: convergencia, con resto ∑k>nkα∼−(α+1)nα+1 por el encuadre (1) (ambas cotas integrales son equivalentes a ese valor).
Ejercicio 6.3★★
Demuestra que Hn=lnn+γ+2n1+o(n1). (Estudia vn=Hn−lnn−γ: prueba que vn−vn+1=2n21+O(n−3) y suma la cola, comparando con ∑k≥n2k21∼2n1, Teorema 6.6 (1).)
Para n≥2, demuestra que xn+x=1 tiene una única solución xn∈(0,1), que xn→1, y establece
xn=1−nlnn+o(nlnn).
(A partir de xnn=1−xn: toma logaritmos y arranca con xn=1−εn.)
Solución
Solución de Ejercicio 6.5.
g(x)=xn+x−1 crece estrictamente sobre [0,1] de −1 a 1: hay una única raíz xn. Como xnn=1−xn∈(0,1): si xn≤c<1 a lo largo de una subsucesión, entonces xnn≤cn→0, luego 1−xn→0, en contradicción con xn≤c. Por tanto xn→1.
Escribamos xn=1−εn con εn→0+. La ecuación se lee (1−εn)n=εn, es decir,
nln(1−εn)=lnεn⟹−nεn(1+o(1))=lnεn.
Así pues, nεn=−lnεn(1+o(1))→+∞ y, tomando logaritmos de nuevo: lnn+lnεn=ln(−lnεn)+o(1). Como ln(−lnεn)=o(ln(1/εn)), esto da lnεn∼−lnn, de donde εn=n−lnεn(1+o(1))∼nlnn:
Luego 1≤n!1∑k!≤1+n2: el límite es 1. Afinando: n!(n−1)!=n1, y la cota tosca ∑k≤n−2k!≤(n−1)! puede refinarse del mismo modo: ∑k≤n−2k!=(n−2)!(1+O(n1))=O(n2n!). Por tanto
k=0∑nk!=n!(1+n1+O(n21)).
Ejercicio 6.8★★★
Sea u0>0 y un+1=un+un1. Demuestra que un→∞, después que un∼2n(estudia un2: sus incrementos son 2+un−2; suma), y refina:
un=2n(1+8nlnn+o(nlnn)).
(De un2=2n+∑k<nuk−2+u02 y uk2∼2k: la suma vale ∼21lnn por el Teorema 6.6.)
Solución
Solución de Ejercicio 6.8.
(un) crece; si estuviera acotada convergería a un ℓ con ℓ=ℓ+ℓ1: absurdo. Luego un→∞.
Cuadrados: un+12=un2+2+un−2, de donde
un2=u02+2n+k=0∑n−1uk21.
La suma es o(n) (los términos tienden a 0; Cesàro), luego un2∼2n y un∼2n.
Refinamiento: uk21∼2k1, así que por comparación (Teorema 6.6, o mediante equivalentes de sumas parciales de series positivas) ∑k<nuk−2∼21lnn. Por tanto
(Una suma de Riemann con truco) Determina el comportamiento asintótico de
Sn=k=1∑nn+klnn1.
(Saca factor n: Sn=n1∑k(1+nklnn)−1; reconoce una suma de tipo Riemann con un parámetro t=lnn que varía lentamente, calcula ∫011+tudu=tln(1+t) y concluye Sn∼lnnlnlnn.)
Solución
Solución de Ejercicio 6.9.
Saquemos factor n y pongamos t=lnn:
Sn=n1k=1∑n1+tnk1.
Para t fijo, la suma es una suma de Riemann de u↦1+tu1 sobre [0,1]; la función es monótona en u, de modo que la suma de Riemann queda encuadrada por la integral desplazada en un paso de malla:
∫011+tudu−n1≤Sn≤∫011+tudu+n1
(comparación de las sumas de Riemann de una función monótona con su integral, válida para cada n con su propio t=lnn). Ahora bien, ∫011+tudu=tln(1+t) y n1=o(tlnt); por tanto
Sn=lnnln(1+lnn)+O(n1)∼lnnlnlnn.
Ejercicio 6.10★
Demuestra la identidad (lnn)lnn=nlnlnn y ordena después las siguientes expresiones en orden creciente por o(⋅) en el infinito, con demostración: n2, (lnn)lnn, 2n, n!, nn.
Solución
Solución de Ejercicio 6.10.
Identidad: (lnn)lnn=elnnlnlnn=(elnn)lnlnn=nlnlnn. Orden: comparemos logaritmos. ln(n2)=2lnn; ln((lnn)lnn)=lnnlnlnn; ln(2n)=nln2; ln(n!)=nlnn−n+O(lnn) (Stirling, o el encuadre más tosco lnn!∼nlnn); ln(nn)=nlnn. Como 2lnn=o(lnnlnlnn), lnnlnlnn=o(n), nln2=o(nlnn−n), y nlnn−n∼nlnn pero n!/nn→0 (la diferencia de logaritmos es −n+O(lnn)→−∞):
n2=o((lnn)lnn),(lnn)lnn=o(2n),2n=o(n!),n!=o(nn).
(En cada paso: la diferencia de logaritmos tiende a +∞, luego el cociente tiende a 0.)
Ejercicio 6.11★★
(Cola de ∑1/k2, dos términos) Usando la suma telescópica exacta ∑k>nk(k+1)1=n+11 y la descomposición k21=k(k+1)1+k2(k+1)1, demuestra que
k>n∑k21=n1−2n21+O(n31).
Solución
Solución de Ejercicio 6.11.
Descompongamos k21=k(k+1)1+k2(k+1)1 y sumemos para k>n:
k>n∑k21=n+11+k>n∑k2(k+1)1,
donde la primera suma es telescópica exacta (k(k+1)1=k1−k+11). Para la segunda: k2(k+1)1=k31+O(k41) (pues k2(k+1)1−k31=k3(k+1)−1), y por comparación integral ∑k>nk31=2n21+O(n31) y ∑k>nk41=O(n31). Por tanto
Sumando primero vk+1−vk=1+O(1) se obtiene vn=n+O(n), de donde vn≥cn a partir de un índice; volviendo a sumar con 2vk1=O(k1) resulta vn=n+O(lnn). Una pasada más: 2vk1=2k1(1+O(klnk)), luego
vn=n+k<n∑2k1+O(1)=n+2lnn+O(1).
Por último, un=lnvn=lnn+ln(1+2nlnn+O(n1))=lnn+2nlnn+O(n1).
6.6 Problema: arranques sucesivos, de Euler–Maclaurin a los primos
Una cantidad implícita o acumulada rara vez entrega su asintótica de una vez; se extrae en pasadas, cada una realimentando la estimación previa en la relación que la define. Este problema de fin de semana entrena ese bucle sobre ecuaciones nuevas, demuestra la fórmula de Euler–Maclaurin de primer orden (la mejora trapezoidal de la comparación serie–integral, con barras de error rigurosas), invierte xlnx=n y cobra el cheque más famoso del método: a partir del teorema de los números primos, admitido, la ley asintótica pn∼nlnn del n-ésimo primo.
Problema 6.1
Problema de fin de semana — la corrección de Euler–Maclaurin y la asintótica del n-ésimo primo
Parte I — El bucle de arranque sobre una ecuación nueva.
Demuestra la afirmación de unicidad de la Definición 6.2: si f=∑i≤kciφi+o(φk)=∑i≤kci′φi+o(φk) sobre la misma escala, entonces ci=ci′ para todo i. Lleva después el ejemplo mixto del curso un peldaño más allá:
x−lnx1=x1+x2lnx+x3(lnx)2+o(x3(lnx)2)(x→+∞),
y explica por qué no aparece ningún término x2c.
Prueba que para todo n≥1 la ecuación ex+x=n tiene exactamente una solución real xn, y que xn→+∞ con xn∼lnn.
Arranca dos veces:
xn=lnn−nlnn−2n2(lnn)2+o(n2(lnn)2).
Comprueba numéricamente en n=1000: compara x1000≈6.90083 con los valores de uno, dos y tres términos de la pregunta 3, con cinco decimales.
Parte II — Euler–Maclaurin, orden uno.
Demuestra la identidad del núcleo trapezoidal: para g de clase C2 sobre [0,1],
∫01g(t)dt=2g(0)+g(1)−21∫01t(1−t)g′′(t)dt
(integra 21t(1−t)g′′ por partes dos veces).
Sea f de clase C2 sobre [1,+∞) con ∫1∞∣f′′∣<∞. Prueba que
En=k=1∑nf(k)−∫1nf−2f(1)+f(n)
converge a una constante E, con la cota de la cola ∣E−En∣≤81∫n∞∣f′′∣: la fórmula de Euler–Maclaurin de primer orden.
Extrae el coeficiente siguiente: prueba que εn=−12n21+o(n21)(los incrementos de En son 21∫01t(1−t)f′′(n+t)dt=121f′′(n)+o(f′′(n)); suma la cola con el Teorema 6.6).
Aplica la pregunta 6 a f=ln: vuelve a deducir en tres líneas la convergencia de dn=lnn!−(n+21)lnn+n (paso 1 del Teorema 6.13), con la velocidad de error de regalo dn=d+O(n1).
Aplica la pregunta 6 a f(t)=t1: prueba que
k=1∑nk1=2n+c+2n1+O(n3/21)
para cierta constante c, y evalúa todos los términos en n=104 (la constante es c≈−1.4604).
Parte III — Inversión: la ecuación xlnx=n.
Prueba que xlnx=n tiene exactamente una solución xn∈[1,+∞) para n≥1, que xn→∞ y que lnxn∼lnn.
Deduce la inversión de un término xn∼lnnn y arranca una vez más:
Comprueba en n=106: la raíz verdadera es x≈87848; compárala con los valores de un término (≈72382) y de dos términos (≈86140), y explica la lentitud de la ganancia (el parámetro del desarrollo es lnnlnlnn, que solo vale ≈0.19 en n=106).
Admitamos ahora el teorema de los números primos: el número π(x) de primos ≤x cumple π(x)∼lnxx cuando x→∞ (demostrado honestamente en el volumen del tercer año). Escribiendo pn para el n-ésimo primo, justifica π(pn)=n y aplica la inversión de las preguntas 11–12 para demostrar que
pn∼nlnn.
Dividendos: (a) prueba que ∑k≤npk∼2n2lnn(compara ∑klnk con ∫tlntdt); (b) calcula la probabilidad aproximada de que un entero elegido al azar uniformemente con 100 cifras sea primo (ln10100≈230.26: alrededor de uno de cada 230).
Parte IV — El método exportado: xtanx=1.
Prueba que para cada n≥1 la ecuación tanx=x1 tiene exactamente una solución xn en (nπ,nπ+2π), y que zn=xn−nπ→0+.
Un término: zn∼nπ1.
Prueba que el desarrollo de znno tiene término n2c: zn=nπ1+O(n31).
Tres términos: usando arctanu=u−3u3+O(u5) y xn1=nπ1−(nπ)2zn+O(n−3⋅zn2), demuestra que
xn=nπ+nπ1−3π3n34+o(n31).
Comprueba en n=3: la raíz verdadera es x3≈9.5293344; compara los valores de uno y de tres términos, y contrasta en una frase con el tanx=x del curso (Ejemplo 6.23): dónde se sitúa cada sucesión dentro de su ventana, y por qué.
Parte V — Un arranque dinámico, reglas del juego, síntesis.
Sea u0∈(0,π) y un+1=sinun. Prueba que un→0 de manera decreciente y calcula el límite de un+121−un21(desarrolla sin−2 mediante sinu=u−6u3+o(u3)).
Deduce, mediante las medias de Cesàro (volumen del primer año), el clásico
un∼n3.
(Numérica certificada) Usando la cota rigurosa de la pregunta 7, prueba que evaluar lnn+γ+2n1 en n=106 da H106 con un error de a lo sumo 1.25⋅10−13: una suma de un millón de términos calculada con trece cifras mediante tres términos.
(Reglas del juego) Demuestra o refuta, con demostraciones o contraejemplos: (a) si un∼vn→+∞, entonces lnun∼lnvn; (b) si un∼vn, entonces eun∼evn; (c) si f∼g en +∞ (f,g derivables), entonces f′∼g′.
(Síntesis) Una frase para cada punto: el bucle de arranque del Método 6.22 tal como se usa en las partes I, III y IV; qué añade la corrección trapezoidal al Teorema 6.6; por qué la inversión de xlnx es exactamente el puente de π(x) a pn; y cuál de las reglas de la pregunta 24 protegió cada paso. Nombra las dos cumbres: la fórmula de Euler–Maclaurin (de primer orden) y la ley asintótica del n-ésimo primo.
Solución
Solución de Problema 6.1.
1. Restando los dos desarrollos: ∑i(ci−ci′)φi=o(φk). Si algún coeficiente difiere, sea i0 el primero: dividiendo por φi0 y usando φj=o(φi0) para j>i0 se obtiene ci0−ci0′=o(1), es decir, cero: contradicción. Para el desarrollo: con u=xlnx→0,
No aparece ningún término x2c porque el desarrollo es una serie geométrica en u=xlnx: cada término lleva tantas potencias de lnx como de x1 más allá de la primera; el peldaño x21 de la escala (coeficiente de (lnx)0) sencillamente no está, con coeficiente 0.
2.f(x)=ex+x es continua, estrictamente creciente y con límites −∞ y +∞: es una biyección R→R, de modo que xn=f−1(n) existe y es único, y xn→+∞ (f−1 crece hacia +∞). De exn=n−xn resulta xn=ln(n−xn)≤lnn, luego xn/n→0 y xn=lnn+ln(1−xn/n)=lnn+o(1)∼lnn.
3. Escribamos un=xn/n. Segunda pasada: un=nlnn+o(1), así que
4. Para n=1000: ln1000≈6.90776 (error 7⋅10−3); dos términos: 6.90085 (error 2⋅10−5); tres términos: 6.90082 (error por debajo de 10−5), frente a x1000≈6.90083. Cada pasada compra aproximadamente el factor previsto nlnn.
5. Dos integraciones por partes, empezando por la derecha: con dtd[21t(1−t)]=21−t y anulándose t(1−t) en ambos extremos,
Como 0≤t(1−t)≤41: ∣En+1−En∣≤81∫nn+1∣f′′∣, cuya suma sobre n converge por hipótesis; luego (En) converge (incrementos absolutamente sumables) a cierto E, con
∣E−En∣≤k≥n∑∣Ek+1−Ek∣≤81∫n∞∣f′′∣.
7.f(t)=t1: f′′(t)=t32 y ∫1∞∣f′′∣=1<∞. Por la pregunta 6:
Hn=lnn+21+n1+E+(En−E)=lnn+(E+21)+2n1+εn,
con ∣εn∣=∣En−E∣≤81∫n∞t32dt=8n21. Comparando con Hn=lnn+γ+o(1) (Ejemplo 6.7) se identifica E+21=γ.
8. A partir de la fórmula del incremento de la pregunta 6,
de modo que dn=1+E+O(n1): convergencia de (dn) —paso 1 del Teorema 6.13— más la velocidad O(1/n). (El valor del límite dado por Stirling da E=ln2π−1.)
10.f(t)=t−1/2: f′′(t)=43t−5/2, absolutamente integrable. La pregunta 6 da
con c=E−23. Para n=104: 2n=200, c≈−1.46035, 2n1=0.005; predicción 198.54465, y en efecto ∑k≤104k−1/2=198.544645…: tres términos, siete cifras.
11.t↦tlnt es continua y estrictamente creciente sobre [1,∞) (derivada lnt+1≥1), de 0 a +∞: existe un único xn, y xn→∞ (en otro caso xnlnxn quedaría acotado). Tomando logaritmos en xnlnxn=n: lnxn+lnlnxn=lnn; como lnlnxn=o(lnxn), dividiendo por lnxn resulta lnxnlnn→1, es decir, lnxn∼lnn.
12. De xn=lnxnn y lnxn∼lnn: xn∼lnnn. Pasada siguiente: lnlnxn=ln(lnn(1+o(1)))=lnlnn+o(1), luego lnxn=lnn−lnlnn+o(1) y
13. Para n=106: lnnn≈72382 (error del 18%); dos términos dan ≈86140 (error del 1.9%), frente al valor verdadero x≈87848. La ganancia por pasada es solo el factor lnnlnlnn≈13.82.63≈0.19: las escalas logarítmicas convergen con desesperante lentitud, hecho de la vida allí donde intervienen los primos.
14. Hay exactamente n primos ≤pn (a saber, p1,…,pn): π(pn)=n. El teorema de los números primos (admitido; volumen del tercer año) da n=π(pn)∼lnpnpn, es decir, pn∼nlnpn: es la ecuación xlnx≈n leída al revés. Tomando logaritmos: lnpn=lnn+lnlnpn+o(1), y lnlnpn=o(lnpn) obliga a lnpn∼lnn como en la pregunta 11. Sustituyendo de vuelta:
pn∼nlnpn=nlnnlnnlnpn∼nlnn.
15. (a) Fijemos ε>0; para k grande, (1−ε)klnk≤pk≤(1+ε)klnk. Por comparación con la función creciente tlnt (encuadre del tipo del Teorema 6.6), ∑k≤nklnk=∫1ntlntdt+O(nlnn)=2n2lnn−4n2+O(nlnn)∼2n2lnn. Por tanto ∑k≤npk=2n2lnn(1+O(ε)+o(1)) para todo ε: ∑k≤npk∼2n2lnn. (b) Por el teorema de los números primos, entre los enteros hasta 10100 una proporción ∼ln101001=230.26…1 son primos: un entero de 100 cifras elegido uniformemente al azar es primo con probabilidad aproximada 2301.
16. Sobre (nπ,nπ+2π), la función g(x)=tanx−x1 es continua y estrictamente creciente (g′=1+tan2x+x21>0), con g→−nπ1<0 en el extremo izquierdo y g→+∞ en el derecho: hay exactamente una raíz xn. Como tanzn=tanxn=xn1→0 con zn∈(0,2π): zn→0+.
17.tanzn∼zn y xn1∼nπ1: zn∼nπ1.
18.zn=arctanxn1 y arctanu=u+O(u3). Con zn=O(n1):
luego zn=nπ1+O(n31): el peldaño n2c tiene coeficiente 0, porque la primera corrección de xn1 es ya de tamaño n2zn=O(n−3).
19. Insertemos zn=nπ1+O(n−3) en la fórmula anterior:
xn1=nπ1−n3π31+O(n51),
y entonces zn=arctanxn1=xn1−31(xn1)3+O(n51)=nπ1−n3π31−3n3π31+O(n51):
xn=nπ+nπ1−3π3n34+O(n51).
20. Para n=3: un término da 9.53088, tres términos 9.52929, y la raíz verdadera es 9.52933; los errores son 1.5⋅10−3 y 5⋅10−5. Contraste: para tanx=x la raíz ha de hacer enorme la tangente, así que se pega al extremo derechonπ+2π de la ventana, a distancia ∼nπ1 antes de la asíntota; para xtanx=1 la raíz ha de hacer diminuta la tangente, así que se sitúa justo después del extremo izquierdonπ, a distancia ∼nπ1 tras el cero. Mismo método, geografía especular.
21.sinu<u sobre (0,π), y sin envía (0,π) dentro de (0,1]⊆(0,π): tras un paso, u1∈(0,1], y entonces (un) decrece y está acotada inferiormente por 0, luego converge a un punto fijo de sin, es decir, a 0. Desarrollo: sinu=u(1−6u2+o(u2)), de modo que
luego un2∼n3 y, al ser todos los términos positivos, un∼3/n.
23. Por la pregunta 7, Hn−lnn−γ−2n1≤8n21. En n=106 esa cota vale 8⋅10121=1.25⋅10−13: tres términos calculados entregan la suma armónica de un millón de sumandos con trece cifras y un certificado de error plenamente riguroso, que es todo el sentido de una fórmula asintótica con resto explícito.
24. (a) Cierto: lnun−lnvn=lnvnun→0 mientras que lnvn→+∞, luego el cociente de los logaritmos tiende a 1. (b) Falso: un=n+1∼vn=n, pero eun/evn=e=1. La equivalencia tolera errores aditivos o(1) en el exponente, no O(1). (c) Falso: f(x)=x+sin(x2)∼g(x)=x en +∞, pero f′(x)=1+2xcos(x2) oscila sin acotación mientras que g′=1: las derivadas de funciones equivalentes no tienen por qué ser comparables en absoluto.
25. El bucle del Método 6.22 se ejecutó idénticamente tres veces: localizar la raíz, extraer un término tosco y realimentarlo para el orden siguiente —sobre ex+x=n (parte I), sobre xlnx=n (parte III) y sobre xtanx=1 (parte IV)—. La corrección trapezoidal mejora la comparación serie–integral desde “la diferencia converge” hasta un término explícito 2f(1)+f(n) con resto certificado O(∫n∞∣f′′∣): constantes y barras de error en lugar de mera convergencia. El puente hacia los primos es pura inversión: el teorema de los números primos dice que π(x)lnx≈x, de modo que pn, definido por π(pn)=n, resuelve una ecuación xlnx=n y hereda su asintótica. La regla (a) de la pregunta 24 legitimó todo paso de un∼vn a lnun∼lnvn (preguntas 11 y 14); y la falsedad de (b) es la razón de que nunca exponenciemos equivalencias. Cumbres: la fórmula de Euler–Maclaurin de primer orden (pregunta 6) y la ley asintótica pn∼nlnn del n-ésimo primo (pregunta 14).