Taller Eje 4 - Integracion y Dif Numerica 2
Taller Eje 4 - Integracion y Dif Numerica 2
Taller Eje 4 - Integracion y Dif Numerica 2
DIFERENCIACIÓN E INTEGRACIÓN
NUMÉRICA
En las aplicaciones prácticas, a menudo sucede que tenemos valores de la variación de una
función en algunos puntos sin tener su expresión analítica, y queremos inferir de ellos la
función derivada en algún otro punto. De la misma forma, podemos necesitar la integral
definida de esa función de la que sólo sabemos sus valores en algunos puntos; en el caso de
la integral, además, puede darse el caso de que, incluso teniendo la expresión analítica de
casos la solución pasa por aproximar la derivada y la integral usando los métodos del
cálculo numérico.
para el problema formulado como una integral definida. Hay varias razones para llevar a
integración de forma analítica. Los métodos que se abordan son: método del trapecio,
a) h=0.1
𝑓(𝑥+ℎ)−𝑓(𝑥)
𝑓′ (𝑥 ) =
ℎ
𝑙𝑛 (1 + 0.1) − ln(1)
𝑓′ (1) =
0.1
𝑓′ (1) = 0.95310
b) h=0.01
𝑓(𝑥+ℎ)−𝑓(𝑥)
𝑓′ (𝑥) =
ℎ
c) h=0.001
𝑓 (𝑥 + ℎ) − 𝑓(𝑥)
𝑓′ (𝑥 ) =
ℎ
𝑙𝑛 (1 + 0.001) − 𝑙𝑛(1)
𝑓′ (1) =
0.001
𝑓′ (1) = 0.99950
2. Regla de la Diferencia centrada de tres puntos para la segunda derivada.
𝜋
Aproximar f´’( ) donde f(x) = cos (x)
4
a) h=0.1
𝑓 (𝑥 + ℎ) − 2𝑓 (𝑥 ) + 𝑓(𝑥 − ℎ)
𝑓′′ (𝑥 ) =
ℎ2
𝜋 𝜋 𝜋
𝜋 𝑐𝑜𝑠 ( + 0.1) − 2 cos ( ) + 𝑐𝑜𝑠( − 0.1)
𝑓′′ ( ) = 4 4 4
4 0.12
𝜋
𝑓 ′′ ( ) = −0.70652
4
b) h=0.01
𝑓 (𝑥 + ℎ) − 2𝑓 (𝑥 ) + 𝑓(𝑥 − ℎ)
𝑓′′ (𝑥 ) =
ℎ2
𝜋 𝜋 𝜋
𝜋 𝑐𝑜𝑠 ( + 0.01) − 2 cos ( ) + 𝑐𝑜𝑠( − 0.01)
𝑓′′ ( ) = 4 4 4
4 0.012
𝜋
𝑓 ′′ ( ) = −0.70710
4
c) h=0.001
𝑓 (𝑥 + ℎ) − 2𝑓 (𝑥 ) + 𝑓(𝑥 − ℎ)
𝑓′′ (𝑥 ) =
ℎ2
𝜋 𝜋 𝜋
𝜋 𝑐𝑜𝑠 ( + 0.001) − 2 cos ( ) + 𝑐𝑜𝑠( − 0.001)
𝑓′′ ( ) = 4 4 4
4 0.0012
𝜋
𝑓 ′′ ( ) = −0.70711
4
3. Regla del Trapecio
0.8
∫ 0.2 + 25𝑥 − 200𝑥 2 + 675𝑥 3 − 900𝑥 4 + 400𝑥 5 𝑑𝑥
0
𝑏
𝑓(𝑎) + 𝑓(𝑏)
∫ 𝑓(𝑥)𝑑𝑥 = (𝑏 − 𝑎) ( )
𝑎 2
𝑏
0.2 + 0.232
∫ 𝑓(𝑥)𝑑𝑥 = (0.8 − 0) ( ) = 0,1728
𝑎 2
Por lo que:
0.8
∫ 0.2 + 25𝑥 − 200𝑥 2 + 675𝑥 3 − 900𝑥 4 + 400𝑥 5 𝑑𝑥 ≅ 0.1728
0
4. Regla del Trapecio Compuesta
Evaluar con n= 2 y n= 4
𝜋/2
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥
0
a) n= 2
x0 x1 x2
0 𝜋 𝜋
4 2
f (x0) f (x1) f (x2)
9 8.12132 6
𝑏 𝑛−1
ℎ
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑓 (𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏))
𝑎 2
𝑖=1
𝜋
2 𝜋
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ (9 + 2(8.12132) + 6)
0 8
𝜋
2
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ 12.26896
0
b) n= 4
x0 x1 x2 x3 x4
0 𝜋 𝜋 3𝜋 𝜋
8 4 8 2
f (x0) f (x1) f (x2) f (x3) f (x4)
9 8.77164 8.12132 7.14805 6
𝑏 𝑛−1
ℎ
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑓 (𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏))
𝑎 2
𝑖=1
𝜋
2 𝜋
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ (9 + 2(8.77164 + 8.12132 + 7.14805) + 6)
0 16
𝜋
2
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ 12.38613
0
5. Regla de Simpson 1/3 con n=4
Evalúe:
3
∫ 𝑥 2 𝑒 𝑥 𝑑𝑥
0
𝑛=1 𝑛=2
𝑏 𝑓(𝑥0 ) + 4 ∑𝑖=1,3,5 𝑓(𝑥𝑖 ) + 2 ∑𝑖=2,4 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑏 − 𝑎)
𝑎 3𝑛
𝑎=0
𝑏=3
𝑛=4
𝑏−𝑎 3−0 3
ℎ= = = = 0.75
𝑛 4 4
𝑥0 = 𝑎 = 0
𝑓(𝑥0 ) = 0
𝑥1 = 𝑥0 + ℎ = 0.75
𝑓(𝑥1 ) = 𝑓 (0.75) = 1.19081
𝑏
0 + 4(1.19081 + 48.03166) + 2(10.08380) + 180.76983
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (3)
𝑎 3∗4
3
∫ 𝑥 2 𝑒 𝑥 𝑑𝑥 ≅ 99.4568275
0
6. Regla de Simpson 1/3 y 3/8 con n=5
Evalúe:
3
∫ 5 + 3𝑐𝑜𝑠𝑥 𝑑𝑥
0
𝑎=0
𝑏=3
𝑛=4 𝑠𝑒 𝑜𝑐𝑢𝑝𝑎𝑟á 4 𝑦𝑎 𝑞𝑢𝑒 𝑛 𝑑𝑒𝑏𝑒 𝑠𝑒𝑟 𝑝𝑎𝑟
𝑏−𝑎 3−0 3
ℎ= = = = 0.75
𝑛 4 4
𝑥0 = 𝑎 = 0
𝑓(𝑥0 ) = 5 + 3 cos(0) = 8
𝑥1 = 𝑥0 + ℎ = 0.75
𝑓(𝑥1 ) = 5 + 3𝑐𝑜𝑠(0.75) = 7.19507
𝑏
8 + 4(7.19507 + 3.11548) + 2(5.21221) + 2.03002
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (3)
𝑎 3∗4
3
∫ 5 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ 15.42416
0
𝑥1 = 𝑥0 + ℎ = 1
𝑓(𝑥1 ) = 5 + 3𝑐𝑜𝑠(1) = 6.62091
𝑥2 = 𝑥1 + ℎ = 1 + 1 = 2
𝑓(𝑥2 ) = 5 + 3𝑐𝑜𝑠(2) = 3.75156
𝑥3 = 𝑥2 + ℎ = 2 + 1 = 3
𝑓(𝑥3 ) = 5 + 3𝑐𝑜𝑠(3) = 2.03002
La integración con métodos numéricos es una herramienta útil cuando se trata de integrar
El método de Derivación centrada es mucho más preciso que el método de diferencias hacia
atrás y el método de diferencias hacia adelante, pero por otro lado tiene más funciones de
evaluar por lo que su algoritmo es más tedioso y echo analíticamente más propenso a fallar
regla de Simpson de 1/3 y la de 3/8 se pueden aplicar juntas sobre una misma curva para
obtener exactitudes de tercer orden sobre todo un intervalo. Además, la regla de Simpson de
1/3 y 3/8 es más exacta que otras formas numéricas de aproximar la integral como la regla
del trapecio.
notar que algunas tienen un grado de dificultad mayor que otras, pero la utilización de ellas
depende del tipo de diferencial a resolver, es decir del tipo de problema a enfrentar, debemos
saber identificar cuál de ellas utilizar para su correcta resolución. Lo mismo sucede son los
https://www.ugr.es/~arobles/Telecos/Tema10.pdf
https://www.academia.edu/32875352/Diferenciaci%C3%B3n_Num%C3%A9rica
https://vdocuments.mx/metodos-de-diferenciacion-e-integracion-numerica.html