NOTA TÉCNICA 2
MÉTODOS NUMÉRICOS
Cálculo de raíces de ecuaciones
El objeto del cálculo de las raíces de una ecuación es determinar los valores de 𝑥 para los que
se cumple:
𝑓(𝑥)
La determinación de las raíces de una ecuación es uno de los problemas más antiguos en
matemáticas y se han realizado un gran número de esfuerzos en este sentido. Su importancia
radica en que si podemos determinar las raíces de una ecuación también podemos determinar
máximos y mínimos, valores propios de matrices, resolver sistemas de ecuaciones lineales y
diferenciales, etc...
La determinación de las soluciones de la ecuación 𝑓(𝑥) puede llegar a ser un problema muy
difícil. Si 𝑓(𝑥) es una función polinómica de grado 1 ó 2, conocemos expresiones simples que
nos permitirán determinar sus raíces. Para polinomios de grado 3 o 4 es necesario emplear
métodos complejos y laboriosos. Sin embargo, si 𝑓(𝑥) es de grado mayor de cuatro o bien
no es polinómica, no hay ninguna fórmula conocida que permita determinar los ceros de la
ecuación (excepto en casos muy particulares).
Existen una serie de reglas que pueden ayudar a determinar las raíces de una ecuación:
El teorema de Bolzano, que establece que si una función continua, 𝑓(𝑥), toma en los extremos
del intervalo [𝑎, 𝑏] valores de signo opuesto, entonces la función admite, al menos, una raíz
en dicho intervalo.
En el caso en que 𝑓(𝑥) sea una función algebraica (polinómica) de grado n y coeficientes
reales, podemos afirmar que tendrá 𝑛 raíces reales o complejas.
La propiedad más importante que verifican las raíces racionales de una ecuación algebraica
establece que si 𝑝/𝑞 es una raíz racional de la ecuación de coeficientes enteros:
𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + ⋯ + 𝑎𝑛 𝑥 𝑛 = 0 (𝑎𝑖 ∈ ℵ)
entonces el denominador 𝑞 divide al coeficientes 𝑎𝑛 y el numerador 𝑝 divide al término
independiente 𝑎0 .
Ejemplo: Pretendemos calcular las raíces racionales de la ecuación:
3𝑥 3 + 3𝑥 2 − 𝑥 − 1 = 0
𝑦
Primero es necesario efectuar un cambio de variable 𝑥 = 3:
𝑦3 𝑦2 𝑦
3 3+3 2− +1=0
3 3 3
y después multiplicamos por 32 :
𝑦 3 + 3𝑦 2 − 3𝑦 − 9 = 0
con lo que los candidatos a raíz del polinomio son:
𝑦 = 9; 𝑦 = −9;
𝑦 = 3; 𝑦 = −3;
𝑦 = 1; 𝑦 = −1;
Este método , que se utiliza para resolver ecuaciones de una variable, está basado en el
“Teorema de los Valores Intermedios” (TVM), en el cual se establece que toda función
continua 𝑓, en un intervalo cerrado [𝑎, 𝑏] , toma todos los valores que se hallan entre 𝑓(𝑎)
y 𝑓(𝑏), de tal forma que la ecuación 𝑓(𝑥) = 0 tiene una sola raíz que verifica 𝑓(𝑎). 𝑓(𝑏) < 0.
Es el método más elemental y antiguo para determinar las raíces de una ecuación. Está basado
directamente en el teorema de Bolzano explicado con anterioridad. Consiste en partir de un
intervalo [𝑥0 , 𝑥1 ] tal que 𝑓(𝑥0 )𝑓(𝑥1 ) < 0, por lo que sabemos que existe, al menos, una raíz
real. A partir de este punto se va reduciendo el intervalo sucesivamente hasta hacerlo tan
pequeño como exija la precisión que hayamos decidido emplear.
El método de bisección es un algoritmo de búsqueda de raíces que trabaja dividiendo el
intervalo a la mitad y seleccionando el subintervalo que tiene la raíz. Esto se logra llevar a cabo
a través de varias interacciones que son aplicadas en un intervalo para por medio de ello
encontrar la raíz de la función.
Este es uno de los métodos más sencillos de fácil intuición para resolver ecuaciones en una
variable, también conocido como método del intervalo medio, este se basa en el teorema del
valor intermedio, el cual establece que toda función continua 𝑓𝑥 en un intervalo cerrado
[𝑎, 𝑐] toma todos los valores que se hallan entre 𝑓(𝑎) y 𝑓(𝑏). Esto es que todo valor entre
𝑓(𝑎) y 𝑓(𝑏) es la imagen de al menos un valor del intervalo [𝑎, 𝑏]. En caso de que 𝑓(𝑎) y 𝑓(𝑏)
tengan signos opuestos, el valor cero sería un valor intermedio entre 𝑓(𝑗) y 𝑓(𝑒), por lo que
con certeza existe un 𝑝 en [𝑎, 𝑏] que cumple 𝑓(𝑝) = 0. De esta forma, se asegura la existencia
de al menos una solución de la ecuación 𝑓(𝑎) = 0.
El método consiste en lo siguiente:
Debe existir seguridad sobre la continuidad de la función f(x) en el intervalo [𝑎, 𝑏]
A continuación se verifica que 𝑓(𝑎). 𝑓(𝑏) < 0.
Se calcula el punto medio m del intervalo [𝑎, 𝑏] y se evalúa 𝑓(𝑚) si ese valor es igual a cero,
ya hemos encontrado la raíz buscada.
En caso de que no lo sea, verificamos si 𝑓(𝑚) tiene signo opuesto con 𝑓(𝑎) o con 𝑓(𝑏).
Se re define el intervalo [𝑎, 𝑏] como [𝑎, 𝑚] ó [𝑚, 𝑏] según se haya determinado en cuál de
estos intervalos ocurre un cambio de signo.
Con este nuevo intervalo se continúa sucesivamente encerrando la solución en un intervalo
cada vez más pequeño, hasta alcanzar la precisión deseada
El método de bisección es menos eficiente que el método de Newton, pero es mucho más
seguro para garantizar la convergencia.
Use las formulas
𝑥𝑎 + 𝑥𝑏
𝑥𝑟 =
2
𝑓(𝑥𝑎 ) × 𝑓(𝑥𝑟 )
𝑥𝑟 (𝐴𝑐𝑡𝑢𝑎𝑙)+𝑥𝑟 (𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟)
Error Porcentual 𝐸𝑝 = | (100)|
𝑥𝑟 (𝐴𝑐𝑡𝑢𝑎𝑙)
Ejemplo: 𝑥 4 + 3𝑥 3 − 2 con un error menor al 15%
Generar valores arbitrarios a X para encontrar el cambio de signo:
𝑥 𝑓(𝑥)
Cambio de Signo:
- -7.6
1.5 𝑥𝑎 = 0
-1 -4 𝑥𝑏 = 1
0 -2
1 2
1.5 13.18
Iteración 1:
0+1
𝑥𝑟 = = 0.5
2
𝑓(𝑥𝑟 ) = (0.5)4 + 3(0.5)3 − 2 = −1.56
𝑓(𝑥𝑎 ) × 𝑓(𝑥𝑟 ) = (−2) ∗ (−1.56) = 3.12 > 0
Iteración 2:
𝑥𝑎 = 0.5
𝑥𝑏 = 1
0.5 + 1
𝑥𝑟 = = 0.75
2
0.75 + 0.5
𝐸𝑝 = | (100) = 33%|
0.75
𝑓(𝑥𝑟 ) = (0.75)4 + 3(0.75)3 − 2 = −0.41
𝑓(𝑥𝑎 ) × 𝑓(𝑥𝑟 ) = (−1.56) ∗ (−0.41) = 0.63 > 0
Iteración 3:
𝑥𝑎 = 0.75
𝑥𝑏 = 1
0.75 + 1
𝑥𝑟 = = 0.87
2
0.87 + 0.75
𝐸𝑝 = | (100) = 13%|
0.87
Como el error ya es menor al solicitado nos indica que es la última iteración.
Resultado
𝑥𝑟 = 0.75
Método de Newton
Este método parte de una aproximación inicial 𝑥0 y obtiene una aproximación mejor, 𝑥1 , dada
por la fórmula:
𝑓(𝑥0 )
𝑥1 = 𝑥0 −
𝑓 ′ (𝑥0 )
La expresión anterior puede derivarse a partir de un desarrollo en serie de Taylor.
Efectivamente, sea 𝑟 r un cero de 𝑓 y sea 𝑥 una aproximación a 𝑟 tal que 𝑟 = 𝑥 + ℎ. Si 𝑓′
existe y es continua, por el teorema de Taylor tenemos:
0 = 𝑓(𝑟) = 𝑓(𝑥 + ℎ) = 𝑓(𝑥 + ℎ𝑓(𝑥) + 𝑂(ℎ2 )
en donde ℎ = 𝑟 − 𝑥. Si 𝑥 está próximo a 𝑟 (es decir es pequeña), es razonable ignorar el
término 𝑂(ℎ2 ):
0 = 𝑓(𝑥) + ℎ𝑓(𝑥)
por lo que obtenemos la siguiente expresión para ℎ:
𝑓(𝑥)
ℎ=−
𝑓′(𝑥)
𝑓(𝑥)
A partir de la ecuación ℎ = − 𝑓′(𝑥) y teniendo en cuenta que 𝑟 = 𝑥 + ℎ es fácil derivar la
𝑓(𝑥 )
ecuación 𝑥1 = 𝑥0 − 𝑓′ (𝑥0 ).
0
El método de Newton tiene una interpretación geométrica sencilla. De hecho, el método de
Newton consiste en una linealización de la función, es decir, 𝑓 se reemplaza por una recta tal
que contiene al punto (𝑥0 , 𝑓(𝑥0 )) y cuya pendiente coincide con la derivada de la función en
el punto, 𝑓′(𝑥0 ). La nueva aproximación a la raíz, 𝑥1 , se obtiene de la intersección de la función
linear con el eje 𝑋 de ordenadas.
𝑓(𝑥 )
Veamos cómo podemos obtener la ecuación 𝑥1 = 𝑥0 − 𝑓′ (𝑥0 ) a partir de lo dicho en el párrafo
0
anterior. La ecuación de la recta que pasa por el punto (𝑥0 , 𝑓(𝑥0 )) y de pendiente 𝑓′(𝑥0 ) es:
𝑦 − 𝑓(𝑥0 ) = 𝑓(𝑥0 )(𝑥 − 𝑥0 )
Dos situaciones en las que el método de Newton no funciona adecuadamente: (a) el método
no alcanza la convergencia y (b) el método converge hacia un punto que no es un cero de la
ecuación.
Ejemplo:
𝑥 3 − 4𝑥 2 − 2 = 0
𝑓(𝑥 )
Solución: 𝑥𝑛+1 = 𝑥𝑛 − 𝑓′ (𝑥𝑛 ) , 𝑓 ′ (𝑥𝑛 ) ≠ 0 (0)
𝑛
Sea 𝑓(𝑥) = 𝑥 3 − 4𝑥 2 − 2 (1)
→ 𝑓 ′ (𝑥) = 3𝑥 2 − 8𝑥 (2)
DE (0), (1) y (2), se deduce que la segunda aproximación está dada por:
𝑥 3 −4𝑥 2 −2
𝑥2 = 𝑥1 − (3)
3𝑥 2 −8𝑥
Gráfica de 𝑓
Fijándose en la gráfica de 𝑓, escogemos a
𝑥1 = 4.0000 (4)
Ahora 𝑓(4.000) = −2.000 (5)
𝑓′(4.000) = 16.000 (6)
Sustituyendo (4), (5) (6) en (3), se obtiene:
−2.0000
𝑥2 = 4.0000 − ,
16.0000
𝑥2 = 4.1250
Para este nuevo valor de 𝑥, calculamos:
𝑓(4.1250) = 0.1270
𝑓 ′ (4.1250) = 18.0469
De tal manera que:
0.1270
𝑥3 = 4.1250 − = 4.1180
18.0469
Para este nuevo valor de 𝑥, calculamos:
𝑓(4.1180) = 0.0010
𝑓 ′ (4.1180) = 17.9298
0.0010
𝑥4 = 4.1180 − = 4.1179
17.0469
Para este nuevo valor de 𝑥, calculamos:
𝑓(4.1179) = 0.0008
𝑓 ′ (4.1179) = 17.9281
0.0008
𝑥5 = 4.1179 − = 4.1179
17.9281
Comparando 𝑥4 = 𝑥5 = 4.1179, la respuesta de la raíz real de 𝑥 3 − 4𝑥 2 − 2 = 0, dada con
una aproximación de 4 cifras decimales es 4.1179.
Método de la secante
El principal inconveniente del método de Newton estriba en que requiere conocer el valor de
la primera derivada de la función en el punto. Sin embargo, la forma funcional de 𝑓(𝑥) dificulta
en ocasiones el cálculo de la derivada. En estos casos es más útil emplear el método de la
secante.
El método de la secante parte de dos puntos (y no sólo uno como el método de Newton) y
estima la tangente (es decir, la pendiente de la recta) por una aproximación de acuerdo con
la expresión:
𝑓(𝑥1 ) − 𝑓(𝑥0 )
𝑓 ′ (𝑥0 ) =
𝑥1 − 𝑥0
𝑓(𝑥 )
Sustituyendo esta expresión en la ecuación 𝑥1 = 𝑥0 − 𝑓′ (𝑥0 ) del método de Newton,
0
obtenemos la expresión del método de la secante que nos proporciona el siguiente punto de
iteración:
𝑥1 − 𝑥0
𝑥2 = 𝑥0 𝑓(𝑥0 )
𝑓(𝑥1 ) − 𝑓(𝑥0 )
CUADRO INFORMATIVO
UNIDAD LINKS BIBLIOGRAFÍA
1 Método de Bisección: Nieves Hurtado, A. (2015). Métodos numéricos:
aplicados a la ingeniería.. Grupo Editorial Patria.
[Link] [Link]
ch?v=mdG6gpzE54k
Flórez Calderón, T. (2016). Métodos numéricos
Método de la secante: que debes saber.. Editorial Universidad Nacional
[Link] de Colombia.
ch?v=kbyWO1i-OcM [Link]
Método de Newton:
García, I. (2013). Métodos numéricos: problemas
[Link] resueltos y prácticas.. Edicions de la Universitat
ch?v=n53t8CtaLrM de Lleida.
[Link]
