Numérico I: Aproximación numérica

Entregable 3

1. En un experimento, la cantidad de calor despedida en relación con el tiempo, fue la siguiente: a

las 2 horas del comienzo, 3 kilocalorías1 (Kcal), a las 4 horas 1 Kcal y a las 6 horas 2 Kcal.

a. Calcular el polinomio que interpola estos datos, aplicando el método de Lagrange:

Los puntos de datos dados son: (2,3), (4,1), (6,2)

El polinomio interpolante de Lagrange está definido por la siguiente ecuación:

n
P ( x )=∑ y i∗Li ( x )
i=0

En donde n es el número de puntos menos 1, yi es la coordenada y del punto i-ésimo y Li(x) es el

término de Lagrange para el punto i-ésimo, dado por la ecuación:
n
x−x j
Li ( x ) = ∏ x i−x j
j=0 , j ≠i

Entonces para los datos del ejercicio, se procede con lo siguiente; x0= 2 horas, x1= 4 horas, x2= 6
horas

- Para i = 0,

(x−4 )(x−6) (x−4)(x−6)

L0 ( x ) = =
(2−4 )(2−6) 8

- Para i = 1
(x−2)(x−6) −(x −2)(x−6)
L1 ( x ) = =
(4−2)(4−6) 4
- Para i = 2
(x−2)(x−4) ( x −2)(x−4)
L2 ( x ) = =
(6−2)(6−4) 8

Entonces, el polinomio que interpola los datos dados P(x) es:

3∗( x−4 ) ( x−6 ) 1∗−( x−2 )( x−6 ) 2∗(x−2)( x−4)

P ( x )= + +
8 4 8
Al simplificar lo anterior, se obtiene el polinomio completo:

3 ( x−4 )( x−6 ) ( x−2 )( x−6 ) (x−2)(x−4)

P ( x )= − +
8 4 4
3 ( x−4 )( x−6 ) ( x−2 )( x−6 ) (x−2)(x−4)
P ( x )= − +
8 4 4
3 ( x−4 )( x−6 ) 2 ( x−2 ) ( x−6 ) 2(x−2)(x−4 )
P ( x )= − +
8 8 8
 3 ( x−4 )( x−6 )−2 ( x−2 ) ( x−6 ) +2(x−2)(x−4)
P ( x )=
8

3 ( x 2−10 x +24 )−2 ( x 2−8 x +12 ) +2(x 2−6 x +8)

P ( x )=
8
2 2 2
3 x −30 x +72−2 x +16 x−24+2 x −12 x +16
P ( x )=
8
2
x −2 X + 64
P ( x )=
8
1 2 1
P ( x )= x − x+ 8
8 4
b. Calcular el polinomio que interpola estos datos, aplicando el método de Newton:

Calcular las diferencias divididas:

f [ x 0 ]= y 0

f [ x 1 ]−f [ x 0 ]
f [ x 0 , x1 ] =
x 1−x 0

f [ x 1 , x 2 ]−f [ x 0 , x 1 ]
f [ x 0 , x1 , x2 ]=
x 2−x 0

Con esto se construye el polinomio interpolante con los datos: x0= (2,3), x1= (4,1), x2= (6,2)

f [ 2 ] =3 , f [ 4 ] =1 , f [ 6 ]=2
1−3
f [2 , 4 ]= = -1
4−2
2−1
f [ 4 , 6 ]= =0.5
6−4
0.5−(−1) 1.5
f [ 2 , 4 ,6 ] = = =0.375
6−2 4
Con lo anterior, el polinomio interpolante de Newton es:

P ( x )=f [ x 0 ] +f [ x 0 , x 1 ] ( x−x 0 ) + f [ x 0 , x 1 , x2 ] ( x−x 0 ) ( x−x 1 )

P ( x )=3−1 ( x−2 ) +0.375 (x−2)(x−4 )

c. Calcular una estimación de la cantidad de calor que se desprenderá a las 3 y a las 10 horas del
comienzo del experimento. ¿Qué resultado es “más fiable”? ¿Por qué?

Con el polinomio de Lagrange:

1 2 1
P ( x )= x − x+ 8
8 4
Para x = 3,

