Mathematics · Libro 3 · Bachelor Year 1

Matemáticas universitarias — Grado 1

Matemáticas universitarias — Grado 1 · Bachelor Year 1

16Fórmulas de Taylor y desarrollos asintóticos

Cerca de un punto, una función regular vale tanto como un polinomio — con un error controlable. Las fórmulas de Taylor lo hacen exacto en tres versiones (resto integral, resto de Lagrange, resto de Young), y los desarrollos asintóticos resultantes, manipulados algebraicamente, se convierten en la herramienta más afilada del análisis elemental: límites, equivalentes, comportamiento local, asíntotas.

16.1 Notación de comparación

Definición 16.1 (Notación de Landau)

Sean f,gf, g definidas cerca de x0x_0 (x0Rx_0 \in \R o ±\pm\infty). Se escribe, cuando xx0x \to x_0:

  • f=o(g)f = o(g) («o pequeña») cuando f=εgf = \varepsilon g con ε(x)0\varepsilon(x) \to 0;
  • f=O(g)f = O(g) («O grande») cuando f=ugf = u g con uu acotada cerca de x0x_0;
  • fgf \sim g («equivalente») cuando f=(1+ε)gf = (1 + \varepsilon) g con ε0\varepsilon \to 0 — equivalentemente, fg=o(g)f - g = o(g).

La misma notación se aplica a las sucesiones (nn \to \infty).

Proposición 16.2 (Reglas)

Cuando xx0x \to x_0:

  1. \sim es una relación de equivalencia; fgf \sim g implica que ff y gg comparten los límites, los signos (cerca de x0x_0) y la ausencia de ceros;
  2. los equivalentes se multiplican y dividen: f1g1f_1 \sim g_1, f2g2f_2 \sim g_2 implican f1f2g1g2f_1 f_2 \sim g_1 g_2 y f1f2g1g2\frac{f_1}{f_2} \sim \frac{g_1}{g_2};
  3. los equivalentes no se suman: x+1xx + 1 \sim x y xx+2-x \sim -x + 2 en ++\infty y, sin embargo, las sumas 11 y 22 no son equivalentes. Para sumar hay que volver a los desarrollos con términos o()o(\cdot) explícitos;
  4. o(g)+o(g)=o(g)o(g) + o(g) = o(g),   uo(g)=o(ug)\;u \cdot o(g) = o(ug),   o(o(g))=o(g)\;o(o(g)) = o(g), y fg    f=g+o(g)f \sim g \iff f = g + o(g).

Demostración. Cada una es una manipulación breve de las definiciones; por ejemplo, f1f2=(1+ε1)(1+ε2)g1g2f_1 f_2 = (1+\varepsilon_1)(1+\varepsilon_2) g_1 g_2 y (1+ε1)(1+ε2)1(1 + \varepsilon_1)(1+\varepsilon_2) \to 1. Dos puntos de (4) merecen su línea. uo(g)=o(ug)u\cdot o(g) = o(ug): si f=εgf = \varepsilon g con ε0\varepsilon \to 0, entonces uf=ε(ug)uf = \varepsilon\,(ug) con el mismo ε\varepsilon. o(o(g))=o(g)o(o(g)) = o(g): si f=ε1hf = \varepsilon_1 h y h=ε2gh = \varepsilon_2 g con los dos εi0\varepsilon_i \to 0, entonces f=(ε1ε2)gf = (\varepsilon_1\varepsilon_2) g, y el producto de los dos infinitésimos lo es. La equivalencia fg    f=g+o(g)f \sim g \iff f = g + o(g) es la definición leída dos veces: fg=εgf - g = \varepsilon g. El contraejemplo de (3) es la demostración de (3).

Ejemplo 16.3 (La escala de comparación)

Cuando x+x \to +\infty, la escala estándar dice, en orden creciente de fuerza:

1  =  o(lnx),lnx=o(x0.01),x0.01=o(x),x=o(x10),x10=o(ex),ex=o(e2x),1 \;=\; o(\ln x), \quad \ln x = o(x^{0.01}), \quad x^{0.01} = o(\sqrt x), \quad \sqrt x = o(x^{10}), \quad x^{10} = o(\eu^{x}), \quad \eu^x = o(\eu^{2x}) ,

siendo cada paso una instancia de las comparaciones de crecimiento de la Proposición 4.6 (las potencias ganan a los logaritmos, las exponenciales ganan a las potencias y, dentro de una misma familia, decide el exponente). Dos hábitos que conviene adquirir: primero, un O()O(\cdot) que cae en una clase menor asciende en silencio (O(lnx)O(\ln x) también es o(x0.01)o(x^{0.01})); segundo, en x0+x \to 0^+ toda la escalera se invierte con la sustitución x1xx \mapsto \frac1x —ahí lnx=o(x0.01)\ln x = o(x^{-0.01})—, de modo que «xαlnx0x^\alpha \ln x \to 0» vale para todo α>0\alpha > 0. Tener clara la escala es la mitad de todo argumento asintótico del Capítulo 17.

Ejemplo 16.4 (Unicidad de los desarrollos, y un dividendo de paridad)

Si una función admite dos desarrollos en 00 al mismo orden,

a0+a1x++anxn+o(xn)=b0+b1x++bnxn+o(xn),a_0 + a_1 x + \dots + a_n x^n + o(x^n) = b_0 + b_1 x + \dots + b_n x^n + o(x^n),

entonces ak=bka_k = b_k para todo kk: restando y poniendo ck=akbkc_k = a_k - b_k, evalúese la identidad c0+c1x++cnxn=o(xn)c_0 + c_1 x + \dots + c_n x^n = o(x^n) cuando x0x \to 0 para obtener c0=0c_0 = 0; divídase por xx y repítase — cada división es legítima porque la expresión restante vuelve a ser o(xnk)o(x^{n-k}). Los coeficientes son, por tanto, intrínsecos, y se pueden calcular por cualquier camino (derivadas de Taylor, álgebra sobre desarrollos conocidos, integración): todos los caminos han de coincidir. Dividendo: una función par solo tiene potencias pares en su desarrollo —sustitúyase xx por x-x e invóquese la unicidad—, e igualmente las impares tienen potencias impares. Por eso cos\cos lleva o(x2p+1)o(x^{2p+1}) y no o(x2p)o(x^{2p}) en la tabla de más abajo: el término impar ausente es información gratis, un orden de precisión por nada.

16.2 Las tres fórmulas de Taylor

Teorema 16.5 (Taylor con resto integral)

Sea ff de clase Cn+1C^{n+1} en un intervalo que contiene a aa y a xx. Entonces

f(x)=k=0nf(k)(a)k!(xa)k+ax(xt)nn!f(n+1)(t) ⁣dt.f(x) = \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}\,(x - a)^k + \int_a^x \frac{(x - t)^n}{n!}\, f^{(n+1)}(t)\, \dd t .

Demostración. Inducción sobre nn. Para n=0n = 0: f(x)=f(a)+axf(t) ⁣dtf(x) = f(a) + \int_a^x f'(t)\dd t es el teorema fundamental (Teorema 15.9). Paso: intégrese el resto por partes,

ax(xt)nn!f(n+1)(t) ⁣dt=[(xt)n+1(n+1)!f(n+1)(t)]ax+ax(xt)n+1(n+1)!f(n+2)(t) ⁣dt,\int_a^x \frac{(x-t)^n}{n!} f^{(n+1)}(t)\,\dd t = \Bigl[-\frac{(x-t)^{n+1}}{(n+1)!} f^{(n+1)}(t)\Bigr]_a^x + \int_a^x \frac{(x-t)^{n+1}}{(n+1)!} f^{(n+2)}(t)\,\dd t ,

contribuyendo el corchete con el término f(n+1)(a)(n+1)!(xa)n+1\frac{f^{(n+1)}(a)}{(n+1)!}(x - a)^{n+1}.

Ejemplo 16.6 (Un desarrollo exacto con su resto)

Para ln(1+x)\ln(1 + x), el resto integral se puede hacer completamente explícito sin derivar nada nn veces: intégrese la identidad geométrica finita 11+t=k=0n1(t)k+(t)n1+t\frac{1}{1+t} = \sum_{k=0}^{n-1}(-t)^k + \frac{(-t)^n}{1+t} de 00 a xx:

ln(1+x)=k=1n(1)k1xkk+(1)n0xtn1+t ⁣dt,\ln(1 + x) = \sum_{k=1}^{n} \frac{(-1)^{k-1}x^k}{k} + (-1)^n \int_0^x \frac{t^n}{1 + t}\,\dd t ,

y para 0x10 \leq x \leq 1 el resto está acotado por 0xtn ⁣dt=xn+1n+1\int_0^x t^n\,\dd t = \frac{x^{n+1}}{n+1}. Esto es más fuerte que Taylor–Young en dos sentidos: es una identidad válida para un xx fijo (no solo cuando x0x \to 0), y la cota del error es numérica. El problema del fin de semana (Problema 16.1) vive de formas exactas así; Taylor–Young, más abajo, es la herramienta ligera para los límites, donde solo importa la forma del error.

Teorema 16.7 (Desigualdad de Taylor–Lagrange)

Sea ff de clase Cn+1C^{n+1} con f(n+1)M\abs{f^{(n+1)}} \leq M entre aa y xx. Entonces

f(x)k=0nf(k)(a)k!(xa)kMxan+1(n+1)!.\Bigl| f(x) - \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}(x-a)^k \Bigr| \leq M\, \frac{\abs{x - a}^{n+1}}{(n+1)!} .

Demostración. Acótese el resto integral: ax(xt)nn!f(n+1)(t) ⁣dtMaxxtnn! ⁣dt=Mxan+1(n+1)!\bigl|\int_a^x \frac{(x-t)^n}{n!} f^{(n+1)}(t)\,\dd t\bigr| \leq M \bigl|\int_a^x \frac{\abs{x-t}^n}{n!}\dd t\bigr| = M\frac{\abs{x-a}^{n+1}}{(n+1)!}.

Ejemplo 16.8 (Cálculo numérico certificado)

¿Cuánto vale 1.02\sqrt{1.02}? Aplíquese Taylor–Lagrange a f(t)=1+tf(t) = \sqrt{1 + t} en a=0a = 0, al orden 22, con x=0.02x = 0.02:

