Problemas Método de bisección. Método de Newton. - CAA EII

Problemas
Método de bisección.
1. La ecuación e x − 3x = 0 tiene por raíz a r = 0,61906129. Comenzando con el intervalo
[0, 1], realizar seis iteraciones por el Método de bisección para encontrar la raíz
aproximada. ¿Cuántos decimales significativos tiene dicha aproximación?. ¿Cuántas iteraciones
son necesarias para que la raíz obtenida tenga un error menor que 10 −4 ?
2. Utilizar el Método de bisección para encontrar una solución aproximada con un error
menor que 10 −2 en el intervalo [4, 4,5] para la ecuación x = tg(x).
3. Sabiendo que existe una raíz de la ecuación x 3 + x = 6 entre 1.55 y 1.75, ¿cuántas
iteraciones son necesarias hasta obtener mediante el método de bisección, un intervalo de
amplitud menor o igual que 10 −3 que contenga a la raíz?. Calcular todas las iteraciones
necesarias.
4. Aplicar el Método de bisección a F(x) = x 3 −17 = 0, a fin de determinar la raíz cúbica
de 17 con un error menor que 0.125.
Método de Newton.
5. Aplicando el Método de Newton, encontrar una raíz próxima a x0 = 0 para la ecuación
f(x) = 3x + senx − e x = 0.
Redondear los cálculos a cinco cifras significativas e iterar hasta que se cumpla | xi −
xi−1 |≤ 0,001.
6. La función f(x) = 4x−7
x−2 tiene una raíz en x=1.75. Utilizar el método de Newton con las
siguientes aproximaciones iniciales, estudiando en cada caso, previamente, si se produce
un proceso convergente o no a la raíz.
a)x0 = 1,6 , b)x0 = 1,5 , c)x0 = 3
7. Mediante el Método de Newton modificado, encontrar una raíz próxima a x0 = 0 de la
ecuación x − 2 −x = 0.
Utilizar tres decimales redondeados en cada iteración hasta que se cumpla | xi −xi−1 |≤
10 −3 .
8. La concentración c de una bacteria contaminante en un lago decrece según la expresión:
c(t) = 80e −2t + 20e −0,5t
siendo t el tiempo en horas. Determinar el tiempo que se necesita para que el número
de bacterias se reduzca a 7. (Utilizar el Método de Newton).
9. Una determinada sustancia se desintegra según la ecuación A = P · e −0,0248t , donde P
es la cantidad inicial en el tiempo t = 0 y A la cantidad resultante después de t anõs. Si
inicialmente se depositan 500 miligramos de dicha sustancia, ¿cuánto tiempo habrá de
transcurrir para que quede el 1 por ciento de ésta? Utilizar el Método de Newton.
Ingeniería Técnica
Forestal
1 Fundamentos Matemáticos
Curso 2004/05

10. Demostrar que para encontrar la raíz r-ésima de un número a, la fórmula iterativa de
Newton se puede expresar como
xn+1 = 1
r [(r − 1) · xn + a
xn r−1]
11. Hallar la raíz cuadrada de 10 usando tres iteraciones mediante el método de Newton
y comenzando con el valor inicial x0 = 3. Utilizar dos decimales redondeados en los
cálculos.
12. Se considera la función F(x) = x 5 +2x. Mediante el Método de Newton, hallar el menor
número positivo x (con tres decimales) para el cual F(x) = 4.
13. En los casos siguientes, aplicar el método de Newton con la estimación inicial propuesta,
y explicar por qué falla el método.
a) y = 2x 3 − 6x 2 + 6x − 1, x1 = 1.
b) y = 4x 3 − 12x 2 + 12x − 3, x1 = 3
2 .
c) y = −x 3 + 3x 2 − x + 1, x1 = 1.
d) y = 3√ x − 1, x1 = 2.
14. Probar, mediante el método de Newton, que la ecuación
xn+1 = xn(2 − axn)
se puede utilizar para aproximar 1
a si x1 es una estimación inicial del recíproco de a.
Nótese que este método de aproximar recíprocos utiliza sólo operaciones de suma y
− a.]
multiplicación. [Ayuda: Considerar f(x) = 1
x
15. Aproximar, con ayuda del resultado del ejercicio anterior, con tres cifras decimales, los
siguientes recíprocos:
a)
b)
1
3 .
1
11 .
16. Una medicina administrada a un paciente produce una concentración en la sangre dada
por c(t) = Ate −t/3 mg/ml, t horas después de que se hayan administrado A unidades.
La máxima concentración sin peligro es de 1 mg/ml, y a esta cantidad se le denomina
concentración de seguridad.
a) ¿Qué cantidad debe ser inyectada para alcanzar como máximo esta concentración
de seguridad?. ¿Cuándo se alcanza este máximo?.
b) Una cantidad adicional se debe administrar al paciente cuando la concentración
baja a 0 ′ 25 mg/ml. Determínese con un error menor de 1 minuto cuándo debe
ponerse esta segunda inyección.
Ingeniería Técnica
Forestal
2 Fundamentos Matemáticos
Curso 2004/05