1 2 1
P ( 3 )= 3 − 3+8
8 4
9 3
P ( 3 )= − + 8
8 4
9−6+64 67
P ( 3 )= =
8 8
P ( 3 )=8.375

Para x = 10

1 2 1
P ( 10 )= 10 − 10+8
8 4
100 10
P ( 10 )= − +8
8 4
100−20+ 64 144
P ( 10 )= = =18
8 8
Con el polinomio de Newton:

Para x = 3,

P ( 3 )=3−1 (3−2 ) +0.375 (3−2)(3−4)

P ( 3 )=2 ( 1 ) +0.375(1)(−1)
P ( 3 )=2−0.375=1 , 25

Para x = 10

P ( 10 )=3−1 (10−2 ) +0.375(10−2)(10−4)

P ( 10 )=3−1 ( 8 ) +0.375(8)(6)
P ( 10 )=2 ( 8 ) +0.375 (8)(6)
P ( 10 )=16+ 18=31

Teniendo en cuenta los resultados, es posible decir que el polinomio de Lagrange parece ser más
confiable para el problema dado, ya que las estimaciones están mas cercanas a los valores
observados reales de la investigación. Es importante mencionar, que dichos polinomios solo
interpolan los datos por lo que no son muy precisos para estimar datos mas allá de los dados en el
problema.

2. Empíricamente se han obtenido los siguientes valores de una función y = f(x) (desconocida):

xi 0 2 4 6
yi 0 1 4 2
 a. Hallar el polinomio de interpolación que aproxima estos datos por el método que
consideres:

Método de interpolación de Lagrange, siguiendo el procedimiento del primer punto

- Para i = 0
(x−2)(x−4 )(x−6) x 3−12 x 2+ 44 x−48
L0 ( x ) = =
(0−2)(0−4)( 0−6) −48

- Para i = 1
(x−0)(x−4)(x−6) x 3−10 x 2 +24 x
L1 ( x ) = =
(2−0)(2−4)(2−6) 16
- Para i = 2
(x−0)(x−2)( x−6) x 3−8 x 2 +12 x
L2 ( x ) = =
(4−0)(4−2)( 4−6) −16
- Para i = 3

(x−0)(x−2)( x−4) x 3−6 x 2 +8 x

L3 ( x ) = =
(6−0)(6−2)(6−4) 48

Entonces, el polinomio de interpolación sería:

3 2 3 2 3 2 3 2
0∗x −12 x + 44 x−48 1∗x −10 x +24 x 4∗x −8 x + 12 x 2∗x −6 x +8 x
P ( x )= + + +
−48 16 −16 48

x −10 x +24 x 4(x −8 x +12 x ) 2(x ¿ ¿ 3−6 x +8 x )

3 2 3 2 2
P ( x )= + + ¿
16 −16 48
3 2
x −10 x +24 x
−4 x −32 x +48 x ¿ ¿ + x −6 x +8 x ¿ ¿
3 2 3 2
( )
P x=
16 16 24
3 2 3 2 3 2
x −10 x +24 x−4 x + 32 x −48 x + x −6 x +8 x
P ( x )=
16
3 2
−2 x +16 x −16 x
P ( x )=
16
−1 3 2
P ( x )= x + x −x
8
b. Si a la tabla anterior se añade una nueva observación (x s; ys) = (8,1), obtener el nuevo
polinomio de interpolación correspondiente a todos los datos.

Con la nueva posición el procedimiento sería:

- Para i = 0
(x−2)(x−4 )(x−6)(x−8) x 4−20 x 3 +140 x2 −192 x +384
L0 ( x ) = =
(0−2)(0−4 )(0−6)(0−8) 384
 - Para i = 1
(x−0)(x−4)(x−6)(x−8) x 4−18 x 3 +104 x 2−192 x
L1 ( x ) = =
(2−0)(2−4)(2−6)(2−8) −96
- Para i = 2
(x−0)(x −2)( x−6)(x−8) x 4−16 x 3 +76 x 2−96 x
L2 ( x ) = =
(4−0)(4−2)( 4−6)(4−8) 64
- Para i = 3

(x−0)(x−2)( x−4)(x−8) x 4−14 x 3 +56 x 2−64 x

