Ingeniería en Sistemas Computacionales

MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel
2 MÉTODOS DE SOLUCIÓN DE ECUACIONES

2.1 Métodos de Intervalo

Los métodos de los intervalos utilizan una propiedad muy importante, consistente en el hecho
del cambio de signo de una función en inmediaciones de una raíz.

Se llaman métodos de los intervalos porque se necesitan como mínimo dos valores que forman
un intervalo que encierra la raíz.

f(x)

+ f(x)
+ f(x)

d e f
x
a c b
- f(x)
- f(x)

Gráfica 2.1

En la gráfica anterior se observa como la función cambia de +f(x) a -f(x), cuando pasa por la
raíz c. Esto ocurre porque f (c)= 0 y necesariamente la función pasa del cuadrante positivo al
negativo de x. En algunos casos, que se verán más adelante esto no ocurre así, por ahora se
asumirá como se ha mostrado.
Los métodos abiertos utilizan estos cambios de signo para poder ubicar la raíz (punto c), pero
es necesario entonces establecer un intervalo (como el [a,b]).
De igual manera sucede cuando la función pasa por el punto e, el cambio ocurre de -f(x) a +f(x),
para hallar la raíz el método necesita un intervalo como el [d,f].
Los métodos de Intervalos que se verán son:
a. Método de Bisección
b. Método de Aproximaciones Sucesivas
c. Método de Interpolación Lineal

2.2 Método de Bisección

Si f es una función continua sobre el intervalo [a,b] y si f(a) f(b)<0, entonces f debe tener un cero
en (a,b). Dado que f(a)f(b)<0, la función cambia de signo en el intervalo [a,b] y por lo tanto tiene
por lo menos un cero en el intervalo.

1
Ingeniería en Sistemas Computacionales
MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel

Esta es una consecuencia del teorema del valor intermedio para funciones continuas, que
establece que si f es continua en [a,b] y si k es un número entre f(a) y f(b) , entonces existe por
lo menos un c (a,b) tal que f(c)=k. (para el caso en que f(a)f(b)<0 se escoge k=0, luego f(c)=0,
c (a,b)).

El método de bisección consiste en dividir el intervalo en 2 subintervalos de igual magnitud,

reteniendo el subintervalo en donde f cambia de signo, para conservar al menos una raíz o
cero, y repetir el proceso varias veces.

Por ejemplo, suponga que f tiene un cero en el intervalo [a,b].

Primero se calcula el punto medio del intervalo ; después se averigua sí f(a)f(c)<0. Si

lo es, entonces f tiene un cero en [a,c].

A continuación se renombra a c como b y se comienza una vez más con el nuevo intervalo [a,b],
cuya longitud es igual a la mitad del intervalo original.

Si f(a)f(c)>0 , entonces f(c)f(b)<0 y en este caso se renombra a c como a.

En ambos casos se ha generado un nuevo intervalo que contiene un cero de f, y el proceso

puede repetirse.

2
Ingeniería en Sistemas Computacionales
MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel

Para detener el método de bisección y dar una aproximación del cero de una función se
pueden usar varios criterios (llamados criterios de parada).
Uno de los criterios de parada consiste en examinar si |f(cn)| < , donde es una tolerancia
previamente establecida (por ejemplo = 10-3). Otro criterio que puede utilizarse es examinar sí

También se puede usar como criterio de parada el error relativo entre dos aproximaciones del

cero de f. .

Teorema. (Error en el método de bisección).

Si f es continua en [a, b] y f(a) f(b) < 0, el método de bisección genera una sucesión que aproxima

un cero c de f con la propiedad que : ,n 1.

Notas:

 El método de bisección tiene la desventaja que es lento en cuanto a convergencia (es

decir que se necesita un n grande para que sea pequeño). Otros métodos
requieren menos iteraciones para alcanzar la misma exactitud, pero entonces no
siempre se conoce una cota para la precisión.

 El método de bisección suele recomendarse para encontrar un valor aproximado del

cero de una función, y luego este valor se refina por medio de métodos más eficaces.
La razón es porque la mayoría de los otros métodos para encontrar ceros de funciones
requieren un valor inicial cerca de un cero; al carecer de dicho valor, pueden fallar por
completo.

3
Ingeniería en Sistemas Computacionales
MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel
 Resolver una ecuación en una variable como por ejemplo: xex=1 es equivalente a
resolver la ecuación xex-1=0, o a encontrar el cero de la función f(x) = xex-1. Para
aproximar el cero de f o la raíz de la ecuación se puede hacer la gráfica de f en una
calculadora o usar MathLab para determinar un intervalo donde f tenga un cero.
También se pueden ensayar números a y b de tal manera que f(a)f(b)<0. Para el caso
de f(x) =xex-1 por ejemplo f(0) = -1, f(1) = e-1 1.71828 entonces f tiene un cero en el
intervalo [0,1].

 Cuando hay raíces múltiples, el método de bisección quizá no sea válido, ya que la
función podría no cambiar de signo en puntos situados a cualquier lado de sus raíces.
Una gráfica es fundamental para aclarar la situación. En este caso sería posible hallar
los ceros o raíces trabajando con la derivada f’(x), que es cero en una raíz múltiple.

2.3 Método de Aproximaciones Sucesivas

Los métodos abiertos se basan en fórmulas que requieren de un solo valor para predecir la
raíz.

La ecuación f(x)= 0 se arregla de tal modo que x quede del lado izquierdo de la ecuación (esto
es, despejar una x de primer grado, o agregar una x en ambos lados de la ecuación para no
alterarla):

X = g (x )

Ejemplo:

f(x) = x2 −2x + 3 = 0

reordenando : (x = x2 + 3)/2