17. El crecimiento de poblaciones grandes puede modelarse en períodos cortos suponiendo
que el crecimiento de la población es una función continua en t mediante una ecuación
diferencial cuya solución es
N(t) = N0e λt + v
e
λ
λt
− 1 ,
donde N(t) es el número de individuos en el tiempo t (medido en años), λ es la razón
de natalidad, N0 es la población inicial y v es un razón constante de inmigración, que
se mide en número de inmigrantes al año.
Supóngase que una población dada tiene un millón de individuos inicialmente y una
inmigración de 400,000 individuos al año. Se observa que al final del primer año la
población es de 1,506,000 individuos. Se pide:
a) Determinar la tasa de natalidad.
b) Hacer una previsión de la población al cabo de tres años.
Ingeniería Técnica
Forestal
3 Fundamentos Matemáticos
Curso 2004/05

Soluciones a algunos ejercicios
1. En efecto, si llamamos f(x) = e x −3x, por el Teorema de Bolzano, al ser f(1) = e−3 < 0,
y f(0) = 1 − 0 = 1 > 0, sabemos que hay al menos una 1 raíz en el intervalo [0, 1].
1.4
1.2
1
0.8
0.6
0.4
0.2
0
−0.2
f(x)=exp(x)−3x
−0.4
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
f(x) = e x − 3x.
Nos piden que hagamos seis iteraciones por el Método de Bisección, esto es, vamos calculando
puntos medios de los intervalos, y el valor de la función en dichos puntos, quedándonos
con aquél donde haya cambio de signo (dicho de otro modo, ci+1 sustituye a extremo cuya
imagen tenga el mismo signo que f(ci+1)), e iteramos de nuevo:
c0 = 1/2 f(c0) = 0 ′ 148 > 0 ⇒ [c0, 1]
c1 = 3/4 f(c1) = −0 ′ 133 < 0 ⇒ [c0, c1]
c2 = 5/8 f(c2) = −0 ′ 006 < 0 ⇒ [c0, c2]
c3 = 9/16 f(c3) = 0 ′ 067 > 0 ⇒ [c3, c2]
c4 = 19/32 f(c4) = 0 ′ 029 > 0 ⇒ [c4, c2]
c5 = 39/64 f(c5) = 0 ′ 011 > 0 ⇒ [c5, c2].
De modo que c6 = 79/128 ≡ 0 ′ 6171875, y la función ahí vale f(c6) ≡ 0 ′ 002. En efecto, el
valor es relativamente pequeño, aunque no una aproximación excelente, y es que c6 dista de
la solución exacta que nos da el enunciado. En realidad, lo único que sabíamos por el método
a priori es que la cota de error entre la raíz exacta y la aproximada era |x ∗ − cn| ≤ b−a
2 n+1 =
1
ln(b−a)−ln(δ)
ln 2
27 ≡ 0 ′ 0078125, por lo que efectivamente no cabía esperar más de dos decimales exactos.
Las iteraciones necesarias para obtener orden δ = 10−4
son n = E
= 13.
Nota sobre su resolución en el ordenador:
Como el orden del método es uno, la aproximación es lenta. Sin embargo, se comprueba que
es fácilmente automatizable (i.e. podemos llevar el método al ordenador, evaluando la función
en el punto medio, y pidiéndole a la computadora que compare el signo y en función de si sale
igual o distinto tomar un intervalo u otro para la nueva iteración). Con ayuda de Microsoft
Excel c○ podemos implementar fácilmente el esquema (buscar la función SI, véase material
1 De hecho, f ′ (x) = e x − 3, con lo que deducimos que en x = ln 3 la derivada se anula y la función tiene un
mínimo; como después crecerá mucho, pasará el cero de nuevo una única vez.
Ingeniería Técnica
Forestal
4 Fundamentos Matemáticos
Curso 2004/05

complementario de este tema en la web www.uhu.es/pedro.marin/docencia/extraT2C.zip),
aunque por contra la precisión de cálculo deja de ser buena muy pronto. MATLAB, sin embargo,
ofrece mayor exactitud, aunque la adaptación a este lenguaje puede ser algo más tediosa
(se puede hacer con un proceso por lotes, o directamente con una función .m, en todo caso
usando el condicional, que en este programa es IF).
4. Para hallar la raíz cúbica de 17, es una elección evidente tomar la función f(x) = x3 −17.
Tanteamos mentalmente para asegurarnos un intervalo inicial en el que ejecutar el Algoritmo
de Bisección (p.ej. [2, 3]). Con error menor que 0 ′
ln 1−ln 0 ′ 125
125 debemos efectuar n = E ln 2 =
0−ln(1/8)
E ln2 = E
3ln 2
ln 2 = 3 iteraciones2 , y la solución aproximada resultante es 2 ′ 5625,
mientras que una aproximación mejor es 3√ 17 ∼ 2,5712815906154.
5. Usando el Teorema de Bolzano sabemos que f(x) = 3x − sin x − e x tiene una raíz en el
intervalo [0, 1].
1.5
1
0.5
0
−0.5
−1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
f(x)=3x−sen(x)−exp(x) .
Como f ′ (x) = 3 + cos(x) − e x , el Método de Newton-Raphson resulta
xn+1 = xn − f(xn)
f ′ (xn) = xn − 3xn + sin(xn) − exn . xn 3 + cos(xn) − e
El primer dato lo podemos elegir a nuestro antojo, tomamos por ejemplo x0 = 0, con lo
que x1 = 1/3, x2 = 0 ′ 360170714 y x3 = 0 ′ 36042168. Ya hemos cumplido la condición |xi −
xi−1| ≤ 0 ′ 001 pues |x3 − x2| = 0 ′ 00025. Además, comprobamos efectivamente que f(x3) =
−5,744246611705250 10 −8 . Vemos gráficamente la aproximación realizada en dos etapas (la
derecha es una ampliación)
2 Observa que la calculadora no es indispensable siempre.
Ingeniería Técnica
Forestal
5 Fundamentos Matemáticos
Curso 2004/05

