MÉTODOS NUMÉRICOS

Métodos numéricos Versión: 22 de septiembre de 2019 La mayor parte

de las matemáticas estudiadas hasta ahora se han dedicado a
desarrollar métodos que nos proporcionen la solución exacta de un
problema. Por ejemplo, calcular la solución de una ecuación del tipo
f(x)=0 realizando operaciones elementales sobre la misma para
conseguir despejar la incógnita x. O bien, calcular el valor de una
integral definida Z b a f(x) dx calculando una primitiva de f(x) y luego
aplicando la Fórmula de Barrow. Desgraciadamente, en la gran mayoría
de los casos que se presentan en la práctica, estos métodos no son de
aplicación. Ello puede deberse a que el método para calcular la solución
exacta sea muy complicado, a que no se conozca un método adecuado,
o incluso a que no exista un método que nos permita, mediante
cálculos elementales, encontrar la solución. En estos casos es necesario
recurrir a métodos numéricos, denominados así porque, usualmente,
consisten en realizar una sucesión más o menos larga de operaciones
numéricas (normalmente mediante la ayuda de un ordenador), al cabo
de las cuales encontramos un valor numérico que, si bien no es la
solución exacta del problema, se le parece mucho, es decir, aproxima
la solución buscada con una precisión razonablemente buena. 4.1
Resolución numérica de ecuaciones Uno de los problemas que más se
presenta en matemáticas es el de calcular la solución de una ecuación.
En algunas (pocas) ocasiones, esto puede hacerse por métodos
analíticos, es decir, se puede “despejar” la incógnita para encontrar el
o los valores que resuelven la ecuación. En la gran mayoría de las
ocasiones con algún interés práctico esto no es posible y es necesario
recurrir a un método numérico que, con la ayuda de un ordenador, nos
permita calcular un valor aproximado de la solución. 4.1.1 Teoremas
del Valor Intermedio y de Bolzano Cuando se plantea el problema de
calcular la solución de una ecuación como f(x)=0 existen dos cuestiones
previas que conviene analizar: ¿Tiene solución esta ecuación? ¿Dónde
está (aunque sea más o menos) la solución? En los casos en que la
solución se puede calcular exactamente por métodos elementales, se
tiene la respuesta a ambas preguntas: existe, puesto que la hemos
encontrado y sabemos dónde está, puesto que sabemos exactamente
su valor. 159 4. Métodos numéricos 160 En muchos de los otros casos,
la respuesta a estas preguntas se obtiene con ayuda de los siguientes
teoremas. Teorema del Valor Intermedio Una función continua en un
intervalo [a, b] toma todos los valores comprendidos entre f(a) y f(b).
Figura 4.1: Teorema del Valor Intermedio: como se puede observar, la
función toma todos los valores comprendidos entre f(a) y f(b). Figura
4.2: Teorema del Valor Intermedio: además de todos los valores
comprendidos entre f(a) y f(b) la función f puede tomar otros valores.
Teorema de Bolzano Sea f una función continua en un intervalo [a, b] y
tal que f(a) y f(b) tienen signos opuestos (es decir f(a)f(b) < 0). Entonces
existe c 2 (a, b) tal que f(c)=0. Figura 4.3: Teorema de Bolzano: como
se puede observar, la función tiene signos opuestos en a y b (f(a) < 0 y
f(b) > 0). En concecuencia, toma el valor 0 en algún punto del intervalo
(a, b) (de hecho lo toma en tres puntos). Matemáticas Aplicadas a la
Biología - Grado en Biología R. Echevarría - Dpto. EDAN - Univ. de Sevilla
4. Métodos numéricos 161 Ejemplo 4.1 Utilizando el Teorema de
Bolzano, probar que la ecuación x = 2x tiene al menos una solución real.
En primer lugar, hay que escribir la ecuación en la forma f(x)=0 y, luego,
encontrar un intervalo [a, b] en el cual se verifiquen las hipótesis del
Teorema, para así poder concluir que existe algún punto en el intervalo
en el que la función se anula, es decir, alguna solución de la ecuación.
Se tiene: x = 2x () f(x) = x 2x = 0 Esta función está definida y es continua
en todo R. Es fácil ver que f(0) = 0 20 = 1 < 0. Por otro lado, teniendo
en cuenta que cuando x tiende a +1, l´ım x!+1 2x = 0, tampoco es difícil
comprender que, para x suficientemente grande, x será mayor que 2x
y por tanto x 2x será positivo. Por ejemplo: f(1) = 1 21 = 1 1 2 = 1 2 > 0
En consecuencia, f verifica las hipótesis del Teorema de Bolzano en el
intervalo [0, 1]: es continua y f(0) y f(1) tienen signos opuestos. Luego
podemos afirmar que f(x) tiene al menos un cero en el intervalo (0, 1).
O, lo que es lo mismo, que la ecuación x = 2x tiene al menos una
solución en dicho intervalo. Ejemplo 4.2 Utilizando el Teorema de
Bolzano, probar que la ecuación x4 = 1+3ex tiene al menos una raíz
real. Razonando como en el ejercicio anterior, se tiene x4 =1+3ex () f(x)
= x4 1 3ex = 0 La función f(x) está definida y es continua en todo R. Se
tiene, por ejemplo, f(0) = 1 3 = 4 < 0. Por otro lado, igual que en el
ejemplo anterior, x4 1 tiende a +1 cuando x ! +1 mientras que l´ım x!+1
3ex = 0, y no resulta difícil comprender que, para x suficientemente
grande, x4 1 será mayor que 3ex y por tanto x4 1 3ex será positivo. Por
ejemplo,

