116 Integraci´on Num

´erica

Soluci´on. (i) Notemos que la integral del enunciado

∫ 1
sin x
I = 0 f (x) dx , f (x) = √ ,
1 − x2
es impropia de segunda especie puesto que la funci´on integrando f (x) no est´a
aco- tada en el intervalo de integraci´on [0, 1]. En particular, la gra´fica de f (x)
tiene una as´ıntota vertical en x = 1. Realizando el cambio de variable x = sin
t en la integral I se obtiene
∫ π/2 ∫ π/2
sin(sin t) = sin(sin )
=
√
I cos t dt t dt .
0 1 − sin2 0
t
Notemos que, tras el cambio de varible, la integral I deja de ser impropia puesto
que la funci´on integrando g(t) = sin(sin t) est´a acotada en el intervalo de
inte- graci´on [0, π/2].
A continuacio´n aproximamos el valor de I mediante la regla de los
trapecios compuesta con n = 4 subintervalos, es decir, con una longitud de
paso h = (π/2 − 0)/4 = π/8.

I ≈ g0 + + g4
(g) = h +g + 2
Th g2 g3
1
2
π g(0) g(π/2)
= 8 2 + g(π/8) + g(π/4) + g(3π/8) + 2

π8 02 0,841471
2
= + 0,373411 + 0,649637 + 0,797946 + = 0,880325 .

∫b
por(ii) El error absoluto Δ en la aproximaci´ona g(t) dt ≈ h (g) viene dado
T
h2 ′′
b− (ξ) , ξ (a, b) .
Δ= (
12 a)g ∈
En nuestro caso se tiene

