← Índice de temas
Módulo II · Funciones derivables · Tema 2.4

Método de Newton. Método de punto fijo

Dos formas de resolver ecuaciones que no se pueden resolver «con fórmulas»: seguir la tangente hasta que corte el eje (Newton) y aplicar una y otra vez la misma función hasta que el resultado deje de moverse (punto fijo). Aprenderás a calcular, a justificar que el método funciona y a decidir cuántas iteraciones hacen falta.

Libro oficial: Cálculo para Ingenieros, §2.4 (pp. 86-93; Ej. 2.18 en p. 94)Ejercicios: E Ej. 2.34–2.38 (pp. 72-76)⏱ ~5 h de estudio
U oficial S Stewart L Larson E ejercicios P propio X exámenes — pasa el ratón (o toca) por cada etiqueta para ver libro y página; en pantallas anchas aparecen en el margen derecho. Decimales con punto, como en E.

0Antes de empezar U p.86

La idea del tema en una frase

Casi ninguna ecuación $f(x)=0$ se puede resolver de forma exacta. Este tema enseña dos algoritmos que fabrican una sucesión $x_0,x_1,x_2,\dots$ cuyos términos se acercan a la solución. Newton es muy rápido pero no tiene garantía; el punto fijo es más lento pero, si se comprueban unas hipótesis, sí la tiene (y además te dice cuántas iteraciones bastan).

U p.86

Qué necesitas saber ya

Mapa del tema

ApartadoQué contieneDónde está
1 Newtonidea de la tangente, fórmula, criterio de parada, cuándo fallaU pp.86-90 · L pp.225-228 · S pp.338-341
2 Punto fijodefinición, iteración, Teoremas 2.2 y 2.3, Ejemplo 2.17U pp.90-94
3 ¿Cuántas iteraciones?tabla para decidir según el casoU pp.54, 89-92 · propio
4 Comparaciónbisección / punto fijo / Newtonpropio
5 EjerciciosE 2.34-2.38 + 6 ejemplos propiosE pp.72-76
6 Chuletafórmulas, lista de comprobación, errores

Qué se pide en los ejercicios

Los cinco ejercicios de E son de este tipo: (a) hacer $n$ iteraciones de Newton para aproximar un número o una raíz (2.34, 2.35, 2.36); (b) comprobar las hipótesis del Teorema 2.2 y hacer iteraciones del punto fijo (2.38); (c) calcular con el Teorema 2.3 cuántas iteraciones garantizan un error dado (2.37). E pp.72-76

Lo que U no dice (y aquí se completa)

Comprobado en el texto: U no da ninguna condición de convergencia para Newton (solo una figura al margen con los fallos, p. 89), no demuestra los Teoremas 2.2 y 2.3 y no cuantifica la velocidad de convergencia (lineal o cuadrática; solo dice cualitativamente que Newton es «rápido», pp. 86-87). Se completa con Larson §3.8 y Stewart §4.8 (marcados L S) y con explicaciones propias (P). Lo que cae en examen es lo que hay en U; el resto es para entender por qué.

1Método de Newton U pp.86-90

Idea intuitiva: sustituir la curva por su tangente

Estás en un punto $x_0$ que no es solución ($f(x_0)\ne0$). No sabes dónde corta la curva al eje, pero sí sabes dónde lo corta una recta. Así que sustituyes la curva por su recta tangente en $(x_0,f(x_0))$, que se le parece mucho cerca de $x_0$, y tomas como nueva aproximación el punto $x_1$ donde esa tangente corta al eje $x$. Repites desde $x_1$. Es como bajar una cuesta siguiendo su pendiente hasta llegar a nivel del suelo: si la cuesta es casi recta, llegas casi exacto.