recordando que ex < 1 8x > 0 y, en consecuencia, que 3ex < 3 8x > 0, se

tiene f(2) = 24 1 3e2 = 15 3e2 > 12 > 0 Luego, por el Teorema de
Bolzano, f(x) tiene, al menos, un cero en el intervalo (0, 2), es decir, la
ecuación dada tiene, al menos, una raíz en dicho intervalo. 4.1.2
Resolución numérica de ecuaciones: método de bisección Se presenta
aquí un método sencillo, basado directamente en el

Teorema de Bolzano, que permite, en determinadas circunstancias,

calcular la solución de una ecuación. Hay que comenzar por decir que
cualquier ecuación en una variable se puede siempre escribir (y no de
manera única) en la forma de una equivalente (es decir, que tiene las
mismas soluciones) pero con segundo miembro nulo f(x)=0
Matemáticas Aplicadas a la Biología - Grado en Biología R.
existe c 2 (a, b) tal que f(c)=0. Figura 4.3: Teorema de Bolzano: como
se puede observar, la función tiene signos opuestos en a y b (f(a) < 0 y
f(b) > 0). En concecuencia, toma el valor 0 en algún punto del intervalo
(a, b) (de hecho lo toma en tres puntos). Matemáticas Aplicadas a la
Biología - Grado en Biología R. Echevarría - Dpto. EDAN - Univ. de Sevilla
4. Métodos numéricos 161 Ejemplo 4.1 Utilizando el Teorema de
Bolzano, probar que la ecuación x = 2x tiene al menos una solución real.
En primer lugar, hay que escribir la ecuación en la forma f(x)=0 y, luego,
encontrar un intervalo [a, b] en el cual se verifiquen las hipótesis del
Teorema, para así poder concluir que existe algún punto en el intervalo
en el que la función se anula, es decir, alguna solución de la ecuación.
Se tiene: x = 2x () f(x) = x 2x = 0 Esta función está definida y es continua
en todo R. Es fácil ver que f(0) = 0 20 = 1 < 0. Por otro lado, teniendo
en cuenta que cuando x tiende a +1, l´ım x!+1 2x = 0, tampoco es difícil
comprender que, para x suficientemente grande, x será mayor que 2x
y por tanto x 2x será positivo. Por ejemplo: f(1) = 1 21 = 1 1 2 = 1 2 > 0
En consecuencia, f verifica las hipótesis del Teorema de Bolzano en el
Los métodos para aproximar raíces de ecuaciones son, en general
iterativos, es decir consisten en construir una sucesión de valores x1,
x2, x3, x4 . . . mediante una relación de recurrencia, esto es, se calcula
cada uno de ellos a partir del anterior: x1 ! x2 ! x3 ! x4, etc. Cuando la
sucesión de valores x1, x2, x3 . . . tiende hacia la raíz ↵ de f (es decir,
se acerca cada vez más a ella, tanto como se quiera: l´ımn!1 xn = ↵),
se dice que el método iterativo es convergente. Método de bisección
Sin mucha precisión, el método de bisección consiste en lo siguiente:
1. Subdividir en dos partes el intervalo en que se sabe que la función
cambia de signo y tiene una sola raíz. 2. Averiguar, utilizando el
Teorema de Bolzano, en cual de las dos mitades se encuentra la raiz y
descartar la otra mitad del intervalo. 3. Reiniciar este proceso con el
subintervalo elegido. 4. Continuar con este proceso hasta que el
subintervalo elegido tenga una longitud lo suficientemente pequeña
como para que cualquiera de sus puntos sea una aproximación
aceptable de la solución. La elección óptima como aproximación es,
entonces, el punto medio del subintervalo. a b ↵ x1 a b ↵ x2 a b ↵ x3
Figura 4.5: Tres etapas del método de dicotomía. En cada iteración se
descarta la mitad del intervalo que no contiene a la raíz (en la que f
no cambia de signo). El intervalo donde se encuentra la raíz es cada
vez más pequeño y, su punto medio se acerca cada vez más a la
solución buscada. Matemáticas Aplicadas a la Biología - Grado en
Biología R. Echevarría - Dpto. EDAN - Univ. de Sevilla 4. Métodos
numéricos 164 Ejemplo 4.5 Utilizando el método de dicotomía,
aproximar la solución de la ecuación x 2x = 0 en el intervalo [0, 1] Sea
f(x) = x 2x. Intervalo Punto medio [0, 1] f(0) < 0 f(1) > 0 x0 = 0+1 2 =
0.5 [0.5, 1] f(0.5) < 0 x1 = 0.5+1 2 = 0.75 [0.5, 0.75] f(0.75) > 0 x2 =
0.5+0.75 2 = 0.625 [0.625, 0.75] f(0.625) < 0 x3 = 0.625 + 0.75 2 =
0.6875 [0.625, 0.6875] f(0.6875) > 0 x4 = 0.625 + 0.6875 2 = 0.65625 .
. . Por lo que una aproximación de la solución es ↵ ⇡ 0.65625,
obtenida aplicando el proceso de subdivisión 4 veces y eligiendo
como aproximación el punto medio del último subintervalo.
Obsérvese que, si se elige como aproximación x0, el error máximo
que se comete es la mitad de la longitud del intervalo inicial e0 = b a 2
. Si se elige como aproximación x1, el error máximo es la mitad del
anterior e1 = e0 2 = b

