Pauta: Evaluación de Recuperación: Cálculo Numérico 2022-1

1. La siguiente es la ecuación de van der Waals que generaliza la ecuación de la ley ideal de los gases.

n2 a
 
P + 2 (x − nb) = nRT, (1)
x

donde P es la presión (en atmósferas), x el volumen (en litros), T la temperatura (en Kelvin), n el
número de moles del gas y R la constante molar.

(a) [0.2 pts.] Se sabe que P = 15, T = 300 R = 0.08, a = 1, b = 0 y n = 1, pero se desconoce el
volumen x. Reemplace estos valores en (1) y escriba la ecuación resultante.
(b) [0.8 pts.] Partiendo con x0 = 2, calcule x1 utilizando el método de Newton-Raphson aplicado
a la ecuación obtenida en (a).

Solución:

(a) Reemplazando los valores en (1) se tiene que

 
1
15 + 2 x = 24,
x

de donde 15x2 − 24x + 1 = 0.

(b) Sea f (x) = 15x2 − 24x + 1. Entonces,f 0 (x) = 30x − 24 y por tanto, como x0 = 2,

f (x0 ) f (2) 13 59
x1 = x0 − 0
=2− 0 =2− = ≈ 1.6389.
f (x0 ) f (2) 36 36

2. Considere la siguiente tabla que muestra la cantidad de personas infectadas y(t) producto de una
pandemia en el dı́a t.

t (dı́a) 0 1 2
y (número de personas) 1 2 16

Se propone el siguiente modelo: y(t) = αeβt .

(a) [0.2 pts.] Linealice el modelo.

(b) [0.6 pts.] Determine los parámetros α y β que mejor ajustan los datos de la tabla en el sentido
de los mı́nimos cuadrados.
(c) [0.2 pts.] ¿Cuántas personas contagiadas habrá el tercer dı́a?.

Solución:

(a) Tenemos que ln(y(t)) = ln(α) + βt = a0 + βt con a0 = ln(α).

(b) Queremos determinar α, β ∈ R de modo que
   
1 t0   ln(y(t0 ))
1 t1  a0 = ln(y(t1 )) .
β
1 t2 ln(y(t2 ))

Es decir, debemos resolver    

1 0   ln(1)
1 1 a0 =  ln(2) 
β
1 2 ln(16)
 en el sentido de los mı́nimos cuadrados. Para ello resolvemos el sistema de ecuaciones normales
   
  1 0     ln(1)
1 1 1  a 1 1 1 
1 1 0 = ln(2) 
0 1 2 β 0 1 2
1 2 ln(16)
      
3 3 a0 ln(2) + ln(16) 5 ln(2)
= = .
3 5 β ln(2) + 2 ln(16) 9 ln(2)
1
De donde β = 2 ln(2) y a0 = − ln(2) = ln(2−1/3 ). Por tanto α = ea0 = 2−1/3 .
3
Ası́, y(t) = 2−1/3 e2 ln(2)t .
(c) y(3) = 2−1/3 e2 ln(2)3 = 2−1/3 26 ≈ 50.80.

3. [1.2 pts.] Considere los datos de siguiente tabla. Utilice los polinomios de Lagrange para determinar

x 0 π/2 π
y 0 3 0

el polinomio que interpola a los datos de la tabla.

Solución: Tenemos que p(x) = y0 l0 (x) + y1 l1 (x) + y2 l2 (x) = 3l1 (x).
x(x − π) 4 12
Además, l1 (x) = = − 2 x(x − π) y por tanto p(x) = − 2 x(x − π).
π/2(π/2 − π) π π
R1
4. Llamemos I(f ) = −1 f (x)dx y sean

Q(f ) = f (s) + f (−s), (2)

reglas de cuadratura para aproximar I(f ) en las que s ∈ [0, 1].

a) [0.6 pts.] Determine para qué valor de s se cumple que esta regla es exacta para polinomios
del mayor grado posible.
b) [0.6 pts.] ¿Cuál es la regla de cuadratura que resulta de escribir
Z 2 Z −1 Z 1 Z 2
f (x)dx = f (x)dx + f (x)dx + f (x)dx
−2 −2 −1 1

y aproximar las integrales en la parte derecha de la igualdad anterior por la regla (2) con el
valor de s determinado en item anterior.
Observación: Tenga en cuenta
R −1 R 2 que (2) no puede aplicarse directamente al cálculo de aproxi-
maciones de −2 f (x)dx y de 1 f (x)dx.

Solución:

a) Como
Z 1
dx = 2, Q(1) = 1 + 1 = 2,
−1