U pp.86-87 S p.339 L p.225
x₀ = 3 x₁ = 7/3 x₂ ≈ 2.2381 raíz √5 ≈ 2.2361 (círculo vacío) f(x) = x² − 5 tangente en x₀ tangente en x₁
Newton para $f(x)=x^2-5$ con $x_0=3$: la tangente en $(x_0,f(x_0))$ corta el eje en $x_1=x_0-\frac{f(x_0)}{f'(x_0)}=\frac73$; repetir desde $x_1$ da $x_2\approx2.2381$, a $2\times10^{-3}$ de $\sqrt5$. Ya en el segundo paso la tangente y la curva son casi indistinguibles.

De la idea a la fórmula

Método de Newton (algoritmo)

Hipótesis. $f$ derivable en un intervalo que contiene a la raíz $\tilde x$ que se busca (y que, en cada paso, $f'(x_n)\neq0$).

Datos. Una aproximación inicial $x_0$ cercana a la raíz.

Iteración. $$x_{n+1}=x_n-\frac{f(x_n)}{f'(x_n)},\qquad n=0,1,2,\dots$$

(La derivada $f'(x_n)$ es la pendiente de la tangente en $(x_n,f(x_n))$, y por eso «da» el paso siguiente.) También se llama método de Newton-Raphson o de Newton-Fourier (U p.86).

U p.88 L p.225 S p.339
Deducción (por qué esa fórmula)
  1. La recta tangente a la gráfica de $f$ en $(x_0,f(x_0))$ tiene pendiente $f'(x_0)$: $\;y-f(x_0)=f'(x_0)(x-x_0)$.
  2. Corta al eje $x$ donde $y=0$: $\;-f(x_0)=f'(x_0)(x_1-x_0)$.
  3. Si $f'(x_0)\neq0$ podemos despejar: $\;x_1=x_0-\dfrac{f(x_0)}{f'(x_0)}$. Si $f'(x_0)=0$ la tangente es horizontal, no corta al eje y hay que cambiar $x_0$.
U pp.87-88
Ejemplo 2.15 de U: $\sqrt[4]5$ como raíz de $x^4-5$
  1. $f(x)=x^4-5$, $f'(x)=4x^3$; iteración $x_{n+1}=x_n-\dfrac{x_n^4-5}{4x_n^3}$.
  2. Con $x_0=2$: $f(2)=11$ (lejos de ser solución). $x_1=2-\frac{11}{32}=\frac{53}{32}\approx1.6563$.
  3. Siguiendo: $x_2\approx1.5173$, $x_3\approx1.4958$. El error $|x_n-\sqrt[4]5|$ va de $0.505$ a $0.161$, $0.0220$ y $0.00045$ (con $\sqrt[4]5=1.495349\ldots$).

Fíjate en el patrón de los errores: $10^{-1}\to10^{-2}\to10^{-4}$. Cada error es más o menos el cuadrado del anterior (ver el apartado de velocidad, más abajo).

U p.88
Error típico: aplicar Newton sin haber escrito «$f(x)=0$»

Newton resuelve ecuaciones de la forma $f(x)=0$. Para $\sqrt5$ hay que decir «$\sqrt5$ es la raíz positiva de $x^2-5=0$» y tomar $f(x)=x^2-5$. Aplicar la fórmula a «$x=\sqrt5$» no significa nada. Y si la ecuación es $e^{-x}=x$, hay que pasarla a $e^{-x}-x=0$ antes de derivar.

Si converge, converge a una raíz

La sucesión de Newton, si converge, converge a una solución

U lo comprueba para $x^4-5$ (Ejemplo 2.16, tomando límites en la iteración). En general: si $f$ y $f'$ son continuas, $x_n\to\tilde x$ y $f'(\tilde x)\neq0$, al tomar límites en $x_{n+1}=x_n-f(x_n)/f'(x_n)$ queda $\tilde x=\tilde x-f(\tilde x)/f'(\tilde x)$, es decir, $f(\tilde x)=0$. Ojo: esto no dice que la sucesión converja; solo que, si lo hace, el límite es una raíz.

U p.89 L p.227 P

Cuándo parar: los errores aproximados

Error absoluto y relativo aproximados; criterios de parada

No conocemos la solución exacta $\tilde x$ (por eso usamos un método numérico), así que medimos cuánto cambia la sucesión en el último paso: $$\varepsilon_a=|x_n-x_{n-1}|,\qquad \varepsilon_r=\frac{|x_n-x_{n-1}|}{|x_n|}.$$ Se para cuando $\varepsilon_r$ (o $\varepsilon_a$) es menor que una tolerancia dada. Otro criterio, muy conveniente al programar, es un número máximo de iteraciones fijado de antemano, para no entrar en un proceso infinito.

Se toman valores absolutos porque normalmente da igual si el error es por defecto o por exceso (no siempre: U cita la velocidad mínima de despegue de un avión, donde un error por defecto significa que no despega y uno por exceso solo gasta más combustible, p.89); se divide entre $|x_n|$ porque un error de $10$ céntimos no pesa igual en una caña que en un coche. Rigurosamente habría que usar el error relativo exacto (con $\tilde x$), pero es imposible, y se supone que $\tilde x$ está cerca de $x_n$.

U pp.89-90 L p.225 S p.340
Error típico: creer que $\varepsilon_r$ es el error verdadero

$\varepsilon_r$ solo mide «cuánto me he movido». Si el método va muy despacio o se ha atascado, $\varepsilon_r$ puede ser pequeño y estar la solución lejos. En Newton, cuando funciona bien, suele sobreestimar el error real (Ejercicio 2.36). Y si se redondea a pocos decimales, puede salir $\varepsilon_r=0$ sin que el error sea nulo: solo es menor de lo que se ve.

¿Es infalible? No

Cuándo falla Newton
  • $f'(x_n)=0$ para algún $n$: no se puede dividir (tangente horizontal). Solución: elegir otro $x_0$.
  • Ciclos: la sucesión salta entre dos puntos y no converge (figura al margen de U).
  • Mal punto de partida: la tangente puede llevarte lejos de la raíz o fuera del dominio de $f$ (S, figura 4).

Ejemplo de ciclo (L ej. 22, p.229). Con $f(x)=x^3-2x+2$ y $x_0=0$ salen $0,1,0,1,\dots$ y no converge. Detalle en el ejemplo propio P3.

Ejemplo de divergencia (L Ejemplo 3, p.227). $f(x)=\sqrt[3]{x}$ (no derivable en su raíz $0$). Aquí $x_{n+1}=x_n-\frac{x_n^{1/3}}{\frac13x_n^{-2/3}}=x_n-3x_n=-2x_n$: con $x_0=0.1$ salen $-0.2,\ 0.4,\ -0.8,\dots$ que se alejan del $0$ en cada paso.

Moraleja: antes de elegir $x_0$, localiza la raíz (Bolzano o una gráfica).

U p.89 L pp.227,229 S pp.340,343
Una condición suficiente de convergencia (no está en U)

Hipótesis. $f$ con derivada segunda y $f'\neq0$ en un intervalo abierto $I$ que contiene al cero $c$, y punto inicial $x_0\in I$ tal que las iteraciones permanecen en $I$ (por ejemplo, si $I$ es simétrico respecto de $c$).

Conclusión (Larson). Si $\left|\dfrac{f(x)\,f''(x)}{[f'(x)]^2}\right|<1$ para todo $x\in I$, el método de Newton converge al cero. Para $f(x)=x^2-2$ el cociente es $\frac{x^2-2}{2x^2}=\frac12-\frac1{x^2}$, que en $(1,3)$ toma valores en $(-\frac12,\frac7{18})$: su valor absoluto es $<1$ y converge. Para $f(x)=\sqrt[3]x$ el cociente vale $-2$ en todo punto ($|{\cdot}|=2$): no se cumple, y en efecto diverge.

Es solo suficiente; U no la pide ni la usa, y Larson no la demuestra. Sirve como criterio de sentido común: Newton va bien si $f$ no se «curva» demasiado respecto a lo empinada que es.

L p.228

Velocidad: por qué Newton es tan rápido (a veces)

Los decimales correctos se duplican

Con una raíz simple ($f'(\tilde x)\ne0$), el error nuevo es aproximadamente el cuadrado del anterior: $10^{-2}\to10^{-4}\to10^{-8}\to10^{-16}$ (convergencia cuadrática). Con una raíz doble ($f'(\tilde x)=0$, la curva «toca» el eje sin cruzarlo) el error solo se reduce a la mitad en cada paso, y Newton no es más rápido que la bisección (Ejercicio 2.35). La razón se ve con álgebra pura en el paso 6 de la solución del Ejercicio 2.34 (explicación propia, no viene en E). U no cuantifica la velocidad.

P
Truco: simplifica $x-\frac{f}{f'}$ antes de iterar

Casi siempre sale una expresión más sencilla que se calcula mejor: para $x^2-a$ es $\tfrac12\!\left(x+\tfrac ax\right)$ (regla de Herón, la media entre $x$ y $a/x$); para $x^m-a$ es $\frac{(m-1)x^m+a}{m\,x^{m-1}}$; para $e^{-x}-x$ es $\frac{x+1}{1+e^{x}}$.

L p.229 S p.342

2Método de punto fijo U pp.90-94

Idea intuitiva: la calculadora que pulsa «cos» sin parar

Pon la calculadora en radianes, escribe cualquier número y pulsa la tecla cos una y otra vez. Los resultados se estabilizan en $0.7390\ldots$ ¿Por qué? Ese número $\tilde x$ cumple $\cos\tilde x=\tilde x$: la función lo deja quieto. Un punto que una función no mueve es un punto fijo. Si repetir la función te acerca cada vez más a un punto, ese punto es la solución de «$x=g(x)$».

P
Definición 2.4 (punto fijo)

Dada una función $g$, $\tilde x$ es un punto fijo de $g$ si $g(\tilde x)=\tilde x$. Gráficamente: los puntos donde la gráfica de $g$ corta a la recta $y=x$. Calcular puntos fijos = resolver la ecuación $x=g(x)$.

U p.90
Método de punto fijo

Datos. $g$ continua y una primera aproximación $x_0$. Iteración (ecuación (2.3) de U): $$x_{n+1}=g(x_n),\qquad n=0,1,2,\dots$$

U p.90
Si converge, converge a un punto fijo

Hipótesis. $g$ continua y $x_n\to a$. Conclusión. $a=g(a)$, pues $\lim_{n\to\infty} x_n=\lim_{n\to\infty} g(x_{n-1})=g(\lim_{n\to\infty} x_{n-1})$ (la continuidad de $g$ permite meter el límite dentro). De nuevo: si converge. Falta una condición que lo garantice.

U p.91

Cuándo funciona: los teoremas 2.2 y 2.3

Teorema 2.2 (existencia, unicidad y convergencia)

Hipótesis (las tres, a la vez):

  • (H0) $g$ continua en $[a,b]$. U la supone desde que presenta el método (p.90: «$g$ es una función continua»); hay que tenerla presente porque (H2) solo habla del abierto $(a,b)$.
  • (H1) $g:[a,b]\to[a,b]$: la función lleva el intervalo cerrado dentro de sí mismo. Es decir, $g$ está definida en $[a,b]$ y $g(x)\in[a,b]$ para todo $x\in[a,b]$.
  • (H2) Existe $k\in[0,1)$ con $|g'(x)|\le k$ para todo $x\in(a,b)$. (Implica que $g$ es derivable en $(a,b)$, pero no dice nada de los extremos ni de la continuidad en ellos.)

Conclusión. (1) $g$ tiene un único punto fijo $\tilde x$ en $[a,b]$; (2) para cualquier $x_0\in[a,b]$ la sucesión $x_{n+1}=g(x_n)$ converge a $\tilde x$.

U p.91
Teorema 2.3 (cota del error)

En las condiciones del Teorema 2.2, si $\tilde x$ es el único punto fijo de $g$ en $[a,b]$ y $k$ es la cota de $|g'|$ (con $k>0$: si $k=0$, $g$ es constante, $x_1=\tilde x$, la desigualdad estricta $0<0$ falla y $\ln k$ no existe): $$|x_n-\tilde x|<\frac{k^n}{1-k}\,(b-a).$$ Número de iteraciones. Para garantizar error menor que $\varepsilon$ hay que tomar el menor natural $n$ tal que $\dfrac{k^n}{1-k}(b-a)<\varepsilon$.

U pp.91-92
Error típico: comprobar solo $|g'|<1$

Hacen falta todas las hipótesis. Si $g$ no lleva $[a,b]$ en sí mismo, una iteración puede salir del intervalo y ya no vale la cota de $g'$ para el paso siguiente. Además, $k$ debe ser un número que de verdad acote a $|g'|$ en todo $(a,b)$: se toma el máximo (o un número mayor), nunca el mínimo. E tiene una errata precisamente en esto (Ejercicio 2.38).

Y (H0) no sobra. Sea $g(0)=1$ y $g(x)=x/2$ para $x\in(0,1]$. Lleva $[0,1]$ en $[0,1]$ y $|g'(x)|=\frac12$ en $(0,1)$, pero no tiene ningún punto fijo ($g(0)=1\neq0$ y $x/2=x$ solo si $x=0$): (H1) y (H2) se cumplen y falla la continuidad en $0$. Casi siempre $g$ es una función elemental continua en su dominio y (H0) se comprueba diciéndolo.

Por qué funciona (explicación propia: «contracciones»)

Si $|g'|\le k<1$, por el teorema del valor medio (U §2.5.3, p. 96, que viene después de este tema) $|g(x)-g(y)|\le k|x-y|$: la función acerca los puntos (es una contracción). Entonces:

  • Existe: $h(x)=g(x)-x$ es continua (H0) y cumple $h(a)\ge0\ge h(b)$ (porque $g(x)\in[a,b]$). Si $h(a)=0$ o $h(b)=0$, ese extremo ya es un punto fijo; si no, $h(a)>0>h(b)$ y Bolzano da un cero.
  • Es único: si hubiera dos, $|\tilde x-\tilde y|=|g(\tilde x)-g(\tilde y)|\le k|\tilde x-\tilde y|$ con $k<1$ obliga a $\tilde x=\tilde y$.
  • Converge: $|x_n-\tilde x|=|g(x_{n-1})-g(\tilde x)|\le k|x_{n-1}-\tilde x|\le\dots\le k^n|x_0-\tilde x|\le k^n(b-a)\to0$. Como $k^n(b-a)<\frac{k^n}{1-k}(b-a)$, sale la cota de U.

Cada paso reduce el error como mucho a una fracción $k$: es una convergencia lineal (comparada con la cuadrática de Newton). Cuanto más cerca de $1$ esté $k$, más lenta.

P
0 1 1 x x₀ x₁ x₂ y = g(x) = e⁻ˣ (curva azul) y = x (recta gris); escalera naranja En vertical hasta la curva, en horizontal hasta la recta y = x y otra vez en vertical: cada escalón es una iteración. Como g′ < 0, es una espiral: las iteraciones oscilan a ambos lados de la solución (círculo, 0.5671).
Punto fijo para $g(x)=e^{-x}$ con $x_0=1$ (Ejercicio 2.38). Cada «escalón» es una iteración $x_{n+1}=g(x_n)$. Si $g'>0$ (y $|g'|<1$) la escalera es monótona; si $g'<0$ es una espiral.
Ejemplo 2.17 de U: $\tfrac12\cos x=x$

a) Comprobar convergencia para todo $x_0\in[0,\frac\pi2]$. Es buscar puntos fijos de $g(x)=\frac12\cos x$ en $[a,b]=[0,\frac\pi2]$.

  1. (H0) $g(x)=\frac12\cos x$ es continua en $\mathbb R$.
  2. (H1) $g:[0,\frac\pi2]\to[0,\frac\pi2]$: como $\cos x\in[0,1]$ para $x\in[-\frac\pi2,\frac\pi2]$ y $[0,\frac\pi2]$ está dentro, $0\le\frac12\cos x\le\frac12<\frac\pi2$.
  3. (H2) $|g'(x)|=\frac12|\operatorname{sen}x|\le\frac12<1$: tomamos $k=\frac12$.
  4. Por el Teorema 2.2: un único punto fijo y convergencia para cualquier $x_0$ del intervalo.

b) Cuatro iteraciones con $x_0=1$ (radianes): $x_1=0.2701$, $x_2=0.4819$, $x_3=0.4430$, $x_4=0.4517$. (U imprime $x_3=0.4430$ por redondear en cada paso; con todas las cifras $x_3=0.443066\ldots$, que redondea a $0.4431$.)

c) Error de $x_4$ (Teorema 2.3): $|x_4-\tilde x|<\dfrac{(1/2)^4}{1-1/2}\cdot\dfrac\pi2=\dfrac\pi{16}\approx0.1963$.

d) Iteraciones para error $<0.01$: $\dfrac{(1/2)^n}{1/2}\cdot\dfrac\pi2<0.01\iff\left(\tfrac12\right)^n<\dfrac{0.01}{\pi}\iff n>\dfrac{\ln(0.01/\pi)}{\ln(1/2)}\approx8.2953$ (se invierte la desigualdad al dividir entre $\ln\frac12<0$). Hacen falta 9 iteraciones.

U pp.92-93

De $f(x)=0$ a $x=g(x)$: la $g$ no es única

Cómo se elige $g$

Para resolver $f(x)=0$ con punto fijo hay que despejar una ecuación equivalente $x=g(x)$; hay infinitas maneras. La más sencilla es $g(x)=f(x)+x$. Ejemplo 2.18 de U con $x^3-x=0$: sumando $x$ sale $g(x)=x^3$; con raíces cúbicas, $g(x)=\sqrt[3]x$; dividiendo por $x^2$ (buscando raíz no nula), $g(x)=\frac1x$. No todas sirven: hay que elegir una con $|g'|<1$ cerca de la raíz (ver Ejemplo propio P6).

U pp.93-94
Newton es un punto fijo

El método de Newton es el método de punto fijo con $g(x)=x-\dfrac{f(x)}{f'(x)}$: los puntos fijos de esa $g$ son exactamente los ceros de $f$ (donde $f'\ne0$). Además $g'(\tilde x)=\dfrac{f(\tilde x)f''(\tilde x)}{[f'(\tilde x)]^2}=0$ en una raíz simple, es decir, «$k=0$»; por eso Newton es tan rápido cerca de la raíz. Fíjate en que es justo el cociente de la condición de Larson.