f(x) = sen x = 0

reordenando: x = sen x + x

Se requiere de un valor inicial xi, que se pueda usar para obtener una aproximación de xi + 1,
expresándolo por la forma iterativa:

El método de las aproximaciones sucesivas es uno de los procedimientos más importantes y

más sencillos de codificar. Se sabe que la ecuación

Donde f(x) es una función continua que se desea determinar sus raíces reales. Se sustituye f(x)
por la ecuación equivalente

Se estima el valor aproximado de la raíz x0, y se sustituye en el segundo miembro de la

ecuación para obtener x1.

4
Ingeniería en Sistemas Computacionales
MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel
Poniendo x1 como argumento de x), se obtendrá un nuevo número x2, y así sucesivamente.
Este proceso se puede sintetizar en la fórmula.

(1)

Si esta secuencia es convergente es decir, tiende hacia un límite, la solución  es:

El método de iteración se explica geométricamente mediante el gráfico de la figura. Se dibuja la

curva y=(x), y la recta y=x, bisectriz del primer cuadrante. La abscisa  del punto de
intersección es la raíz buscada.

Un ejemplo típico es la de encontrar la raíz de la ecuación:

Para encontrar la raíz, se comienza en el punto cualquiera de abscisa x0 dentro del intervalo (0,
/2), y se traza la línea vertical hasta que interseca la curva, luego, desde este punto, se traza
una línea horizontal hasta que se alcanza la recta bisectriz, este punto tendrá por abscisa x1. Se
traza de nuevo, una línea vertical hasta encontrar a la curva, y otra línea horizontal hasta
encontrar la línea recta, el punto de intersección tiene de abscisa x2, y así sucesivamente.
Como se puede apreciar en la figura, la sucesión x1, x2, x3... tiende hacia la raíz  de la ecuación
buscada.

La condición de finalización

Primero, se introduce el valor inicial x, la primera aproximación, se calcula el valor del coseno
de x, el valor devuelto (segunda aproximación), se guarda de nuevo en la variable x, y se repite
el proceso indefinidamente. El método aunque correcto, necesita terminarse en algún
momento, cumpliendo una determinada condición. Cuando el valor absoluto del cociente entre
la diferencia de dos términos consecutivos de la sucesión y uno de los términos, sea menor
que cierta cantidad 

5
Ingeniería en Sistemas Computacionales
MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel
Este criterio, no es completamente riguroso, pero es un buen punto de partida para el estudio
de este método.

El criterio de convergencia

No todas las ecuaciones pueden resolverse por este método, solamente si el valor absoluto de
la derivada de la función (x) en la vecindad de la raíz  es menor que la unidad (la pendiente
de la recta bisectriz del primer cuadrante es uno). En la figura, se aprecia como es imposible
encontrar la solución marcada por un puntito negro en la intersección entre la curva y la recta
bisectriz del primer cuadrante, ya que la sucesión xi diverge.

2.4 Método de Interpolación

Uno de los métodos de interpolación más sencillos es el lineal. En general, en la interpolación

lineal se utilizan dos puntos, (xa,ya) y (xb,yb), para obtener un tercer punto interpolado (x,y) a
partir de la siguiente fórmula:

La interpolación lineal es rápida y sencilla, pero no muy precisa.

La interpolación lineal es un caso particular de la interpolación general de Newton.

Con el polinomio de interpolación de Newton se logra aproximar un valor de la función f(x) en
un valor desconocido de x. El caso particular, para que una interpolación sea lineal es en el que
se utiliza un polinomio de interpolación de grado 1, y se denota de la siguiente manera:

En la gráfica se representan algunos valores de la función, pero no todos, en ocasiones nos

interesa el valor de la función para un valor de la variable independiente distinto de los que
figuran en la gráfica, en este caso podemos tomar el más próximo al buscado, o aproximarnos
un poco más por interpolación, la interpolación casi siempre nos dará un pequeño error
respecto al valor de la función verdadero, pero siempre será menor que tomar el valor más
próximo de los que figuran en la tabla, veamos como se calcula al valor de la función para un
6
Ingeniería en Sistemas Computacionales
MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel
valor de la variable independiente que se encuentre entre dos valores de la tabla por
interpolación lineal.

Por la gráfica se sabe que:

y
Lo que se quiere saber es:

Siendo:

La interpolación lineal consiste en trazar una recta que pasa por (x1,y1) y (x2,y2), y = r(x) y
calcular los valores intermedios según esta recta en lugar de la función y = f(x)

Para ello se usa la semejanza de triángulos y

esto es:

despejando, se tiene:

o lo que es lo mismo:

El valor buscado es:

7
Ingeniería en Sistemas Computacionales
MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel
ésto es:

Considerando que en una raíz f(x)=0, la línea atraviesa el eje x

Considerando a=xi , fa=f(xi), b=xu, fb=f(xu), x=xr fx=f(xr) por triángulos semejantes:

b-a = fb-fa
b-x fb-fx Para encontrar la raíz fx=0;

b-a = fb-fa Hay que despejar x; se pasa fb al primer término;

b-x fb

fb(b-a) = fb-fa Se intercambia (fb-fa) con (b-x);

b-x

fb(b-a) = b-x Se pasa b al primer término, para dejar a -x;

fb-fa

-b + fb(b-a) = -x Se multiplica por (-1) ambos términos.

fb-fa

x = b - fb(a-b)
fa-fb

Verificar la convergencia (fa*fb)<0, asegura los puntos de unión de la hipotenusa cruzan al eje
y.

8
Ingeniería en Sistemas Computacionales
MÉTODOS NUMÉERICOS
Unidad 2. Métodos de Solución de Ecuaciones
MTI. Ana Francisca Lule Rangel