L3 ( x ) = =
(6−0)(6−2)(6−4)(6−8) −96

- Para i = 4

(x−0)(x −2)( x −4)(x−6) x 4 −12 x 3 +44 x 2−48 x

L4 ( x ) = =
(8−0)(8−2)( 8−4)(8−6) 384

Entonces, el polinomio de interpolación sería:

4 3 2 4 3 2 4 3 2 4 3
0∗x −20 x +140 x −192 x+ 384 1∗x −18 x +104 x −192 x 4∗x −16 x +76 x −96 x 2∗x −14 x +
P ( x )= + + +
384 −96 64 −9

(x ¿ ¿ 4−16 x +76 x −96 x ) 2(x ¿ ¿ 4−14 x +56 x −64 x ) x 4−12 x 3+ 44 x 2−

4 3 2 3 2 3 2
x −18 x +104 x −192 x
P ( x )= +4 + +
−96 64 −96 384
4 3 2 4 3 2 4 3 2 4 3 2
−x +18 x −104 x +192 x x −16 x +76 x −96 x −x +14 x −56 x +64 x x −12 x + 44 x −48 x
P ( x )= + + +
96 16 48 384

P ( x )=−x 4 +18 x 3−104 x 2 +192 x + x 4−16 x 3 +76 x 2−96 x −x 4 +14 x 3−56 x 2+ 64 x ¿+ x 4−12 x 3+ 44 x 2−48 x ¿
384
3 2 4 3 2
16 x −84 x +160 x + x −12 x +44 x −48 x
P ( x )=
384
4 3 2
x + 4 x −40 x +112 x
P ( x )=
384
c. Hallar el valor que, aproximadamente, debería corresponder a x=7 mediante ambos
polinomios. ¿Cuál de los valores obtenidos crees que se debería tomar como mejor
aproximación de f(7)? ¿Por qué?

Polinomio 1:
−1 3 2 −1 3 2 −343 −7
P ( x )= x + x −x = 7 +7 −7= +49−7=
8 8 8 8
Polinomio 2:
4 3 2
7 + 4(7) −40 ( 7 ) +112 (7) 2401+1372−1960+ 784 2417
P ( x )= = = =6.294
384 384 384
De lo anterior, el valor del polinomio dos sería el mejor para tomar como aproximación, ya que, al
tener más datos en la tabla, más preciso será el polinomio.

3. Se han obtenido los siguientes tiempos (medidos en segundos) y posiciones (medidas en

metros) para predecir la velocidad de un tren de alta velocidad para cada instante incluido en la
tabla:

Tiempo 0 3 5 8 10 13
Distancia 0 225 38 623 74 993
3 2

a. Halla la velocidad haciendo uso de la fórmula del punto medio (revisad bien los datos para
hacer uso de los casos particulares):

Para trabajar con la fórmula del punto medio, se emplea la siguiente fórmula

t 1 +t 2
V promedio =f ( )
2
En donde V promedio es la velocidad promedio en un intervalo de tiempo específico, f(t) es la función
que describe la velocidad del objeto; y, t 1 y t 2 son los tiempos inicial y final del intervalo de tiempo.

Primero para hacer uso de la fórmula anterior, primero se debe calcular los puntos medio del
intervalo de tiempo dado y así luego poder calcular la velocidad:

1. Intervalo de tiempo t=0 y t=3

0+ 3
t medio = =1.5 s
2
El intervalo para calcular la velocidad de este intervalo de tiempo se encuentra entonces entre 0 m
y 225 m.

2. Intervalo de tiempo t=3 y t=5

3+5
t medio = =4 s
2
El intervalo para calcular la velocidad de este intervalo de tiempo se encuentra entonces entre 225
m y 383 m.

3. Intervalo de tiempo t=5 y t=8

5+ 8
t medio = =6.5 s
2
El intervalo para calcular la velocidad de este intervalo de tiempo se encuentra entonces entre 383
m y 623 m.
4. Intervalo de tiempo t=8 y t=10
8+ 10
t medio = =9 s
2
El intervalo para calcular la velocidad de este intervalo de tiempo se encuentra entonces entre 623
m y 742 m.

5. Intervalo de tiempo t=10 y t=13