P

3¿Cuántas iteraciones? Tabla para decidir P U pp.54,89-92

Lo primero es leer bien qué te piden: casi todos los errores de este tema vienen de mezclar los casos. Esta tabla los separa.

Si el enunciado dice…MétodoCuántas iteracionesFórmula / criterioEjemplo
«realiza $n$ iteraciones»Newton o punto fijoexactamente $n$tabla $x_0,x_1,\dots,x_n$; el resultado es $x_n$E 2.34 (4), 2.35 (3), 2.36 (3), 2.38 (4)
«hasta que el error relativo aproximado sea menor que…»Newton (o punto fijo)no se sabe a priori: se para cuando se cumple$\varepsilon_r=\dfrac{|x_n-x_{n-1}|}{|x_n|}<\text{tol}$ (más un tope máximo)P2 (4 pasos para $10^{-6}$)
«¿cuántas iteraciones garantizan error $<\varepsilon$?» con $g:[a,b]\to[a,b]$ y $|g'|\le k$punto fijoel menor natural $n$ con $\dfrac{k^n}{1-k}(b-a)<\varepsilon$$n>\dfrac{\ln\big(\varepsilon(1-k)/(b-a)\big)}{\ln k}$ (al dividir entre $\ln k<0$ se invierte la desigualdad)E 2.37: 21 con $k=0.8$; 3 con $k=0.2$. U 2.17 d: 9
«¿cuántas iteraciones para error $<\varepsilon$?» con cambio de signo en $[a,b]$bisecciónel menor $n$ con $\dfrac{b-a}{2^n}<\varepsilon$$n>\log_2\dfrac{b-a}{\varepsilon}$U Ej. 1.49: $n>3$, luego 4
«¿cuántas iteraciones de Newton garantizan…?»NewtonU no da fórmulausar $\varepsilon_r$ como criterio de parada práctico, no como garantía—

Lista de pasos para el caso «garantizar error $<\varepsilon$» (punto fijo)

  1. Identifica $g$, $[a,b]$ y $k$. Comprueba (H1) y (H2) antes de usar la fórmula; si no se cumplen, el Teorema 2.3 no dice nada.
  2. Sustituye en $\dfrac{k^n}{1-k}(b-a)<\varepsilon$ y aísla $k^n<\dfrac{\varepsilon(1-k)}{b-a}$.
  3. Toma $\ln$ (conserva el orden, es creciente). Divide entre $\ln k$, que es negativo: la desigualdad se invierte.
  4. El resultado es «$n>$ número». Toma el primer natural estrictamente mayor (si sale $20.64$, es $21$; si saliera exactamente $2$, sería $3$).
  5. Comprueba sustituyendo $n-1$ (no basta) y $n$ (basta).

Valores ya calculados (para comparar tu resultado)

Situación$k$$b-a$$\varepsilon$$n>$$n$
E 2.37 (a)0.820.120.6421
E 2.37 (b)0.220.12 (justo)3
U Ej. 2.17 d0.5$\pi/2$0.018.309
E 2.38 (cota pesimista)$e^{-0.1}\approx0.9048$1.90.0175.9976
P40.51$10^{-3}$10.9711
Cota garantizada ≠ iteraciones necesarias

El Teorema 2.3 da una garantía, no una estimación ajustada. En el Ejercicio 2.38 la garantía exige $76$ iteraciones para error $<0.01$, pero la tabla real alcanza ese error con $7$ (con $6$ el error real es $0.0125$; con $7$, $0.0070$). Cuando $k$ está cerca de $1$ la cota es muy pesimista (con $4$ iteraciones vale $13.4$, mayor que el propio intervalo). Y en Newton, $\varepsilon_r$ es solo una estimación.

4Comparación de los tres métodos U pp.86-87 P

Bisección (U §1.4)Punto fijoNewton
Qué necesita$f$ continua y cambio de signo en $[a,b]$ (Bolzano)$g$ continua, $g:[a,b]\to[a,b]$ y $|g'|\le k<1$ en $(a,b)$$f$ derivable, $f'(x_n)\ne0$ y un buen $x_0$
Garantía a priorisiempre converge; error $<\frac{b-a}{2^n}$sí, Teorema 2.3ninguna en U (solo condición suficiente en L)
Velocidadlenta: el error se divide entre 2lineal: el error se multiplica por $\approx|g'(\tilde x)|$cuadrática (raíz simple); lineal si la raíz es doble
Punto débilno sirve si $f$ no cambia de signo (p. ej. $x^2=0$)elegir bien $g$ (no todas valen)$f'=0$, ciclos, mal $x_0$

Ventajas de Newton sobre bisección según U: no hace falta conocer un intervalo con la solución y, cuando converge, lo hace mucho más rápido. (U p.87.) La comparación de velocidades es explicación propia.

5Ejercicios del libro, resueltos E pp.72-76

Orden recomendado

Primero el bloque de Newton, 2.34 → 2.35 → 2.36; luego punto fijo, 2.38 → 2.37. El 2.38 enseña a comprobar las hipótesis del teorema y el 2.37 solo usa la cota del error, que viene justo después. Aquí aparecen en ese orden (la numeración no cambia). Abre cada solución con su desplegable o con el botón «Mostrar soluciones».

Bloque Newton

Ejercicio 2.34 — aproximar $\sqrt5$ con cuatro iteraciones de Newton E p.72

Usando el método de Newton, calcula un valor aproximado de $\sqrt5$, realizando cuatro iteraciones.

  1. Convertir «calcular $\sqrt5$» en «resolver $f(x)=0$». $\sqrt5$ es el número positivo cuyo cuadrado vale $5$: solución de $x^2-5=0$. Tomamos $f(x)=x^2-5$, derivable en todo $\mathbb R$, con $f'(x)=2x$.
  2. Elegir $x_0$ con criterio. $f(2)=-1<0<f(3)=4$ y $f$ es continua: por Bolzano hay un cero en $(2,3)$. Tomamos $x_0=2$; además $f'(2)=4\ne0$, luego se puede aplicar la fórmula.
  3. Simplificar la iteración. $$x_{n+1}=x_n-\frac{x_n^2-5}{2x_n}=\frac{2x_n^2-x_n^2+5}{2x_n}=\frac12\Big(x_n+\frac5{x_n}\Big).$$ Es la media entre $x_n$ y $5/x_n$: si uno se queda corto, el otro se pasa.
  4. Iterar con fracciones exactas: $x_1=\frac94$, $x_2=\frac{161}{72}$, $x_3=\frac{51841}{23184}$, $x_4=\frac{5374978561}{2403763488}$.