5
4
3
2
1
0
−1
−2
−3
−4
−5
−1 −0.5 0 0.5 1 1.5 2
1
0.5
0
−0.5
−1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
usando las rectas tangentes: y = f ′ (x0)(x − x0) + f(x0), e y = f ′ (x1)(x − x1) + f(x1).
6. Claramente la única raíz de la función f(x) = 4x−7
x−2 es x∗ = 1 ′ 75. Si implementamos el
método de Newton-Raphson para aproximar la solución
f ′ (x) = −1
(x − 2) 2,
⇒ xn+1 = xn − f(xn)
f ′ (xn) = xn + (4xn − 7)(x − 2),
comprobamos que con los tres valores que nos dan se obtienen las siguientes secuencias:
1 ′ 6, 1 ′ 84, 1 ′ 7824, 1 ′ 75419904, 1 ′ 750070528, 1 ′ 75000002, 1 ′ 75, 1 ′ 75, 1 ′ 75 . . .
3, 8, 158, 97658, 38146972658, 5 ′ 82077E + 21, 1 ′ 35525E + 44, 7 ′ 34684E + 88, . . .
1 ′ 5, 2, 2, 2, . . .
es decir que la primera converge a la verdadera solución, la segunda diverge y la tercera
converge a una falsa solución. Ello es debido a que el método es local, y hay que empezar
suficientemente cerca de la solución para tener garantía de convergencia, lo que nos impele a
hacer siempre la comprobación (cuando obtengamos una convergencia hacia cierto valor) de
que la imagen del valor aproximado obtenido está cerca de cero, para evitar falsas soluciones.
[En general los códigos suelen tener dos condiciones para parar: que dos iteraciones consecutivas
estén cerca entre sí y que la imagen de una de ellas esté cerca de cero.]
8. Debemos hallar la raíz de c(t) = 7, o lo que es lo mismo, definiendo f(t) = c(t) − 7,
tenemos que hallar un cero para f. Ambas funciones, c y f tienen por derivada a la función
f ′ (t) = −160e −2t −10e −t/2 . Al ser negativa, sabemos que las funciones c y f son decrecientes
estrictamente, luego si existe solución, es única.
Como c(0) = 100 y límt→∞ c(t) = 0, en tiempo positivo la función c (que es continua)
tomará todos los valores del intervalo (0, 100].
Una vez que hemos concluido que existe una única raíz de c(t) = 7, y por tanto un único
cero de f(t), aplicamos el Método de Newton-Raphson. Como c(0) = 100 dista bastante
del objetivo, conviene, para ahorrar cálculos empezar con un dato inicial del tiempo algo
mayor (no mucho, porque la exponencial decae rápidamente), por tanteo parece conveniente
comenzar con x1 = 2. En cuatro iteraciones conseguimos:
x2 = 2 ′ 27579938312569
x3 = 2 ′ 32768095901589
Ingeniería Técnica
Forestal
6 Fundamentos Matemáticos
Curso 2004/05
.