10+13
t medio = =11.5 s
2
El intervalo para calcular la velocidad de este intervalo de tiempo se encuentra entonces entre 742
m y 993 m.

Ahora que ya se tienen los puntos medios en cada intervalo de tiempo, se procede a calcular la
velocidad promedio. Empleando la siguiente fórmula:

∆ Distancia
V promedio =
∆ Tiempo

1. Intervalo de tiempo t=0 y t=3

225 m−0 m 225 m
V promedio = = =75 m/s
3 s−0 s 3s

2. Intervalo de tiempo t=3 y t=5

383 m−225 m 158 m
V promedio = = =79 m/s
5 s−3 s 2s

3. Intervalo de tiempo t=5 y t=8

623 m−383 m 240 m
V promedio = = =80 m/s
8 s−5 s 3s

4. Intervalo de tiempo t=8 y t=10

742m−623 m 119 m
V promedio = = =59.5 m/s
10 s−8 s 2s

5. Intervalo de tiempo t=10 y t=13

993 m−742 m 251 m
V promedio = = =83.67 m/s
13 s−10 s 3s

Lo anterior son los valores aproximados de la velocidad promedio en cada intervalo de tiempo,
empleando la fórmula del punto medio.
b. ¿Se podría aplicar la derivada de cinco puntos vista en el tema 4?

Si es posible aplicar dicha derivada para estimar la velocidad del tren. Esto debido a que se tienen
los datos de distancia vs el tiempo, se puede calcular la velocidad en cada tiempo. Para esto
primero hay que convertir los puntos de tiempo y distancia en arreglos de datos. Luego, elegir los
pesos apropiados para cada punto dado. Y por último aplicar la formula de la derivada de cinco
puntos. El cálculo exacto y la aplicación de la derivada requiere conocer la fórmula específica y los
espacios entre los puntos de datos, por lo que se necesitaría más información sobre el fenómeno.

4. Aproxima las siguientes integrales con la Regla del punto medio, la Regla del trapecio y con la
Regla de Simpson:
3

a) ∫ x∗e dx
2x

Hay que dividir el intervalo en 4 subintervalos, ya que la regla de Simpson necesita un número par
de subintervalos para aplicarse correctamente. Los puntos de división: 0, 0.75, 1.5, 2.25, 3.

1. Regla del punto medio:

b n
x i−1 + x i
∫ f ( x ) dx ≈ h ∑ f ( 2
)
a i=l

Donde h es la longitud de los subintervalos y xi son los puntos de división. H =(b-a)/n =

(3-0)/4=0.75

≈ 0.75 [ f ( 0.375 )+ f ( 1.125 ) + f ( 1.875 )+ f ( 2.625 ) ]

≈ 0.75 [ ( 0.375 )(e 0.75)+ ( 1.125 ) (e2.25 )+ ( 1.875 ) (e 3.75)+ ( 2.625 ) (e5.25 ) ]

≈ 0.75 [ ( 0.375 ) (2.117)+ ( 1.125 ) (9.487)+ ( 1.875 ) (42.212)+ ( 2.625 ) (188.301) ]

≈ 0.75 [0.796+ 10.684+79.276+ 494.735]

≈ 0.75∗585.491
≈ 439.118
2. Regla del trapecio:
b

∫ f ( x ) dx ≈ h2 ¿ ¿
a

∫ x∗e2 x dx ≈ 0.75
2
¿
0

0.75
≈ ¿
2
0.75
≈ ¿
2
 0.75
≈ ¿
2
0.75
≈ [0+ 2 ( 52.82 ) +20.085 ]
2
0.75
≈ [105.64 +20.085]
2
0.75
≈ ∗125.725 ≈ 47.171
2
3. Regla de Simpson:
b

∫ f ( x ) dx ≈ 3h ¿ ¿
a
3

∫ x∗e2 x dx ≈ 0.75
3
¿
0
0.75
≈ ¿
3

0.75
≈ ¿
3
0.75
≈ [0+ 4 ( 0.796 +10.684+79.276 +494.735 )+2 (1.121+ 9.487+ 42.212)]
3
0.75
≈ [0+ 4 ( 585.491 ) +2 ( 52.82 ) ]
3
0.75
≈ [0+ 2341.964+105.64 ]
3
0.75
≈ ∗2447.604 ≈ 611.901
3