$n$$x_n$$f(x_n)$error real $|x_n-\sqrt5|$$\varepsilon_r$
02$-1$$2.36\times10^{-1}$—
12.25$0.0625$$1.39\times10^{-2}$$1.11\times10^{-1}$
22.236111111111$1.93\times10^{-4}$$4.31\times10^{-5}$$6.21\times10^{-3}$
32.236067977916$1.86\times10^{-9}$$4.16\times10^{-10}$$1.93\times10^{-5}$
42.236067977500$1.7\times10^{-19}$$3.9\times10^{-20}$$1.86\times10^{-10}$
  1. Conclusión. $\sqrt5\approx x_4=2.2360679775$ ($\sqrt5=2.2360679774997\ldots$). Con 4 decimales, como en E, $x_2=x_3=x_4=2.2361$; ojo, esa «estabilización» es solo a 4 decimales: el error real de $x_2$ es $4.3\times10^{-5}$.
  2. Por qué converge tan rápido (propio, álgebra pura). Restando $\sqrt5$ en la fórmula del paso 3: $$x_{n+1}-\sqrt5=\frac{x_n^2-2\sqrt5\,x_n+5}{2x_n}=\frac{(x_n-\sqrt5)^2}{2x_n}.$$ Si $x_n>0$: (i) $x_{n+1}\ge\sqrt5$ (desde $x_1$ todas son por exceso, como se ve); (ii) el error nuevo es aproximadamente el cuadrado del anterior dividido entre $2\sqrt5\approx4.47$. Los errores $10^{-2}\to10^{-5}\to10^{-10}\to10^{-20}$ duplican los decimales correctos: convergencia cuadrática.
Receta: para aproximar $\sqrt[m]a$, toma $f(x)=x^m-a$, localiza la raíz con Bolzano, parte de un entero cercano, simplifica $x_{n+1}=x_n-\frac{x_n^m-a}{m\,x_n^{m-1}}$ y rellena una tabla $n\,|\,x_n\,|\,f(x_n)\,|\,f'(x_n)$.

Error típico: aplicar Newton a «$x=\sqrt5$», o redondear demasiado pronto (con 2 decimales la tabla «se estabiliza» en $2.24$, falso en el tercer decimal).

Dónde está la teoría

U §2.4.2 pp.86-89 (algoritmo p.88; Ejemplo 2.15 es el mismo esquema con $x^4-5$). Comparar con L ej. 19 (regla de Herón), p.229.

U p.88
Ejercicio 2.35 — Newton para $x^2=0$ desde $x_0=0.4$; ¿se puede usar bisección? E p.73

Utiliza el método de Newton para aproximar la solución de $x^2=0$, tomando como dato inicial $x_0=0.4$ y realizando tres iteraciones. ¿Es posible utilizar directamente el método de bisección para aproximar la solución de la ecuación $x^2=0$?

  1. Función y derivada. $f(x)=x^2$, $f'(x)=2x$, derivable en $\mathbb R$. La solución exacta es $\tilde x=0$; el ejercicio sirve para ver cómo se comporta el método.
  2. Iteración simplificada. Si $x_n\ne0$ (para poder dividir entre $2x_n$): $x_{n+1}=x_n-\dfrac{x_n^2}{2x_n}=\dfrac{x_n}2$. Con $x_0=0.4$ todas las $x_n=0.4/2^n$ son no nulas: el método nunca se bloquea.
  3. Tres iteraciones. $x_1=0.2$, $x_2=0.1$, $x_3=0.05$.
$n$$x_n$$f(x_n)$$f'(x_n)$error $|x_n-0|$
00.40.160.80.4
10.20.040.40.2
20.10.010.20.1
30.050.00250.10.05
  1. Observación (no está en E). El error solo se divide entre 2 en cada paso; no se eleva al cuadrado como en 2.34. Motivo: $0$ es una raíz doble, $f(0)=0$ y también $f'(0)=0$. La tangente se vuelve casi horizontal al acercarse y Newton pierde su rapidez.
  2. ¿Bisección? Necesita un intervalo $[a,b]$ con $f$ continua y $f(a)f(b)<0$ (Teorema de Bolzano). Pero $x^2\ge0$ siempre, así que nunca hay cambio de signo. Respuesta: no, no se puede aplicar directamente. Bolzano da una condición suficiente, no necesaria: aquí hay raíz, pero la gráfica «toca» el eje sin cruzarlo.
