INSTRUCCIONES: Puede utilizar todo tipo de material auxiliar y calculadora programable.

No
es suficiente poner la solución, en cada ejercicio deben especificar y exponer los pasos intermedios
realizados.

1. (2.5 puntos) Un método iterativo para encontrar las raíces de una ecuación de una
variable f (x) = 0 propone realizar la búsqueda
2f (xn )f ′ (xn )
xn+1 = xn −
2[f ′ (xn )]2 − f (xn )f ′′ (xn )
partiendo de una semilla inicial x0 y entendiendo que la notación ′
representa la
derivada con respecto a x.
Calcule
√ el orden de convergencia de este método cuando se aplica para encontrar el
valor 3. Para ello, haga uso del error cometido en cada iteración expresado como
√
ei = |xi − 3|, para i = 0, 1, 2, . . .

Ayuda:
n
n
X n!
(a + b) = an−k bk
k=0
k!(n − k)!

√
a) Con el algoritmo se trata de obtener el valor de 3 = 1.73205080. Dado que el mismo resuelve una
ecuación de la forma f (x) = 0, podemos plantear como función, por ejemplo, la siguiente f (x) = x2 − 3.
Dado que f ′ (x) = 2x y f ′′ (x) = 2 el algoritmo iterativo queda

2f (xn )f ′ (xn ) 4x (x2n − 3)

xn+1 = xn − = xn −
2[f ′ (xn )]2 − f (xn )f ′′ (xn ) 8x2n − 2 (x2 − 3)
Por lo tanto, podemos ver este método como uno de punto fijo xn+1 = g(xn ) donde la función g(x) es

4x (x2 − 3) (x3 + 9x)

g(x) = x − =
8x2 − 2 (x2 − 3) 3(x2 + 1)

Utilizando la expresión (a − b)3 = a3 − 3a2 b + 3ab2 − b3 , el error cometido en la iteración n + 1 se puede

expresar en términos del error en la iteración anterior n
√
√ (x3n + 9xn ) √ (x3n + 9xn ) − 33(x2n + 1)
en+1 = xn+1 − 3 = − 3=
3(x2n + 1) 3(x2n + 1)
√
3
xn − 3 e3n
= =
3(x2n + 1) 3(x2n + 1)

De aquí se sigue el orden de convergencia del método es cúbico en+1 ∝ e3n porque se puede calcular
|en+1 | 1 1
lı́m 3 = n→∞
lı́m 2
=
n→∞ |en | 3xn + 3 12
 √
dado que en el proceso de convergencia se tiene lı́mn→∞ xn = 3, demostrándose que el método iterativo
presenta un orden de convergencia cúbico.

2. (2.5 puntos) Consideremos la siguiente ecuación integral en x ∈ [0, 1]

Z 1

u (x) = x + xt + t2 u (t) dt
0

Utilice el método de cuadratura de Gauss con n = 5 para transformar la ecuación

anterior en un sistema de ecuaciones con 5 incógnitas. ¿Se trata de un sistema lineal
o no lineal? Discuta qué método utilizaría para resolver dicho sistema.

Solución:

Para aplicar la cuadratura de Gauss-Legendre debemos transformar el intervalo [0, 1] en [−1, 1]. Para
ello hacemos el cambio de variable
1 + t′ 1
t= , dt = dt′
2 2
En la siguiente tabla se muestran las raíces del polinomio de Legrende de grado 5 y los pesos de la
cuadratura de Gauss-Legendre

ti raíces en [0, 1] t′i raíces en [−1, 1] wi , pesos de ponderación

t1 = 0.046910 t′1 = −0.906180 w1 = 0.236927
t2 = 0.230765 t′2 = −0.538469 w2 = 0.478629
t3 = 0.5 t′3 = 0. w3 = 0.568889
t4 = 0.769235 t′4 = 0.538469 w4 = 0.478629
t5 = 0.953090 t′5 = 0.906180 w5 = 0.236927

(Vea Tabla 5.14 del libro de [Gerald] pag 392, donde aparecen los factores de ponderación así como los
valores de las raíces del polinomio de Legrendre de grado n = 5).