a 22 . Reiterando este razonamiento, si se elige como aproximación xn,

el error máximo es en = b a 2n+1 . Esto permite saber, a priori, cuantas
iteraciones hay que realizar para conseguir una aproximación con un
error tan pequeño como se quiera. En efecto, si en el intervalo [a, b]

hay una solución ↵, ¿qué número n de veces hay que aplicar el proceso
de subdivisión para conseguir que el error cometido no sea mayor que
una cantidad dada "? Se ha visto que, si se aplica n veces, el error
máximo que se comete tomando xn como aproximación es en = b a
2n+1 En consecuencia habrá que elegir n de forma que se tenga b a
2n+1 < " , b a " < 2n+1 , ln ✓b a " ◆ < (n + 1) ln(2) , n + 1 > ln ✓b a " ◆
ln(2) Ejemplo 4.6 ¿Cuántas iteraciones del método de bisección hay
que realizar para aproximar la solución de la ecuación x 2x = 0,
partiendo del intervalo [0, 1] , con un error menor que una centésima?
Se desea que el error sea e < 0.01. Por la fórmula anterior,
hay que tomar n + 1 > ln ✓b a " ◆ ln(2) = ln ✓ 1 0.01◆ ln(2) = ln(100)
ln(2) ⇡ 6.64 () n > 6.64 1=5.64 Luego hay que realizar 6 iteraciones. Si
se hicieran estas 6 iteraciones se obtendría como aproximación x6 =
0.6484, cuyas dos primeras cifras decimales son exactas, y por lo tanto,

el error de aproximación es menor que 0.01. Matemáticas Aplicadas a

la Biología - Grado en Biología R. Echevarría - Dpto. EDAN - Univ. de
Sevilla 4. Métodos numéricos 165 Ejemplo 4.7 Utilizando el método de
dicotomía, aproximar la solución de la ecuación del Ejercicio 4.2, x4 =1+3ex, en el
interval