Práctico 1
Ejercicio 9: Una recurrencia inestable
Apellido: Arriola
Nombre: Iván
Cédula de Identidad: 5.536.679-6
1
Enunciado
R1
Los números pn = 0
xn ex dx satisfacen p1 > p2 > . . . > 0 y la relación de recurrencia
pn+1 = e − (n + 1)pn , p1 = 1.
(a) Demostrar la relación de recurrencia.
(b) Escribir un programa en Octave que genere los primeros 20 valores de pn y explicar
por qué los resultados obtenidos no satisfacen la desigualdad p1 > p2 > . . . > 0.
(c) Tomar p20 = 18 y usar la relación de recurrencia para calcular p19 , . . . , p1 . ¿Los
números obtenidos satisfacen las desigualdades 1 = p1 > p2 > . . . > 0? Explicar la
diferencia entre los dos procedimientos en términos de estabilidad numérica.
Repetir con p20 = 20 y p20 = 100, y explicar qué ocurre y por qué.
Soluciones
Demostración de la relación de recurrencia
Hipótesis: Z 1
+
n∈Z ⇒ pn = xn ex dx
0
Tesis: Cumple la siguiente relación de recurrencia:
p1 = 1,
pn+1 = e − (n + 1)pn .
Para demostrar la relación de recurrencia, procederemos por inducción.
Caso base:
Tesis: Se cumple que p1 = 1.
1 ∈ Z+ , entonces sabemos que: Z 1
p1 = x1 ex dx.
0
La primitiva de xex se puede resolver por integración por partes:
(Aplicando integración por partes)
u=x ⇒ du = dx
dv = ex dx ⇒ v = ex
Z Z
⇒ xe dx = xe − ex dx = xex − ex + C = ex (x − 1) + C
x x
Aplicando el teorema de Barrow, sabemos que:
Z 1
p1 = xex dx = [ex (x − 1)]10
0
Calculamos:
= e1 (1 − 1) − e0 (0 − 1) = e(0) − (1)(−1) = 0 + 1 = 1
Por lo tanto, p1 = 1 y eso cumple con la relación de recurrencia
2
Caso inductivo:
Hipótesis: La siguiente igualdad se cumple:
pk = e − kpk−1
Tesis: La siguiente igualdad se cumple:
pk+1 = e − (k + 1)pk
Demostración:
Usamos la definición de pk+1 : Z 1
pk+1 = xk+1 ex dx
0
Aplicando integración por partes, elegimos:
u = xk+1 ⇒ du = (k + 1)xk dx
dv = ex dx ⇒ v = ex
Entonces tenemos:
1
Z 1
xk+1 ex 0 ex (k + 1)xk dx
pk+1 = −
0
Ahora evaluamos en x = 1 y x = 0 y sacamos la constante de la integral:
Z 1
k+1 1 k+1 0
xk ex dx
pk+1 = 1 e − 0 e − (k + 1)
0
Ası́ que, la evaluación se convierte en:
Z 1
pk+1 = e − (k + 1) xk ex dx
0
R1
Como sabemos que pk = 0
xk ex dx, sustituimos esto en la ecuación:
pk+1 = e − (k + 1)pk
Esto concluye la demostración de que la relación es cierta para n = k + 1.
Conclusión:
Por lo tanto, como se cumple para el caso base y para el caso inductivo, por el principio
de inducción, la relación se cumple para todo n ≥ 1.
(b) Programa en Octave
La relación de recurrencia pn se imlementa en Octave de la siguiente manera:
3
1 % p_ { n } = e - ( n ) p_ {n -1}.
2 function res = pn ( n )
3 if n < 1
4 error ( ’n debe ser mayor o igual a 1 ’) ;
5 end
6 if n == 1
7 res = 1;
8 else
9 res = exp (1) - ( n ) * pn ( n - 1) ;
10 end
11 end
Listing 1: Implementación de la relación de recurrencia en pn.m
A continuación, se presenta el código en Octave que genera los primeros 20 valores de pn :
1 diary ( ’ salida . txt ’) ; % Iniciar el registro de la salida en el archivo
’ salida . txt ’
2
3 k = 20; % N m e r o de t r m i n o s a calcular
4 p = zeros (1 , k ) ;
5
6 % Calcular los valores de p_n usando la r e l a c i n de recurrencia
7 for n = 1: k
8 p ( n ) = pn ( n ) ;
9 end
10
11 % Mostrar los resultados
12 disp ( ’ Valores de p_n : ’) ;
13 disp ( p ) ;
14
15 diary off ; % Finalizar el registro
Listing 2: Código principal en entregable1.m
La salida del programa es:
1 Valores de p_n :
2 Columns 1 through 7:
3
4 1.0000 0.7183 0.5634 0.4645 0.3956 0.3447
0.3055
5
6 Columns 8 through 14:
7
8 0.2744 0.2490 0.2280 0.2103 0.1951 0.1820
0.1705
9
10 Columns 15 through 20:
11
12 0.1605 0.1503 0.1624 -0.2043 6.5991 -129.2637
Listing 3: Salida del programa
Explicación de los resultados:
Los valores obtenidos no satisfacen la desigualdad p1 > p2 > . . . > 0 a partir de n = 17.
Esto se debe a la inestabilidad numérica introducida por el método de recurrencia, que
amplifica los errores de cálculo a medida que se computan valores sucesivos.
4
Cálculo de la propagación del error
La propagación del error en la relación de recurrencia para pn se puede expresar de la
siguiente manera:
P n+1 = e − (n + 1) · P n
Sustituyendo P n :
P n = (1 + ϵpn )Pn
Por lo tanto, podemos reescribir la primera ecuación como:
P n+1 = e − (n + 1) · (1 + ϵpn )Pn