Aplicamos el cambio de variable a la ecuación integral:

Z "    2 #  
1 1 1 + t′ 1 + t′ 1 + t′
u (x) = x + x + u dt′ ,
2 −1 2 2 2

y mediante la cuadratura de Gauss-Legendre

1



u (x) ≈ x + w1 xt1 + t21 u (t1 ) + w2 xt2 + t22 u (t2 ) + ... + w5 xt5 + t25 u (t5 )
2
En esta ecuación, debemos resolver la función u(x), y en el lado de la derecha aparecen como incógnitas
5 valores de la función u en 5 puntos determinados u1 ≡ u(t1 ), u2 ≡ u(t2 ), · · · , u5 ≡ u(t5 ), es decir,
identificamos u(ti ) por ui , que va a representar el valor aproximado de la solución en el punto ti , i =
1, 2, · · · , 5. Si sustituimos en la ecuación x por los diferentes nodos t1 , t2 , t3 , t4 y t5 obtenemos el siguiente
sistema lineal de 5 ecuaciones con 5 incógnitas:
1


u1 = t1 + 2w1 t21 u1 + w2 t1 t2 + t22 u2 + ... + w5 t1 t5 + t25 u5
2
1


u2 = t2 + w1 t2 t1 + t21 u1 + 2w2 t22 u2 + ... + w5 t2 t5 + t25 u5
2
..
.
1

u5 = t5 + w1 t5 t1 + t21 u1 + w2 (t5 t2 + t22 )u2 + ... + 2w5t25 u5
2
Podemos expresar este sistema lineal de 5 ecuaciones en la forma Au = b donde
 1     
w1 t21 − 1 w (t t + t22 ) · · · 21 w5 (t1 t5 + t25 )
2 2 1 2
u1 −t1
 12 w1 (t2 t1 + t21 ) w2 t22 − 1 · · · 21 w5 (t2 t5 + t25 )   u   −t 
A= .. .. ..  , u =  .2  , b =  . 2 
 . . ··· .   ..   .. 
1
w (t t + t21 ) 21 w2 (t5 t2 + t22 ) · · ·
2 1 5 1
w5 t25 − 1 u5 −t5

El método de resolución podría ser, por ejemplo, el de Crout o el de Gaus-Jordan.

Recuerde que si ha realizado alguna actividad de evaluación continua no debe

contestar a ninguno de los siguientes problemas.

3. (2.5 puntos) Para la modelización del desarrollo de una enfermedad contagiosa en el

tiempo t se utiliza una función de crecimiento logístico
1
C(t) =
a + be−2t
Plantee el sistema de ecuaciones que necesita resolver para calcular el aproximante de
Padé de la función C(t) alrededor de t0 6= 0 con grado uno en el numerador y también
grado uno en el denominador. Se recomienda escribir el aproximante como
a0 + a1 (t − t0 )
1 + b1 (t − t0 )

Solución:
1
a) Nuestra función es C(t) = a+be−2t
tiene como desarrollo de Taylor alrededor del punto t0
∞
1 X f (n) (t0 )
C(t) ≡ = (t − t0 )n
a + be−2t n=0
n!

Buscamos la aproximación de Padé a C(t) que sea de la forma

a0 + a1 (t − t0 ) + · · · + an (t − t0 )n
C(t) ≈ RN (t) = , N =n+m
1 + b1 (t − t0 ) + · · · + bm (t − t0 )m
En este caso nos piden encontrar el aproximante R2 (t) con n = 1 y m = 1.

a0 + a1 (t − t0 )
R2 (t) =
1 + b1 (t − t0 )
Se empieza utilizando el desarrollo de la serie de Taylor alrededor de t0 hasta términos de orden N, en
nuestro caso N = 2, corresponde a la expresión
1 2be−2t0 2 (be2t0 (ae2t0 − b))
CT2 (t) = + 2 (t − t0 ) − 3 (t − t0 )2
a + be−2t0
| {z } |(a + be
−2t 0 ) 2t
(ae + b)
0
{z } | {z }
c0 c1 c2