Receta: simplifica $x-\frac{f}{f'}$ antes de iterar. Para saber si puedes usar bisección, busca dos puntos donde $f$ tenga signos opuestos; si no cambia de signo, no se puede.

Error típico: pensar que «si Bolzano no se cumple, no hay raíz»; o creer que Newton siempre converge cuadráticamente (falla si $f'(\tilde x)=0$).

Dónde está la teoría

U §2.4.2 p.88; U §1.4 pp.52-55 (Bolzano y bisección).

U pp.52-55,88
Ejercicio 2.36 — Newton para $e^{-x}-x=0$ desde $x_0=0.5$ y $x_0=5$; error relativo E pp.73-74

Utilizando el método de Newton, encuentra una solución aproximada de la ecuación $e^{-x}-x=0$ realizando tres iteraciones y partiendo de las condiciones iniciales $x_0=0.5$ y $x_0=5$. Estima el error relativo aproximado en cada caso.

  1. Hay una única solución y sabemos dónde. $f(x)=e^{-x}-x$, continua y derivable en $\mathbb R$. $f(0)=1>0$ y $f(1)=e^{-1}-1\approx-0.632<0$: por Bolzano hay una raíz en $(0,1)$. Es única porque $e^{-x}$ y $-x$ son estrictamente decrecientes, luego $f$ también y no puede anularse dos veces. (Es la ecuación del Ejemplo 1.49 de U, p.55, resuelta allí por bisección.)
  2. Derivada que nunca se anula. $f'(x)=-e^{-x}-1<0$ para todo $x$, porque $e^{-x}>0$: la división de Newton siempre es posible.
  3. Iteración simplificada. $$x_{n+1}=x_n-\frac{e^{-x_n}-x_n}{-e^{-x_n}-1}=x_n+\frac{e^{-x_n}-x_n}{e^{-x_n}+1}=\frac{x_n+1}{1+e^{x_n}}$$ (en el último paso se multiplica numerador y denominador por $e^{x_n}$).
  4. Caso $x_0=0.5$.
$n$$x_n$$f(x_n)$$\varepsilon_r$error real
00.5$1.065\times10^{-1}$—$6.7\times10^{-2}$
10.566311003197$1.305\times10^{-3}$$1.17\times10^{-1}$$8.3\times10^{-4}$
20.567143165035$1.96\times10^{-7}$$1.47\times10^{-3}$$1.3\times10^{-7}$
30.567143290410$4.5\times10^{-15}$$2.21\times10^{-7}$$2.8\times10^{-15}$

Solución exacta: $\tilde x=0.567143290409784\ldots$ Error relativo aproximado tras 3 iteraciones: $\varepsilon_r=\frac{|x_3-x_2|}{|x_3|}\approx2.2\times10^{-7}$. Con 4 decimales (como hace E): $x_1\approx0.5663$, $x_2\approx x_3\approx0.5671$ y sale $\varepsilon_r=0$. Ese «0» no significa error nulo: el error es menor de lo que 4 decimales pueden mostrar.

  1. Caso $x_0=5$.
$n$$x_n$$f(x_n)$$\varepsilon_r$error real
05$-4.99326$—4.43
10.0401571055460.92048123.5$5.3\times10^{-1}$
20.5096375311610.09108$9.21\times10^{-1}$$5.8\times10^{-2}$
30.566534509041$9.54\times10^{-4}$$1.004\times10^{-1}$$6.1\times10^{-4}$

$\varepsilon_r=\frac{|x_3-x_2|}{|x_3|}=\frac{|0.566535-0.509638|}{0.566535}\approx0.1004$, igual que en E.

  1. Por qué el primer salto es tan grande. En $x_0=5$ la función es casi la recta $y=-x$ ($e^{-5}\approx0.0067$), de pendiente casi $-1$; su tangente corta el eje cerca de $0$, al otro lado de la raíz. Desde ahí ya está cerca y converge deprisa. El resultado tras un número fijo de pasos depende mucho de $x_0$: conviene localizar antes la raíz.
  2. Cómo interpretar $\varepsilon_r$. No es el error verdadero: mide cuánto ha cambiado la última iteración y se usa como estimación. Con $x_0=5$, $\varepsilon_r\approx0.1$ sobreestima el error relativo real ($\approx10^{-3}$): en Newton suele ocurrir, porque el error real ya es mucho menor que el último salto.
Receta: (1) Bolzano para localizar la raíz. (2) Comprobar $f'\ne0$. (3) Simplificar $x-f/f'$. (4) Tabla con $x_n$ y $\varepsilon_r=|x_n-x_{n-1}|/|x_n|$. (5) Parar cuando $\varepsilon_r$ sea menor que la tolerancia o se agoten las iteraciones pedidas.

Error típico: perder el signo de $f'$: el denominador es $-e^{-x_n}-1$. Si se usa $+e^{-x_n}+1$, el paso va en sentido contrario (desde $0.5$ saldría $0.434$, que se aleja). Otro: dividir $\varepsilon_r$ entre $|x_{n-1}|$ en vez de $|x_n|$ (U usa $|x_n|$).

Dónde está la teoría

U §2.4.2 pp.88-90 (algoritmo p.88; errores aproximados y criterios de parada pp.89-90).

U pp.88-90

Bloque punto fijo

Ejercicio 2.38 — convergencia del punto fijo para $e^{-x}=x$ en $[0.1,2]$ y cuatro iteraciones E pp.75-76

Dada la ecuación $e^{-x}=x$: (1) Comprueba que el método de punto fijo converge para cualquier dato inicial $0.1\le x_0\le2$. (2) Utiliza el método de punto fijo para aproximar la solución de la ecuación partiendo del dato inicial $x_0=1$ y realizando cuatro iteraciones.

Apartado 1.

  1. Forma $x=g(x)$. Ya lo está: $g(x)=e^{-x}$. Sus puntos fijos son exactamente las soluciones de la ecuación (Definición 2.4).
  2. Qué hay que comprobar (Teorema 2.2) en $[a,b]=[0.1,2]$: (H0) $g$ continua (inmediato: $e^{-x}$ es continua en $\mathbb R$); (H1) $g([a,b])\subset[a,b]$; (H2) existe $k\in[0,1)$ con $|g'(x)|\le k$ en $(a,b)$.
  3. (H1). $g'(x)=-e^{-x}<0$: $g$ es estrictamente decreciente, e invierte las desigualdades. Para $0.1\le x\le2$: $$g(2)\le g(x)\le g(0.1)\iff 0.1353\approx e^{-2}\le e^{-x}\le e^{-0.1}\approx0.9048.$$ Como $0.1<0.1353$ y $0.9048<2$, se cumple $g(x)\in[0.1,2]$. Por qué importa: si una iteración saliera de $[a,b]$, ya no podríamos usar la cota de $g'$ en el paso siguiente.
  4. (H2). $|g'(x)|=e^{-x}$ es decreciente, así que su mayor valor en $(0.1,2)$ está cerca de $x=0.1$: $|g'(x)|<e^{-0.1}\approx0.9048<1$. Tomamos $k=e^{-0.1}\approx0.905$. (E escribe «$k=0.091$», que es una errata: sirve cualquier $k\in[e^{-0.1},1)$ (con $k=0.905$ vale), pero $0.091$ no acota a $|g'|$.)
  5. Conclusión. Por el Teorema 2.2, $g$ tiene un único punto fijo en $[0.1,2]$ y $x_{n+1}=e^{-x_n}$ converge a él para cualquier $x_0\in[0.1,2]$.

Apartado 2. Iteraciones con $x_0=1$ (y cuatro más para ver la tendencia; $\tilde x=0.567143\ldots$):

$n$$x_n=e^{-x_{n-1}}$error real $|x_n-\tilde x|$$\varepsilon_r$
010.4329—
10.3678794411710.19931.718
20.6922006275550.12510.4685
30.5004735005640.06670.3831
40.6062435350860.03910.1745
50.5453957859750.02170.1116
60.5796123355030.01250.0590
70.5601154613610.00700.0348
80.5711431150800.00400.0193
  1. Respuesta. $x_4\approx0.6062$. Coincide con E aunque se redondee a 4 decimales en cada paso: $0.3679,\ 0.6922,\ 0.5005,\ 0.6062$.
  2. Lectura de la tabla (propio). Las iteraciones oscilan alrededor de $\tilde x$ (una por defecto, otra por exceso), lo que ocurre cuando $g'<0$. El error se multiplica más o menos por $|g'(\tilde x)|=\tilde x\approx0.567$ en cada paso: convergencia lineal. Newton sobre la misma ecuación (Ej. 2.36) lograba 15 decimales exactos en 3 pasos.
  3. ¿Qué garantiza la teoría con 4 pasos? Teorema 2.3: $\frac{k^4}{1-k}(b-a)=\frac{0.9048^4}{0.0952}\cdot1.9\approx13.4$: inútil (el intervalo mide $1.9$). Para garantizar error $<0.01$ harían falta $n>\ln\!\big(0.01(1-k)/1.9\big)/\ln k\approx75.99$, o sea $76$ iteraciones; en la práctica bastan $7$ (según la tabla, $|x_6-\tilde x|=0.0125$ no basta y $|x_7-\tilde x|=0.0070$ sí).
Receta (convergencia del punto fijo en $[a,b]$): (1) pasa la ecuación a $x=g(x)$; (2) estudia la monotonía de $g$ y calcula $g(a)$, $g(b)$ para ver que $g([a,b])\subset[a,b]$; (3) acota $|g'|$ en $(a,b)$ por su valor mayor $k$ y comprueba $k<1$; (4) cita el Teorema 2.2.

Error típico: comprobar solo $|g'|<1$ y olvidar que $g$ debe llevar $[a,b]$ en sí mismo; o dar como $k$ un número que no acota realmente a $|g'|$ (la errata de E).

Dónde está la teoría

U §2.4.3 pp.90-93 (Definición 2.4, ec. (2.3) p.90; Teoremas 2.2 y 2.3 p.91; el Ejemplo 2.17 es el modelo exacto de este ejercicio).

U pp.90-93
Ejercicio 2.37 — número de iteraciones del punto fijo para error $<0.1$ ($k=0.8$ y $k=0.2$) E pp.74-75

Sea $f:[-1,1]\to[-1,1]$. Sabiendo que la única solución de la ecuación $f(x)=x$ en el intervalo $[-1,1]$ puede ser aproximada utilizando el método de punto fijo, ¿cuántas iteraciones necesitaremos para garantizar que el error cometido sea menor que $0.1$ si $|f'(x)|<0.8$ en $[-1,1]$? ¿Y si $|f'(x)|<0.2$ en $[-1,1]$?

  1. Datos del Teorema 2.3. Aquí la función de iteración se llama $f$ (en U es $g$). $[a,b]=[-1,1]$, $f([-1,1])\subset[-1,1]$ y $|f'|<k$ con $k=0.8$ o $k=0.2$, ambos en $[0,1)$. Se cumplen las hipótesis del Teorema 2.2, y por el Teorema 2.3 $$|x_n-\tilde x|<\frac{k^n}{1-k}(b-a)=\frac{2k^n}{1-k}.$$
  2. Planteamiento. Buscamos el menor $n\in\mathbb N$ con $\dfrac{2k^n}{1-k}<0.1$.
  3. Caso $k=0.8$. $\dfrac{2\cdot0.8^n}{0.2}<0.1\iff10\cdot0.8^n<0.1\iff0.8^n<0.01$. Tomamos $\ln$ (creciente, conserva el sentido): $n\ln0.8<\ln0.01$. Dividimos entre $\ln0.8<0$, que invierte la desigualdad: $$n>\frac{\ln0.01}{\ln0.8}\approx20.6377.$$ El primer natural es $n=21$. Comprobación: $n=20$ da cota $0.1153$ (no basta); $n=21$ da $0.0922<0.1$.
  4. Caso $k=0.2$. $\dfrac{2\cdot0.2^n}{0.8}<0.1\iff0.2^n<0.04$. Como $0.2^2=0.04$ no es $<0.04$ y $0.2^3=0.008<0.04$, sale $n=3$. (Cotas: $n=2\to0.1$, $n=3\to0.02$.)
  5. Moraleja. Con el mismo intervalo, pasar de $k=0.8$ a $k=0.2$ reduce la garantía de $21$ a $3$ iteraciones: cuanto menor es $k$, más rápida es la convergencia.
  6. Matiz de rigor (omisión en E). El Teorema 2.3 de U tiene desigualdad estricta. Con $k=0.2$ y $n=2$ la cota vale exactamente $0.1$, y como el teorema dice «error $<$ cota», en rigor ya daría error $<0.1$. Pero la receta de U (p.92) pide el natural más pequeño para el que la cota $\frac{k^n}{1-k}(b-a)$ sea $<\varepsilon$, y con $n=2$ la cota es $=0.1$, no menor. La respuesta esperada es 3, que es la de E y la que se da en el examen. (Con la cota más fina $k^n(b-a)$ de la demostración bastarían $14$ y $2$, pero no es la fórmula del libro: no la uses salvo que la justifiques.)
Receta: resuelve $\frac{k^n}{1-k}(b-a)<\varepsilon\Rightarrow n>\dfrac{\ln\big(\varepsilon(1-k)/(b-a)\big)}{\ln k}$ y toma el primer natural estrictamente mayor. Recuerda que $\ln k<0$ invierte la desigualdad.

Error típico: no invertir la desigualdad al dividir entre $\ln k$ (sale $n<20.6$, absurdo); redondear $20.64$ a $20$ en vez de a $21$; olvidar el factor $(b-a)=2$ o el $1/(1-k)$.

Dónde está la teoría

U §2.4.3, Teorema 2.3 p.91 y procedimiento p.92 (Ejemplo 2.17 d, p.93).

U pp.91-93
Erratas y omisiones de E en este tema
  • Ej. 2.38 (p.76): «$k=0.091$» es errata; debe ser $k=e^{-0.1}\approx0.905$ (vale cualquier $k\in[e^{-0.1},1)$). La conclusión (convergencia) sigue siendo correcta.
  • Ej. 2.36 (p.74): con $x_0=5$, $x_1=0.0401571\ldots$ redondea a $0.0402$; E trunca. No afecta a $x_2$, $x_3$ ni a $\varepsilon_r\approx0.1004$. El «$\varepsilon_r=0$» con $x_0=0.5$ es solo redondeo a 4 decimales.
  • Errata de U, Ej. 2.15 (p.88): la tabla imprime $f(1{,}5173)=-6{,}5241\cdot10^{-4}$, pero el valor correcto es $f(1.5173)=1.5173^4-5\approx+0.3001$ (confirmado en el PDF). El resto de la tabla es correcto; por eso aquí solo se reproducen $x_i$ y errores.
  • Ej. 2.37 (p.75): con $k=0.2$ y $n=2$ la cota vale exactamente $0.1$; E da $3$ aplicando la receta de U. Se mantiene $3$.
  • Omisiones: E no explica por qué $x_0=2$ en 2.34 (Bolzano), que $0$ es raíz doble en 2.35, ni el valor del error de $x_4$ en 2.38 ($0.039$).
E pp.72-76

Ejemplos propios (de fácil a difícil)

P1 (fácil) — $\sqrt[3]{10}$ con tres iteraciones de Newton P
  1. $f(x)=x^3-10$, $f'(x)=3x^2$. Como $f(2)=-2<0<f(3)=17$, por Bolzano la raíz está en $(2,3)$. Tomamos $x_0=2$.
  2. $x_{n+1}=x_n-\dfrac{x_n^3-10}{3x_n^2}=\dfrac{2x_n^3+10}{3x_n^2}$.
$n$$x_n$error real
02$1.5\times10^{-1}$
1$13/6=2.166666666667$$1.2\times10^{-2}$
22.154503616042$6.9\times10^{-5}$
32.154434692237$2.2\times10^{-9}$

$\sqrt[3]{10}=2.154434690031\ldots$: $x_3$ tiene 8 decimales exactos.

Receta: la misma que en 2.34, con $m=3$.
P2 (medio) — resolver $x^3-x-1=0$ con Newton hasta $\varepsilon_r<10^{-6}$ P
  1. $f(1)=-1<0$, $f(2)=5>0$: hay raíz en $(1,2)$. $f'(x)=3x^2-1\ge2$ en $[1,2]$, luego no se anula.
  2. $x_{n+1}=x_n-\dfrac{x_n^3-x_n-1}{3x_n^2-1}$, con $x_0=1.5$.
$n$$x_n$$\varepsilon_r$
01.5—
11.347826086957$1.13\times10^{-1}$
21.325200398951$1.71\times10^{-2}$
31.324718173999$3.64\times10^{-4}$
41.324717957245$1.64\times10^{-7}$

Paramos en $n=4$ porque $\varepsilon_r<10^{-6}$: $\tilde x\approx1.3247180$ (exacto $1.32471795724\ldots$).

Receta: criterio de parada por $\varepsilon_r$: itera y calcula $\varepsilon_r$ en cada fila hasta que baje de la tolerancia.
P3 (difícil) — Newton puede no converger: $x^3-2x+2$ desde $x_0=0$ P
  1. $f'(x)=3x^2-2$. Desde $x_0=0$: $x_1=0-\frac{2}{-2}=1$; desde $x_1=1$: $x_2=1-\frac11=0$. La sucesión es $0,1,0,1,\dots$: un ciclo, sin límite.
  2. ¿Dónde está la raíz? $f(-2)=-2<0$, $f(-1)=3>0$: en $(-2,-1)$.
  3. Desde $x_0=-2$: $-1.8,\ -1.769948186529,\ -1.769292662906,\ -1.769292354239$, con errores $3\times10^{-2},\ 6.6\times10^{-4},\ 3.1\times10^{-7},\ 6.8\times10^{-14}$.
Moraleja: localiza la raíz con Bolzano antes de elegir $x_0$.
P4 (fácil) — iteraciones de punto fijo con $k=0.5$ en $[0,1]$ P

$g:[0,1]\to[0,1]$ con $|g'(x)|\le0.5$ en $(0,1)$. ¿Cuántas iteraciones garantizan error $<10^{-3}$?

$\dfrac{0.5^n}{0.5}\cdot1<10^{-3}\iff0.5^n<5\cdot10^{-4}\iff n>\dfrac{\ln(5\cdot10^{-4})}{\ln0.5}\approx10.97$, luego $n=11$. (Cotas: $n=10\to1.95\times10^{-3}$; $n=11\to9.8\times10^{-4}$.)

P5 (medio) — $x=\cos x$ en $[0,1]$: hipótesis, iteraciones y garantía P
  1. (H1). $g(x)=\cos x$ es decreciente en $[0,1]\subset[0,\frac\pi2]$, luego $g([0,1])=[\cos1,\cos0]=[0.5403,1]\subset[0,1]$.
  2. (H2). $|g'(x)|=\operatorname{sen}x\le\operatorname{sen}1\approx0.8415<1$ en $(0,1)$ (el seno crece en $[0,\frac\pi2]$). Tomamos $k=\operatorname{sen}1$. Por el Teorema 2.2: un único punto fijo y convergencia desde cualquier $x_0\in[0,1]$.
  3. Con $x_0=1$: $x_1=0.540302305868$, $x_2=0.857553215846$, $x_3=0.654289790498$, $x_4=0.793480358743$. Valor exacto $\tilde x=0.739085133\ldots$; error de $x_4$: $0.054$ (oscila, $g'<0$).
  4. Garantía para error $<0.01$: $\dfrac{k^n}{1-k}<0.01\iff n>\dfrac{\ln(0.01(1-k))}{\ln k}\approx37.35$, luego $n=38$ (cota con $37$: $0.0106$; con $38$: $0.0089$). Newton resuelve la misma ecuación con seis decimales en unos $4$ pasos (S Ej. 3, p.341).
P6 (difícil) — elegir bien $g$: $x^3-x-1=0$ en $[1,2]$ P
  1. Opción mala: $x=x^3-1$, $g_1(x)=x^3-1$. Entonces $g_1'(x)=3x^2\ge3$ en $[1,2]$: (H2) falla, y $g_1([1,2])=[0,7]\not\subset[1,2]$: (H1) también. Desde $x_0=1.5$: $2.375,\ 12.40,\ 1904,\ 6.9\times10^{9}$. Diverge.
  2. Opción buena: $x^3=x+1\iff x=\sqrt[3]{x+1}$, $g_2(x)=\sqrt[3]{x+1}$. Es creciente: $g_2([1,2])=[\sqrt[3]2,\sqrt[3]3]=[1.2599,1.4422]\subset[1,2]$ (H1). $g_2'(x)=\frac{1}{3(x+1)^{2/3}}$ es positiva y decreciente, luego $|g_2'|\le g_2'(1)=\frac1{3\sqrt[3]4}\approx0.2100=k$ (H2).
  3. Desde $x_0=1$: $1.259921049895,\ 1.312293836683,\ 1.322353819139,\ 1.324268744552,\ 1.324632625251$ (errores $6.5\times10^{-2}\to8.5\times10^{-5}$).
  4. Error $<10^{-6}$ garantizado si $\frac{k^n}{1-k}<10^{-6}\iff n>9.003$, o sea $n=10$ (con $n=9$ la cota es $1.005\times10^{-6}$, que por muy poco no basta).
Moraleja: la misma ecuación admite muchas $g$ (U pp.93-94). Elige una con $|g'|<1$ cerca de la raíz.

6Así lo preguntan en el examen X

Estos son los tipos reales que caen de este tema en los exámenes resueltos (X = Examenes-resueltos-Calculo.pdf), ordenados de fácil a difícil. Las preguntas cortas (1 punto) piden una cuenta concreta: iteraciones de Newton o de punto fijo, o el número de pasos de bisección. Los ejercicios largos (3 puntos) añaden justificar las hipótesis del Teorema 2.2 y acotar el error con el Teorema 2.3. La dificultad no está en la teoría, sino en no saltarse hipótesis y en el orden de las cuentas.

X1 (fácil) — bisección: ¿cuántas iteraciones para error $<0.15$? X p.73

Función continua $f:[-1,1]\to\mathbb R$ con $f(-1)=1$ y $f(1)=-2$. Se busca una solución aproximada de $f(x)=0$ por bisección. ¿Cuántas iteraciones aseguran un error menor que $0.15$?

La idea en una frase

Cada paso de bisección parte el intervalo por la mitad, así que su longitud se divide entre $2$. Solo hay que contar cuántas mitades hacen falta para bajar de $0.15$.

  1. ¿Se puede aplicar? $f$ es continua y cambia de signo ($f(-1)=1>0$, $f(1)=-2<0$): por Bolzano hay una raíz entre $-1$ y $1$.
  2. Longitud de partida: de $-1$ a $1$ hay $2$.
  3. Cuenta de mitades (el error es menor que la longitud, U p.54): tras $1$ paso queda $1$; tras $2$, $0.5$; tras $3$, $0.25$; tras $4$, $0.125$.
  4. Primera que baja de $0.15$: la de $4$ pasos ($0.125<0.15$; con $3$ pasos son $0.25$, no vale).

Resultado: $n=4$ iteraciones.

Lo que escribes en el examen

«Por Bolzano hay raíz en $(-1,1)$. Tras $n$ pasos el error es $<\frac{2}{2^n}$. Con $n=4$: $\frac{2}{16}=0.125<0.15$; con $n=3$: $0.25$. Luego $n=4$.»

Receta: menor $n$ con $\frac{b-a}{2^n}<\varepsilon$; se tantea con potencias de $2$. No hace falta conocer $f$, solo la longitud del intervalo.

Parte del tema: §3 (fila de bisección); U §1.4 p.54.

X2 (fácil) — Newton, dos iteraciones: $4-2x=e^x$ desde $x_0=0$ X p.133

Aproxime una solución de $4-2x=e^x$ por el método de Newton con $x_0=0$ y dos iteraciones.

La idea en una frase

Newton: desde el punto actual se traza la tangente y el sitio donde corta al eje $x$ es el siguiente punto. En fórmula, «punto menos $f$ entre $f'$». Aquí solo hay que aplicarla dos veces.

  1. Pasar a $f(x)=0$. Todo a un lado: $f(x)=e^x+2x-4$.
  2. Derivar. $f'(x)=e^x+2$. Es siempre positiva, así que nunca se divide entre $0$ (U p.88).
  3. Fórmula. $x_{n+1}=x_n-\dfrac{e^{x_n}+2x_n-4}{e^{x_n}+2}$.
  4. Primera iteración ($x_0=0$): $f(0)=1+0-4=-3$ y $f'(0)=3$. Así $x_1=0-\dfrac{-3}{3}=1$.
  5. Segunda iteración ($x_1=1$): $f(1)=e-2$ y $f'(1)=e+2$. Así $x_2=1-\dfrac{e-2}{e+2}=\dfrac{4}{e+2}\approx0.8477$.

Resultado: $x_2=\dfrac{4}{e+2}\approx0.8477$ (comprobación: la raíz real es $0.8408\ldots$).

Lo que escribes en el examen

«$f(x)=e^x+2x-4$, $f'(x)=e^x+2\ne0$. $x_{n+1}=x_n-f(x_n)/f'(x_n)$. $x_1=1$, $x_2=\frac{4}{e+2}\approx0.8477$.»

Receta: escribe $f$ y $f'$, comprueba $f'\ne0$ y rellena una tabla $n\,|\,x_n\,|\,f(x_n)\,|\,f'(x_n)$ fila a fila.

Parte del tema: §1 (U pp.86-90).

X3 (fácil) — punto fijo, cinco iteraciones: $f(x)=1-\frac1{x^2+2}$ desde $x_0=0$ X p.1

Dada $f(x)=1-\dfrac1{x^2+2}$, encuentre un punto fijo con el método del punto fijo: 5 iteraciones, $x_0=0$.

La idea en una frase

Un punto fijo es un $x$ que la función deja igual: $f(x)=x$. Se encuentra repitiendo: metes un número en $f$, el resultado lo vuelves a meter, y así cinco veces.

  1. Regla. $x_{n+1}=f(x_n)=1-\dfrac1{x_n^2+2}$, empezando en $x_0=0$ (U p.90).
  2. $x_1=1-\dfrac1{0+2}=\dfrac12=0.5$.
  3. $x_2=1-\dfrac1{0.25+2}=\dfrac59\approx0.55556$.
  4. $x_3\approx0.56684$.
  5. $x_4\approx0.56921$.
  6. $x_5\approx0.56971$.

Resultado: $x_5\approx0.5697$ (el punto fijo real es $0.56984\ldots$).

Lo que escribes en el examen

La lista $x_1,\dots,x_5$ con la regla $x_{n+1}=f(x_n)$ y el resultado final. Con 1 punto basta eso. Si piden justificar la convergencia, mira el recuadro de abajo.

Extra (solo si piden justificar): ¿por qué converge?

En $[0,1]$: (H0) $f$ es continua; (H1) $f$ es creciente y $f(0)=\frac12$, $f(1)=\frac23$, así que $f$ lleva $[0,1]$ dentro de $[0,1]$; (H2) $|f'(x)|=\frac{2x}{(x^2+2)^2}\le\frac12=k<1$ (numerador $\le2$, denominador $\ge4$). Teorema 2.2 (U p.91): hay un único punto fijo y la iteración converge. Cota del Teorema 2.3: $|x_5-\tilde x|<\frac{(1/2)^5}{1/2}\cdot1=0.0625$ (pesimista: el error real es $\approx1.3\cdot10^{-4}$).

Receta: punto fijo con $g$ dada: itera $x_{n+1}=g(x_n)$ guardando todas las cifras; si piden justificar, comprueba H0-H2 en un intervalo que contenga a $x_0$.
Aviso sobre la solución impresa

X redondea $x_2$ a $0.5556$ y por eso imprime $x_3\approx0.56685$; con todas las cifras es $0.56684$. No cambia nada.

Parte del tema: §2, Teoremas 2.2 y 2.3 (U pp.90-92).

X4 (medio) — Newton con criterio de parada: $x-2^{-x}=0$ desde $x_0=0$ X p.20

Mediante el método de Newton, encontrar una raíz próxima a $x_0=0$ de $x-2^{-x}=0$. Utilizar tres decimales redondeados en cada iteración hasta que se cumpla $|x_{n-1}-x_n|<10^{-3}$.

La idea en una frase

Es Newton de siempre; lo único nuevo es cuándo parar: cuando dos iterados seguidos se diferencian en menos de $10^{-3}$. Se hace una tabla y se mira la última columna.

  1. Función. $f(x)=x-2^{-x}$. Como $f(0)=-1<0$ y $f(1)=\frac12>0$, hay raíz entre $0$ y $1$ (Bolzano).
  2. Derivada. $f'(x)=1+\ln2\cdot2^{-x}$ (la derivada de $2^{-x}$ es $-\ln2\cdot2^{-x}$). Siempre es positiva: nunca se divide entre $0$.
  3. Fórmula. $x_{n+1}=x_n-\dfrac{x_n-2^{-x_n}}{1+\ln2\cdot2^{-x_n}}$.
  4. Tabla (tres decimales en cada fila) y parada en la primera fila con $|x_{n+1}-x_n|<10^{-3}$:
$n$$x_n$$2^{-x_n}$$f(x_n)$$f'(x_n)$$x_{n+1}$$|x_{n+1}-x_n|$
001$-1$1.6930.5910.591
10.5910.664$-0.073$1.4600.6410.050
20.6410.641$\approx-0.0001$1.4440.641$0<10^{-3}$

Resultado: raíz $\approx0.641$ (se cumple el criterio en la tercera fila). Sin redondear sale lo mismo: $x_3=0.64119$.

Lo que escribes en el examen

«$f$ continua con cambio de signo en $[0,1]$; $f'\ne0$; fórmula de Newton; tabla; como $|x_2-x_1|=0.050$ y $|x_3-x_2|=0<10^{-3}$, la raíz es $0.641$.»

Receta: Newton «hasta $|x_n-x_{n-1}|<$ tol»: tabla con una columna para la diferencia entre iterados y parar en la primera fila que cumpla la tolerancia.
Ojo

Este criterio de parada es práctico: no garantiza que el error real sea $<10^{-3}$. La solución impresa tiene una errata de escritura que no afecta al resultado.

Parte del tema: §1, «Cuándo parar» (U pp.89-90).

X5 (medio) — punto fijo completo: $\frac12\operatorname{sen}x-x+1=0$ (convergencia, 4 iteraciones, error) X p.141

Dada $\frac12\operatorname{sen}x-x+1=0$: (a) compruebe que el punto fijo converge para todo $0\le x_0\le\pi$ (1 p.); (b) aproxime la solución con $x_0=2$ y cuatro iteraciones (1 p.); (c) estime el error absoluto cometido (1 p.).

La idea en una frase

El Teorema 2.2 funciona como una lista de tres comprobaciones sobre $g$ en un intervalo; si pasan, la iteración converge. Después se calcula (b) y se pone la cota de error (c).

  1. Despejar $x=g(x)$. $g(x)=\frac12\operatorname{sen}x+1$, intervalo $[a,b]=[0,\pi]$.
  2. (H0) $g$ es continua en todo $\mathbb R$.
  3. (H1) $g$ no se sale del intervalo. Para $x$ en $[0,\pi]$, $\operatorname{sen}x$ está entre $0$ y $1$, así que $g(x)$ está entre $1$ y $\frac32$. Eso cabe dentro de $[0,\pi]$ (porque $\pi\approx3.14$).
  4. (H2) $g$ es «contractiva». $g'(x)=\frac12\cos x$ y $|\cos x|\le1$, así que $|g'|\le\frac12=k<1$.
  5. (a) Conclusión. Por el Teorema 2.2 (U p.91) hay un único punto fijo $\tilde x$ en $[0,\pi]$ y $x_{n+1}=g(x_n)$ converge a él desde cualquier $x_0\in[0,\pi]$.
  6. (b) Iterar (calculadora en radianes, $x_0=2$): $x_1=g(2)\approx1.454649$, $x_2\approx1.496631$, $x_3\approx1.498626$, $x_4\approx1.498698$.
  7. (c) Error. Teorema 2.3: $|x_4-\tilde x|<\dfrac{k^4}{1-k}(b-a)=\dfrac{(1/2)^4}{1/2}\cdot\pi=\dfrac\pi8\approx0.3927$. (Es pesimista: el error real es $2.7\cdot10^{-6}$.)

Resultado: (a) converge; (b) $x_4\approx1.4987$; (c) error $<\pi/8\approx0.39$.

Lo que escribes en el examen

Las tres líneas (H0), (H1), (H2) con su cuenta, «Teorema 2.2», la lista de iterados y la cota $\frac{k^n}{1-k}(b-a)$ con los números sustituidos.

Receta: despeja $x=g(x)$, comprueba las tres hipótesis (continuidad, $g([a,b])\subset[a,b]$, $|g'|\le k<1$) y cita el Teorema 2.2; para el error, sustituye en $\frac{k^n}{1-k}(b-a)$.
Aviso sobre la solución impresa

X omite (H1), que en el examen sí hay que escribir; además tiene dos erratas menores ($g(1)$ por $g(2)$ y la página del teorema).

Parte del tema: §2, Teoremas 2.2 y 2.3; mismo esquema que el Ejemplo 2.17 de U (pp.92-93) y el Ej. 2.38 de E.

X6 (difícil) — punto fijo: número de iteraciones para error $\le0.1$ en $\cos(x^2)\operatorname{sen}(x+1)-2x=0$ X pp.29-30

Aplíquese el método del punto fijo para encontrar una solución aproximada, con 5 cifras significativas y error $\le0.1$, de $\cos x^2\operatorname{sen}(x+1)-2x=0$ con $x_0=0$. ¿Cuántas iteraciones hay que hacer? (La función $f(x)=\frac12\cos x^2\operatorname{sen}(x+1)$ es contractiva en $[0,1]$.) Aquí $\cos x^2$ significa $\cos(x^2)$.

La idea en una frase

La cota del Teorema 2.3 es $\frac{k^n}{1-k}(b-a)$. Cuanto más pequeño sea el intervalo y más pequeña la pendiente $k$, menos iteraciones hacen falta. Así que se busca un intervalo pequeño donde $g$ se cumpla, se saca $k$ y se despeja $n$.

  1. Despejar $x=g(x)$. $g(x)=\frac12\cos(x^2)\operatorname{sen}(x+1)$.
  2. Elegir el intervalo. Para $x\in[0,1]$ se cumple $0<g(x)\le\frac12$. Entonces $g$ lleva $[0,\frac12]$ dentro de sí mismo, y contiene a $x_0=0$. Tomamos $[a,b]=[0,\frac12]$.
  3. (H0) y (H1). $g$ es continua y $g([0,\frac12])\subset[0,\frac12]$ (paso anterior).
  4. Derivar. $g'(x)=-x\operatorname{sen}(x^2)\operatorname{sen}(x+1)+\frac12\cos(x^2)\cos(x+1)$.
  5. (H2) Acotar cada sumando en $(0,\frac12)$. El primero es $\le\frac12\operatorname{sen}\frac14\approx0.124$ y el segundo $\le\frac12\cos1\approx0.270$. Luego $|g'|\le0.394<\frac12=:k$.
  6. Despejar $n$ (Teorema 2.3): $\frac{(1/2)^n}{1/2}\cdot\frac12=\frac1{2^n}\le0.1\iff2^n\ge10$. Como $2^3=8<10\le16$: $n=4$ (de hecho error $<0.0625$).
  7. Iterar (5 cifras): $x_1=0.42074$, $x_2=0.48666$, $x_3=0.48432$, $x_4=0.48449$ ($x_5=0.48448$).

Resultado: $n=4$ iteraciones; solución aproximada $0.48449$ (punto fijo real $0.4844793\ldots$).

Lo que escribes en el examen

$g$, el intervalo $[0,\frac12]$ con (H0), (H1), (H2) y $k=\frac12$; la desigualdad $\frac1{2^n}\le0.1$ con $n=4$; y la lista $x_1,\dots,x_4$.

Receta: para «¿cuántas iteraciones garantizan error $<\varepsilon$?», elige el intervalo $[a,b]$ más pequeño que contenga a $x_0$ y que $g$ lleve en sí mismo, acota $|g'|$ por $k<1$ y toma el menor natural con $\frac{k^n}{1-k}(b-a)<\varepsilon$.
Aviso sobre la solución impresa
  • X trabaja en $[0,1]$ con $k=\frac12$ y obtiene $n=5$ (calcula hasta $x_5=0.48448$). El $5$ también garantiza el error, pero la justificación es incorrecta.
  • Falta el factor $2x$ de la regla de la cadena en $(\cos x^2)'=-2x\operatorname{sen}x^2$. Con la derivada bien, $|g'(1)|\approx0.878>\frac12$: $k=\frac12$ no vale en $[0,1]$ (haría falta $k=0.88$ y $35$ iteraciones).
  • X afirma $g([0,1])\subset[0,1]$ sin comprobarlo; hay que escribirlo (paso 2).

Parte del tema: §2 y §3, «Lista de pasos para garantizar error $<\varepsilon$» (U pp.91-92).

Tipos que también caen y no están en el PDF resuelto
  • Teoría (pregunta corta): «describa el método de Newton: para qué sirve y cómo se usa» (X p.129); «condiciones para que $g:[a,b]\to[a,b]$ tenga un único punto fijo» (Febrero 2014, Modelo B, P1); «describa el método del punto fijo y sus condiciones» e «importancia de la cota de error» (X p.9, sin solución impresa). Se responden con los Teoremas 2.2-2.3 y la deducción geométrica de Newton (recta tangente).
  • Newton para raíces cúbicas: $\sqrt[3]9$ con $x_0=1$ (Septiembre 2018, Modelo A, X pp.96-97; $x_1=3.6667$, $x_2=2.6676$) y $\sqrt[3]5$ con $x_0=2$ (Febrero 2013, X p.32): como X2 y el Ej. 2.34.
  • Newton dentro de un ejercicio largo: función a trozos (continuidad y derivabilidad) y después Newton para $f(x)=1$ (Febrero 2017, Modelo B, X pp.82-83); una iteración para $x^5-x^4+2x^3-1=0$ (Febrero 2016, Modelo A).
  • Punto fijo dentro de un estudio de funciones: $g(x)=1/(x^2+4)$, 3 iteraciones en $[0,1]$ (Febrero 2018, Modelo A, X pp.89-90).
  • Bisección con la ecuación explícita: $x^3+2x^2-1=0$ en $(0,1)$ con error $<0.3$ o $<0.2$ (X p.10, pp.17-18, pp.116-117): como X1, escribiendo además los intervalos sucesivos.
  • No aparece ninguna pregunta de número de iteraciones de Newton garantizado (U no da esa cota) ni de elegir, entre varias $g$, la que converge (Ejemplo 2.18 de U); esta última podría caer como teoría.

7Chuleta y lista de comprobación

Chuleta
  • Newton (para $f(x)=0$, $f$ derivable, $f'(x_n)\ne0$): $\;x_{n+1}=x_n-\dfrac{f(x_n)}{f'(x_n)}$. Raíz de $x^2-a$: $\tfrac12(x_n+a/x_n)$.
  • Errores aproximados: $\varepsilon_a=|x_n-x_{n-1}|$, $\;\varepsilon_r=\dfrac{|x_n-x_{n-1}|}{|x_n|}$. Parar cuando sea menor que la tolerancia, o por número máximo de iteraciones.
  • Punto fijo: $g(\tilde x)=\tilde x$; iteración $x_{n+1}=g(x_n)$ (2.3).
  • Teorema 2.2: $g$ continua, $g:[a,b]\to[a,b]$ y $|g'|\le k<1$ en $(a,b)$ $\Rightarrow$ punto fijo único en $[a,b]$ y convergencia desde todo $x_0\in[a,b]$.
  • Teorema 2.3: $|x_n-\tilde x|<\dfrac{k^n}{1-k}(b-a)$. Número de iteraciones: menor $n$ con $\dfrac{k^n}{1-k}(b-a)<\varepsilon$, es decir $n>\dfrac{\ln(\varepsilon(1-k)/(b-a))}{\ln k}$.
  • Bisección: error $<\dfrac{b-a}{2^n}$; $n>\log_2\dfrac{b-a}{\varepsilon}$.
  • Fallos de Newton: $f'(x_n)=0$, ciclos, mal $x_0$; raíz doble $\Rightarrow$ convergencia lenta.
U pp.54,88-92

Lista de comprobación: «sé hacer…»

Errores frecuentes
  • Olvidar que $g$ debe llevar $[a,b]$ en sí mismo (solo mirar $|g'|<1$).
  • Perder el signo de $f'$ (por ejemplo con $e^{-x}$).
  • No invertir la desigualdad al dividir entre $\ln k<0$; redondear el $n$ hacia abajo.
  • Redondear demasiado pronto y creer que el método «ya converge».
  • Confundir la garantía del Teorema 2.3 con las iteraciones que realmente hacen falta.
  • Empezar Newton en un $x_0$ con $f'(x_0)=0$ o sin haber localizado la raíz.
  • Calcular en grados en vez de en radianes (ejemplos con $\cos x$).

★Fuentes de este tema

Mostrar lista de fuentes

Verificadas contra el texto de los libros; todas las páginas son las impresas. U §2.4 ocupa las pp. 86-94 (§2.5 empieza en la p. 94). Se han usado L §3.8 (pp. 225-229) y S §4.8 (pp. 338-343).