Original 1P

Ejercicio 1 Sea A una matriz cuadrada de dimensión n cuyos coeficientes

son aij = min{i, j} para i, j = 1, · · · , n. Calcular el determinante de A
calculando para ello la factorización de Cholesky.

Solución: En la primera columna se obtiene l11 = 1, li,1 = 1 para i =

2, · · · , n. Por inducción en j se supone que li,j−1√= 1 si i ≥ j − 1 y se
probará que li,j = 1 para i ≥ j. En efecto ljj = j − (j − 1) = 1 y lji =
i−(i−1)
1
= 1. Consecuentemente, lji = 1 si j ≤ i para todo i, j = 1, · · · , n.
Finalmente puesto que A = Lt L se tiene que det(A) = det(L)2 = 1.

Ejercicio 2 Sea e un vector unitario de IR 3 y A = I − 23 eet . ¿Si se utiliza el

método de la potencia iterada (sin normalización) qué valor se obtiene como
lı́mite, partiendo de x(0) tal que a = x(0) · e ̸= 0.?

Solución: La sucesión de Krylov generada partiendo de x(0) es la siguiente:

2
Ax(0) = x(0) − ae,
3
8
A2 x(0) = x(0) − ae,
9
26
A3 x(0) = x(0) − ae,
27
···
3k − 1
Ak x(0) = x(0) − ae
3k
como fácilmente puede probarse por inducción. Consecuentemente
( k+1 (0) ) ( )
x1 − 1 − 3−k e1
(0)
A x 1
lim = lim (0) = 1.
k→∞ (Ak x(0) ) k→∞ x − (1 − 3−k+1 ) e
1 1 1

Ejercicio 3 Se consideran los puntos a = x0 < x1 < · · · < xn−1 < xn = b

en el intervalo [a, b]. Probar que si p es un polinomio de grado p, con p > n,
la diferencia dividida p[x0 , x1 , · · · , xn , x] es un polinomio de grado p − n − 1.
Calcular p[0, 1/3, 1/2, 1, x] si p(x) = x2 .

Solución: p[0, 1/3, 1/2, 1, x] = 0.

Ejercicio 4 Se considera el producto escalar en C[0, 1] definido por

∫ 1
⟨f, g⟩ = f (x)g(x)w(x)dx
0

donde w(x) = − ln(x). Calcular los tres primeros términos de la sucesión de

polinomios ortogonales con el coeficiente principal 1.

Solución: Si se usa la fórmula de recurrencia de tres términos se obtiene

p0 (x) = 1
a1 = 14 b1 = 0 p1 (x) = x − 14
a2 = 13
28
7
b2 = 144 p2 (x) = x2 − 57 x − 17
252
Ejercicio 5 ¿Qué propiedades destacables tienen los polinomios estandari-
zados de Chebyshev (coeficiente principal igual a 1) respecto a la norma
uniforme en [−1, 1].? Justificar la respuesta. Calcular los cuatro primeros
polinomios estandarizados de Chebyshev.

Solución: De acuerdo con la fórmula de recurrencia se tiene que

T̂0 = 1,
T̂1 = x,
1
T̂2 = x2 − ,
2
3
T̂3 = x3 − x.
4
Ejercicio 6 Calcular el polinomio de interpolación de Hermite de la función
f (x) = cos x utilizando {0, π4 , π2 } como puntos de interpolación

Solución:

Ejercicio 7 Para aproximar el valor de la integral

∫ 1
dx
Iε = 1
2 2
0 (ε + x ) 2

donde ε representa un número positivo destinado a tender a cero, se desea

utilizar el método de los trapecios compuesto. Se pide acotar el error de cua-
dratura. ¿Es la estimación encontrada independiente de ε?.

Solución: Sea
1
f (x) = 1 .
(ε2 + x2 ) 2
Sus derivadas segunda y tercera son

d2 f 2x2 − ε2 d3 f 3x(3ε2 − 2x2 )

= − 5 , = − 7 .
dx2 (ε2 + x2 ) 2 dx3 (ε2 + x2 ) 2
La acotación
d2 f
(η) 2 2η 2 − ε2 1
|E1,m | = dx2
h = 5 h2 ≤ 5 h2
12 12(ε2 + η2) 2 6(5/2) 2 ε3

es dependiente de ε. No hay garantı́a de que el error de cuadratura no crezca

cuando ε → 0 salvo que h3 → 0.
ε2
 Original 2P
Ejercicio[ 8 La] función f (x) = xex transforma[ el] intervalo [−1, 0] en el
intervalo − 1e , 0 . La función inversa de f en − 1e , 0 se conoce como función
de Lambert (rama -1). Para evaluar la función de Lambert se[ puede] utilizar el
método de Newton a la ecuación xex = a para cualquier a ∈ − 1e , 0 . Analizar
la convergencia del método cuando se toma como punto inicial x(0) = 0 y
cómo se comporta cuando a tiene un valor muy próximo a − 1e .

Solución: La iteración de Newton es x(k+1) = g(x(k) ) donde

xex − a x2 + ae−x
g(x) = x − = .
(1 + x)ex 1+x

Por otra parte las derivadas segunda y tercera de f son

df d2 f
(x) = (1 + x)ex , (x) = (2 + x)ex ,
dx dx2
de donde se deduce que (f es es] creciente y estrictamente convexa en [−1, 0]
y f (−1) < 0 para a ∈ − 1e , 0 , y la sucesión es monótona decreciente a la
raı́z de acuerdo con el teorema 33. Para a = − 1e , la raı́z es x = −1.
Por otra parte

f (x)f ′′ (x) (xex − a)(2 + x)

g ′ (x) = =
f ′ (x)2 (1 + x)2 ex

