Nota: Conviene revisar el capítulo 3, página 62 del Libro de Shoichiro Nakamura
para complementar los conceptos. El libro está en la plataforma Paideia, en la
carpeta Bibliografía.
Capítulo 2 (Continuación)
2.2 Método de Bisección
Si f es una función continua sobre el intervalo a0 , b0 y si f ( a0 ) . f ( b0 ) 0 , la función
f cambia de signo en el intervalo a0 , b0 y, por lo tanto, tiene por lo menos un cero en el
intervalo. Esto como consecuencia del Teorema del Valor Intermedio para funciones
continuas.
El método de bisección hace uso de esta idea de la siguiente forma: si f ( a0 ) . f ( b0 ) 0,
entonces calculamos x0 = ( a0 + b0 ) / 2 y averiguamos si f ( a0 ) . f ( x0 ) 0 . Si lo es
entonces f ( x) tiene un cero en a0 , x0 . A continuación, redefinimos:
a1 = a0 y b1 = x0 y comenzamos con el nuevo intervalo a1 , b1 , cuya longitud es
igual a la mitad del intervalo original. Si f ( a0 ) f ( x1 ) 0 entonces f ( x1 ) f ( b0 ) 0,
y en este caso redefinimos:
a1 = x0 , y b1 = b0 .
En ambos casos se ha generado un nuevo intervalo que contiene un cero de f, y el proceso
puede repetirse. Claro está que si f ( a0 ) f ( x0 ) = 0 , entonces f ( x0 ) = 0 y con ello se
ha encontrado un cero de f ( x ) . Sin embargo, por los errores de redondeo o por el
proceso mismo, es poco factible que en la computadora f ( x0 ) sea cero. Así, el criterio
para concluir no deberá depender de que f ( x0 ) = 0 . Se debe permitir una tolerancia
razonable, tal como f ( x0 ) 10− m , para algún entero m positivo El método de bisección
se conoce también como método de Bipartición.
A continuación, se ilustra el método, para la función:
f ( x ) = e x − sen x , x −8, − 4
Página 1 de 5
Elizabeth Doig
Entonces al calcular al punto medio, se tiene:
x0 = (−8 − 4) / 2 = −6
f ( x ) = e x − sen x , x −8, − 4
f (−8) f (−4) = -0.730876
f (−8) = 0.989694
f (−6) = -0.276937
f (−4) = -0.738487
Notar que :
f (a0 ) f ( x0 ) = f ( −8) f ( −6) 0
Luego, se define:
a1 = −8; b1 = −6
Observar los valores de las imágenes de la función en la gráfica:
1.0
0.5
7 6 5 4
0.5
1.0
Gráfico 1
Es decir, como f(a).f(x0) < 0 entonces: a1 = a, b1 = x0 y como f(x1)f(b1) < 0 entonces a2 =
x1 y b2 = b1, con x1 = (a1 + b1 ) / 2 punto medio del nuevo intervalo.
Teorema
Si 𝑓 es una función continua definida sobre [a, b] y tal que f(a).f(b) < 0, entonces la
sucesión generada por el método de Bisección {xn}(n0) converge a un cero, , de f en [a,
b]; esto es lim x n = y, además:
n→
xn - < (𝑏 − 𝑎) / 2n+1, n 0
Página 2 de 5
Elizabeth Doig
Algoritmo de Bisección
Para encontrar una solución de f(x) = 0 dada la función continua en el intervalo [a, b]
donde f(a) y f(b) tienen signos opuestos, hacer:
ENTRADA: extremos a, b : tolerancia TOL ; máximo número de iteraciones N0.
SALIDA: solución aproximada x* o mensaje de fracaso.
Paso 1 : Tomar i=1.
Paso 2 : Mientras que i N0 seguir pasos 3-6
Paso 3 : Tomar x* = a + (b-a)/2. (Calcular xi )
Paso 4 : Si f(x*) = 0 ó (b-a)/2< TOL entonces
SALIDA(x*); (procedimiento satisfactoriamente completado)
PARAR.
Paso 5 : Tomar i=i+1
Paso 6 : Si f(a)f(x*) > 0 entonces tomar a = x* (calcular ai , bi ).
si no tomar b = x*
Paso 7 : SALIDA (‘El método fracasó después de N0 iteraciones,
N0 =‘, N0 );
(Procedimiento completado sin éxito).
PARAR.
Ejemplo 2:
En la ecuación:
f ( x) = 1 + xe − x − x = 0 ,
hallaremos la raíz ubicada en el intervalo 1.2, 1.4 1, 2 , haciendo uso del Método de
Bisección.
Página 3 de 5
Elizabeth Doig
i ai bi xi f(xi)
0 1.2 1.4 1.3 0.0543
1 1.3 1.4 1.35 -2.56x10-5
2 1.3 1.35 1.325 0.027
3 1.325 1.35 1.3375 0.013595
4 1.3375 1.35 1.34375 6.78x10-3
5 1.34375 1.35 1.346875 3.38x10-3
6 1.346875 1.35 1.3484375 1.67x10-3
7 1.3484375 1.35 1.34921875 8.16x10-4
8 1.34921875 1.35 1.349609375 4.004x10-4
9 1.349609375 1.35 1.349804688 1.87x10-4
10 1.349804688 1.35 1.349902344 8.086x10-5
11 1.349902344 1.35 1.349951172 2.761x10-5
12 1.349951172 1.35 1.349975586 9.81x10-7
13 1.349975586 1.35 1.349987793 -1.233x10-5
14 1.349975586 1.349987793 1.34998169 -5.68x10-6
15 1.349975586 1.34998169 1.349978638 -2.35x10-6
16 1.349975586 1.349978638 1.349977112 -6.8x10-7
17 1.349975586 1.349977112 1.349976349 1.48x10-7
18 1.349976349 1.349977112 1.349976731 -2.68x10-7
19 1.349976349 1.349976731 1.34997654 -5.95x10-8
20 1.349976349 1.34997654 1.349976445 4.406x10-8
21 1.349976349 1.349976540 1.349976444 4.5x10-8
22 1.349976444 1.349976540 1.349976492 -7x10-9
23 1.349976444 1.349976492 1.349976468 1.9x10-8
24 1.349976468 1.349976492 1.349976480 6x10-9
25 1.349976480 1.349976492 1.349976486 -10-9
Se observa una precisión aproximada de ocho decimales exactos en el resultado obtenido
para la raíz x25 = 1.349976492 , con una TOL = 10-9.
Observar la gráfica de esta función:
2 1 1 2
Gráfico 2
En donde se ubican las raíces.
Página 4 de 5
Elizabeth Doig
*Nota: Criterios de Tolerancia
Para detener el algoritmo se pueden utilizar cualesquiera de los siguientes criterios:
1. xn +1 − x n
xn+1− xn
2.
xn+1
3. f ( xn )
4. Fijar un número máximo de iteraciones n = N, para interrumpir el algoritmo.
2.3 Orden de Convergencia
Con el fin de medir la rapidez de convergencia se pueden dar los siguientes conceptos:
Definición:
Sea lim x n = y en = xn - , n 0.
n →
Si existen constantes positivas y k, tal que
e n +1 x n +1 −
lim = lim =
n → n → k ,
k xn −
en
se dice entonces que la convergencia de {xn} a es de orden k, con una constante de error
asintótico .
Casos Particulares:
La constante afecta a la rapidez de convergencia, pero no es tan importante como la
constante k, orden de convergencia,
• Si k=1 la convergencia es de orden 1 ó lineal y
en + 1
lim = ,
n→
en
lo que implica que para n suficientemente grande se satisface:
en+1 en.
• Si k=2 la convergencia es de orden 2 ó cuadrática y se cumple que:
en +1
lim 2 = , y para n suficientemente grande se satisface:
n → e
n
en+1 en2
por lo tanto, la convergencia es mucho más rápida.
Página 5 de 5
Elizabeth Doig