Definición 2.2.
Sea f(x) una función continua en un intervalo [a, b], si existe a,b se debe tener f()=0
para un número real. Gráficamente se tendría. Ver figura ([Link]):
Figura [Link]
Nota.
Para que exista debe haber un cambio de signo en la funcion en el intervalo, o su
equivalente seria cuando corta el eje “x”
Teorema [Link]
Sea f(x) una función continua en un intervalo cerrado [a, b] si f (a) * f (b) 0
Entonces existe a, b tal que f()=0.
Este teorema dice que para que una función tenga por lo menos una raíz, se debe tener que
f(a) y f(b) debe tener signos contrarios (signos diferentes). Como se muestra en la figuras
([Link].1 a,b).
Figura ([Link].1 a).
Figura ([Link].1b).
Este teorema dice que si se cumple éste, se encuentra una raíz, pero lo que puede suceder el
siguiente caso. Ver figura [Link].2 en la cual se tiene que f (a) * f (b) 0 , según el teorema
no tendría ninguna raíz pero gráficamente se tienen dos.
Figura [Link].2
2.1 Método de la bisección.
Uno de los métodos más antiguos utilizado para el cálculo de raíces es el método de la
bisección desarrollado por Bolzano-Weierstrass. Este se puede desarrollar de la siguiente
manera:
Usando el teorema ([Link]), como su nombre lo indica bisección “es partir en dos, en este
caso partir en dos el intervalo”, esto quiere decir que se deben realizar las siguientes
operaciones.
1. Verificar que se cumplen las condiciones del teorema ([Link])
2. c=(a+b)/2 dividir el intervalo en dos sub intervalos [a,c] y [c,b]
3. Preguntar en qué sub-intervalo se encuentra, esto se hace realizando la siguiente pregunta:
4.
si f(a)*f(c)< 0 entonces
hacer b=c
de lo contrario
hacer a=c
(de esta manera se ha encontrado una sucesión de intervalos encajados de la forma [an ,
bn]).
5. Esta operación se sigue hasta que se tenga que an -bn donde es el error dado por
el usuario. Mientras no se cumpla ésto regresar al paso (2)
6. Por último hacer c=(a+b)/2 ( en este paso ya se tendría cn=(an+bn)/2)
7. Así que la aproximación a la raíz sería c.
Observe la figura (2.2.1)
Ejemplo 1.
Sea f(x)=x2-2 en el intervalo [0,3] y =0.1
Para resolver este problema primero se debe verificar que se cumple que f (a) * f (b) 0 . En
este caso a = 0 y b = 0 sustituyendo en la función.
f(a) * f(b) 0 se tiene f (0) * f (3) = (−2) * (7) = −14 0 . Esto indica que tiene una raíz en
el intervalo. Entonces sí se puede utiliza el método.
c=
(a + b ) 0+3
= =1.5 de esta manera se obtienen los siguientes sub-intervalos
2 2
[0,1.5] y [1.5, 3]
Ahora se pregunta en qué sub-intervalo se encuentra la raíz
si f(0)*f(1.5)=(-2)*(0.25)<0 y este producto si es menor que cero, entonces el sub- intervalo
en el que se encuentra la raíz es [0, 1.5]
Ahora se debe checar si se cumple ya el criterio de paro.
a1 − b1 Sustituyendo valores 0 − 1.5 = 1.5 0.1 no se cumple, así que se regresa al paso
(2)
c=
(a1 + b1 ) = 0 + 1.5 =0.5 aquí a1=0 y b1=1.5
2 2
en este caso los sub-intervalos son [0,0.75], [0.75,1.5]
Preguntaremos ahora en qué sub-intervalo se encuentra la raíz
si f(0)*f(0.57)<0 evaluando la función f(0)=-2, f(0.75)=-1.43
f(0)*f(0.57)=(-2)(-1.43)>0 este producto es mayor que cero, por lo tanto el sub-intervalo
donde se encuentra la raíz es [0.75,1.5 ]=[a2,b2] se probará si ya se cumple el criterio de paro.
a 2 − b2 Sustituyendo valores 0.75 − 1.5 = 0.75 0.1 no se cumple así que regresar
al paso (2).
c=
(a 2 + b 2 ) 0.75 + 1.5
= =1.125
2 2
los sub-intervalos son:
[0.75, 1.125], [1.125, 1.5]
se verifica en que sub-intervalo está
si f(0.75)*f(1.125)=(-1.43)*(10.7343)<0 no se cumple, por lo tanto está en el sub-intervalo
[1.125, 1.5] se verifica si ya se llegó a la tolerancia.
1.125 − 1.5 = 0.375 0.1 Aquí se ve que todavía no se cumple el criterio de paro.
Así que se regresa al punto (2).
Una manera de saber cuántas iteraciones hay que realizar para llegar a tener la tolerancia
dada es el siguiente:
b-a
lo que necesitamos conocer de esta ecuación es n
2n
b-a
2n Como >0 se puede hacer lo siguiente.
b-a
2n
Aplicando logaritmo natural a los dos miembros
b-a
ln ln 2n
b-a
ln (n ) ln 2
b-a
ln
n
ln 2
pero los valores de a=0, b=3, =0.1
3-0
ln
0.1 n
ln 2
ln 30
n
ln 2
3.40119738 2
n
0.69314718
4.906890596 n
realizando una operación de redondeo:
5n
Tomando n=6 se debe tener el resultado adecuado, para la tolerancia dada.
2.2 MÉTODO DE LA SECANTE.
Este método se deduce apartar del de la figura 2.4.1 de donde se tienen dos triángulos, los
cuales están formados por las siguientes coordenadas:
Para el triángulo (1) 1 son (x0,f(x0)), (x1,f(x0)), (x1,f(x1)) de la misma manera para el
triángulo (2) y 2 (x2,0), (x1,0), (x1,f(x1)), estas coordenadas son los vértices de los
triángulos.
Estos dos triángulos son semejantes (usando la congruencia ángulo, ángulo, ángulo).
Relacionando la bases y las alturas de estos triángulos, para encontrar la incógnita que es x2
:
f(x1 ) − f(x0 ) f(x1 )
= (2.4.1)
x1 - x 0 x1 - x 2
de la ecuación 2.4.1 se despejará la incógnita que es x2 :
f(x1 )(x1 - x 0 )
x1 - x 2 = (2.4.2)
f(x1 ) − f(x0 )
f(x1 )(x1 - x 0 )
x 2 = x1 − (2.5.3)
f(x1 ) − f(x0 )
se generalizará este procedimiento, tomando para x0 el valor del sub-índice de n=0, en x1
el valor del sub-índice es n=n+1, por último x2 el sub-índice será n=n+2, usando estos
cambios la ecuación (2.4.3) se expresará de la siguiente manera:
f(xn +1 )(x n +1 - x n )
x n + 2 = x n +1 − (2.4.4)
f(xn +1 ) − f(xn )
el criterio de paro para este método es el siguiente:
x n + 2 - x n +1 o f(xn + 2 ) (2.4.5)
Ejemplo:
Se realizará este procedimiento usando la siguiente función:
f(x)=x2-2 en el intervalo [0,3] ,=0.1
los valores de x0=0 , x1=3 para n=0 se tiene
x1 − x0 = 3 − 0 = 3
f ( x1 )( x1 − x0 )
x 2 = x1 −
f ( x1 ) − f ( x0 )
f (3)(3 − 0)
x2 = 3 −
f (3) − f (0)
(7)(3) 21
x2 = 3 − = 3− = .66
7 − (−2) 9
x2 − x1 = 0.66 − 3 = 2.33 =0.1
No se cumple la tolerancia se continúa el proceso.
Para n=1 se tiene:
f ( x2 )( x2 − x1 )
x3 = x2 −
f ( x2 ) − f ( x1 )
Los valores ahora son: x1= 3 x2=0.66
f (0.66)(0.66 − 3) (−1.5644 )(−2.34)
x3 = 0.66 − = 0.66 − = 1.08
f (0.66) − f (3) (−1.5644 ) − (7)
x3 − x2 = 1.08 − 0.66 = 0.42 =0.1
Aquí no se cumple la tolerancia así que se continúa el proceso
para n=2 se tiene:
f ( x3 )( x3 − x2 )
x4 = x3 −
f ( x3 ) − f ( x2 )
Los valores ahora son: x2= 0.66 x3=1.08
f (1.08)(1.08 − 0.66) (−0.8336)(0.42)
x4 = 1.08 − = 1.08 − = 1.55
f (1.08) − f (0.66) (−0.8336) − (−1.5644)
x4 − x3 = 1.55 − 1.08 = 0.47 =0.1
aquí no se cumple la tolerancia así que se continúa el proceso
Para n=3 se tiene:
Los valores ahora son: x3= 1.08 x4=1.55
f ( x 4 )( x 4 − x3 )
x5 = x 4 −
f ( x 4 ) − f ( x3 )
f (1.55)(1.55 − 1.08) (0.4025)(0.47)
x5 = 1.55 − = 1.55 − = 1.396
f (1.55) − f (1.08) (0.4025) − (−0.8336 )
x5 − x4 = 1.396 − 1.55 = 0.154
aquí no se cumple la tolerancia así que se continúa el proceso
Para n=4 se tiene:
f ( x5 )( x5 − x4 )
x6 = x5 −
f ( x5 ) − f ( x4 )
Los valores ahora son: x4= 1.55 x5=1.396
f (1.396)(1.396 − 1.55) (−0.051)(−0.154.)
x6 = 1.396 − = 1.396 − = 1.413
f (1.396) − f (1.55) (−0.051) − (0.4025)
x6 − x5 = 1.413 − 1.396 = 0.017 = 0.1
aquí ya se cumplió el criterio de paro
por lo tanto, la raíz aproximada es:
x6 = 1.413
MÉTODO DE LA FALSA POSICIÓN.
Este método es una combinación del método de la secante con el método de la bisección.
Para utilizar este método sólo se pide que f(x) sea una función continua en a, b , ahora
usando la ecuación (2.5.3):
f ( x1 )( x1 − x0 )
c = x1 −
f ( x1 ) − f ( x0 )
Donde x0 =a x1 =b, sustituyendo estos valores en la ecuación anterior se tiene:
f (b)(b − a)
c =b− (2.6.1)
f (b) − f (a)
A partir de la ecuación (2.6.1) el intervalo [a,b], se divide en dos sub- intervalos de la forma:
[a,c] y [c,b]
En esta parte hay que elegir el sub- intervalo en la cual se encuentra la raíz, esto se realiza
usando la condición:
si f(a)*f(c)<0 entonces
c=b
de lo contrario
a=c;
el criterio para detener el procedimiento, uno de ellos puede ser f( c )
Si no se cumple, éste regresa a la ecuación (2.6.1). En caso contrario, la aproximación a la
raíz es c.
Ejemplo.
Calcular la raíz de la función f(x)=x2-2 en el intervalo [0,3] y la tolerancia =0.1.
Solución.
En este caso a=0 b=3 usando la ecuación (2.6.1).
f (b)(b − a) f (3)(3 − 0)
c =b− = 3−
f (b) − f (a) f (3) − f (0)
7 * (3)
c = 3− = 0.6666
7 − (−2)
Los sub-intervalos formados son:
[0,0.666] y [0.666,3]
Para saber en que intervalo se encuentra la raíz se utiliza el Teorema [Link].
f(a)*f(c)<0
f(0) * f(0.666) =(-2)*(-1.555)<0 la cual no se cumple, por lo tanto el sub- intervalo donde
se encuentra la raíz es:
[0.666,3]
Se probará si ya se alcanzó la tolerancia.
f(c) = f(0.666) = 1.556 = 0.1 Esta condición no se cumple.
Ahora a=0.666 b=3
f (b)(b − a) f (3)(3 − 0.666)
c =b− = 3−
f (b) − f (a) f (3) − f (0.666)
7 * (2.334)
c = 3− = 1.0909
7 − (−1.556)
Los sub-intervalos formados son:
[0.666,1.0909] y [1.0909,3]
Para saber en qué intervalo se encuentra la raíz se utiliza el Teorema [Link].
f(a)*f(c)<0
f(0.666) * f(1.0909) =(-1.555)*(-0.8099)<0 la cual no se cumple, por lo tanto el sub-
intervalo donde se encuentra la raíz es:
[1.0909,3]
Se probará si ya se alcanzó la tolerancia.
f(c) = f(1.0909) = 0.8099 = 0.1 Esta condición no se cumple.
Ahora a=1.0909 b=3
f (b)(b − a) f (3)(3 − 1.0909)
c =b− = 3−
f (b) − f (a) f (3) − f (1.0909)
7 * (1.9091)
c = 3− = 1.1127
7 − (−0.8099)
Los sub-intervalos formados son:
[0.666,1.1127] y [1.1127,3]
Para saber en qué intervalo se encuentra la raíz, se utiliza el Teorema [Link].
f(a)*f(c)<0
f(0.666) * f(1.1127) =(-1.555)*(-0.7618)<0 la cual no se cumple, por lo tanto el sub-
intervalo donde se encuentra la raíz es:
[1.1127,3]
Se probará si ya se alcanzó la tolerancia.
f(c) = f(1.1127) = 0.7618 = 0.1 Esta condición no se cumple.
Ahora a=1.1127 b=3
f (b)(b − a) f (3)(3 − 1.1127 )
c =b− = 3−
f (b) − f (a) f (3) − f (1.1127 )
7 * (1.8873)
c = 3− = 1.2979
7 − (−0.7618)
Los sub-intervalos formados son:
[1.1127,1.2979] y [1.2979,3]
Para saber en qué intervalo se encuentra la raíz se utiliza el Teorema [Link].
f(a)*f(c)<0
f(1.1127) * f(1.2979) =(-0.7618)*(-0.3154)<0 la cual no se cumple, por lo tanto el sub-
intervalo donde se encuentra la raíz es:
[1.2979,3]
Se probará si ya se alcanzó la tolerancia.
f(c) = f(1.2979) = 0.3154 = 0.1 Esta condición no se cumple.
Ahora a=1.2979 b=3
f (b)(b − a) f (3)(3 − 1.2979)
c =b− = 3−
f (b) − f (a) f (3) − f (1.2979)
7 * (1.7021)
c = 3− = 1.3712
7 − (−0.3154)
Los sub-intervalos formados son:
[1.2979, 1.3712] y [1.3712,3]
Para saber en qué intervalo se encuentra la raíz se utiliza el Teorema [Link].
f(a)*f(c)<0
f(1.2979) * f(1.3712) =(-0.3154) *(-0.1198)<0 la cual no se cumple, por lo tanto el sub-
intervalo donde se encuentra la raíz es:
[1.3712,3]
Se probará si ya se alcanzó la tolerancia.
f(c) = f(1.3712) = 0.1198 = 0.1 Esta condición no se cumple.
Ahora a=1.3712 b=3
f (b)(b − a) f (3)(3 − 1.3712)
c =b− = 3−
f (b) − f (a) f (3) − f (1.3712)
7 * (1.6288)
c = 3− = 1.3986
7 − (−0.1198)
Los sub-intervalos formados son:
[ 1.3712, 1.3986] y [1.3986,3]
Para saber en qué intervalo se encuentra la raíz se utiliza el Teorema [Link].
f(a)*f(c)<0
f(1.3712) * f(1.3986) =(-0.1198) *(-0.04391)<0 la cual no se cumple, por lo tanto el sub-
intervalo donde se encuentra la raíz es:
[1.3986,3]
Se probará si ya se alcanzó la tolerancia.
f(c) = f(1.3986) = 0.0439 = 0.1 Esta condición se cumple.
Así la aproximación a la solución es c=1.39861.4
Método del punto fijo (aproximaciones sucesivas).
Este método usa sólo un valor x0 arbitrario, que se da como primera aproximación, para
generar una sucesión de puntos x0, x1,...xn que cuando n → se obtiene un punto fijo (ver
teorema 2.04). El es sencillo, sólo consiste en dado una función f(x)=0 continua en [a,b], se
deja de un solo miembro una "x" de f(x) quedando de la siguiente manera:
x=G(x) (2.12.1)
Así que cualquier solución de x=G(x), es también solución de f(x), en general podemos
realizar diferentes despejes de G(x), pero no quiere decir que todos deben funcionar.
Para nuestro cálculo.
Ejemplo 1
Sea f(x) = x2-a =0 para a > 0
Se encontrará G(x).
1) x = x2 + x -a donde G(x)= x2 + x –a
2) x = x + c ( x2 - a ) para c 0 donde G(x)= x + c ( x2 - a )
a a
3) x = donde G(x)=
x x
1 a 1 a
4) x = x + donde G(x)= x +
2 x 2 x
Ejemplo 2
Sea f(x) = x2 - x -2 =0 encontrar G(x)
1) x = x2 – 2 donde G(x)= x2 – 2
2) x= 2 + x donde G(x)= 2 + x
2 2
3) x= 1 + donde G(x)= 1 +
x x
4)
x2 − x − 2 x2 − x − 2
5) x= x - con m 0 donde G(x)= x - con m 0
m m
Ahora la forma recurrencia de la ecuación (2.12.1).
xn+1=G(xn) (2.12.2)
Que es la que se utiliza para generar la sucesión, como criterio de paro se utiliza.
xn+1-xn
Como una observación, se graficará por separado la ecuación (2.12.2),
Donde.
1. Y1(x)=x
2. Y2(x)=g(x)
Aquí la función (1) es una línea recta a 450 que pasa por el origen, la función (2) es cualquiera
ver figura (2.12.1).
De la ecuación (a) sé esta afirmando que la gráfica (1) es igual que la (2), lo cual no es posible,
lo que sí es que las dos gráficas son iguales en ciertos puntos, ¿cuáles son estos?, Son aquellos
en los que se cortan. Las dos gráficas (ver figura). A estos puntos se les suele llamar puntos
fijos i , son a su vez solución de f(x) (f( i )=0).
Ejemplo 1.
Sea la función f(x)= x2 - 2=0 x0 =3, = 0 .01 calcular la raíz.
Se usará como despeje la siguiente función para G(x):
2
G(x)=
x
La función de recurrencia para este caso:
x n +1 = G (x n ) =
1
xn
Sí n=0
2 2
x1=G(x0)= = 0.66
x0 3
x1 − xo = 0.1
0.66 − 3 = − 2.33 = 2.33 = 0.1
no se cumple el criterio de paro, se continúa el proceso
Sí n=1
2 2
X 2=G(x1)= = 3
x1 0.66
el cual no se cumple
Sí n=2
x 3 = G (x 2 ) =
1 2
= 0.66
x2 3
x 3 - x 2 = 0.66 − 3 = − 2 − 33 = 2.33 = 0.01
El cual no se sigue cumpliendo.
El proceso continúa hasta que se cumpla la tolerancia
.
Nota. Observe que el proceso nunca converge, de las iteraciones se tiene que se entró en
un ciclo.
Ejemplo.
Se resolverá el mismo problema, pero en este caso se utilizará como G(x) a:
1 2
G(x) = x +
2 x
para n=0
1 2
x n +1 = G(x n ) = xn +
2 xn
1 2 1 2 1 9 + 2 1 11
x1 = G(x 0 ) = x0 + = 3 + = = 1.83
2 x0 2 3 2 3 2 3
x1 - x 0 = 1.81 − 3 = − 1.16 = 1.16 = 0.01
No se cumple el criterio de paro, se continúa el proceso.
Para n=1
1 2 1 2 1
= (1.83 + 1.09 ) = (2.92 ) 1.46
1
x 2 = G(x 1 ) = x1 + = 1.83 +
2 x1 2 1.83 2 2
x 2 − x1 = 1.46 − 1.81 = − 0.36 = 0.36 = 0.01
No se cumple el criterio de paro, se continúa con el proceso.
para n=2
1 2 1 2 1
= (1.46 + 1.37 ) = (2.92 ) 1.414
1
x 3 = G(x 2 ) = x2 + = 1.46 +
2 x2 2 1.46 2 2
x3 − x 2 = 1.414 − 1.46 = − 0.046 = 0.046 = 0.01
Aquí se cumple el criterio de paro:
La solución al problema es:
x 3 1.414
Una manera de asegurar que el método de las aproximaciones sucesivas va a converger es la
utilización del siguiente teorema.
Una condición necesaria para que converja es:
Dada G(x) una función continua con primera derivada continúa si G(x0 ) 1 el método del
punto fijo converge.
Método de Newton (o de la tangente).
Uno de los primeros métodos que trabajan a partir de un valor inicial arbitrario x0, el cual da
una buena aproximación, este se puede desarrollar de varias maneras una de ellas es por
medio de un desarrollo de Taylor alrededor de x0, una forma más sencilla de formularlo es
de la siguiente manera:
Dada una función continua, con un valor arbitrario x0, de la figura (2.8.1)
Se tienen las coordenadas de dos puntos que definen a la recta, estos son:
(x0 , f(x0)), (x1,0) recordando un poco de geometría analítica, se sabe que la pendiente de la
recta que pasa por estos dos puntos es:
y -y
m= 1 0 (2.8.1)
x1 - x 0
donde y0=f(x0), y1=0 sustituyendo en la ecuación (2.8.1)
0 - f(x0 )
m= (2.8.2)
x1 - x 0
Ahora, usando el cálculo diferencial se calcula la derivad, evaluada en el punto x0, ésta se le
denota por: f (x 0 ) que a su vez debe ser la misma que se calculó por medio de la pendiente
de una recta..
m = f (x 0 ) (2.8.3)
Igualando las ecuaciones (2.8.2) y (2.8.3).
= f (x 0 )
0 - f(x0 )
m= (2.8.4)
x1 - x 0
De la ecuación (2.8.4) el único valor desconocido es x1, realizando un despeje:
- f(x0 ) = f (x 0 )(x1 - x 0 ) (2.8.5)
- f(x0 )
= x1 - x 0 (2.8.6)
f (x 0 )
f(x )
x1 = x 0 − 0 (2.8.7)
f (x 0 )
Generalizando la ecuación (2.8.7):
f(x )
x n +1 = x n − n (2.8.8)
f (x n )
A esta ecuación se le llama forma de recurrencia del método de Newton.
El criterio de paro para este es:
x n +1 - x n = Tolerancia (2.8.9)
Ejemplo:
Calcular la raíz de la siguiente función:
f(x) = x 2 − 2 , f (x ) = 2x para x0 =3 y =0.1
Sustituyendo estos valores en la ecuación (2.8.8):
x2 − 2
x n +1 = x n − n
2 * xn
para n=0
x2 − 2 32 − 2 11
x1 = x 0 − 0 = 3 − = = 1.83
2 * x0 2*3 6
x1 - x 0 = 1.83 − 3 = 1.17 = 0.1, No se cumple la tolerancia así que se continúa el
proceso.
Para n=1
x12 − 2 1.832 − 2
x 2 = x1 − = 1.83 − = 1.83 − 0.36 = 1.46
2 * x1 2 * 1.83
x 2 - x1 = 1.46 − 1.83 = 0.37 = 0.1 , No se cumple la tolerancia así que se continúa el
proceso.
x 22 − 2 1.46 2 − 2
x3 = x 2 − = 1.46 − = 1.46 − 0.45 = 1.415
2 * x2 2 * 1.46
x 3 - x 2 = 1.415 − 1.46 = 0.045 = 0.1 , aquí se cumple la tolerancia, el resultado a este
problema es:
x3=1.415