de donde hemos calculado los coeficiente del desarrollo de Taylor de grado 2 haciendo las derivadas de
la función

′ 2be−2t
C (t) =
(a + be−2t )2
8b2 e−4t 4be−2t
C ′′ (t) = −
(a + be−2t )3 (a + be−2t )2
Por otra parte
a0 + a1 (t − t0 )
R2 (t) =
1 + b1 (t − t0 )
a1 − a0 b1
R2′ (t) =
[1 + b1 (t − t0 )]2
−2b1 (a1 − a0 b1 )
R2′′ (t) =
[1 + b1 (t − t0 )]3

La construcción de Padé requiere de la obtención de 3 parámetros a0 , a1 , b1 , lo cual implica obtener 3

ecuaciones. Se impone que coincidan C(t) = R2 (t) para t0 y también las N = 2 primeras derivadas de
ambas funciones coincidan en t0 .
1
C(t0 ) − R2 (t0 ) = 0 →a0 =
a + be−2t0
2be−2t0
C ′ (t0 ) − R2′ (t0 ) = 0 → = a1 − a0 b1
(a + be−2t0 )2
8b2 e−4t0 4be−2t0
C ′′ (t0 ) − R2′′ (t0 ) = 0 → − = −2b1 (a1 − a0 b1 )
(a + be−2t0 )3 (a + be−2t0 )2
El primer coeficiente a0 lo dicta la primera ecuación. Insertando el valor de a1 − a0 b1 de la segunda
ecuación en la tercera, se obtiene el coeficiente b1
a − be−2t0
b1 =
a + be−2t0
de donde, se obtiene que el tercer coeficiente del Padé es
1
a1 =
a + be−2t0
 4. (2.5 puntos) El agua almacenada en un tanque cilíndrico, cuya base reposa en el suelo,
fluye por un orificio circular de radio r0 a una velocidad
√
dy 2 2gy
= −πr0
dt A(y)
donde r0 = 0.1 m es el radio del orificio, y es la altura de la superficie del líquido con
respecto a la base del cilindro, g = 10 m/s2 es la aceleración de la gravedad y A(y) es
el área transversal del depósito a una altura y. Supongamos que la altura inicial del
agua es de h0 = 10 m y el volumen inicial de agua es V0 = 100π m3 .
Halle la altura del nivel del agua y al cabo de un tiempo ∆t utilizando el método de
Euler modificado con un paso temporal de ∆t = 0.1 segundos. Calcule el error cometido
con la aproximación numérica en relación a la solución exacta del problema.
Nota: No es necesario sustituir los valores numéricos.

Solución:
Nos dicen que inicialmente el líquido tiene un volumen V0 a una altura h0 . Como se trata
√ de un cilindro
de radio R, V0 = πR2 h0 , se tiene que el radio del cilindro del que hablamos es R = 10 m.
Recordemos que se miden las alturas respecto de la base del cilindro de radio R. El área transversal a
la altura y en un cilindro es constante, e independiente de y

A(y) = πR2 = 10π

Por lo que la ecuación diferencial a resolver es

√
dy 2gy √
=− = −c y
dt 1000
√
20
con c = 1000
.
√
La solución exacta de la ecuación diferencial ẏ = −c y es
√ c 2
y(t) = y0 − t
2
donde se ha hecho uso de la condición inicial dada en el enunciado y0 = h0 = 10.
El método de Euler modificado nos dicta
∆t 1/2
yn+1/2 = yn − cy
2 n
1/2 ∆t 1/2 1/2
yn+1 = yn − ∆tcyn+1/2 = yn − ∆tc(yn − cy )
2 n
y por lo tanto, al primer paso ∆t = 0.1
∆t 1/2
y0+1/2 = y0 + cy
2 0
∆t 1/2 1/2
y1 = y0 + ∆tc(y0 + cy )
2 0
El error cometido para este ∆t es
∆t 1/2 1/2  c 1/2
2
E∆t = |y1 − y(∆t)| = |y0 + ∆tc(y0 + cy0 ) − ∆t + y0 |
2 2