x4 = 2 ′ 32908663665517
x5 = 2 ′ 32908761684562.
9. Es propio del material radioactivo desintegrarse con respecto a esa ley (solución de cierta
ecuación diferencial que veremos más adelante en el Tema 4). Obsérvese que resolver el
problema equivale a hallar el cero de f(t) = 500e −0′ 00248t −5, o lo que es lo mismo (pero mejor
para hacer los cálculos), de la función g(t) = 100e −0′ 00248t − 1. Antes de aplicar el Método
de Newton-Raphson es adecuado pensar en qué dato inicial tomar (de no ser adecuado,
tardaremos mucho o puede, como se ve en algunos ejercicios del tema, que no lleguemos a la
solución).
Es claro geométricamente que, por la forma de la exponencial, el método va a ser efectivo
con cualquier dato inicial.
100
90
80
70
60
50
40
30
20
10
0
10 20 30 40 50 60 70 80 90 100
g(t) = 100e −0′ 00248t − 1.
Pero por el bajo coeficiente que afecta al exponente, ha de pasar mucho tiempo para que la
exponencial comience a decrecer de forma notable. Por ello, comenzamos con t1 = 100. En 6
iteraciones llegamos a la solución:
t2 = 135 ′ 507554670222
t3 = 164 ′ 214770256907
t4 = 180 ′ 865946787540
t5 = 185 ′ 414686713111
t6 = 185 ′ 691392424226
t7 = 185 ′ 692346197917 años han de pasar para quedar 3 el 1 %.
10. Se obtiene tras transformaciones aritméticas inmediatas a
siendo f(x) = x r − a.
xn+1 = xx − f(xn)
f ′ (xn)
11. El ejercicio anterior se aplica con r = 2 y a = 10 generando la sucesión recurrente
Ingeniería Técnica
Forestal
xn+1 = 1
xn +
2
10
.
xn
7 Fundamentos Matemáticos
Curso 2004/05