1.021+0.0220.0228=1.00995,f(t)=38(1+t)5/238,\sqrt{1.02} \approx 1 + \frac{0.02}{2} - \frac{0.02^2}{8} = 1.00995 , \qquad \abs{f'''(t)} = \frac{3}{8}(1+t)^{-5/2} \leq \frac 38 ,

luego el error es a lo sumo 380.0236=5107\frac38 \cdot \frac{0.02^3}{6} = 5\cdot10^{-7}: 1.02=1.00995\sqrt{1.02} = 1.00995 con seis decimales certificados (valor verdadero 1.00995049381.0099504938\dots — la cota es casi óptima). La idea de cierre: Taylor–Young dice solamente con qué rapidez se anula el error; Taylor–Lagrange convierte el mismo polinomio en un certificado, un número más una barra de error demostrada. Cada vez que en este libro se afirma un decimal, hay detrás una cota de tipo Lagrange; el problema del fin de semana (Problema 16.1) industrializa la idea.

Teorema 16.9 (Taylor–Young)

Sea ff nn veces derivable en aa. Entonces, cuando xax \to a:

f(x)=k=0nf(k)(a)k!(xa)k+o((xa)n).f(x) = \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}\,(x-a)^k + o\bigl((x-a)^n\bigr) .

Demostración. Inducción sobre nn. Para n=1n = 1, esto es la definición de la derivada (Definición 14.1). Supóngase el enunciado al orden n1n - 1 y sea ff nn veces derivable en aa. Aplíquese la hipótesis de inducción a ff' (que es n1n-1 veces derivable en aa):

f(t)=k=0n1f(k+1)(a)k!(ta)k+r(t),r(t)=o((ta)n1).f'(t) = \sum_{k=0}^{n-1} \frac{f^{(k+1)}(a)}{k!}(t-a)^k + r(t), \qquad r(t) = o\bigl((t-a)^{n-1}\bigr).

Sea g(x)=f(x)k=0nf(k)(a)k!(xa)kg(x) = f(x) - \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}(x-a)^k; entonces g=rg' = r y g(a)=0g(a) = 0. Dado ε>0\varepsilon > 0, elíjase δ\delta con r(t)εtan1\abs{r(t)} \leq \varepsilon\abs{t - a}^{n-1} para taδ\abs{t-a} \leq \delta; la desigualdad del valor medio (Teorema 14.9) aplicada en el segmento de aa a xx (donde gεxan1\abs{g'} \leq \varepsilon\abs{x-a}^{n-1}) da g(x)εxan\abs{g(x)} \leq \varepsilon\abs{x - a}^n: exactamente g(x)=o((xa)n)g(x) = o((x-a)^n).

Observación 16.10 (Tres fórmulas, tres precios, tres productos)

Las hipótesis se gradúan exactamente con las conclusiones. Taylor–Young pide lo mínimo (nn derivadas en el punto, nada más) y da lo mínimo: un o((xa)n)o((x-a)^n) cualitativo, perfecto para los límites e inútil para cifras certificadas. La desigualdad de Lagrange pide Cn+1C^{n+1} en el intervalo y una cota MM ahí, y devuelve una barra numérica de error. La forma integral pide la misma regularidad y devuelve lo máximo: el error como objeto explícito que se puede transformar (integrar por partes, acotar a trozos, cambiar de variable) — es la forma que movía la máquina de irracionalidad del Problema 15.1. Elegir la fórmula más débil que sostiene el objetivo no es pedantería: la función plana del Problema 16.1 cumple Taylor–Young a todos los órdenes, mientras que toda conclusión más fuerte sobre ella es falsa lejos de 00.

Proposición 16.11 (Desarrollos usuales en 00)

Cuando x0x \to 0, para todo orden nn fijo:

ex=1+x+x22!++xnn!+o(xn),cosx=1x22!+x44!+(1)px2p(2p)!+o(x2p+1),sinx=xx33!++(1)px2p+1(2p+1)!+o(x2p+2),11x=1+x+x2++xn+o(xn),ln(1+x)=xx22+x33+(1)n1xnn+o(xn),(1+x)α=1+αx+α(α1)2!x2++(αn)xn+o(xn),\begin{align*} \eu^x &= 1 + x + \frac{x^2}{2!} + \dots + \frac{x^n}{n!} + o(x^n),\\ \cos x &= 1 - \frac{x^2}{2!} + \frac{x^4}{4!} - \dots + \frac{(-1)^p x^{2p}}{(2p)!} + o(x^{2p+1}),\\ \sin x &= x - \frac{x^3}{3!} + \dots + \frac{(-1)^p x^{2p+1}}{(2p+1)!} + o(x^{2p+2}),\\ \frac{1}{1 - x} &= 1 + x + x^2 + \dots + x^n + o(x^n),\\ \ln(1 + x) &= x - \frac{x^2}{2} + \frac{x^3}{3} - \dots + \frac{(-1)^{n-1} x^n}{n} + o(x^n),\\ (1 + x)^\alpha &= 1 + \alpha x + \frac{\alpha(\alpha-1)}{2!}x^2 + \dots + \binom{\alpha}{n} x^n + o(x^n), \end{align*}

donde (αn)=α(α1)(αn+1)n!\binom{\alpha}{n} = \frac{\alpha(\alpha - 1)\cdots(\alpha - n + 1)}{n!} para α\alpha real. (cosh\cosh y sinh\sinh: igual que cos\cos y sin\sin, pero sin los signos alternados.)

Demostración. Cada función es regular cerca de 00 con derivadas fáciles de evaluar: (ex)(k)=ex(\eu^x)^{(k)} = \eu^x; las derivadas de sin\sin y de cos\cos son cíclicas de período 44; ((1+x)α)(k)=α(α1)(αk+1)(1+x)αk\bigl((1+x)^\alpha\bigr)^{(k)} = \alpha(\alpha - 1)\cdots(\alpha - k + 1)(1+x)^{\alpha - k}; (ln(1+x))(k)=(1)k1(k1)!(1+x)k\bigl(\ln(1+x)\bigr)^{(k)} = \frac{(-1)^{k-1}(k-1)!}{(1+x)^k}. Aplíquese Taylor–Young en a=0a = 0. (El geométrico es exacto: 11x0nxk=xn+11x=o(xn)\frac{1}{1-x} - \sum_0^n x^k = \frac{x^{n+1}}{1 - x} = o(x^n).)

Método 16.12 (Calcular con desarrollos)

  1. Fíjese primero el orden objetivo nn y trúnquese ahí todo resultado intermedio — arrastrar términos de orden superior es trabajo perdido, y dejar caer los inferiores es un error.
  2. Sumas y productos: desarróllese cada factor al orden nn y multiplíquese, descartando más allá de xnx^n.
  3. Composición f(u(x))f(u(x)) con u(x)0u(x) \to 0: sustitúyase el desarrollo de uu en el de ff, orden a orden.
  4. Cocientes: escríbase 11+v\frac{1}{1 + v} con v0v \to 0 y úsese el desarrollo geométrico.
  5. Intégrese un desarrollo término a término (derivar exige más cuidado — justificación: la integral de o(tn)o(t^n) de 00 a xx es o(xn+1)o(x^{n+1}), por acotación directa).

Ejemplo 16.13 (Composición, con la contabilidad a la vista)

Desarróllese esinx\eu^{\sin x} al orden 33. Desarrollo interior: u=sinx=xx36+o(x3)u = \sin x = x - \frac{x^3}{6} + o(x^3), que en efecto tiende a 00. Exterior: eu=1+u+u22+u36+o(u3)\eu^u = 1 + u + \frac{u^2}{2} + \frac{u^3}{6} + o(u^3), y o(u3)=o(x3)o(u^3) = o(x^3) puesto que uxu \sim x. Potencias de uu, truncadas en x3x^3:

u2=x2+o(x3),u3=x3+o(x3)u^2 = x^2 + o(x^3), \qquad u^3 = x^3 + o(x^3)

(el término cruzado 2x(x36)2x\cdot(-\frac{x^3}{6}) ya es x4x^4). Ensámblese:

esinx=1+(xx36)+x22+x36+o(x3)=1+x+x22+o(x3):\eu^{\sin x} = 1 + \Bigl(x - \frac{x^3}{6}\Bigr) + \frac{x^2}{2} + \frac{x^3}{6} + o(x^3) = 1 + x + \frac{x^2}{2} + o(x^3) :

las dos contribuciones en x3x^3 se cancelan exactamente. La idea de cierre: esinx\eu^{\sin x} y ex\eu^x coinciden hasta el orden 33 — no porque sinxx\sin x \approx x groseramente, sino porque el primer desacuerdo de los exponentes (x36-\frac{x^3}{6}) entra multiplicado por e0=1\eu^0 = 1 y se topa después con el término cúbico de la exponencial exterior; la contabilidad orden a orden detecta esas conspiraciones, y el ojo nunca. (El término siguiente es x48-\frac{x^4}{8}: la tregua acaba en el orden 44.)

Ejemplo 16.14

Desarrollo de tan\tan al orden 55. Escríbase tanx=sinx1cosx\tan x = \sin x \cdot \frac{1}{\cos x}:

1cosx=11(x22x424+o(x5))=1+(x22x424)+(x22) ⁣2+o(x5)=1+x22+5x424+o(x5),\frac{1}{\cos x} = \frac{1}{1 - \bigl(\frac{x^2}{2} - \frac{x^4}{24} + o(x^5)\bigr)} = 1 + \Bigl(\frac{x^2}{2} - \frac{x^4}{24}\Bigr) + \Bigl(\frac{x^2}{2}\Bigr)^{\!2} + o(x^5) = 1 + \frac{x^2}{2} + \frac{5x^4}{24} + o(x^5),

y después

tanx=(xx36+x5120)(1+x22+5x424)+o(x5)=x+x33+2x515+o(x5).\tan x = \Bigl(x - \frac{x^3}{6} + \frac{x^5}{120}\Bigr) \Bigl(1 + \frac{x^2}{2} + \frac{5x^4}{24}\Bigr) + o(x^5) = x + \frac{x^3}{3} + \frac{2x^5}{15} + o(x^5) .
El seno (oscuro) frente a sus polinomios de Taylor en 0: T_1 = x, T_3 = x - x3/6 (a trazos) y T_5 = x - x3/6 + x5/120 (punteado). Cada nuevo par de términos abraza la curva en una ventana visiblemente más ancha, pero todo polinomio acaba despegándose: un desarrollo de Taylor es un contrato local, afinado en 0 y mudo a lo lejos. La cota de Taylor–Lagrange |x|n+1(n+1)! cuantifica la ventana; y la función plana del  muestra que el contrato puede incluso ser vacío más allá del propio punto.
El seno (oscuro) frente a sus polinomios de Taylor en 00: T1=xT_1 = x, T3=xx36T_3 = x - \frac{x^3}{6} (a trazos) y T5=xx36+x5120T_5 = x - \frac{x^3}{6} + \frac{x^5}{120} (punteado). Cada nuevo par de términos abraza la curva en una ventana visiblemente más ancha, pero todo polinomio acaba despegándose: un desarrollo de Taylor es un contrato local, afinado en 00 y mudo a lo lejos. La cota de Taylor–Lagrange xn+1(n+1)!\frac{\abs{x}^{n+1}}{(n+1)!} cuantifica la ventana; y la función plana del Problema 16.1 muestra que el contrato puede incluso ser vacío más allá del propio punto.

16.3 Aplicaciones

Ejemplo 16.15 (Límites)

limx0xsinxx3:xsinx=x36+o(x3)x36,luego el lıˊmite es 16\lim_{x \to 0} \frac{x - \sin x}{x^3}: \qquad x - \sin x = \frac{x^3}{6} + o(x^3) \sim \frac{x^3}{6}, \qquad\text{luego el límite es } \frac16

— lo que zanja la cuestión planteada en el Ejercicio 4.9. Del mismo modo, limx0(sinxx)1/x2\displaystyle\lim_{x\to0}\Bigl(\frac{\sin x}{x}\Bigr)^{1/x^2}: el logaritmo es

1x2ln(1x26+o(x2))=1x2(x26+o(x2))16,lıˊmite e1/6.\frac{1}{x^2}\ln\Bigl(1 - \frac{x^2}{6} + o(x^2)\Bigr) = \frac{1}{x^2}\Bigl(-\frac{x^2}{6} + o(x^2)\Bigr) \longrightarrow -\frac16, \qquad\text{límite } \eu^{-1/6}.

Observación 16.16 (Errores frecuentes con los desarrollos)

(i) No se suman ni se restan nunca equivalentes: de tanxx\tan x \sim x y sinxx\sin x \sim x no se puede concluir tanxsinx0\tan x - \sin x \sim 0 (que no significa nada) — el camino honesto son los desarrollos:

tanxsinx=(x+x33)(xx36)+o(x3)=x32+o(x3)x32.\tan x - \sin x = \Bigl(x + \frac{x^3}{3}\Bigr) - \Bigl(x - \frac{x^3}{6}\Bigr) + o(x^3) = \frac{x^3}{2} + o(x^3) \sim \frac{x^3}{2} .

(ii) Desarróllese más allá de la masacre: en ese mismo cálculo, el orden 11 solo ve xx=0x - x = 0; siempre que los términos dominantes se cancelen, súbase el orden hasta que sobreviva un coeficiente no nulo y solo entonces vuélvase a un equivalente. (iii) Los equivalentes no pasan por las exponenciales: n2+nn2n^2 + n \sim n^2 y, sin embargo, en2+n=enen2\eu^{n^2+n} = \eu^{n}\,\eu^{n^2} no es equivalente a en2\eu^{n^2} — exponéncense solo desarrollos del exponente cuyo error tienda a 00, nunca equivalentes del exponente. (Los logaritmos son más benignos: si unvn1u_n \sim v_n \to \ell \neq 1, con >0\ell > 0, entonces lnunlnvn\ln u_n \sim \ln v_n.) (iv) El cálculo con o()o(\cdot) es de un solo sentido: o(x2)+o(x2)=o(x2)o(x^2) + o(x^2) = o(x^2), 5o(x2)=o(x2)5\,o(x^2) = o(x^2), xo(x2)=o(x3)x\cdot o(x^2) = o(x^3) — pero un o(x2)o(x^2) no es una función concreta, así que no se cancelen nunca dos de ellos entre sí: o(x2)o(x2)o(x^2) - o(x^2) es o(x2)o(x^2), no 00.

Proposición 16.17 (Comportamiento local)

Supóngase f(x)=f(a)+c(xa)p+o((xa)p)f(x) = f(a) + c\,(x - a)^p + o\bigl((x-a)^p\bigr) con c0c \neq 0 (primer término no nulo tras la constante; p2p \geq 2 en un punto crítico).

  • Si pp es par: ff tiene un mínimo local en aa si c>0c > 0, y un máximo local si c<0c < 0.
  • Si pp es impar: no hay extremo (ff(a)f - f(a) cambia de signo); y si además el desarrollo empieza tras un término lineal f(a)(xa)f'(a)(x - a), la gráfica cruza su tangente: es un punto de inflexión.

Demostración. Cerca de aa, f(x)f(a)=(xa)p(c+o(1))f(x) - f(a) = (x-a)^p\bigl(c + o(1)\bigr) tiene el signo de c(xa)pc\,(x-a)^p: constante si pp es par, cambiante si pp es impar.

Ejemplo 16.18 (Los exponentes hay que desarrollarlos hasta el o(1)o(1))

Hállese un equivalente de un=(1+1n)n2u_n = \bigl(1 + \frac1n\bigr)^{n^2}. Desarróllese el exponente hasta que su error tienda a 00:

n2ln(1+1n)=n2(1n12n2+O(1n3))=n12+O(1n),n^2 \ln\Bigl(1 + \frac1n\Bigr) = n^2\Bigl(\frac1n - \frac{1}{2n^2} + O\Bigl(\frac{1}{n^3}\Bigr)\Bigr) = n - \frac12 + O\Bigl(\frac1n\Bigr),

luego un=en1/2eO(1/n)u_n = \eu^{\,n - 1/2}\,\eu^{O(1/n)} con eO(1/n)1\eu^{O(1/n)} \to 1:

un    en12.u_n \;\sim\; \eu^{\,n - \frac12} .

Obsérvese qué habría salido mal con menos cuidado: detener el exponente en n21n=n+O(1)n^2\cdot\frac1n = n + O(1) deja un factor eO(1)\eu^{O(1)} —acotado, pero que no tiende a 11— y no se puede afirmar ningún equivalente. La regla de los errores frecuentes de más arriba, en versión positiva: un equivalente de ean\eu^{a_n} exige el desarrollo de ana_n hasta un término que tienda a cero, conservando exactamente todos los coeficientes anteriores.

Ejemplo 16.19 (Clasificar un punto crítico plano)

Estúdiese f(x)=cosx+x22f(x) = \cos x + \frac{x^2}{2} cerca de 00. Se tiene f(0)=0f'(0) = 0 y f(0)=cos0+1=0f''(0) = -\cos 0 + 1 = 0: el criterio de la segunda derivada enmudece. Desarróllese en su lugar:

f(x)=(1x22+x424+o(x4))+x22=1+x424+o(x4):f(x) = \Bigl(1 - \frac{x^2}{2} + \frac{x^4}{24} + o(x^4)\Bigr) + \frac{x^2}{2} = 1 + \frac{x^4}{24} + o(x^4) :

primer término no nulo cxpc\,x^p con p=4p = 4 par y c=124>0c = \frac{1}{24} > 0: un mínimo local, de una planitud inusual (la gráfica abandona su valor mínimo como x4x^4, no como x2x^2). La idea de cierre: el desarrollo ve en una línea lo que la derivación iterada oscurece — y la Proposición 16.17 es el diccionario sistemático que va del «primer término superviviente» a la «forma local».

Ejemplo 16.20 (Desarrollos en el infinito)

Dos cálculos en los que la variable se va a ++\infty y la sustitución h=1x0+h = \frac1x \to 0^+ importa todo el utillaje. Primero,

xx2ln(1+1x)=xx2(1x12x2+O(1x3))=12+O(1x)12.x - x^2\ln\Bigl(1 + \frac1x\Bigr) = x - x^2\Bigl(\frac1x - \frac{1}{2x^2} + O\Bigl(\frac{1}{x^3}\Bigr)\Bigr) = \frac12 + O\Bigl(\frac1x\Bigr) \longrightarrow \frac12 .

Segundo, el arcotangente en el infinito: de arctanx+arctan1x=π2\arctan x + \arctan\frac1x = \frac\pi2 para x>0x > 0 (Proposición 4.12) y del desarrollo de arctan\arctan en 00 (Ejercicio 16.3),

arctanx=π2arctan1x=π21x+13x3+o(1x3):\arctan x = \frac\pi2 - \arctan\frac1x = \frac\pi2 - \frac1x + \frac{1}{3x^3} + o\Bigl(\frac{1}{x^3}\Bigr) :

la gráfica se acerca a su asíntota y=π2y = \frac\pi2 por debajo, a velocidad 1x\frac1x. La idea de cierre: no hay una teoría aparte de los desarrollos en el infinito — una sola sustitución recíproca los reduce a desarrollos en 00, siempre que todo OO y todo oo intermedio se arrastre con honestidad.

Ejemplo 16.21 (Asíntota por desarrollo)

Cuando x+x \to +\infty,

x2+x=x1+1x=x(1+12x18x2+o(1x2))=x+1218x+o(1x):\sqrt{x^2 + x} = x\sqrt{1 + \tfrac1x} = x\Bigl(1 + \frac{1}{2x} - \frac{1}{8x^2} + o\bigl(\tfrac{1}{x^2}\bigr)\Bigr) = x + \frac12 - \frac{1}{8x} + o\bigl(\tfrac 1x\bigr):

la recta y=x+12y = x + \frac12 es una asíntota, alcanzada por debajo (el término siguiente 18x-\frac{1}{8x} es negativo).

Observación 16.22 (Dónde funcionan los desarrollos a continuación)

Los desarrollos asintóticos son el idioma corriente del resto del libro: en el Capítulo 17 deciden la convergencia (los equivalentes alimentan los criterios de comparación, y el estudio de 1nα\sum \frac{1}{n^\alpha} es un desarrollo disfrazado); en el volumen del segundo año se convierten en series de potencias, donde el polinomio de Taylor adquiere infinitos términos y un radio de convergencia; y toda linealización en física —el péndulo, la perturbación de primer orden— es un enunciado de Taylor–Young con el o()o(\cdot) suprimido en silencio. El único aviso que merece grabarse: un desarrollo describe una función solo cerca de un punto — véase la función plana del problema del fin de semana, cuyo desarrollo en 00 es idénticamente nulo sin que la función lo sea.

Observación 16.23 (Perspectivas dentro de este volumen)

Los desarrollos son el idioma de trabajo del análisis restante y de la geometría que viene. Capítulo 17 los convierte en veredictos de convergencia: un equivalente del término general es un desarrollo truncado en su primer término, y los criterios más finos (los alternados con control del error) consumen también el segundo. Capítulo 24 lee la geometría local en los desarrollos de las dos funciones coordenadas: que una curva parametrizada se cruce, sea tangente o presente un punto de retroceso se decide por qué potencias de tt sobreviven en x(t)x(t) e y(t)y(t) — la versión plana de la Proposición 16.17. Y el Capítulo 25 se detiene en el orden uno a propósito: el plano tangente es un enunciado de Taylor–Young en dos variables, y toda la teoría de segundo orden (hessianas, puntos de silla) queda aplazada al volumen del segundo año. El hilo común: toda pregunta «local» de este libro se responde escribiendo el primer término superviviente de un desarrollo.

16.4 Ejercicios

Ejercicio 16.1

Dense los desarrollos en 00: de e2x\eu^{2x} al orden 33; de   ln(1x)\;\ln(1 - x) al orden 44; de   1+x\;\sqrt{1 + x} al orden 33; de   11+x2\;\dfrac{1}{1 + x^2} al orden 66.

Solución

Solución de Ejercicio 16.1.

e2x=1+2x+2x2+4x33+o(x3);ln(1x)=xx22x33x44+o(x4);\eu^{2x} = 1 + 2x + 2x^2 + \frac{4x^3}{3} + o(x^3); \qquad \ln(1 - x) = -x - \frac{x^2}{2} - \frac{x^3}{3} - \frac{x^4}{4} + o(x^4);
1+x=1+x2x28+x316+o(x3);11+x2=1x2+x4x6+o(x6),\sqrt{1+x} = 1 + \frac x2 - \frac{x^2}{8} + \frac{x^3}{16} + o(x^3); \qquad \frac{1}{1 + x^2} = 1 - x^2 + x^4 - x^6 + o(x^6),

el último sustituyendo x2-x^2 en el desarrollo geométrico.

Ejercicio 16.2

Calcúlense los límites:

limx0ex1xx2,limx0cosx1x2x4,limx0ln(1+x)sinxx2.\lim_{x\to 0} \frac{\eu^x - 1 - x}{x^2}, \qquad \lim_{x\to 0} \frac{\cos x - \sqrt{1 - x^2}}{x^4}, \qquad \lim_{x\to 0} \frac{\ln(1+x) - \sin x}{x^2}.
Solución

Solución de Ejercicio 16.2.

ex1x=x22+o(x2)\eu^x - 1 - x = \frac{x^2}{2} + o(x^2): límite 12\dfrac12.

cosx=1x22+x424+o(x4)\cos x = 1 - \frac{x^2}{2} + \frac{x^4}{24} + o(x^4) y 1x2=1x22x48+o(x4)\sqrt{1 - x^2} = 1 - \frac{x^2}{2} - \frac{x^4}{8} + o(x^4): diferencia x424+x48=x46+o(x4)\frac{x^4}{24} + \frac{x^4}{8} = \frac{x^4}{6} + o(x^4): límite 16\dfrac16.

ln(1+x)sinx=(xx22)x+o(x2)=x22+o(x2)\ln(1+x) - \sin x = \bigl(x - \frac{x^2}{2}\bigr) - x + o(x^2) = -\frac{x^2}{2} + o(x^2): límite 12-\dfrac12.

Ejercicio 16.3

Desarróllese arctanx\arctan x en 00 al orden 55 integrando el desarrollo de 11+x2\frac{1}{1 + x^2}, y arcsinx\arcsin x al orden 55 integrando el de (1x2)1/2(1 - x^2)^{-1/2}.

Solución

Solución de Ejercicio 16.3.

11+t2=1t2+t4+o(t5)\frac{1}{1+t^2} = 1 - t^2 + t^4 + o(t^5); integrando de 00 a xx (Método 16.12 (5)):

arctanx=xx33+x55+o(x5)  (incluso o(x6), por imparidad).\arctan x = x - \frac{x^3}{3} + \frac{x^5}{5} + o(x^5)\ \ (\text{incluso }o(x^6)\text{, por imparidad}).

(1t2)1/2=1+t22+38t4+o(t5)(1 - t^2)^{-1/2} = 1 + \frac{t^2}{2} + \frac38 t^4 + o(t^5) (desarrollo binomial con α=12\alpha = -\frac12, x=t2x = -t^2: (1/22)=(12)(32)2=38\binom{-1/2}{2} = \frac{(-\frac12)(-\frac32)}{2} = \frac38); integrando:

arcsinx=x+x36+3x540+o(x5).\arcsin x = x + \frac{x^3}{6} + \frac{3x^5}{40} + o(x^5) .

Ejercicio 16.4

Usando Taylor–Lagrange para exp\exp en [0,1]\intcc{0}{1}, demuéstrese que

ek=0n1k!3(n+1)!,\Bigl| \eu - \sum_{k=0}^{n} \frac{1}{k!} \Bigr| \leq \frac{3}{(n+1)!},

y determínese un nn que garantice 66 decimales exactos de e\eu.

Solución

Solución de Ejercicio 16.4.

Taylor–Lagrange (Teorema 16.7) para exp\exp en a=0a = 0, x=1x = 1: la derivada de orden (n+1)(n+1) es ete<3\eu^t \leq \eu < 3 en [0,1]\intcc{0}{1}, luego

ek=0n1k!3(n+1)!.\Bigl|\eu - \sum_{k=0}^{n} \frac{1}{k!}\Bigr| \leq \frac{3}{(n+1)!} .

Para 66 decimales exactos se quiere 3(n+1)!<5×107\frac{3}{(n+1)!} < 5\times 10^{-7}, es decir, (n+1)!>6×106(n+1)! > 6\times 10^{6}: como 10!=362880010! = 3\,628\,800 y 11!=3991680011! = 39\,916\,800, basta n+1=11n + 1 = 11, es decir, n=10n = 10.

Ejercicio 16.5 ★★

Desarróllese al orden 22 en 1n\frac1n y dedúzcanse el límite y la velocidad de convergencia:

(1+1n) ⁣n=e(112n+1124n2+o(1n2)).\Bigl(1 + \frac 1n\Bigr)^{\!n} = \eu\Bigl(1 - \frac{1}{2n} + \frac{11}{24n^2} + o\Bigl(\frac{1}{n^2}\Bigr)\Bigr).
Solución

Solución de Ejercicio 16.5.

nln(1+1n)=n(1n12n2+13n3+o(1n3))=112n+13n2+o(1n2)n\ln\bigl(1 + \frac1n\bigr) = n\Bigl(\frac1n - \frac{1}{2n^2} + \frac{1}{3n^3} + o\bigl(\frac{1}{n^3}\bigr)\Bigr) = 1 - \frac{1}{2n} + \frac{1}{3n^2} + o\bigl(\frac{1}{n^2}\bigr). Exponenciando, con u=12n+13n2u = -\frac{1}{2n} + \frac{1}{3n^2} y eu=1+u+u22+o(u2)\eu^u = 1 + u + \frac{u^2}2 + o(u^2):

(1+1n)n=eeu=e(112n+13n2+18n2+o(1n2))=e(112n+1124n2+o(1n2)).\Bigl(1 + \frac1n\Bigr)^n = \eu\cdot \eu^{u} = \eu\Bigl(1 - \frac{1}{2n} + \frac{1}{3n^2} + \frac{1}{8n^2} + o\Bigl(\frac{1}{n^2}\Bigr)\Bigr) = \eu\Bigl(1 - \frac{1}{2n} + \frac{11}{24n^2} + o\Bigl(\frac{1}{n^2}\Bigr)\Bigr).

Límite e\eu; el error es e2n\sim \dfrac{\eu}{2n}: lento (una cifra por cada multiplicación de nn por diez).

Ejercicio 16.6 ★★

Estúdiese el comportamiento local en 00 de f(x)=x2x4f(x) = x^2 - x^4 y de g(x)=x3+x5g(x) = x^3 + x^5; y hállese la posición de la gráfica de h(x)=exh(x) = \eu^x respecto de su tangente en a=1a = 1, primero localmente y después globalmente.

Solución

Solución de Ejercicio 16.6.

f(x)=x2x4=x2(1+o(1))f(x) = x^2 - x^4 = x^2(1 + o(1)): primer término x2x^2, p=2p = 2 par, coeficiente >0> 0: mínimo local en 00 (no global: f(2)=12f(2) = -12).

g(x)=x3+x5g(x) = x^3 + x^5: primer término x3x^3, con pp impar: no hay extremo; gg cruza su tangente (horizontal): punto de inflexión en 00.

h=exph = \exp en a=1a = 1: h(x)=e+e(x1)+e2(x1)2+o((x1)2)h(x) = \eu + \eu(x-1) + \frac{\eu}{2}(x-1)^2 + o((x-1)^2); la diferencia con la tangente es e2(x1)2+o()>0\frac{\eu}{2}(x-1)^2 + o(\cdot) > 0 localmente. Globalmente: exex0\eu^x - \eu x \geq 0 para todo xx por convexidad (Teorema 14.19 (3)): la gráfica queda por encima de toda tangente, con igualdad solo en el punto de contacto.

Ejercicio 16.7 ★★

Determínense las asíntotas en ±\pm\infty de f(x)=x3+x23f(x) = \sqrt[3]{x^3 + x^2} y la posición de la curva respecto de ellas.

Solución

Solución de Ejercicio 16.7.

Para x+x \to +\infty:

f(x)=x(1+1x)1/3=x(1+13x19x2+o(1x2))=x+1319x+o(1x):f(x) = x\Bigl(1 + \frac1x\Bigr)^{1/3} = x\Bigl(1 + \frac{1}{3x} - \frac{1}{9x^2} + o\Bigl(\frac{1}{x^2}\Bigr)\Bigr) = x + \frac13 - \frac{1}{9x} + o\Bigl(\frac1x\Bigr):

asíntota y=x+13y = x + \frac13, con la curva por debajo cerca de ++\infty. Cuando xx \to -\infty, el mismo cálculo vale (la raíz cúbica está definida para todos los reales y 1x0\frac1x \to 0): la misma asíntota y=x+13y = x + \frac13, pero ahora 19x>0-\frac{1}{9x} > 0: la curva queda por encima de la recta.

Ejercicio 16.8 ★★

Hállese el equivalente, cuando nn \to \infty, de

un=n+1n,vn=ln(n+1)lnn,wn=sin1ntan1n,u_n = \sqrt{n+1} - \sqrt n, \qquad v_n = \ln(n+1) - \ln n, \qquad w_n = \sin\frac{1}{n} - \tan\frac{1}{n},

cada uno como una potencia de nn por una constante.

Solución

Solución de Ejercicio 16.8.

un=n(1+1n1)=n(12n+o(1n))12nu_n = \sqrt n\bigl(\sqrt{1 + \tfrac1n} - 1\bigr) = \sqrt n\bigl(\frac{1}{2n} + o(\frac1n)\bigr) \sim \dfrac{1}{2\sqrt n}.

vn=ln(1+1n)1nv_n = \ln\bigl(1 + \frac1n\bigr) \sim \dfrac 1n.

wnw_n: con h=1n0h = \frac1n \to 0, sinhtanh=(hh36)(h+h33)+o(h3)=h32+o(h3)\sin h - \tan h = \bigl(h - \frac{h^3}{6}\bigr) - \bigl(h + \frac{h^3}{3}\bigr) + o(h^3) = -\frac{h^3}{2} + o(h^3), luego wn12n3w_n \sim -\dfrac{1}{2n^3}.

Ejercicio 16.9 ★★★

Sea ff de clase C2C^2 en R\R. Demuéstrese que para todos xx y h>0h > 0:

f(x)f(x+h)f(xh)2h+h2sup[xh,x+h]f,\abs{f'(x)} \leq \frac{\abs{f(x+h) - f(x-h)}}{2h} + \frac{h}{2}\sup_{\intcc{x-h}{x+h}}\abs{f''} ,

y dedúzcase la desigualdad de tipo Landau–Kolmogórov: si fM0\abs f \leq M_0 y fM2\abs{f''} \leq M_2 en R\R, entonces f2M0M2\abs{f'} \leq \sqrt{2 M_0 M_2} en todas partes. (Optimícese en hh.)

Solución

Solución de Ejercicio 16.9.

Taylor–Lagrange al orden 11 en torno a xx, por los dos lados:

f(x+h)=f(x)+hf(x)+R+,f(xh)=f(x)hf(x)+R,R±h22supf.f(x + h) = f(x) + h f'(x) + R_+,\quad f(x - h) = f(x) - h f'(x) + R_-, \qquad \abs{R_\pm} \leq \frac{h^2}{2} \sup \abs{f''} .

Restando: f(x+h)f(xh)=2hf(x)+(R+R)f(x+h) - f(x-h) = 2h f'(x) + (R_+ - R_-), luego

f(x)f(x+h)f(xh)2h+h2sup[xh,x+h]f.\abs{f'(x)} \leq \frac{\abs{f(x+h) - f(x-h)}}{2h} + \frac{h}{2}\sup_{\intcc{x-h}{x+h}}\abs{f''} .

Con cotas globales: f(x)M0h+M2h2\abs{f'(x)} \leq \frac{M_0}{h} + \frac{M_2 h}{2} para todo h>0h > 0. El miembro derecho se minimiza en h=2M0/M2h = \sqrt{2M_0/M_2} (derivada nula), con valor 2M0M2\sqrt{2M_0M_2} — luego f2M0M2\abs{f'} \leq \sqrt{2M_0M_2}. (Si M2=0M_2 = 0, hágase hh \to \infty: f=0f' = 0, coherente.)

Ejercicio 16.10 ★★★

La sucesión u0(0,π)u_0 \in \intoo{0}{\pi}, un+1=sinunu_{n+1} = \sin u_n decrece hacia 00 (justifíquese brevemente). Para hallar su velocidad, considérese vn=1un2v_n = \frac{1}{u_n^2}:

  1. usando el desarrollo de sin\sin, demuéstrese vn+1vn13v_{n+1} - v_n \to \frac13;
  2. con Cesàro (Ejercicio 11.10), dedúzcase vnn13\frac{v_n}{n} \to \frac13 y después el equivalente un3nu_n \sim \sqrt{\dfrac{3}{n}}.
Solución

Solución de Ejercicio 16.10.

En (0,π)\intoo{0}{\pi}: 0<sinu<u0 < \sin u < u, luego (un)(u_n) es estrictamente decreciente y positiva, y por tanto convergente; el límite es un punto fijo de sin\sin en [0,π]\intcc{0}{\pi}, y sin=\sin \ell = \ell obliga a =0\ell = 0 (ya que sinx<x\sin x < x para x>0x > 0).

  1. Usando sinu=uu36+o(u3)\sin u = u - \frac{u^3}{6} + o(u^3) cuando u0u \to 0:

    vn+1vn=1sin2un1un2=1un2((1un26+o(un2)) ⁣21)=1un2(un23+o(un2))13.v_{n+1} - v_n = \frac{1}{\sin^2 u_n} - \frac{1}{u_n^2} = \frac{1}{u_n^2}\Bigl(\Bigl(1 - \frac{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 .
  2. Por el Ejercicio 11.10 (3) (Cesàro para diferencias), vnn13\frac{v_n}{n} \to \frac13, es decir, vnn3v_n \sim \frac n3, es decir, un23nu_n^2 \sim \frac 3n: como un>0u_n > 0,

    un3n.u_n \sim \sqrt{\frac{3}{n}} .

Ejercicio 16.11 ★★

(Diferencias de infinitos) Calcúlense

limx0(1x21sin2x)ylimx0+(1x1ln(1+x))\lim_{x \to 0} \Bigl(\frac{1}{x^2} - \frac{1}{\sin^2 x}\Bigr) \qquad\text{y}\qquad \lim_{x \to 0^+} \Bigl(\frac 1x - \frac{1}{\ln(1 + x)}\Bigr)

reduciendo a común denominador y desarrollando numerador y denominador por separado.

Solución

Solución de Ejercicio 16.11.

Denominadores comunes. Primer límite:

1x21sin2x=sin2xx2x2sin2x,sin2x=(xx36+o(x4)) ⁣2=x2x43+o(x5),\frac{1}{x^2} - \frac{1}{\sin^2 x} = \frac{\sin^2 x - x^2}{x^2\sin^2 x}, \qquad \sin^2 x = \Bigl(x - \frac{x^3}{6} + o(x^4)\Bigr)^{\!2} = x^2 - \frac{x^4}{3} + o(x^5) ,

de modo que el numerador es x43+o(x4)-\frac{x^4}{3} + o(x^4) mientras que el denominador es x4\sim x^4: el límite es 13-\dfrac13.

Segundo: 1x1ln(1+x)=ln(1+x)xxln(1+x)\dfrac1x - \dfrac{1}{\ln(1+x)} = \dfrac{\ln(1+x) - x}{x\ln(1+x)}; el numerador es x22+o(x2)-\frac{x^2}{2} + o(x^2) y el denominador, x(x+o(x))x2x\bigl(x + o(x)\bigr) \sim x^2: el límite es 12-\dfrac12.

Ejercicio 16.12 ★★★

(Asintótica de raíces implícitas) Véase que para todo kNk \in \N^* la ecuación tanx=x\tan x = x tiene exactamente una solución xkx_k en (kππ2,kπ+π2)\intoo{k\pi - \frac\pi2}{k\pi + \frac\pi2}, que xk=kπ+π2εkx_k = k\pi + \frac\pi2 - \varepsilon_k con εk=arctan1xk\varepsilon_k = \arctan\frac{1}{x_k}, y dedúzcase el desarrollo

xk=kπ+π21kπ+o(1k)(k).x_k = k\pi + \frac\pi2 - \frac{1}{k\pi} + o\Bigl(\frac 1k\Bigr) \qquad (k \to \infty).
Solución

Solución de Ejercicio 16.12.

En Ik=(kππ2,kπ+π2)I_k = \intoo{k\pi - \frac\pi2}{k\pi + \frac\pi2}, la función g(x)=tanxxg(x) = \tan x - x tiene derivada tan2x0\tan^2 x \geq 0, que se anula solo en el único punto kπk\pi: gg es estrictamente creciente en IkI_k (Corolario 14.12 (2)), con límites -\infty y ++\infty en los extremos: exactamente un cero xkx_k. Para k1k \geq 1, g(kπ)=kπ<0g(k\pi) = -k\pi < 0, luego xk(kπ,kπ+π2)x_k \in \intoo{k\pi}{k\pi + \frac\pi2}: escríbase xk=kπ+π2εkx_k = k\pi + \frac\pi2 - \varepsilon_k con εk(0,π2)\varepsilon_k \in \intoo{0}{\frac\pi2}. Entonces

xk=tanxk=tan(π2εk)=1tanεkεk=arctan1xk,x_k = \tan x_k = \tan\Bigl(\frac\pi2 - \varepsilon_k\Bigr) = \frac{1}{\tan\varepsilon_k} \quad\Longrightarrow\quad \varepsilon_k = \arctan\frac{1}{x_k} ,

usando tanεk=1xk\tan\varepsilon_k = \frac{1}{x_k} y εk(0,π2)\varepsilon_k \in \intoo{0}{\frac\pi2}. Como xkkπx_k \geq k\pi \to \infty: εk0\varepsilon_k \to 0, y

εk=arctan1xk=1xk+O(1xk3)=1kπ+O(1)+O(1k3)=1kπ+O(1k2),\varepsilon_k = \arctan\frac{1}{x_k} = \frac{1}{x_k} + O\Bigl(\frac{1}{x_k^3}\Bigr) = \frac{1}{k\pi + O(1)} + O\Bigl(\frac{1}{k^3}\Bigr) = \frac{1}{k\pi} + O\Bigl(\frac{1}{k^2}\Bigr) ,

de donde xk=kπ+π21kπ+o(1k)x_k = k\pi + \frac\pi2 - \frac{1}{k\pi} + o\bigl(\frac1k\bigr).

16.5 Problema: sumas alternadas, cifras certificadas y la irracionalidad de cos1\cos 1

Problema 16.1

Problema del fin de semana — la estimación alternada SSnan+1\abs{S - S_n} \leq a_{n+1}: ln2\ln 2 y π\pi con decimales demostrados, la fórmula de Machin y cos1Q\cos 1 \notin \Q

Una suma alternada de términos decrecientes es el objeto más amable del análisis numérico: su error está acotado por el primer término omitido, con signo conocido. Este problema demuestra ese principio con el teorema de las sucesiones adyacentes y después lo gasta de tres maneras: decimales certificados para ln2\ln 2 (tres caminos rivales) y para π\pi (Leibniz y después la fórmula de Machin de 1706, todavía la idea que hay detrás de los cálculos récord durante siglos), la irracionalidad de cos1\cos 1, sin1\sin 1 y cosh1\cosh 1 y, como contrapeso, la igualdad de Taylor–Lagrange y la función plana cuyo desarrollo de Taylor miente. En todo el problema, el lenguaje de «series» es informal: toda suma de aquí es una sucesión de sumas parciales, como en el Ejemplo 11.12; la teoría propiamente dicha se abre en el Capítulo 17.

Parte I — La estimación alternada. Sea (ak)k0(a_k)_{k \geq 0} decreciente hacia 00 y Sn=k=0n(1)kakS_n = \sum_{k=0}^{n} (-1)^k a_k.

  1. Véase que (S2n+1)(S_{2n+1}) es no decreciente, que (S2n)(S_{2n}) es no creciente y que son adyacentes (Teorema 11.11): las dos convergen a un mismo SS con, para todo nn,

    S2n+1SS2n,SSnan+1,S_{2n+1} \leq S \leq S_{2n}, \qquad \abs{S - S_n} \leq a_{n+1} ,

    teniendo el error el signo del primer término omitido. Véase además que si el decrecimiento es estricto, todas estas desigualdades son estrictas.

  2. Primer dividendo: para x=1x = 1 en la serie exponencial, compárese con el Teorema 16.7 en a=0a = 0: véase que Tn=k=0n(1)kk!T_n = \sum_{k=0}^{n} \frac{(-1)^k}{k!} converge a e1\eu^{-1} con e1Tn<1(n+1)!\abs{\eu^{-1} - T_n} < \frac{1}{(n+1)!}.
  3. (Leibniz, 1674) A partir de la identidad finita exacta

    11+t2=k=0n(1)kt2k+(1)n+1t2n+21+t2,\frac{1}{1 + t^2} = \sum_{k=0}^{n} (-1)^k t^{2k} + \frac{(-1)^{n+1} t^{2n+2}}{1 + t^2} ,

    integrada en [0,1]\intcc{0}{1}, demuéstrese

    π4=k=0n(1)k2k+1+ρn,ρn12n+3.\frac\pi4 = \sum_{k=0}^{n} \frac{(-1)^k}{2k+1} + \rho_n, \qquad \abs{\rho_n} \leq \frac{1}{2n+3} .
  4. Lentitud: ¿cuántos términos de Leibniz garantizan seis decimales exactos de π\pi? (Unos dos millones.) Evalúese 4S4=4(113+1517+19)4 S_4 = 4\bigl(1 - \frac13 + \frac15 - \frac17 + \frac19\bigr) y su distancia a π\pi, para sentir el dolor.

Parte II — ln2\ln 2 de tres maneras.

  1. (Camino 1: la armónica alternada) De 11+t=k=0n1(1)ktk+(1)ntn1+t\frac{1}{1+t} = \sum_{k=0}^{n-1}(-1)^k t^k + \frac{(-1)^n t^n}{1+t} integrada en [0,1]\intcc{0}{1}:

    ln2=k=1n(1)k1k+(1)nRn,12(n+1)Rn1n+1:\ln 2 = \sum_{k=1}^{n} \frac{(-1)^{k-1}}{k} + (-1)^n R_n, \qquad \frac{1}{2(n+1)} \leq R_n \leq \frac{1}{n+1} :

    el error es de orden exactamente 1n\frac 1n — un millón de términos para seis decimales.

  2. (Camino 2: la serie rápida) Intégrese 11t2=k=0nt2k+t2n+21t2\frac{1}{1 - t^2} = \sum_{k=0}^{n} t^{2k} + \frac{t^{2n+2}}{1-t^2} de 00 a x(0,1)x \in \intoo{0}{1} y evalúese en x=13x = \frac13 (nótese 1+1/311/3=2\frac{1 + 1/3}{1 - 1/3} = 2):

    ln2=2k=0n(1/3)2k+12k+1+ρ~n,0<ρ~n94(1/3)2n+32n+3:\ln 2 = 2\sum_{k=0}^{n} \frac{(1/3)^{2k+1}}{2k+1} + \tilde\rho_n, \qquad 0 < \tilde\rho_n \leq \frac{9}{4}\cdot \frac{(1/3)^{2n+3}}{2n+3} :

    convergencia geométrica, aproximadamente una cifra por término.

  3. (Camino 3: sumas de Riemann y una identidad oculta) Demuéstrese por inducción la identidad

    k=12n(1)k1k=H2nHn=k=1n1n+k,\sum_{k=1}^{2n} \frac{(-1)^{k-1}}{k} = H_{2n} - H_n = \sum_{k=1}^{n} \frac{1}{n + k} ,

    y recupérese ln2\ln 2 como el límite de sumas de Riemann del Ejemplo 15.21: los caminos 1 y 3 son en secreto el mismo número visto dos veces.

  4. Duelo con seis términos: compárese k=16(1)k1k=0.6167\sum_{k=1}^{6} \frac{(-1)^{k-1}}{k} = 0.6167 con el camino 2 en n=5n = 5, que ya da ln2=0.693147\ln 2 = 0.693147 con error 1.1107\leq 1.1\cdot10^{-7}. Explíquese la razón estructural (punto de evaluación bien dentro del intervalo de convergencia frente a en su frontera).
  5. ¿Cuántos términos del camino 2 certifican diez decimales de ln2\ln 2? Véase que basta n=10n = 10.

Parte III — La fórmula de Machin.

  1. Calcúlese (5+i)4(5 + \iu)^4 y compruébese la identidad compleja

    (5+i)4=2(1+i)(239+i).(5 + \iu)^4 = 2\,(1 + \iu)\,(239 + \iu) .

    Tomando argumentos (con los convenios del Capítulo 3), dedúzcase la fórmula de Machin

    π4=4arctan15arctan1239.\frac\pi4 = 4\arctan\frac15 - \arctan\frac{1}{239} .

    (Compruébese que ningún argumento sale de (π2,π2)\intoo{-\frac\pi2}{\frac\pi2}.)

  2. Como en la pregunta 3, establézcase, para 0<x<10 < x < 1:

    arctanx=k=0n(1)kx2k+12k+1+rn(x),rn(x)x2n+32n+3.\arctan x = \sum_{k=0}^{n} \frac{(-1)^k x^{2k+1}}{2k+1} + r_n(x), \qquad \abs{r_n(x)} \leq \frac{x^{2n+3}}{2n+3} .
  3. Certifíquense siete decimales de π\pi con seis términos: acótese el error total de

    π16k=04(1)k(1/5)2k+12k+14(1239132393)\pi \approx 16\sum_{k=0}^{4} \frac{(-1)^k (1/5)^{2k+1}}{2k+1} - 4\Bigl(\frac{1}{239} - \frac{1}{3\cdot239^3}\Bigr)

    por 16(1/5)1111+4(1/239)55<510816\,\frac{(1/5)^{11}}{11} + 4\,\frac{(1/239)^5}{5} < 5\cdot10^{-8}, y dese el valor resultante 3.14159263.1415926\dots

  4. Compárense los tres caminos hacia π\pi ya disponibles —Leibniz (pregunta 4), las integrales de Dalzell del Problema 15.1 (error 415m4^{1-5m}) y Machin (error 165(2n+3)\approx 16\cdot 5^{-(2n+3)})— en cifras por término, y explíquese por qué encoger el punto de evaluación gana a todo lo demás.

Parte IV — La trampa del entero, edición alternada.

  1. Supóngase cos1=pq\cos 1 = \frac pq. Multiplíquese el encuadre alternado estricto de k(1)k(2k)!\sum_{k} \frac{(-1)^k}{(2k)!} (pregunta 1) por (2n)!(2n)! con 2nq2n \geq q y dedúzcase una contradicción: cos1\cos 1 es irracional.
  2. Adáptese a sin1=k(1)k(2k+1)!\sin 1 = \sum_k \frac{(-1)^k}{(2k+1)!} (multiplíquese por (2n+1)!(2n+1)!): sin1Q\sin 1 \notin \Q. Los dos irracionales y, sin embargo, cos21+sin21=1\cos^2 1 + \sin^2 1 = 1: la irracionalidad no es estable por operaciones algebraicas.
  3. El primo no alternado: cosh1=k1(2k)!\cosh 1 = \sum_k \frac{1}{(2k)!} (en el sentido de sumas parciales, con la cota bilateral de la cola 0<cosh1kn1(2k)!<2(2n+2)!0 < \cosh 1 - \sum_{k \leq n} \frac{1}{(2k)!} < \frac{2}{(2n+2)!}, que hay que demostrar). Conclúyase cosh1Q\cosh 1 \notin \Q con la misma trampa.
  4. Empújese hasta cos1m\cos\frac 1m para todo entero m1m \geq 1: multiplíquese por m2n(2n)!m^{2n}(2n)! y conclúyase cos1mQ\cos\frac1m \notin \Q. ¿Dónde se rompe el mismo intento para cosab\cos\frac ab con b>1b > 1 general? (Identifíquese el denominador que ya no se elimina.)

Parte V — Más fino y más oscuro: la forma con igualdad, y una función que engaña a Taylor.

  1. (Taylor–Lagrange, forma con igualdad) Sea ff n+1n + 1 veces derivable entre aa y xx. Defínase g(t)=f(x)k=0nf(k)(t)k!(xt)kA(xt)n+1(n+1)!g(t) = f(x) - \sum_{k=0}^{n} \frac{f^{(k)}(t)}{k!}(x - t)^k - A\,\frac{(x-t)^{n+1}}{(n+1)!} con la constante AA elegida de modo que g(a)=0g(a) = 0. Calcúlese gg' (la suma telescopa), aplíquese Rolle en [a,x]\intcc{a}{x} y conclúyase que existe cc estrictamente entre aa y xx con

    f(x)=k=0nf(k)(a)k!(xa)k+f(n+1)(c)(n+1)!(xa)n+1.f(x) = \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}(x-a)^k + \frac{f^{(n+1)}(c)}{(n+1)!}\,(x-a)^{n+1} .
  2. Dividendo de la igualdad: para x>0x > 0 véase que

    ex>1+x+x22!++xnn!\eu^x > 1 + x + \frac{x^2}{2!} + \dots + \frac{x^n}{n!}

    (estrictamente, para todo nn), y localícese dónde se invierte la desigualdad para x<0x < 0 según la paridad de nn.

  3. Compárense los restos sobre sin(0.5)\sin(0.5) al orden 33: Young da solo o(x3)o(x^3) (ningún número); Lagrange da sin0.5(0.50.536)0.55120=2.61104\abs{\sin 0.5 - (0.5 - \frac{0.5^3}{6})} \leq \frac{0.5^5}{120} = 2.61\cdot10^{-4}; y la estimación alternada da la misma cota más la información del signo sin0.5>0.50.536\sin 0.5 > 0.5 - \frac{0.5^3}{6}. Compárese con el error verdadero 2.591042.59\cdot10^{-4}: la cota casi se alcanza. ¿A qué herramienta se echaría mano, y cuándo?
  4. (La función plana) Sea f(x)=e1/x2f(x) = \eu^{-1/x^2} para x0x \neq 0 y f(0)=0f(0) = 0. Véase que ff es continua en 00, que f(0)=0f'(0) = 0 y, más en general —demostrando que toda derivada tiene la forma f(k)(x)=Pk(1x)e1/x2f^{(k)}(x) = P_k\bigl(\frac1x\bigr)\eu^{-1/x^2} para un polinomio PkP_k (inducción)—, que f(k)(0)=0f^{(k)}(0) = 0 para todo kk (comparación de crecimientos, Proposición 4.6). Conclúyase: todos los polinomios de Taylor de ff en 00 se anulan y, sin embargo, f(x)>0f(x) > 0 para x0x \neq 0: Taylor–Young vale a todos los órdenes y no dice nada de ff lejos de 00. Los desarrollos describen gérmenes, no funciones.

Parte VI — Síntesis.

  1. Hágase correr la trampa una vez más, sobre e1=k(1)kk!\eu^{-1} = \sum_k \frac{(-1)^k}{k!}: multiplíquese el encuadre alternado estricto por n!n! y conclúyase e1Q\eu^{-1} \notin \Q, y por tanto eQ\eu \notin \Q — la tercera demostración de este hecho en el volumen. Enumérense las tres (sucesiones adyacentes, Ejercicio 11.9; integrales, Problema 15.1; sumas alternadas, aquí) y lo que necesitó cada una.
  2. La letra pequeña: la monotonía no es decorativa. Sea bk=1kb_k = \frac1k para kk impar y bk=1k2b_k = \frac{1}{k^2} para kk par: los bkb_k son positivos y tienden a 00 y, sin embargo, las sumas parciales de (1)kbk\sum (-1)^k b_k divergen a -\infty. Demuéstrese (pártase la suma parcial en la parte par, acotada mediante el Ejemplo 11.22, y la parte impar, que domina la mitad de la serie armónica, Ejercicio 11.5), y dígase exactamente qué paso de la pregunta 1 usó la monotonía.
  3. Compruébese la identidad más sencilla de Euler arctan12+arctan13=π4\arctan\frac12 + \arctan\frac13 = \frac\pi4 mediante (2+i)(3+i)=5(1+i)(2 + \iu)(3 + \iu) = 5(1 + \iu), estímense los términos necesarios para seis decimales de π\pi por este camino (basta n=10n = 10) y colóquese entre Leibniz y Machin en la clasificación de la pregunta 13.
  4. Síntesis, una frase para cada punto: (i) enúnciese la estimación alternada y sus dos productos (cota y signo); (ii) por qué las identidades finitas exactas con restos explícitos ganan a los enunciados con límites para el cálculo numérico certificado; (iii) inventario del problema (π\pi con error 10710^{-7} a mano, ln2\ln 2 con diez decimales, cuatro demostraciones de irracionalidad, un teorema con igualdad, un ejemplo de aviso); (iv) cuáles de estos hilos recogerá el Capítulo 17 (el criterio de las series alternadas, la convergencia absoluta frente a la condicional y el drama de las reordenaciones de su problema del fin de semana).
Solución

Solución de Problema 16.1.

1. S2n+1S2n1=a2na2n+10S_{2n+1} - S_{2n-1} = a_{2n} - a_{2n+1} \geq 0 y S2n+2S2n=a2n+2a2n+10S_{2n+2} - S_{2n} = a_{2n+2} - a_{2n+1} \leq 0, mientras que S2nS2n+1=a2n+10S_{2n} - S_{2n+1} = a_{2n+1} \to 0: las sucesiones (S2n+1)(S_{2n+1}) y (S2n)(S_{2n}) son adyacentes y convergen a un mismo SS (Teorema 11.11) con S2n+1SS2nS_{2n+1} \leq S \leq S_{2n}. Para nn par: Sn+1SSnS_{n+1} \leq S \leq S_n da an+1SSn0-a_{n+1} \leq S - S_n \leq 0; y para nn impar: 0SSnan+10 \leq S - S_n \leq a_{n+1}. En los dos casos, SSnan+1\abs{S - S_n} \leq a_{n+1} y SSnS - S_n tiene el signo de (1)n+1(-1)^{n+1}, el primer término omitido. El decrecimiento estricto hace estricta toda desigualdad anterior, en particular 0<SSn<an+10 < \abs{S - S_n} < a_{n+1}.

2. ak=1k!a_k = \frac{1}{k!} decrece estrictamente hacia 00: se aplica la pregunta 1. Taylor–Lagrange (Teorema 16.7) para exp\exp entre 1-1 y 00: e1Tn1(n+1)!\abs{\eu^{-1} - T_n} \leq \frac{1}{(n+1)!} (la derivada et\eu^t es 1\leq 1 ahí), luego Tne1T_n \to \eu^{-1}, y el límite SS de la pregunta 1 es e1\eu^{-1}, con las cotas estrictas 0<e1Tn<1(n+1)!0 < \abs{\eu^{-1} - T_n} < \frac{1}{(n+1)!}.

3. Integrando la identidad en [0,1]\intcc{0}{1}: el miembro izquierdo es arctan1=π4\arctan 1 = \frac\pi4 (teorema fundamental), el término kk-ésimo da (1)k2k+1\frac{(-1)^k}{2k+1}, y

ρn=(1)n+101t2n+21+t2 ⁣dt,ρn01t2n+2 ⁣dt=12n+3.\rho_n = (-1)^{n+1}\int_0^1 \frac{t^{2n+2}}{1+t^2}\dd t, \qquad \abs{\rho_n} \leq \int_0^1 t^{2n+2}\dd t = \frac{1}{2n+3} .

4. El error sobre π\pi es 4ρn42n+34\abs{\rho_n} \leq \frac{4}{2n+3}: bajar de 10610^{-6} exige 2n+3>41062n + 3 > 4\cdot10^6, unos dos millones de términos. Mientras tanto, 4S4=4(113+1517+19)=4×0.834921=3.3396834S_4 = 4\bigl(1 - \frac13 + \frac15 - \frac17 + \frac19\bigr) = 4 \times 0.834921 = 3.339683, a casi 0.20.2 de π\pi: cinco términos y ni siquiera una cifra.

5. Integrando en [0,1]\intcc{0}{1}: ln2=k=1n(1)k1k+(1)nRn\ln 2 = \sum_{k=1}^{n} \frac{(-1)^{k-1}}{k} + (-1)^n R_n con Rn=01tn1+t ⁣dtR_n = \int_0^1 \frac{t^n}{1+t}\dd t; y de 1211+t1\frac12 \leq \frac{1}{1+t} \leq 1: 12(n+1)Rn1n+1\frac{1}{2(n+1)} \leq R_n \leq \frac{1}{n+1}. El error queda atrapado entre dos múltiplos de 1n\frac1n: seis decimales cuestan alrededor de un millón de términos.

6. Integrando de 00 a xx: 12ln1+x1x=k=0nx2k+12k+1+0xt2n+21t2 ⁣dt\frac12\ln\frac{1 + x}{1 - x} = \sum_{k=0}^{n} \frac{x^{2k+1}}{2k+1} + \int_0^x \frac{t^{2n+2}}{1-t^2}\dd t. En x=13x = \frac13: 1+1/311/3=2\frac{1 + 1/3}{1 - 1/3} = 2, y en [0,13]\intcc{0}{\frac13} se tiene 11t298\frac{1}{1 - t^2} \leq \frac98:

ln2=2k=0n(1/3)2k+12k+1+ρ~n,0<ρ~n94(1/3)2n+32n+3:\ln 2 = 2\sum_{k=0}^{n} \frac{(1/3)^{2k+1}}{2k+1} + \tilde\rho_n, \qquad 0 < \tilde\rho_n \leq \frac{9}{4}\cdot \frac{(1/3)^{2n+3}}{2n+3} :

cada término adicional divide el error por unos 99.

7. Inducción: para n=1n = 1: 112=12=H2H11 - \frac12 = \frac12 = H_2 - H_1. Paso:

H2n+2Hn+1=(H2nHn)+12n+1+12n+21n+1=(H2nHn)+12n+112n+2,H_{2n+2} - H_{n+1} = (H_{2n} - H_n) + \frac{1}{2n+1} + \frac{1}{2n+2} - \frac{1}{n+1} = (H_{2n} - H_n) + \frac{1}{2n+1} - \frac{1}{2n+2},

que es exactamente el incremento de la suma alternada. Y H2nHn=k=1n1n+kH_{2n} - H_n = \sum_{k=1}^{n}\frac{1}{n+k} es la suma de Riemann del Ejemplo 15.21, que converge a ln2\ln 2: las sumas parciales pares del camino 1 son las sumas de Riemann del camino 3.

8. k=16(1)k1k=0.61667\sum_{k=1}^{6}\frac{(-1)^{k-1}}{k} = 0.61667, con error 0.07650.0765; el camino 2 en n=5n = 5 da 0.69314710.6931471 con error 94(1/3)1313=1.1107\leq \frac94\cdot\frac{(1/3)^{13}}{13} = 1.1\cdot10^{-7}. La razón: el camino 1 evalúa la serie del logaritmo en el punto de la frontera x=1x = 1, donde los términos decrecen como 1k\frac1k; el camino 2 evalúa en x=13x = \frac13, bien dentro, donde cada término lleva un factor fresco 19\frac19.

9. Diez decimales: se quiere ρ~n51011\tilde\rho_n \leq 5\cdot10^{-11}. En n=10n = 10: 94(1/3)2323=941.061011231.01012<51011\frac94 \cdot \frac{(1/3)^{23}}{23} = \frac94\cdot\frac{1.06\cdot10^{-11}}{23} \approx 1.0\cdot10^{-12} < 5\cdot10^{-11}: bastan once términos.

10. (5+i)2=24+10i(5+\iu)^2 = 24 + 10\iu, y después (5+i)4=(24+10i)2=476+480i(5+\iu)^4 = (24 + 10\iu)^2 = 476 + 480\iu; y 2(1+i)(239+i)=2(238+240i)=476+480i2(1+\iu)(239+\iu) = 2(238 + 240\iu) = 476 + 480\iu: iguales. Argumentos: arg(5+i)=arctan15\arg(5 + \iu) = \arctan\frac15, luego el miembro izquierdo tiene argumento 4arctan150.79(0,π)4\arctan\frac15 \approx 0.79 \in \intoo{0}{\pi}; y el miembro derecho tiene argumento π4+arctan1239(0,π)\frac\pi4 + \arctan\frac{1}{239} \in \intoo{0}{\pi}. Dos complejos iguales con argumentos en un mismo intervalo de longitud <2π< 2\pi:

4arctan15=π4+arctan1239,4\arctan\frac15 = \frac\pi4 + \arctan\frac{1}{239} ,

que es la fórmula de Machin.

11. Intégrese 11+t2=k=0n(1)kt2k+(1)n+1t2n+21+t2\frac{1}{1+t^2} = \sum_{k=0}^n (-1)^k t^{2k} + \frac{(-1)^{n+1}t^{2n+2}}{1+t^2} de 00 a xx:

arctanx=k=0n(1)kx2k+12k+1+rn(x),rn(x)0xt2n+2 ⁣dt=x2n+32n+3.\arctan x = \sum_{k=0}^{n}\frac{(-1)^k x^{2k+1}}{2k+1} + r_n(x), \qquad \abs{r_n(x)} \leq \int_0^x t^{2n+2}\dd t = \frac{x^{2n+3}}{2n+3} .

12. Errores: 16(1/5)1111=3.010816\,\frac{(1/5)^{11}}{11} = 3.0\cdot 10^{-8} y 4(1/239)55110124\,\frac{(1/239)^5}{5} \approx 1\cdot10^{-12}: en total <5108< 5\cdot10^{-8}. La suma anterior vale 3.141592683.14159268\dots, luego π=3.1415926\pi = 3.1415926\dots certificado con error 51085\cdot10^{-8}: siete decimales a partir de seis términos (cinco en 15\frac15 y dos en 1239\frac1{239}, contando con generosidad).

13. Leibniz: error 1n\sim \frac1n, de modo que cada cifra nueva multiplica el trabajo por diez. Dalzell (Problema 15.1, pregunta 22): error 415m4^{1-5m}, unas tres cifras por paso, y cada paso con un polinomio más pesado. Machin: razón de error 125\frac{1}{25} por término, unas 1.41.4 cifras por término, y cada término una división. La moraleja: el resto de un desarrollo de tipo geométrico escala como x2nx^{2n}, así que hacer xx pequeño compra cifras a un coste fijo por término — la identidad compleja de Machin es precisamente una máquina de encoger xx.

14. ak=1(2k)!a_k = \frac{1}{(2k)!} decrece estrictamente hacia 00; por la pregunta 1 y la Proposición 16.11 (cota de Lagrange como en la pregunta 2), kn(1)k(2k)!cos1\sum_{k \leq n}\frac{(-1)^k}{(2k)!} \to \cos 1 con el encuadre estricto 0<cos1Sn<1(2n+2)!0 < \bigl|\cos 1 - S'_n\bigr| < \frac{1}{(2n+2)!}. Supóngase cos1=pq\cos 1 = \frac pq y tómese 2nq2n \geq q: entonces (2n)!Sn=kn(1)k(2n)!(2k)!Z(2n)!\,S'_n = \sum_{k\leq n} (-1)^k \frac{(2n)!}{(2k)!} \in \Z y (2n)!pqZ(2n)!\,\frac pq \in \Z, mientras que

0<(2n)!cos1(2n)!Sn<(2n)!(2n+2)!=1(2n+1)(2n+2)<1:0 < \Bigl|(2n)!\cos 1 - (2n)!S'_n\Bigr| < \frac{(2n)!}{(2n+2)!} = \frac{1}{(2n+1)(2n+2)} < 1 :

un entero no nulo de valor absoluto <1< 1. Contradicción: cos1Q\cos 1 \notin \Q.

15. Idénticamente con ak=1(2k+1)!a_k = \frac{1}{(2k+1)!}, multiplicando por (2n+1)!(2n+1)! con 2n+1q2n + 1 \geq q: sin1Q\sin 1 \notin \Q. Y, sin embargo, cos21+sin21=1Q\cos^2 1 + \sin^2 1 = 1 \in \Q: los productos y las sumas de irracionales pueden ser racionales — la irracionalidad no pasa gratis por ninguna operación algebraica.

16. Cota de la cola: para m>nm > n,

k=n+1m1(2k)!1(2n+2)!(1+12+14+)2(2n+2)!,\sum_{k=n+1}^{m} \frac{1}{(2k)!} \leq \frac{1}{(2n+2)!}\Bigl(1 + \frac12 + \frac14 + \dots\Bigr) \leq \frac{2}{(2n+2)!} ,

puesto que cada razón sucesiva es 1(2k+1)(2k+2)12\frac{1}{(2k+1)(2k+2)} \leq \frac12; y la cola es positiva (lo es su primer término). Así, 0<cosh1kn1(2k)!<2(2n+2)!0 < \cosh 1 - \sum_{k\leq n}\frac{1}{(2k)!} < \frac{2}{(2n+2)!}, y multiplicando por (2n)!(2n)! con 2nq2n \geq q se vuelve a atrapar un entero no nulo en (0,1)\intoo{0}{1}: cosh1Q\cosh 1 \notin \Q.

17. cos1m=k(1)km2k(2k)!\cos\frac1m = \sum_k \frac{(-1)^k}{m^{2k}(2k)!}: los términos decrecen estrictamente hacia cero, y m2n(2n)!1m2k(2k)!=m2(nk)(2n)!(2k)!Zm^{2n}(2n)!\cdot\frac{1}{m^{2k}(2k)!} = m^{2(n-k)}\frac{(2n)!}{(2k)!} \in \Z para knk \leq n. Si cos1m=pq\cos\frac1m = \frac pq, multiplíquese el encuadre estricto por qm2n(2n)!q\,m^{2n}(2n)!: el error queda acotado por qm2(2n+1)(2n+2)<1\frac{q}{m^2(2n+1)(2n+2)} < 1 para nn grande: contradicción. Para ab\frac ab con a2a \geq 2: quitar denominadores multiplica la cola por b2n(2n)!b^{2n}(2n)!, pero el primer término omitido es a2n+2b2n+2(2n+2)!\frac{a^{2n+2}}{b^{2n+2}(2n+2)!}, y el producto a2n+2b2(2n+1)(2n+2)\frac{a^{2n+2}}{b^2(2n+1)(2n+2)} explota: el numerador a2n+2a^{2n+2} ya no se elimina y la trampa se atasca. (El resultado sigue siendo cierto — por maquinaria al estilo de Niven, no por esta.)

18. En t=xt = x todos los términos de gg se anulan salvo f(x)f(x)=0f(x) - f(x) = 0: g(x)=0g(x) = 0; y AA se elige de modo que g(a)=0g(a) = 0. Derivando, la suma telescopa:

g(t)=f(n+1)(t)n!(xt)n+A(xt)nn!=(xt)nn!(Af(n+1)(t)).g'(t) = -\frac{f^{(n+1)}(t)}{n!}(x - t)^n + A\,\frac{(x-t)^n}{n!} = \frac{(x-t)^n}{n!}\bigl(A - f^{(n+1)}(t)\bigr) .

Rolle en el segmento de aa a xx da cc estrictamente entre ambos con g(c)=0g'(c) = 0; y como (xc)n0(x - c)^n \neq 0: A=f(n+1)(c)A = f^{(n+1)}(c). Desplegando g(a)=0g(a) = 0 se obtiene la igualdad de Taylor con resto f(n+1)(c)(n+1)!(xa)n+1\frac{f^{(n+1)}(c)}{(n+1)!}(x-a)^{n+1}.

19. Para x>0x > 0, el resto es ec(n+1)!xn+1>0\frac{\eu^{c}} {(n+1)!}x^{n+1} > 0: la exponencial supera estrictamente a cada uno de sus polinomios de Taylor, a todo orden. Para x<0x < 0, el signo del resto es el de xn+1x^{n+1}: ex\eu^x queda por encima del polinomio si nn es impar y por debajo si nn es par — alternando de lado, como ya muestran las gráficas de 1+x1 + x y de 1+x+x221 + x + \frac{x^2}{2} frente a ex\eu^x.

20. Error verdadero: sin0.50.4791667=2.59104\sin 0.5 - 0.4791667 = 2.59\cdot10^{-4}, frente a la cota 0.55120=2.60104\frac{0.5^5}{120} = 2.60\cdot10^{-4}: casi se alcanza (el término siguiente domina la cola). Young: para límites y análisis local, donde no hace falta ninguna constante. Lagrange: para decimales certificados. Alternada: cuando es aplicable, la misma cota más la dirección del error — la mejor de las tres, pero la más rara.

21. Continuidad en 00: con u=1x2+u = \frac{1}{x^2} \to +\infty, f(x)=eu0=f(0)f(x) = \eu^{-u} \to 0 = f(0). Derivada en 00: f(h)h=ueu0\bigl|\frac{f(h)}{h}\bigr| = \sqrt u\,\eu^{-u} \to 0 (Proposición 4.6): f(0)=0f'(0) = 0. Para x0x \neq 0, f(x)=2x3e1/x2f'(x) = \frac{2}{x^3}\eu^{-1/x^2}: la forma P1(1x)e1/x2P_1\bigl(\frac1x\bigr)\eu^{-1/x^2} con P1(X)=2X3P_1(X) = 2X^3; y, por inducción, derivar Pk(1x)e1/x2P_k(\frac1x)\eu^{-1/x^2} da Pk+1(X)=2X3Pk(X)X2Pk(X)P_{k+1}(X) = 2X^3 P_k(X) - X^2 P_k'(X), un polinomio. Entonces

f(k)(h)0h=vPk(v)ev2v=1/h0\frac{f^{(k)}(h) - 0}{h} = v\,P_k(v)\,\eu^{-v^2} \Big|_{v = 1/h} \longrightarrow 0

(polinomio contra ev2\eu^{-v^2}, comparación de crecimientos en ±\pm\infty): por inducción, f(k)(0)=0f^{(k)}(0) = 0 para todo kk. Todos los polinomios de Taylor de ff en 00 se anulan y, sin embargo, f>0f > 0 fuera de 00: Taylor–Young es exacto a todos los órdenes y ciego más allá del germen. Un desarrollo es información local, nada más.

22. Por la pregunta 2, 0<e1Tn<1(n+1)!0 < \abs{\eu^{-1} - T_n} < \frac{1}{(n+1)!}, estrictamente. Si e1=pq\eu^{-1} = \frac pq, tómese nqn \geq q y multiplíquese por n!n!: n!TnZn!\,T_n \in \Z y n!pqZn!\frac pq \in \Z, luego un entero no nulo tiene valor absoluto <n!(n+1)!=1n+1<1< \frac{n!}{(n+1)!} = \frac{1}{n+1} < 1: contradicción. Por tanto, e1Q\eu^{-1} \notin \Q, y e=1e1\eu = \frac{1}{\eu^{-1}} también es irracional. Las tres demostraciones: sucesiones adyacentes apretando q!eq!\,\eu (Ejercicio 11.9); la recurrencia integral An=enAn1A_n = \eu - nA_{n-1} (Problema 15.1); y el encuadre alternado (aquí). Una trampa, tres certificados de pequeñez.

23. Agrúpense las sumas parciales por parejas: k=12n(1)kbk=EnOn\sum_{k=1}^{2n} (-1)^k b_k = E_n - O_n con En=j=1n14j2E_n = \sum_{j=1}^{n} \frac{1}{4j^2}, acotada (por la cota telescópica del Ejemplo 11.22, En12E_n \leq \frac12), y On=j=1n12j112Hn+O_n = \sum_{j=1}^{n}\frac{1}{2j-1} \geq \frac12 H_n \to +\infty (Ejercicio 11.5): las sumas parciales tienden a -\infty. La monotonía se usó en la pregunta 1 exactamente donde S2n+1S2n1=a2na2n+1S_{2n+1} - S_{2n-1} = a_{2n} - a_{2n+1} necesitaba un signo: sin decrecimiento, las subsucesiones par e impar no tienen por qué ser monótonas, y la adyacencia se derrumba.

24. (2+i)(3+i)=5+5i=5(1+i)(2+\iu)(3+\iu) = 5 + 5\iu = 5(1+\iu); tomando argumentos (todos en (0,π2)\intoo{0}{\frac\pi2}): arctan12+arctan13=π4\arctan\frac12 + \arctan\frac13 = \frac\pi4. Coste en términos para seis decimales: error 4((1/2)2n+32n+3+(1/3)2n+32n+3)\leq 4\bigl(\frac{(1/2)^{2n+3}}{2n+3} + \frac{(1/3)^{2n+3}}{2n+3}\bigr), que en n=10n = 10 es 2108<5107\approx 2\cdot10^{-8} < 5\cdot10^{-7}: once términos. Clasificación: mejor que Leibniz por un margen exponencial, y por detrás de Machin (cuyo punto dominante 15\frac15 es menor que 12\frac12): unas 0.60.6 cifras por término frente a las 1.41.4 de Machin.

25. (i) Para ak0a_k \to 0 decreciente, las sumas parciales alternadas convergen con SSnan+1\abs{S - S_n} \leq a_{n+1} y el error lleva el signo del primer término omitido. (ii) Una identidad finita con resto explícito se puede evaluar y acotar en un punto elegido, mientras que un enunciado con límite solo promete proximidad a la larga — certificar exige lo primero. (iii) Extraído: π\pi con error 51085\cdot10^{-8} por Machin, ln2\ln 2 con diez decimales por la serie en 13\frac13, la irracionalidad de cos1\cos 1, sin1\sin 1, cosh1\cosh 1, cos1m\cos\frac1m y e1\eu^{-1}, la igualdad de Taylor–Lagrange y el aviso de la función plana. (iv) Capítulo 17 eleva la pregunta 1 al criterio de las series alternadas, separa la convergencia absoluta de la condicional y, en su problema del fin de semana, escenifica el drama de las reordenaciones, del que la serie armónica alternada del camino 1 es el testigo estrella.