(π/8)2 3
Δ= π/2 − 0)|g ′′(ξ)| = π | ′′(ξ)| ,
( 1536
12 g

siendo ξ (0, π/2). Derivando 2 veces la funci´on g(t) = sin(sin t) se tiene g′′
∈
(t) = cos(sin x) sin x cos2 x sin(sin x). Para hallar el valor ma´ximo (es decir,
− −
la menor de las cotas superiores) de Δ se ha de obtener el valor
 M = m´a x ′′
0≤t≤π/2{|g (t)|} .

Sin embargo, el problema de extremos absolutos que se plantea para hallar

M , pasa por obtener los extremos relativos de f ′′(t), es decir, por la ecuaci
´on f ′′′(t) = 0 que resulta ser una ecuaci´on trascendente. Entonces, optamos
por hallar una cota superior de |g′′(t)| en el intervalo [0, π/2] por simple
inspecci´on.
Cap´ıtulo 8

Solucio´n de Ecuaciones no
Lineales

8.1. Resumen de teor´ıa

Son muchos los problemas en ciencia e ingenier´ıa que se pueden modelar
matem´aticamente como una ecuaci´on f (x) = 0, siendo f una funci´on de
la variable x. Los valores de x soluciones de dicha ecuaci´on son llamados ceros
de la funci´on f o raices de la ecuaci´on. Existen muchas de estas ecuaciones que
no admiten que su soluci´on pueda ser expresada a trav´es de funciones
elementales. El objetivo del presente cap´ıtulo es el estudio de algunas de
las t´ecnicas num´ericas que nos permitan abordar este tipo de problemas.
Es importante destacar el hecho de que las t´ecnicas que estudiaremos en
este cap´ıtulo son siempre iterativas, es decir, partiremos de una aproximaci´on
inicial x0 del cero exacto x∗ de f y posteriormente construiremos una sucesi´on
de nu´meros reales
{x n } n∞= 0 que converja hacia x∗ cuando n → ∞.
Teorema 8.1 (Bolzano) Sea f una funci´on continua en [a, b] tal que cambia
de signo en los extremos, es decir, f (a)f (b) < 0. Entonces, existe al menos un
ξ ∈ (a, b) tal que f (ξ) = 0.
El algoritmo de bisecci´on se basa en el teorema de Bolzano y consiste en
ir estrechando de manera sistem´atica el intervalo en el cual una funci´on
con- tinua cambia de signo. De esta forma conseguiremos obtener un
intervalo ar- bitrariamente pequen¯o que contenga el cero de la funci´on. El
m´etodo de bisec- ci´on para una funci´on f ∈ C[a0, b0] procede de la forma
siguiente: (a0, b0) ⊃ (a1, b1) ⊃ · · · ⊃ (an, bn) ⊃ · · ·. La sucesi´on de puntos
medios {c n } n∞= 0 siendo cn = (an + bn)/2, verifica l´ımn→∞ cn = x∗ con |x∗
− cn| ≤ (b − a)/2n+1.
M´etodo de Newton-Raphson:
f(xn)
xn+1 = xn − .
f ′(x n)
153
154 Soluci´on de Ecuaciones no
Lineales
Teorema 8.2 (Convergencia de la sucesi´on de Newton-Raphson)
Supongamos que f ∈ C2[a, b] con f (a)f (b) < 0 y f ′(x) y f ′′(x) son
signo constante ∀x ∈ [a, b]. Entonces, tomando x0 ∈ {a, b} que
de
signf (x0) = signf ′′(x0), la sucesi´on de Newton-Raphson {xn}n∈N
verifique
hacia la u´nica ra´ız de f en [a,
converge
b].
M´etodo de la Secante:

xn+1 xn − xn−1
= xn — f (xn )
f (xn) − f (xn−1) .

Consideremos un sistema de n ecuaciones con n inc´ognitas x1, x 2 ,..., xn.

Dicho sistema puede ser escrito de la forma f (x) = 0, siendo

x = (x1, x 2 ,..., xn)T .

La aplicaci´on f : Rn → Rn tiene las funciones componentes fi : Rn → R para

i = 1, 2 , . . . , n, es decir

f (x) = (f1(x), f 2 (x),..., fn(x))T .

Sea x∗ un cero de la funci´on f , es decir f (x∗) = 0. El m´etodo de Newton-Raphson

para aproximar las soluciones del sistema f (x) = 0 es

xi+1 = xi − (Df (xi))−1f (xi) con i = 0, 1,... , (8.1)

siendo Df (xi) la matriz Jacobiana de la funci´on f en el punto xi, es decir

Df (x ) = ⎞
∂f1(xi) ∂f1(xi)
⎛ 1
···
∂x ∂x n
⎟
⎠
⎜ i ⎝
∂fn
..(x ) ..
∂fn(xi) .
i

∂x1
··· ∂xn

Si la matriz Jacobiana Df (xi) es no singular, es decir si det Df (xi) = 0, entonces

se puede realizar la iteraci´on. /

8.2. Problemas resueltos

Problema 8.1 Consideremos el intervalo I = [0, 1/2] y sea g : I → I una
funci´on derivable en I cuya derivada est´a acotada de la forma |g′(x)| ≤
para todo x ∈ I. Partiendo de un punto arbitrario x0 ∈ I, se realiza la iteraci
1/2
xn+1 = g(xn).
´on
(i) Explicar si dicha iteracio´n ser ´a convergente o no. Si lo fuera, es decir
si
l´ımn→∞ xn = x∗, averiguar el valor g(x∗).
(ii) Calcular una cota superior de |x8 − x∗| en funci´on de la distancia entre los
dos primeros t´erminos de la sucesio´n |x1 − x0|.
8.2 Problemas resueltos 155

Soluci´on. La soluci´on del problema viene dada por el siguiente teorema:

g C[a, b] y g([a, b])
Si ∈ ⊆ [a, b], entonces existe almenos un punto fijo x [a,
∗

b] de la funci´on g. Si∈adem´as
∈ g D(a, b) y es contractiva, es decir
| g|(x)
′
≤ L<
1 para todo ∈x (a, b) entonces g tiene un u´ nico punto fijo x∗ [a, b]. Adem´as
la sucesi´on generada mediante la iteracio´n xn+1 = g(xn) converge hacia x∗
para
cualquier vlor inicial x ∈ [a, b]. Adem´as |x∗ − | ≤ Ln | — x |.
x x
0 n 1−L 1 0

(i) Como la funci´on del enunciado es g : I → I, entonces si se toma I = [a, b]

es evidente que g([a, b])⊆[a, b]. Adem´as, como g es una funci´on derivable en I
entonces se tiene que tambi´en es continua en I, es decir g∈ C[a, b]. Finalmente,
del enunciado se desprende que la constante L = 1/2 < 1. Se concluye que
la sucesi´on generada con la recurrencia xn+1 = g(xn) converge hacia el u´ nico
punto fijo x∗ de g y por lo tanto g(x∗) = x∗.

(ii) La cota superior es

∗ (1/2)8 1
|x − x8| ≤ |x1 − x0| = |x1 − x0| = 0,0078125|x1 − x0|.
1 − 1/2 128

Problema 8.2 Considerar el polinomio P (x) = x5 + 5x + 8.

(i) Encontrar dos nu´meros enteros negativos y consecutivos a y a +1 de manera
que P (x) tenga una u´ nica ra´ız x∗ ∈ I = [a, a + 1].
(ii) Averiguar cual de los dos extremos del intervalo I es el adecuado para ini-
ciar la iteracio´n de Newton-Raphson si se pretende asegurar la
convergen- cia hacia x∗ de la sucesi´on obtenida. Calcular la
aproximacio´n x∗ ≈ x2.
Soluci´on. (i) Como P (x) es un polinomio entonces P es continuo en todo R.
Entonces, si encontramos dos nu´meros enteros consecutivos en los cuales
cambie el signo de P se podr´a aplicar el Teorema de Bolzano y concluir que
entre esos dos puntos existe al menos una ra´ız de P . En concreto se tiene
que P ( 1) = 2 > 0 y P ( 2) = 34 < 0 de manera que existe al menos — un x∗
[ 2, −1]. − ∈ − −
Para demostrar que x∗ es u´ nico vamos a demostrar que P (x) tiene derivada
de signo constante en [− 2,− 1] de manera que s´olo puede cortar una vez al eje
de abscisas. En concreto se tiene que P ′(x) = 5x4 +5 = 5(x4 + 1) > 0 para todo
x ∈ R.

(ii) Para resolver la segunda cuesti´on utilizaremos el siguiente teorema: Sea

∈2[a, b] con f (a)f (b) < 0 tal que f ′(x) y f ′′(x) tienen signo constante
f C
para todo∈ x [a, b]. Sea
∈ {c }a, b tal que signf (c) = signf ′′(c). Entonces el m
´etodo de Newton-Raphson para f con valor inicial x0 = c converge hacia la u
´nica ra´ız de f en [a, b].
Ya se ha visto que P ′(x) tiene signo constante. Como P ′′ (x) = 20x3 < 0
para todo x < 0 se tiene que P ′′ es negativa en [− 2,− 1]. Finalmente, como
signP (− 2) = signP ′′( −2), el anterior teorema asegura la convergencia de la
sucesi´on generada con el m´etodo de Newton-Raphson para P tomando
156 Soluci´on de Ecuaciones no
como Lineales
8.2 Problemas resueltos 157

valor inicial x0 = −2.

Realizando la iteraci´on de Newton-Raphson

P ( xn )
x5 + 5xn + 8
n
xn+1 = xn − = xn −
P ′(xn) 5(xn4 + 1)

tomando x0 = −2 se obtiene

n xn
0 -2
1 -1.58025
2 -1.27709

Se concluye que x∗ ≈ x2 = −1,27709.

Problema 8.3 Considerar la funcio´n g(x) = sinh x + sin x − 3x.

(i) Demostrar que g(x) tiene un extremo relativo en el intervalo [1, 2].
Indicacio´ n: Las curvas y = sinh x, y = sin x no tienen ningu´n punto
de corte con abscisa en el intervalo [1, 2].
(ii) Aproximar la abscisa x∗ de dicho extremo mediante dos iteraciones del
m´etodo de Newton-Raphson. Comenzar la recurrencia en el valor m a´ s ade-
cuado x0 = 1 o bien x0 = 2. Indicacio´ n: Ver la figura adjunta.

cosh x - cos x

10

x
0.5 1 1.5 2 2.5 3

Soluci´on. (i) Una condici´on necesaria para que la funci´on g(x) tenga un
ex- tremo relativo en x∗ es que la funci´on f (x) := g′(x) se anule para x = x∗.
Note- mos que f (x) = cosh —x+cos x 3 es una funci´on continua en todo R y
por lo tanto en ∈
particular
C f [1, 2]. Como adem´as− f (1)f (2) =× ( 0,916617)
0,346049 < 0 se puede aplicar el Teorema de Bolzano de manera que se
∈ al menos un valor x∗ (1, 2) tal que f (x∗) = 0.
concluye con la existen- cia de
Una condici´on suficiente para que x∗ sea un extremo/relativo de g es que g′′
(x∗) = 0. Pero esta condici´on est´a garantizada puesto que si las curvas y =
sinh x, y = sin x no tiene ningu´n punto de corte en el intervalo [1, 2]
entonces la funci´on g′′(x) = sinh x − sin x no
158 Soluci´on de Ecuaciones no
Lineales
tiene ningu´ n zero en dicho intervalo.

(ii) El extremo relativo x∗ de g(x) viene dado por la ra´ız de la funci´on f (x)
en el intervalo [1, 2]. Recordemos que ∈ si f 2[1, 2] y adem´as f ′′(x) tiene
signo con- stante ∈ para toda x [1, C2] entonces la recta tangente trazada
desde el extremo del intervalo [1, 2] adecuada para hallar la ra´ız de f por el
m´etodo de Newton- Raphson es el extremo p (con p = 1 o p = 2) en el cual
sigf (p) = sigf ′′(p). A partir de la gr´afica adjunta en el enunciado se
comprueba que f ′′(x) tiene signo constante en [1, 2]. Adem´as, como f (2) =
0,3460 > 0 y f ′′(2) = 4,1783 > 0
resulta p = 2.

La iteraci´on de Newton-Raphson viene dada por

x
x
= − f(xn) = x — cosh xn + cos xn − 3
, x =2
n+1 n
f ′(xn) n
sinh x − sin x 0

de donde se obtiene la siguiente tabla

n xn f (xn)
0 2 0.346049
1 1.87266 0.0323507
2 1.8581 0.000382472

La aproximaci´on buscada es x∗ ≈ x2 = 1,8581.

Problema 8.4 La ecuaci´on de ondas bidimensional c2(uxx + uyy) = utt, donde

c es la velocidad constante de propagaci´on de la onda, se utiliza en ingenier
´ıa mec´anica para modelar las vibraciones u(x, y, t) de una placa el´astica
rectangu- lar en las coordenadas espaciales (x, y) para tiempo t. Tomando el
referencia cartesiano x − y con origen en el v´ertice inferior izquierdo de la placa
y sabiendo que, para un cierto instante de tiempo, el despazamiento de la placa
respecto de su posici´on de equilibrio viene dado por la funci´on

u(x, y) = sin x + sin y + cos(x + y) ,

se pretende hallar en qu´e punto (x∗, y∗) de la placa cuadrada de lado 2π se

tiene un m´ınimo relativo de u(x, y).
(i) Aproximar dicho punto mediante dos iteraciones del m´etodo de Newton,
sabiendo que (x∗, y∗) ∈ [4π/3, 7π/4] ×[4π/3, 7π/4]. Comenzar la iteracio
en
´n el extremo inferior izquierdo del anterior cuadrado.
(ii) Hallar los errores absolutos cometido en el anterior apartado sabiendo que
la soluci´on exacta es (x∗, y ∗ ) = (3π/2, 3π/2).
(iii) Demostrar que la soluci´on hallada es la u´nica.

Soluci´on. (i) Puesto que se ha de hallar un m´ınimo relativo de la funci´on

u(x, y), en primer lugar obtendremos el sistema de ecuaciones que debe verificar
8.2 Problemas resueltos 159

dicho punto, es decir,∇u(x, y) := (ux(x, y), uy(x, y)) = (0, 0). Calculando las
derivadas paciales de u se tiene

ux(x, y) = cos x − sin(x + y) = 0 , uy(x, y) = cos y − sin(x + y) = 0 . (8.2)

La matriz Jacobiana de la aplicaci´on ∇u calculada en el punto (xn, yn) es

uxx(xn, yn) uxy(xn, yn)

D∇u(xn, yn) =
uyx(xn, yn) uyy(xn, yn)
— sin xn − cos(xn + yn) − cos(xn + yn)
— n n − n − n n
= ,
siendo su matriz inversa cos(x + y ) sin y cos(x + y )

D−1
∇ n 1 − sin yn − cos(xn + yn) cos(xn + yn)
(xu n ,y )= ,
δ cos(xn + yn) − sin xn − cos(xn + yn)

donde δ := det[D∇u(xn, yn)] = sin xn cos(xn + y n ) + sin yn cos(xn + y n ) +

sin xn sin yn.
∇ u(x,
La iteraci´on de Newton para aproximar una soluci´on del sistema
y) = (0, 0) viene dada por

xn+1 xn ux(xn, yn)

yn+1 = − D−1 (xn, yn) . (8.3)
yn uy(xn, yn)
∇
u
En el primer paso, tomando el valor inicial (x0, y0) = (4π/3, 4π/3) se tiene
2 √
D−1 (x ,y ) = √
1 + 3 −1 √
0 0 ,
∇u 3+2 −1 1+ 3
3
√ √
(ux(x0, y0), uy(x0, y0)) = (−[1 + 3]/2, −[1 + 3]/2) ,

de manera que el primer iterado vale

4π √ $ √ %
2 1+ 3 −1 − 1+ 3
x1 = 3 — √ √ 1+2√
4π −1 1+ 3
y1 3 3+2 3 3
— 2
√
9+12π+(3+8π) 3
√ 4,92084
⎛ ⎞
9+6 3 √
9+12π+(3+8π) 3
.
√ 4,92084
= ⎝ 9+6 3 ⎠
≈

Los c´alculos en el segundo paso son los siguientes.

D−1 (x1, y1) = D−1 (4,92084, 4,92084)
∇u ∇u
0,689186 −0,33294
−0,33294 0,689186
= ,
160 Soluci´on de Ecuaciones no
(ux(x1, y1), uy(x1, y1)) Lineales
= (ux(4,92084, 4,92084), uy(4,92084, 4,92084))
= (0,611877, 0,611877) .