Puesto que limx→−1 g ′ (x) = ∞, la velocidad de convergencia disminuye cuan-

do a → − 1e .

Ejercicio 9 Usar una sucesión de polinomios de Sturm para determinar el

número de autovalores de la matriz
 
0 1 0
A =  0 0 −1 
α 1 0

que contiene el intervalo (0, 1) dependiendo del valor del parámetro α > 0.

Solución: El polinomio caracterı́stico de A es

p0 (x) = −x3 + x + α
 Con el algoritmo de Euclides se obtiene la siguiente sucesión de Sturm

p0 (x) = −x3 + x + α,
p1 (x) = −3x2 + 1,
2
p2 (x) = − x + α,
3
27
p3 (x) = 1 − α.
4
Del teorema de Cauchy se deduce que las raı́ces reales están en el intervalo
[−4, 4]. Los signos de la sucesión de Sturm en los valores enteros del intervalo
están dados en la tabla
x p0 p1 p2 p3 no de cambios
0 + + + signo(1 − 27
4
α)
1 + − signo(− 23 + α) signo(1 − 27
4
α)

Consecuentemente, si α ≤ 27
4
el número de raı́ces en [0, 1] es 2. Si 4
27
≤α≤ 2
3
no hay ninguna, y en otro caso hay 1.

Ejercicio 10 Invertir la siguiente matriz

 
1 0 ··· 1
 0 1 ··· 1 
 
B =  .. .. . . .. 
 . . . . 
0 0 ··· 1

usando la fórmula de Sherman-Morrison.

Solución: Se aplica la fórmula de Sherman-Morrison con

A = I, u = (1, 1, · · · , 0)t , v = (0, 0, · · · , 1)t .

De este modo, se obtiene

   
1 0 ··· 1 1 0 · · · −1
 0 1 ··· 1   0 1 · · · −1 
−1    
B =I− .... . . ..  =  .. .. . . .. 
 . . . .   . . . . 
0 0 ··· 1 0 0 ··· 1

Ejercicio 11 Hallar la solución de la ecuación en diferencias

xn+1 − 2xn + xn−1 = 0

tal que x0 = 1 x101 − x100 = 1.

Solución: La ecuación caracterı́stica

λ2 − 2λ + 1 = 0

tiene como única raı́z λ = 1 doble. La solución general de la ecuación en

diferencias es
xn = c1 + c2 n.
Si se imponen las condiciones adicionales, se obtiene

c1 = 1, 1 + 101c2 − (1 + 100c2 ) = 1

de donde se deduce que la solución buscada es

xn = 1 + n.

Ejercicio 12 Se considera el problema de valor inicial

dx
= f (t, x(t)), x(t0 ) = x0
dt
que se quiere aproximar un método multipaso de la forma

xn+1 + bxn−1 + cxn−2 = hfn ,

donde a y b son constantes. Determinar a y b para que el esquema sea con-

sistente. Con los valores obtenidos ¿el esquema es convergente?.

Solución: Suli 358

Ejercicio 13 Si g es una función de clase C p+1 en un intervalo I tal que
di g
(α) = 0, para i = 1, 2, · · · , p,
dxi
dp+1 g
(α) ̸= 0,
dxp+1
para un punto fijo α de g, probar toda sucesión convergente generada el méto-
do de punto fijo asociado a g, tiene p+1 como orden de convergencia. Calcular
el orden de convergencia en el caso g(x) =.

Ejercicio 14 Se desea calcular el mı́nimo de la función

f (x, y) = −2xy − 2x + x2 + 2y 2 .

Determinar el punto que se alcanza después de la primera iteración del méto-

do del gradiente con máximo descenso si se parte del punto (−1, 1).
Solución: En este caso, la dirección de descenso está dada por
d(0) = −∇f (x(0) , y (0) ) = (6, −6)t .
Para encontrar el tamaño del paso ρ0 , se busca el mı́nimo de la función
(0) (0)
g(t) = f (x(0) + td1 , y (0) + td2 ) = f (−1 + 6t, 1 − 6t) = 36t2 − 24t + 3.
Puesto que el mı́nimo de g se alcanza en t = 13 . Finalmente se obtiene que
1 1
(x(1) , y (1) ) = (−1 + 6 , 1 − 6 ) = (1, −1).
3 3
Ejercicio 15 Aproximar las soluciones del problema de contorno
−u′′ = 0, u′ (0) = u′ (1) = 0
mediante un método de diferencias finitas que use una diferencia dividida
progresiva para aproximar la condición de contorno en el extremo lateral
izquierdo y una retrógrada para el derecho. Determinar si la solución del
problema aproximado es única para h = 13 .
Solución: Si se utiliza el esquema usual para discretizar la derivada segunda,
se obtiene el siguiente sistema lineal
ui+1 − 2ui + ui−1
− = 0, i = 1, 2, · · · , N − 1, u0 = 1, uN = uN −1 .
h2
Para analizar el problema discreto, se expresa en forma matricial Ax = b
donde
 
2 −1 0 · · · 0 0  
 −1 2 −1 · · · 0  1
 0   0 
 0 −1 2 · · · 0 0   
   .. 
A =  .. .. .. . . .. ..  , b =  . .
 . . . . . .   
   0 
 0 0 0 · · · 2 −1 
0
0 0 0 · · · −1 1
Si se desarrolla el cálculo del determinante de A por la última fila, se com-
prueba que está dado por la fórmula
det(A) = det(AN −2 ) − det(AN −3 ),
donde Aj representa la matriz formada por las j primeras filas y columnas
de la matriz A para j = 1, 2, · · · , N − 2. De acuerdo con el ejercicio ??, se
tiene que
det(A) = N − 1 − (N − 2) = 1 > 0,
lo que prueba que el sistema tiene una única solución. ⋄