INSTRUCCIONES: Puede utilizar todo tipo de material auxiliar y calculadora pro-

gramable. No es suficiente poner la solución, en cada ejercicio debe especificar y exponer

los pasos intermedios realizados.

1. (2.5 puntos) Uno de los fenómenos atmosféricos naturales que atrae bas-
tante atracción pública regularmente es el evento de El Niño en el Pacífico
ecuatorial, el cual ocurre en intervalos de 2-7 años acarreando, la mayoría
de las veces, consecuencias dramáticas.
Un modelo pedagógico muy sencillo, pero válido para introducir el fenó-
meno dinámico de El Niño es el siguiente
dT
= T (t) − [T (t)]3 − αT (t − δ)
dt
En esta ecuación diferencial T (t) representa la anomalía de la tempe-
ratura de la superficie del mar medida en un cierto tiempo t, es decir,
la desviación de la temperatura respecto de una temperatura promedio
definida “a largo plazo” (a escala de tiempos largos). Además aparecen
dos parámetros, δ representa un tiempo de retardo y α es una constante.
Realice la integración numérica de la ecuación diferencial con el méto-
do de Runge-Kutta de orden dos que corresponde al método de Euler
modificado, también llamado método de predicción-corrección de Euler.
Utilice un paso temporal de dt = 0.1 para el conjunto de parámetros
α = 0.75, δ = 2 y con la condición T (t) = 0.55 para t ≤ 0. Observe que
esta condición implica que T (t) es constante e igual a 0.55 en tiempos
anteriores o iguales a t = 0. Describa detalladamente la implementación
del algoritmo para las dos primeras iteraciones.
Solución:
La ecuación diferencial a resolver es
dT
= T − T 3 − αT (t − δ) T (t) = 0.55 t≤0 (1)
dt
Podemos definir la función

g(t, T (t), T (t − δ)) = T (t) − [T (t)]3 − αT (t − δ)

Obsérvese que esta función g en realidad no depende explícitamente del tiempo pero sí
de las funciones T (t) y de T (t − δ).
En cada iteración i, las fórmulas generales del método de Euler mejorado (predicción-
corrección de Euler), también llamado Runge-Kutta de orden 2 para los coeficientes
a2 = b21 = 1 y c1 = c2 = 12 , son
1
Ti+1 = Ti + (k1 + k2 ) (2)
2
donde se hace uso de las siguientes definiciones de los coeficientes

k1 = hf (ti, Ti, T (ti − δ))

k2 = hf (ti + h, Ti + k1, T (ti + h − δ))

Partimos de T0 = T (0) con dt = h. La primera evaluación de los coeficientes será

k1 = h[Ti − Ti3 − αT (ti − δ)]

k2 = h[(Ti + k1 ) − (Ti + k1 )3 − αT (ti+1 − δ)]

Obsérvese que para calcular k1 debemos conocer el valor de T (ti − δ), para k2 se necesita
conocer T (ti+1 − δ). Evidentemente estos valores habrán de tenerse o aproximarse con
los datos de la solución. Para t < 0 se supone conocida la serie T (t − δ).

Existe un programa “el-nigno-RK2.f” en el que se han implementado los resultados para

los parámetro α = 0.75 y δ = 2 con T (t) = 0.55 para t ≤ 0 para calcular numéricamente
la solución con un paso de tiempo de dt = 0.1.

δ
Como δ = 2 tenemos δi = dt = 20, t − δ = (i − δi)dt.

La primera iteración será

Tj = 0.55 j = −∞, · · · , −2, −1, 0

k1 = h[T0 − T03 − α T (t0 − δ)] = 0.1[0.55 − 0.553 − 0.75 · 0.55] = −0.002887500826
| {z }
T0−δi =T−20

k2 = h[(T0 + k1 ) − (T0 + k1)3 − α T (t1 − δ)]

| {z }
T1−δi =T−19

= 0.1[(0.55 + k1 ) − (0.55 + k1 )3 − 0.75 · 0.55] = −0.0029155

1
T1 = T0 + (k1 + k2) = 0.55 + 0.5(−0.0028875 + −0.0029155) = 0.54709846 (3)
2
Para la iteración segunda

k1 = h[T1 − T13 − α T (t1 − δ)] = −0.0029157267

| {z }
T1−δi =T−19

k2 = h[(T1 + k1 ) − (T1 + k1 )3 − α T (t2 − δ)] = −0.00294687455

| {z }
T2−δi =T−18
1
T1 = T0 + (k1 + k2 ) = 0.55 + 0.5(−0.0029157267 − 0.00294687455) = 0.54416716
2
(4)

2. (1.25 puntos) Sea el sistema de ecuaciones no lineales escrito de la forma

x2 − y + 0.5
x =
2
−x − 4y 2 + 8y + 4
2
y =
8
Describa una estrategia iterativa del tipo de Jacobi y otras dos del tipo
de Gauss-Seidel (distintas entre sí) para resolver el sistema de ecuaciones
partiendo de una semilla cualquiera (x0, y0 ).