2
2 x+1
b) ∫ 2
dx
1 x +x

El cambio ∆ x y luego dividir el intervalo [1,2] en subintervalos de igual tamaño.

2−1
∆ x=
n
Se toma un valor de n=2, para mejor aproximación y simplicidad.

2−1 1
∆ x= = =0.5
2 2
 1. Regla del punto medio:
2 2

∫ 2x x+1
2
+x
dx ≈ ∆ x ∑ f (x i ¿ )=( 0.5 ) [f ( 1.25 ) +f (1.75)]¿
¿

1 i=1

2 (1.25 )+ 1 2 (1.75 )+ 1
≈ ( 0.5 ) [ 2
+ 2
]
1.25 +1 1. 75 +1
≈ ( 0.5 ) ¿]
≈ 1.2375
2. Regla del trapecio:
2

∫ 2x x+1
2
+x
dx ≈
∆x
2
¿
1

A y b son los límites inferior y super de la integral respectivamente y ∆ x es el ancho de cada

subintervalo.

∆x
≈ ¿
2
0.5
≈ [f ( 1 )+ 2 f ( 1.5 ) +f ( 2 ) ]
2
0.5 2 ( 1 )+ 1 2 ( 1.5 ) +1 2 ( 2 ) +1
≈ [ +2 + 2 ]
2 12 +1 2
1.5 +1 2 +1
0.5
≈ [ 1.5+1,231+1 ]
2
0.5
≈ ∗3.731 ≈ 0.933
2

3. Regla de Simpson:
2

∫ 2x x+1
2
+x
dx ≈
∆x
3
¿
1

A y b son los límites inferior y superior de la integral respectivamente, y ∆ x es el ancho de cada

subintervalo.

0.5
≈ [f ( 1 )+ 4 f ( 1.25 ) +2 f ( 1.5 ) +f ( 2 ) ]
3
0.5 2 ( 1 )+ 1 2 ( 1. 25 )+1 2 ( 1.5 ) +1 2 ( 2 ) +1
≈ [ 2 +4 +2 + 2 ]
3 1 +1 2
1.2 5 +1
2
1.5 +1 2 +1

≈
0.5 3
3 2 [+4 (1.757 ) +2 (1.231 )+1
]
 ≈
0.5 3
3 2 [+7.028+2.462+1
]
0.5
≈ ∗11.99≈ 1.998
3
1
1
c) ∫ 2
dx
0 x +1

Primero, calculemos los valores exactos de la integral:

|
1

∫ x 21+1 dx=arctan ( x ) ¿ arctan ( 1 )−arctan ( 0 )= π4

0

1−0
Con n= 10 para las aproximaciones, con esto calculamos el tamaño de h = =0.1
10
1. Regla del punto medio:
n
Aproximación= ≈ h ∑ f ( x i ¿ )¿
❑

i=1

b−a 1−0 1
h= = =
n n n
1
a+ h+ (i−1 ) h
2 1 i−1
x i= = +
n 2n n
1 n

∫ x 21+1 dx ≈ h ∑ f ( x ❑i ¿ )¿
0 i =1

10
1
≈ 0.1 ∑
i=1 1 i−1 2
( + ) +1
2n n
Al realizar los cálculos de igual manera en los dos puntos anterior en esta integral utilizando la
regla del punto medio ≈ 0.7854

2. Regla del trapecio:

Aproximación =
1

∫ x 21+1 dx ≈ h2 ¿
0

0.1
≈ ¿
2
Al realizar los cálculos de igual manera en los dos puntos anterior en esta integral utilizando la
regla del trapecio ≈ 0.7854

3. Regla de Simpson:
Aproximación =
1

∫ x 21+1 dx ≈ h3 ¿
0

0.1
≈ ¿
3
Al realizar los cálculos de igual manera en los dos puntos anterior en esta integral utilizando la
regla Simpson ≈ 0.7854

π
Al comparar las tres aproximaciones con el valor exacto calculado al inicio de , es notorio que
4
todas son cercanas al valor real, demostrando la eficacia de las 3 reglas de aproximación numérica.