Elegir x0 = 3 es una buena elección: simple y aproximada, esto último es importante ya
que el Método de Newton-Raphson es local, por lo que cuánto más cerca se comience de la
solución, más garantías de éxito se tendrá al aplicarlo. x0 = 3, x1 = 3 ′ 16, x2 = 3 ′ 1622777,
x3 = 3 ′ 1622777 (que ya coincide con el resultado que podemos obtener usando la calculadora).
13. a) En este caso el Método de Newton-Raphson no es aplicable porque el punto inicial
anula el denominador f ′ (x) que aparece en la expresión recurrente. En cambio, si tomamos
un valor un poco más alejado, por ejemplo, 1 ′ 2, conseguimos aproximar en 9 iteraciones la
solución (con un error/tolerancia de 10 −8 ):
x1 = −3 ′ 03333333333329
x2 = −1 ′ 69913409071934
x3 = −0 ′ 82229976699472
x4 = −0 ′ 26505557694095
x5 = 0 ′ 05248674618149
x6 = 0 ′ 18268167510526
x7 = 0 ′ 20562357143852
x8 = 0 ′ 20629889908101
x9 = 0 ′ 20629947401548.
Empezando en x0 = 0 conseguimos la convergencia (con igual tolerancia) en 4 iteraciones:
x1 = 0 ′ 16666666666667
x2 = 0 ′ 20444444444444
x3 = 0 ′ 20629515192918
x4 = 0 ′ 20629947399236.
b) La derivada de f se anula en x = 1, con lo que si los cálculos pasan por algún xn = 1,
el método no será válido. Y justamente ése es el caso si x1 = 3/2, ya que x2 = 1. Sin embargo,
eso no significa que el método en si sea malo para hallar el cero de esta función, sólo que el
dato inicial no es el adecuado. En efecto, si empezamos por otro valor, x1 = 0, el método
converge en 5 iteraciones:
x2 = 0 ′ 25
x3 = 0 ′ 35185185185185
x4 = 0 ′ 36953388762913
x5 = 0 ′ 37003906971704
x6 = 0 ′ 37003947505230.
c) Con el dato inicial la sucesión que se obtiene es de 0 y 1 alternante, con lo que nunca
convergerá a ningún valor. Aunque no hayamos visto en el desarrollo de teoría las condiciones
de convergencia de los métodos de punto fijo (el MNR lo es para la función g(x) = x −
f(x)
f ′ (x) ), intuimos con este ejemplo que no deben darse condiciones oscilantes en el entorno que
tomemos para comenzar la construcción de la sucesión recurrente, sino que debe ser monótona
en cierto sentido hacia el cero de la función f. Estos ejemplos muestran simplemente malas
elecciones del dato inicial que hacen que en el camino se tope uno con dificultades, que se
evitarían (una vez más) empezando por otro dato, más próximo y adecuado: con dato inicial
3 se consigue la aproximación en 4 iteraciones
x1 = 2 ′ 8
Ingeniería Técnica
Forestal
8 Fundamentos Matemáticos
Curso 2004/05

x2 = 2 ′ 76994818652850
x3 = 2 ′ 76929266290594
x4 = 2 ′ 76929235423870.
d) En este caso la explicación es clara si tenemos en cuenta la gráfica de la función:
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
El pico que, en torno a x = 1, forma la gráfica de la función f(x) = 3√ x − 1, unido a la obvia
raíz que tiene, es una mala condición para poder apoyarnos en el uso de rectas tangentes.
Toda la regularidad que necesitamos, aquí falta, tras acercarse inicialmente, cambia de signo,
lo alterna, y se va alejando (diverge).
14. En este ejercicio se pide probar que una forma de aproximar el valor 1/a es a través de
la sucesión xn+1 = xn(2 −axn), esto es trivial con la indicación que dan: la raíz de la función
f(x) = 1
x − a es justamente x∗ = 1
a , y el Método de Newton-Raphson genera justamente la
anterior relación de recurrencia:
f ′ (x) = −1
x 2
⇒ xn+1 = xn − f(xn)
f ′ (xn) = xn −
1 − a xn
−1
x2 = xn +
n
1 − a xn
1
x2 = xn + xn − ax
n
2 n.
Lo reseñable del ejercicio es que para obtener el cociente 1/a no es necesario dividir, sino
aproximar con sumas y productos.
16. Vemos como la realidad no se reduce siempre al uso de una función de una única variable.
a) Fijada la dosis inicial A, el problema, con la función dada c(t), consiste en hallar el
máximo de concentración, y el momento. Podemos contestar ya a lo segundo: buscamos los
extremos relativos a través de la derivada de la función:
c ′ (t) = Ae −t/3 − 1
3 Ate−t/3 = Ae −t/3 (1 − t/3) .
Por tanto, los candidatos a extremos de la función c (recuérdese que el dominio en que tiene
sentido el problema es [0, ∞)) son {0, 3, ∞}. Sin embargo, c(0) = 0, y límt→∞ c(t) = 0, de
modo que la cantidad c(3) = 3e −1 A > 0 es la concentración máxima.
Ingeniería Técnica
Forestal
9 Fundamentos Matemáticos
Curso 2004/05