ocurre que para todo s ∈]0, 1] la regla es exacta para polinomios de grado menor o igual que
0.
Dado que Z 1
xdx = 0, Q(x) = s − s = 0,
−1

ocurre que para todo s ∈]0, 1] la regla es exacta para polinomios de grado menor o igual que
1.
 Z 1
2
x2 dx = , Q(x2 ) = 2s2 ,
−1 3
ocurre que la regla es exacta para√polinomios de grado menor o igual que 2 si y solo si s ∈]0, 1]
es tal que 2s2 = 32 , es decir, s = 33 .
Con ese valor para s se cumple además que
Z 1
x3 dx = 0, Q(x3 ) = 0,
−1

por tanto, la regla es exacta para polinomios de grado menor o igual que 3.
Por último,
Z 1 √ !4
2 3 1
x4 dx = , Q(x4 ) = 2 = ,
−1 5 3 9
solo podemos garantizar que la regla es exacta para polinomios de grado menor o igual que 3.
b) Como
b 1  
b−a a+b b−a
Z Z
f (x)dx = f + t dt,
a 2 −1 2 2
entonces
Z 2 Z −1 Z 1 Z 2
f (x)dx = f (x)dx + f (x)dx + f (x)dx,
−2 −2 −1 1
1 1
Z 1
1 1
Z   Z  
3 1 3 1
= f − + t dt + f (x)dx + f + t dt,
2 −1 2 2 −1 2 −1 2 2
entonces, aplicando la regla de cuadratura obtenida antes llegamos a
Z 2 √ ! √ !!
1 3 1 3 3 1 3
f (x)dx ≈ f − + f − −
−2 2 2 2 3 2 2 3
√ ! √ !
3 3
+f f −
3 3
√ ! √ !!
1 3 1 3 3 1 3
+ f + f −
2 2 2 3 2 2 3

5. [1.2 pts.] Considere el siguiente problema de valores iniciales

x y 0 (x) = 2x − y(x), x ∈ [1, 2], y(1) = 0.

Calcule una aproximación a y(2) aplicando Euler implı́cito con h = 1.

Solución:
Primero debemos re-escribir el problema en la forma y prime (x) = f (x, y(x)), lo cual es sencillo pues
solo debemos dividir por x la primera ecuación, obteniendo el siguiente PVI
y(x)
y 0 (x) = 2 − , x ∈ [1, 2], y(1) = 0.
x
[0.3 pts.]
Para obtener una aproximación a y(2) con Euler implı́cito y h = 1 solo debemos realizar un paso
del método:

y1 = y0 + hf (x1 , y1 ) = 0 + f (2, y1 ).
 [0.3 pts.]
Por tanto,
y1 4
y1 = 2 − ⇒ y1 = .
2 3
[0.3 pts.]
La aproximación a y(2) que se obtiene con Euler implı́cito y h = 1 es 43 .
[0.3 pts.]

6. [1.2 pts.] Determine las matrices P , L y U que forman la descomposición LU con pivoteo parcial
de la matriz  
1 0 1 0
0 1 0 1
 .
1 0 2 0
0 1 0 2

Solución:
Como el elemento con mayor valor absoluto en la columna 1 está en la fila 1, no es necesario realizar
intercambios de fila. Realizando la operación elemental f3 ← f3 − f1 obtenemos
     
1 0 0 0 1 0 0 0 1 0 1 0
0 1 0 0 0 1 0 0 0 1 0 1
P = 0 0 1 0 ,
 L= 1
, A1 = 
0 0 1 0 .

1 0 
0 0 0 1 0 1 0 1 0 2

La matriz A1 no es aún una matriz triangular superior. Debemos seguir aplicando el algoritmo de
descomposición LU con pivoteo parcial.
Nuevamente el elemento con mayor valor absoluto en la columna 2 (filas 2, 3 o 4) está en la
diagonal principal, no es necesario realizar intercambios de filas. Aplicando la operación elemental
f4 ← f4 − f2 obtenemos
     
1 0 0 0 1 0 0 0 1 0 1 0
0 1 0 0 0 1 0 0 0 1 0 1
P = 0 0 1 0
, L = 
1 0 1 0
, A 2 = 0 0 1 0 .
 

0 0 0 1 0 1 1 0 0 0 1

A2 es una matriz triangular superior. Entonces,

   
1 0 0 0 1 0 0 0
0 1 0 0 0 1 0 0
U = A2 , P = 0 0 1
, L= .
0 1 0 1 0
0 0 0 1 0 1 0 1