Solución: Definamos un par de funciones f (x, y) y g(x, y)

x2 − y + 0.5
x = f (x, y) ≡
2
−x − 4y 2 + 8y + 4
2
y = g(x, y) ≡
8
Para resolver este sistema de dos ecuaciones podemos intentar un método de resolución
iterativo basado en la estrategia de Jacobi y otros dos distintos basados en la estrategia
de Gauss-Seidel.
a) Para la iteración de Jacobi en el cálculo de las dos nuevas variables se utilizan los
valores de la iteración previa, por ello es indistinto si primero se calcula la variable xk+1
o la yk+1, pues el resultado es independiente
x2k − yk + 0.5
xk+1 = f (xk , yk ) ≡
2
−xk − 4yk2 + 8yk + 4
2
yk+1 = g(xk , yk ) ≡
8

b1) Para la iteración de Gauss-Seidel una posibilidad será calcular primero la xk+1 y
pasar a calcular yk+1 usando el valor más actualizado de xk+1 recientemente calculado
en la instrucción anterior
x2k − yk + 0.5
xk+1 = f (xk , yk ) ≡
2
−xk+1 − 4yk2 + 8yk + 4
2
yk+1 = g(xk+1, yk ) ≡
8

b2) Una segunda manera de realizar una iteración Gauss-Seidel sería primero calcular
la yk+1 para luego calcular xk+1 usando el valor más actualizado de yk+1 recientemente
calculado en la instrucción anterior
−x2k − 4yk2 + 8yk + 4
yk+1 = g(xk , yk ) ≡
8
2
x − yk+1 + 0.5
xk+1 = f (xk , yk+1) ≡ k
2

3. (1.25 puntos) Indique cuáles deben ser los pesos de una regla de cuadra-
tura que integre de forma exacta polinomios de orden 2 o menor en el
intervalo [0, 1] y que se apoya en los puntos x0 = 0.1, x1 = 0.5 y x2 = 0.9.

Solución:

Buscamos una regla de cuadratura de la forma

Z 1
dxf (x) = af (x0) + bf (x1) + cf (x2)
0

Imponemos que se calcule de forma exacta la integral de los monomios 1, x y x2 en el

intervalo apoyándonos en los puntos que se indican en el enunciado:
Z 1
1 dx = 1 = a · 1 + b · 1 + c · 1 → a + b + c = 1
0
Z 1
1 1
xdx = = a · 0.1 + b · 0.5 + c · 0.9 → 0.1a + 0.5b + 0.9c =
2 2
Z 01
1 1
x2dx = = a · 0.12 + b · 0.52 + c · 0, 92 → 0.01a + 0.25b + 0.81c =
0 3 3
Lo que nos dice que debemos resolver el siguiente sistema de ecuaciones:

a+b+c=1
1
0.1a + 0.5b + 0.9c =
2
1
0.01a + 0.25b + 0.81c =
3

Resolvemos el sistema usando el método de sustitución, despejamos a de la primera de

las ecuaciones: a = 1−b−c y la sustituimos en la segunda y en la tercera. Tras simplificar
estas quedan como:

b + 2c = 1
97
24b + 80c =
3
despejando c de la primera y sustituyendo en la segunda obtenemos: 24(1 − 2c) + 80c =
97/3 → 32c = 25/3. Es decir que c = 25/96. Sustituyendo hacia atrás encontramos que
los valores de b y a son 46/96 y 25/96 respectivamente.
Comprobamos que la solución es correcta, por ejemplo, calculando la integral de f (x) =
1 + x2 .
Z 1
x3 1 4
dx(1 + x2) = (x + )|10 = 1 + =
0 3 3 3
Usando la regla de integración que hemos propuesto:
Z 1
dx(1 + x2) = af (x0) + bf (x1) + cf (x2)
0
25 46 25
= (1 + 0, 12) + (1 + 0, 52) + (1 + 0, 92)
96 96 96
25 46 25
= (1, 01) + (1, 25) + (1, 81) = 0, 2630208 + 0, 5989580 + 0, 4713541
96 96 96
= 1, 3333332

Expresión que es exacta y que solo difiere con la solución exacta el octavo decimal debido
a los redondeos hechos con la calculadora.

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

no debe contestar a ninguno de los siguientes problemas.
 4. (2.5 puntos) Obtenga una fórmula O(h4 ) de cinco puntos para aproximar
f ′(x0) que utilice f (x0 − h), f (x0), f (x0 + h), f (x0 + 2h) y f (x0 + 3h). Para
ello, parta de la expresión

Af (x0 − h) + Bf (x0 + h) + Cf (x0 + 2h) + Df (x0 + 3h)

desarrolle en serie de Taylor alrededor de x0 hasta cuarto orden y elija

adecuadamente los coeficientes A, B, C y D para poder obtener la ex-
presión válida para la derivada primera en términos de los cinco puntos
x0 − h, x0, x0 + h, x0 + 2h y x0 + 3h.

Solución:

Partamos de la expresión Af (x0 − h) + Bf (x0 + h) + Cf (x0 + 2h) + Df (x0 + 3h) =

Ef (x0) + f ′(x0) y desarrollemos en serie de Taylor en un entorno del punto x0 todas las
funciones hasta términos de orden cuatro:
1 1 1
f (x0 − h) = f (x0) − f ′ (x0)h + f ′′(x0)h2 − f ′′′(x0)h3 + f iv (x0)h4
2 6 24
1 1 1
f (x0 + h) = f (x0) + f ′ (x0)h + f ′′(x0)h2 + f ′′′(x0)h3 + f iv (x0)h4
2 6 24
4 2
f (x0 + 2h) = f (x0) + 2f ′(x0)h + 2f ′′(x0)h2 + f ′′′(x0)h3 + f iv (x0)h4
3 3
9 9 27
f (x0 + 3h) = f (x0) + 3f ′(x0)h + f ′′(x0)h2 + f ′′′(x0)h3 + f iv (x0)h4
2 2 8
Entonces la expresión propuesta queda

Af (x0 − h) + Bf (x0 + h) + Cf (x0 + 2h) + Df (x0 + 3h) =

(A + B + C + D)f (x0)
+(−A + B + 2C + 3D)f ′(x0)h
1
+ (A + B + 4C + 9D)f ′′(x0)h2
2
1
+ (−A + B + 8C + 27D)f ′′′(x0)h3
6
1
+ (A + B + 16C + 81D)f iv (x0)h4 + O(h5 )
24

Si esta expresión debe responder a una aproximación de cinco puntos válida para f ′ (x0)
debemos imponer que sean nulos los sumandos que implican a f ′′(x0), f ′′′(x0), f iv (x0).
Eso se puede obtener si los 4 coeficientes A, B, C, D satisfacen el sistema de 3 ecuaciones

A + B + 4C + 9D = 0
−A + B + 8C + 27D = 0
A + B + 16C + 81D = 0

Nos faltaría una ecuación más para resolver el sistema. Elegimos que en la expresión
propuesta el término (−A + B + 2C + 3D)f ′(x0)h sea realmente f ′ (x0), es decir, que se
debe imponer (−A + B + 2C + 3D)h = 1.

El sistema de cuatro ecuaciones con cuatro incógnitas a resolver es

−A + B + 2C + 3D = 1/h
A + B + 4C + 9D = 0
−A + B + 8C + 27D = 0
A + B + 16C + 81D = 0

con solución
1 3 1 1
A=− B= C=− D=
4h 2h 2h 12h
5
Con este resultado, el término que acompaña a f (x0) con A + B + C + D = 6h y nuestra
expresión queda
5
Af (x0 − h) + Bf (x0 + h) + Cf (x0 + 2h) + Df (x0 + 3h) = f (x0) + f ′(x0) + O(h5 )
6h

5. (2.5 puntos) Aproxime la función f (x) = 3xex − e2x en el punto x = 1.03

mediante un polinomio interpolante de Hermite de grado menor o igual
que 3, utilizando los valores de f (x) y f ′(x) en x0 = 1 y x1 = 1.05.

Solución:

El polinomio interpolante de Hermite más general para 2 puntos es:

1 h
X i
′
H5 (x) = f (xj )H1,j (x) + f (xj )Ĥ1,j (x)
j=0

En primer lugar necesitamos calcular los valores numéricos de la función, f (x), y de su

derivada, f ′(x) = (1 + x)3ex − 2e2x , que necesitaremos más adelante:

f (x0) = 0.765789 f (x1) = 0.835431

f ′(x0) = 1.53158 f ′(x1) = 1.24221
Los correspondientes coeficientes de Lagrange, así como sus derivadas tienen la siguiente
forma;
L0(x) = xx−x 1
0 −x1
= 21 − 20x L1 (x) = xx−x 0
1 −x0
= 20x − 20
′ 1 ′ 1
L0(x) = x0 −x1 = −20 L1 (x) = x1 −x0 = 20
Sus valores numéricos correspondientes son L0 (x0) = 1, L0 (x1) = 0, L1 (x0) = 0,
L1(x1) = 1. Necesitaremos los siguientes coeficientes de Hermite:

H1,0 (x) = [1 − 2(x − 1)L′0(x0)]L20(x) = −17199. + 50400.x − 49200.x2 + 16000.x3

H1,1 (x) = [1 − 2(x − 1)L′1(x1)]L21(x) = 17200. − 50400.x + 49200.x2 − 16000.x3
Ĥ1,0 (x) = (x − x0)L20(x) = −441. + 1281.x − 1240.x2 + 400.x3
Ĥ1,1 (x) = (x − x1)L21(x) = −420. + 1240.x − 1220.x2 + 400.x3

Finalmente, sustituyendo en la expresión de H5(x) y agrupando por términos obtenemos:

H5 (x) = 1198.6 − 3509.94x + 3426.37x2 − 1114.27x3

El valor en el punto deseado es H5(1.03) = 0.810917, frente al valor real que es f (1.03) =
0.809324