De hecho, cualquiera que sea la cantidad inicial A, el máximo de la concentración se
alcanza en t = 3. Para que la cantidad máxima no supere 1 mg/ml, simplemente hay que
poner un A menor del que resuelva 3e −1 A = 1, esto es, A ≤ e/3.
b) Supuesta suministrada esta cantidad inicial, A = e/3, nos preguntan cuándo ocurrirá
que
c(t) = e
3 te−t/3 = 0 ′ 25.
Tenemos que resolver pues la ecuación
f(t) = e
3 te−t/3 − 0 ′ 25 = 0.
Pero del análisis anterior sacamos que al máximo llega desde cero y que después tiende a
cero, es decir, cualquier cantidad en (0, A), en particular 0 ′ 25, es alcanzada dos veces y nos
interesa cuando lo alcanza por segunda vez.
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 5 10 15
c(t) = e
3 te−t/3 .
Usamos el método de Newton-Raphson con algún valor numérico que resulte adecuado
(por lo dicho antes, más bien largo que corto, para que la convergencia sea hacia la solución
por encima de las tres horas). Por la forma de campana, basta tomar cualquier valor inicial
superior a tres horas, p. ej. x0 = 4 genera una sucesión que converge en 4 iteraciones (con
una tolerancia de 10 −8 ; nos pedían precisición de un minuto, que pasado a horas es 0 ′ 01666,
por tanto bastaba con 10 −3 ):
x1 = 10 ′ 30879920381349
x2 = 11 ′ 02140006874606
x3 = 11 ′ 07757133060878
x4 = 11 ′ 07790357510362 horas= 11 horas, 4 minutos.
17. En este problema nos dan el modelo
N(t) = N0e λt + v
e
λ
λt
− 1 ,
los datos N0 = 10 6 , v = 4 10 5 , N(1) = 1506 10 3 , y nos piden simplemente que despejemos el
valor de λ :
Ingeniería Técnica
Forestal
10 Fundamentos Matemáticos
Curso 2004/05

1506 10 3 = 10 6 e λ 4
105
+ e
λ
λ
− 1 .
Simplificamos antes de resolver el cero de cierta función asociada:
f(λ) = 1000e λ + 400
λ
e λ
− 1 − 1506 = 0.
Con un dato inicial bajo x1 = 1 (lo esperable para la natalidad de una población que en
un año sólo ha pasado de 1400000 a algo más de millón y medio de habitantes) la respuesta
se obtiene en 4 iteraciones (con precisión 10 −8 ):
x2 = 0 ′ 390820116865
x3 = 0 ′ 12534036450639
x4 = 0 ′ 08560613200322
x5 = 0 ′ 08483943705533
x6 = 0 ′ 08483915873218 ∼ λ.
b) Consiste en sustituir t = 3 en la expresión (ahora totalmente conocida) de N(t) :
N(3) = 2656373 ′ 589004676, que significa una población en torno a los 2.656373 habitantes.
Ingeniería Técnica
Forestal
11 Fundamentos Matemáticos
Curso 2004/05