ING. EDGAR H.

MEDINA TAPIA METODOS NUMERICOS I

Tema 4.- Integración y Derivación Numéricas

4.1.5.- La Integral Doble. - Dada la integral

� 𝑓𝑓 (𝑥𝑥, 𝑦𝑦)𝑑𝑑𝑑𝑑 (4.1.5.1)

𝑅𝑅
Donde 𝑅𝑅 es la región de Integración

Para aproximar numéricamente una integral doble es necesario considerarla como una combinación de dos
integrales simples y a continuación aplicar las fórmulas de cuadratura numérica para ambas integrales.

(a) Integral Doble sobre una Región No Rectangular

Dependiendo del tipo de la región 𝑅𝑅, si es YS (Y - simple), XS (X - simple) la integral doble puede ser expresada como

𝑏𝑏 𝑑𝑑(𝑥𝑥) 𝑑𝑑 𝑏𝑏(𝑦𝑦)
𝑌𝑌𝑌𝑌: � �� 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) 𝑑𝑑𝑑𝑑� 𝑑𝑑𝑑𝑑 , 𝑋𝑋𝑋𝑋: � �� 𝑓𝑓(𝑥𝑥, 𝑦𝑦) 𝑑𝑑𝑑𝑑� 𝑑𝑑𝑑𝑑 (4.1.5.2)
𝑎𝑎 𝑐𝑐(𝑥𝑥) 𝑐𝑐 𝑎𝑎(𝑦𝑦)

Región 𝑅𝑅 de Integración Y- simple Región 𝑅𝑅 de Integración X- simple

Si se aplica el método de Simpson en ambas direcciones con 𝑀𝑀 intervalos en la dirección de las 𝑥𝑥, 𝑁𝑁 intervalos en
la dirección de las 𝑦𝑦; para una región YS, se realizarán los siguientes pasos:

1º Se aplicará para la integral externa las fórmulas de Simpson con 𝑀𝑀intervalos, esto es

𝑏𝑏 𝑑𝑑(𝑥𝑥𝑖𝑖) 𝑚𝑚 𝑑𝑑(𝑥𝑥𝑖𝑖) 𝑚𝑚
𝑏𝑏 − 𝑎𝑎 𝑏𝑏 − 𝑎𝑎
� � 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 = � � �[𝜔𝜔𝑖𝑖 ] � 𝑓𝑓 (𝑥𝑥𝑖𝑖 , 𝑦𝑦) 𝑑𝑑𝑑𝑑 + 𝐸𝐸 = � � �[𝜔𝜔𝑖𝑖 ]𝑆𝑆(𝑥𝑥𝑖𝑖 ) (4.1.5.3)
𝑎𝑎 𝑐𝑐(𝑥𝑥𝑖𝑖) 3𝑀𝑀 𝑐𝑐(𝑥𝑥𝑖𝑖) 3𝑀𝑀
𝑖𝑖=1 𝑖𝑖=1

2º Luego, se empleará las fórmulas de Simpson con n de puntos a cada una de las integrales internas 𝑆𝑆(𝑥𝑥𝑖𝑖 )

𝑑𝑑(𝑥𝑥𝑖𝑖) 𝑛𝑛
𝑑𝑑 (𝑥𝑥𝑖𝑖 ) − 𝑐𝑐(𝑥𝑥𝑖𝑖 )
𝑆𝑆 (𝑥𝑥𝑖𝑖 ) = � 𝑓𝑓(𝑥𝑥𝑖𝑖 , 𝑦𝑦) 𝑑𝑑𝑑𝑑 = � � ��𝜔𝜔𝑗𝑗 �𝑓𝑓�𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑗𝑗 � (4.1.5.4)
𝑐𝑐(𝑥𝑥𝑖𝑖) 3𝑁𝑁
𝑗𝑗=1

3º La aproximación de la integral doble se obtendrá sustituyendo los resultados obtenidos para las integrales
internas 𝑆𝑆(𝑥𝑥𝑖𝑖 ) en la fórmula de la integral externa.

41
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

Nota. - El número de puntos en la dirección del eje X es igual a 𝑚𝑚 = 𝑀𝑀 + 1 mientras que, el número de puntos en
la dirección eje Y está dado por 𝑛𝑛 = 𝑁𝑁 + 1

Ejemplo

Empleando las Fórmulas de Newton Cotes y la Cuadratura de Gauss – Legendre, aproxime la Integral Doble

𝑑𝑑𝑑𝑑
𝐼𝐼2 = �
𝑅𝑅
�8 − 𝑥𝑥 2 − 𝑦𝑦 2

Donde 𝑅𝑅 es la región ubicada en el primer cuadrante,

limitada por la curva 𝑥𝑥 2 + 𝑦𝑦 2 = 4 y los ejes coordenados

Considerando 𝑅𝑅 como una Región 𝑋𝑋 − 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆, se tiene

2 �4−𝑦𝑦 2
𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝐼𝐼2 = � =� �
𝑅𝑅
�8 − 𝑥𝑥 2 − 𝑦𝑦 2 0 0 �8 − 𝑥𝑥 2 − 𝑦𝑦 2

Resolución.

(1) Aproximación mediante las Fórmulas de Cuadratura de Newton-Cotes

Se empleará la Regla Compuesta de Simpson con 3 puntos (2 intervalos) en ambas direcciones

Obtención de la Fórmula de Integración

𝑑𝑑 − 𝑐𝑐 2 − 0 𝑗𝑗 0 1 2
𝑘𝑘 = = =1
𝑁𝑁 2 𝑦𝑦𝑗𝑗 0 1 2
𝜔𝜔𝑗𝑗 1 4 1

2
�4−𝑦𝑦0 �4−𝑦𝑦1 2 �4−𝑦𝑦2 2
𝑘𝑘 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝐼𝐼2 = �1 � +4� + 1� � (1)
3 0 �8 − 𝑥𝑥 2 − 𝑦𝑦02 0 �8 − 𝑥𝑥 2 − 𝑦𝑦12 0 �8 − 𝑥𝑥 2 − 𝑦𝑦22

1
𝐼𝐼2 = [1𝑆𝑆0 + 4𝑆𝑆1 + 1𝑆𝑆2 ] (2)
3

Cálculo de las Integrales Internas 𝑆𝑆𝑖𝑖 , 𝑖𝑖 = 1,2,3

𝑆𝑆0 𝑦𝑦0 = 0 �4−𝑦𝑦02 2

𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑆𝑆0 = � =�
0 �8 − 𝑥𝑥𝑖𝑖 2 − 𝑦𝑦02 0 �8 − 𝑥𝑥𝑖𝑖 2

𝑏𝑏 − 𝑎𝑎 2 − 0 1
ℎ= = =1 ⇒ 𝑥𝑥𝑖𝑖 = 0 , 1, 2 ⇒ 𝑓𝑓 (𝑥𝑥𝑖𝑖 ) =
𝑀𝑀 2 �8 − 𝑥𝑥𝑖𝑖 2

42
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

ℎ
𝑆𝑆0 = [1𝑓𝑓(𝑥𝑥0 ) + 4𝑓𝑓 (𝑥𝑥1 ) + 1𝑓𝑓 (𝑥𝑥2 )]
3

1
𝑓𝑓 (𝑥𝑥0 ) = 𝑓𝑓 (0) = = 0.353 553 390
√8 − 02

1
𝑓𝑓 (𝑥𝑥1 ) = 𝑓𝑓(1) = = 0.377 964 473
√8 − 12

1
𝑓𝑓 (𝑥𝑥2 ) = 𝑓𝑓 (2) = = 0.500 000 000
√8 − 22

1
𝑆𝑆0 = [1(0.353 553 390) + 4(0.377 964 473) + 1(0.5)] = 0.788 470 427 (𝑎𝑎)
3

𝑆𝑆1 𝑦𝑦1 = 1 �4−𝑦𝑦12 √3

𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑆𝑆1 = � =�
0 �8 − 𝑥𝑥𝑖𝑖 − 2 𝑦𝑦12 0 �7 − 𝑥𝑥𝑖𝑖 2

𝑏𝑏 − 𝑎𝑎 √3 − 0 √3 √3 1
ℎ= = = ⇒ 𝑥𝑥𝑖𝑖 = 0 , , √3 ⇒ 𝑓𝑓 (𝑥𝑥𝑖𝑖 ) =
𝑀𝑀 2 2 2 �7 − 𝑥𝑥𝑖𝑖 2

ℎ
𝑆𝑆0 = [1𝑓𝑓(𝑥𝑥0 ) + 4𝑓𝑓 (𝑥𝑥1 ) + 1𝑓𝑓 (𝑥𝑥2 )]
3

1
𝑓𝑓 (𝑥𝑥0 ) = 𝑓𝑓 (0) = = 0.377 964 473
√7 − 02

1
𝑓𝑓 (𝑥𝑥1 ) = 𝑓𝑓 (1) = 2
= 0.400 000
�7 − �√3
2�

1
𝑓𝑓 (𝑥𝑥2 ) = 𝑓𝑓 (2) = = 0.500 000
2
�7 − �√3�

1
𝑆𝑆1 = [1(0.377 964 473) + 4(0.4) + 1(0.5)] = 0.715 326 727 (𝑏𝑏)
3

𝑆𝑆2 𝑦𝑦2 = 2 �4−𝑦𝑦22 0

𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑆𝑆2 = � =� =0
0 �8 − 𝑥𝑥𝑖𝑖 − 2 𝑦𝑦22 0 �4 − 𝑥𝑥𝑖𝑖 2

𝑆𝑆2 = 0 (𝑐𝑐 )

Obtención de la Integral Doble. - Sustituyendo (𝑎𝑎), (𝑏𝑏), (𝑐𝑐 ) en (2), resulta,

1
𝐼𝐼2 = [1(0.788 470 427) + 4(0.715 326 727) + 1(0)] = 1.216 592 445 ⇽ 𝑆𝑆𝑆𝑆𝑆𝑆.
3

43
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

(2) Aproximación mediante las Fórmulas de Cuadratura de Gauss - Legendre

Se empleará G.L. Con 2 puntos en ambas direcciones

𝜆𝜆1 = 1 𝜆𝜆2 = 1
1 1
𝑡𝑡1 = − 𝑡𝑡2 =
√3 √3

Obtención de la Fórmula de Integración

𝑑𝑑 − 𝑐𝑐 𝑑𝑑 + 𝑐𝑐 2−0 2+0
𝑦𝑦𝑗𝑗 = � � 𝑡𝑡𝑗𝑗 + � �=� � 𝑡𝑡𝑗𝑗 + � � = 𝑡𝑡𝑗𝑗 + 1 , 𝑗𝑗 = 1,2
2 2 2 2

�4−𝑦𝑦12 �4−𝑦𝑦2 2
𝑑𝑑 − 𝑐𝑐 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝐼𝐼2 = � � �𝜆𝜆1 � + 𝜆𝜆2 � � = 1[1𝑆𝑆1 + 1𝑆𝑆2 ] (1)
2 0
2
�8 − 𝑥𝑥 − 𝑦𝑦1
2 0 �8 − 𝑥𝑥 2 − 𝑦𝑦22

𝑦𝑦𝑗𝑗 = 𝑡𝑡𝑗𝑗 + 1

1
𝑦𝑦1 = 𝑡𝑡1 + 1 = − + 1 = 0.422 649 730
√3
1
𝑦𝑦2 = 𝑡𝑡2 + 1 = + 1 = 1.577 350 269
√3

Cálculo de las Integrales Internas 𝑆𝑆𝑖𝑖 , 𝑖𝑖 = 1, 2

𝑆𝑆1 𝑦𝑦1 = 0.422 649 730 �4−𝑦𝑦12 1.954 831 759

𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑆𝑆1 = � =�
0 �8 − 𝑥𝑥𝑖𝑖 − 2 𝑦𝑦12 0 �7.821 367 206 − 𝑥𝑥𝑖𝑖 2

𝑏𝑏 − 𝑎𝑎 𝑏𝑏 + 𝑎𝑎 1.954 831 759 − 0 1.954 831 759 + 0

𝑥𝑥𝑖𝑖 = � � 𝑡𝑡𝑖𝑖 + � �=� � 𝑡𝑡𝑖𝑖 + � �
2 2 2 2

𝑥𝑥𝑖𝑖 = 0. 977 415 879𝑡𝑡𝑖𝑖 + 0. 977 415 879 , 𝑖𝑖 = 1, 2

1
𝑥𝑥1 = 0. 977 415 879 �− � + 0. 977 415 879 = 0.413 104 558
√3
1
𝑥𝑥2 = 0. 977 415 879 � � + 0. 977 415 879 = 1.541 727 200
√3

𝑏𝑏 − 𝑎𝑎
𝑆𝑆1 = � � [𝜆𝜆1 𝑓𝑓1 + 𝜆𝜆2 𝑓𝑓2 ] = 0. 977 415 879[1𝑓𝑓1 + 1𝑓𝑓2 ]
2

1 𝑓𝑓 = 0.361 533 943

𝑓𝑓𝑖𝑖 = ⇒� 1
�7.821 367 206 − 𝑥𝑥𝑖𝑖 2 𝑓𝑓2 = 0.428 571 428

44
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

𝑆𝑆1 = (0. 977 415 879)[1(0.361 533 943) + 1(0.428 571 428)] = 0.772 261 535 (𝑎𝑎)

𝑆𝑆2 𝑦𝑦2 = 1.577 350 269 �4−𝑦𝑦22 1.229 620 319

𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑆𝑆2 = � =�
0 �8 − 𝑥𝑥𝑖𝑖 −
2 𝑦𝑦22 0 �5.511 966 129 − 𝑥𝑥𝑖𝑖 2

𝑏𝑏 − 𝑎𝑎 𝑏𝑏 + 𝑎𝑎
𝑥𝑥𝑖𝑖 = � � 𝑡𝑡𝑖𝑖 + � � = 0.614 810 159𝑡𝑡𝑖𝑖 + 0.614 810 159 , 𝑖𝑖 = 1, 2
2 2

1
𝑥𝑥1 = 0.614 810 159 �− � + 0.614 810 159 = 0.259 849 348
√3
1
𝑥𝑥2 = 0.614 810 159 � � + 0.614 810 159 = 0.969 770 969
√3

𝑏𝑏 − 𝑎𝑎
𝑆𝑆1 = � � [𝜆𝜆1 𝑓𝑓1 + 𝜆𝜆2 𝑓𝑓2 ] = 0.614 810 159[1𝑓𝑓1 + 1𝑓𝑓2 ]
2

1 𝑓𝑓 = 0.428 571 428

𝑓𝑓𝑖𝑖 = ⇒� 1
�5.511 966 129 − 𝑥𝑥𝑖𝑖 2 𝑓𝑓2 = 0.467 702 987

𝑆𝑆2 = (0.614 810 159)[1(0.428 571 428) + 1(0.467 702 987)] = 0.551 038 615 (𝑏𝑏)

Obtención de la Integral Doble

Sustituyendo (𝑎𝑎), (𝑏𝑏) en (1), resulta,

𝐼𝐼2 = (1)[1𝑆𝑆0 + 1𝑆𝑆2 ] = (1)[1(0.772 261 535) + 1(0.551 038 615)] = 1.323 300 151 ⇽ 𝑆𝑆𝑆𝑆𝑆𝑆.

(b) Integral Doble sobre una Región Rectangular

Región 𝑅𝑅 de Integración Rectangular

Si la región 𝑅𝑅 de integración es rectangular, la integral doble estará expresada como

𝑏𝑏 𝑑𝑑 𝑑𝑑 𝑏𝑏
� 𝑓𝑓 (𝑥𝑥, 𝑦𝑦)𝑑𝑑𝑑𝑑 = � �� 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) 𝑑𝑑𝑑𝑑� 𝑑𝑑𝑑𝑑 = � �� 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) 𝑑𝑑𝑑𝑑� 𝑑𝑑𝑑𝑑 (4.1.5.5)
𝑎𝑎 𝑐𝑐 𝑐𝑐 𝑎𝑎
𝑅𝑅
En este caso el procedimiento de cálculo de la integral doble se reduce a una doble suma ponderada de valores de
la función 𝑓𝑓 (𝑥𝑥, 𝑦𝑦), de la siguiente forma:

45
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

1º Empleando las fórmulas de Newton-Cotes

(a) Si se aplica la Regla del Trapecio en ambas direcciones, la fórmula que se empleara es

𝑏𝑏 𝑑𝑑 𝑀𝑀 𝑁𝑁
(𝑏𝑏 − 𝑎𝑎)(𝑑𝑑 − 𝑐𝑐 )
� � 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 ≈ � � � �[𝜔𝜔𝑖𝑖 ]�𝜔𝜔𝑗𝑗 �𝑓𝑓�𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑗𝑗 � (4.1.5.6)
𝑎𝑎 𝑐𝑐 4𝑀𝑀𝑀𝑀
𝑖𝑖=0 𝑗𝑗=0

(b) Si se aplica la Regla de Simpson en ambas direcciones, la fórmula que se empleara es

𝑏𝑏 𝑑𝑑 𝑀𝑀 𝑁𝑁
(𝑏𝑏 − 𝑎𝑎)(𝑑𝑑 − 𝑐𝑐 )
� � 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 ≈ � � � �[𝜔𝜔𝑖𝑖 ]�𝜔𝜔𝑗𝑗 �𝑓𝑓�𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑗𝑗 � (4.1.5.7)
𝑎𝑎 𝑐𝑐 9𝑀𝑀𝑀𝑀
𝑖𝑖=0 𝑗𝑗=0

2º Empleando las fórmulas de Gauss-Legendre

𝑏𝑏 𝑑𝑑 𝑚𝑚 𝑛𝑛
(𝑏𝑏 − 𝑎𝑎)(𝑑𝑑 − 𝑐𝑐 )
� � 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 ≈ � � � �[𝜆𝜆𝑖𝑖 ]�𝜆𝜆𝑗𝑗 �𝑓𝑓�𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑗𝑗 �
𝑎𝑎 𝑐𝑐 4 (4.1.5.8)
𝑖𝑖=1 𝑗𝑗=1
Donde
𝑑𝑑 − 𝑐𝑐 𝑑𝑑 + 𝑐𝑐
𝑦𝑦𝑗𝑗 = � � 𝑡𝑡𝑗𝑗 + � �
� 2 2
𝑏𝑏 − 𝑎𝑎 𝑏𝑏 + 𝑎𝑎
𝑥𝑥𝑖𝑖 = � � 𝑡𝑡𝑖𝑖 + � �
2 2
𝑚𝑚 𝑛𝑛 𝑚𝑚 𝑛𝑛

� � 𝛼𝛼𝑖𝑖 𝛽𝛽𝑗𝑗 = � 𝛼𝛼𝑖𝑖 � 𝛽𝛽𝑗𝑗 = (𝛼𝛼0 + 𝛼𝛼1 + ⋯ + 𝛼𝛼𝑚𝑚 )(𝛽𝛽0 + 𝛽𝛽1 + ⋯ + 𝛽𝛽𝑛𝑛 )
𝑗𝑗=0 𝑖𝑖=0 𝑖𝑖=0 𝑗𝑗=0

Ejemplo

Empleando las Fórmulas de Newton Cotes y la Cuadratura de Gauss – Legendre, aproxime la Integral Doble

𝐼𝐼 = � 𝑒𝑒 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 𝑐𝑐𝑐𝑐𝑐𝑐 𝑦𝑦 𝑑𝑑𝑑𝑑 1 ≤ 𝑥𝑥 ≤ 2

𝑅𝑅: �
3 ≤ 𝑦𝑦 ≤ 4
𝑅𝑅
Resolución. - La Integral doble puede ser escrita como,

4 2 2 4
𝐼𝐼2 = � 𝑒𝑒 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 𝑐𝑐𝑐𝑐𝑐𝑐 𝑦𝑦 𝑑𝑑𝑑𝑑 = � � 𝑒𝑒 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 𝑐𝑐𝑐𝑐𝑐𝑐 𝑦𝑦 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 = � � 𝑒𝑒 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 𝑐𝑐𝑐𝑐𝑐𝑐 𝑦𝑦 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑
3 1 1 3
𝑅𝑅

(1) Fórmulas de Cuadratura de Newton-Cotes

Se empleará la Regla Compuesta de Simpson con 3 puntos (2 intervalos) en ambas direcciones

𝐸𝐸� 𝐼𝐼̂ 𝑃𝑃� 2 2

𝐸𝐸� 1 4 2 �� 𝛼𝛼𝑖𝑖 � �� 𝛽𝛽𝑗𝑗 � = (𝛼𝛼0 + 𝛼𝛼1 + 𝛼𝛼2 )(𝛽𝛽0 + 𝛽𝛽1 + 𝛽𝛽2 )

𝐼𝐼̂ 4 16 8 𝑖𝑖=0 𝑗𝑗=0
𝑃𝑃� 2 8 4 = (𝛼𝛼0 𝛽𝛽0 + 𝛼𝛼0 𝛽𝛽1 + 𝛼𝛼0 𝛽𝛽2 ) + (𝛼𝛼1 𝛽𝛽0 + 𝛼𝛼1 𝛽𝛽1 + 𝛼𝛼1 𝛽𝛽2 ) + (𝛼𝛼2 𝛽𝛽0 + 𝛼𝛼2 𝛽𝛽1 + 𝛼𝛼2 𝛽𝛽2 )

46
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

� 𝜔𝜔𝑖𝑖𝑖𝑖 = (1𝛼𝛼0 𝛽𝛽0 + 4𝛼𝛼0 𝛽𝛽1 + 1𝛼𝛼0 𝛽𝛽2 ) + (4𝛼𝛼1 𝛽𝛽0 + 16𝛼𝛼1 𝛽𝛽1 + 4𝛼𝛼1 𝛽𝛽2 ) + (1𝛼𝛼2 𝛽𝛽0 + 4𝛼𝛼2 𝛽𝛽1 + 1𝛼𝛼2 𝛽𝛽2 )
𝑖𝑖,𝑗𝑗

𝑏𝑏 − 𝑎𝑎 2 − 1 𝑖𝑖 0 1 2
ℎ= = = 0.5
2 2 𝑥𝑥𝑖𝑖 1 1.5 2 𝑓𝑓�𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑗𝑗 � = 𝑒𝑒 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥𝑖𝑖 𝑐𝑐𝑐𝑐𝑐𝑐 𝑦𝑦𝑗𝑗
𝑑𝑑 − 𝑐𝑐 4 − 3
𝑘𝑘 = = = 0.5
2 2 𝑗𝑗 0 1 2 𝑓𝑓 (𝑥𝑥1 , 𝑦𝑦2 ) = 𝑓𝑓 (1.5, 4) = 𝑒𝑒 − 𝑠𝑠𝑠𝑠𝑠𝑠 1.5 𝑐𝑐𝑐𝑐𝑐𝑐 4
𝑦𝑦𝑗𝑗 3 3.5 4

[1𝑓𝑓(𝑥𝑥0 , 𝑦𝑦0 ) + 4𝑓𝑓(𝑥𝑥0 , 𝑦𝑦1 ) + 1𝑓𝑓(𝑥𝑥0 , 𝑦𝑦2 )] +

ℎ 𝑘𝑘 ℎ𝑘𝑘
𝐼𝐼2 = � � � � = � 𝜔𝜔𝑖𝑖𝑖𝑖 𝑓𝑓𝑖𝑖𝑖𝑖 = � +[ 4𝑓𝑓 (𝑥𝑥1 , 𝑦𝑦0 ) + 16𝑓𝑓(𝑥𝑥1 , 𝑦𝑦1 ) + 4𝑓𝑓(𝑥𝑥1 , 𝑦𝑦2 )] +�
3 3 9
+[1𝑓𝑓(𝑥𝑥2 , 𝑦𝑦0 ) + 4𝑓𝑓(𝑥𝑥2 , 𝑦𝑦1 ) + 1𝑓𝑓(𝑥𝑥2 , 𝑦𝑦2 )]

(0.5)(0.5) (1)(2.300 323 951) + (4)(2.198 996 524) + (1)(1.733 291 395) +
𝐼𝐼2 = = �+(4)(2.684 548 486) + (16)(2.544 949 611) + (4)(1.919 378 711) +�
9
+(1)(2.460 089 250) + (4)(2. 343 200 327) + (1)(1.811 864 567)

𝐼𝐼2 = 2.378 035 976 ⇽ 𝑆𝑆𝑆𝑆𝑆𝑆.

(2) Fórmulas de Cuadratura de Gauss - Legendre

Se empleará la Cuadratura de Gauss Legendre con 2 puntos en ambas direcciones

2 2 2

� 𝛼𝛼𝑖𝑖 𝛽𝛽𝑗𝑗 = �� 𝛼𝛼𝑖𝑖 � �� 𝛽𝛽𝑗𝑗 � = (𝛼𝛼1 + 𝛼𝛼2 )(𝛽𝛽1 + 𝛽𝛽2 ) = (𝛼𝛼1 𝛽𝛽1 + 𝛼𝛼1 𝛽𝛽2 ) + (𝛼𝛼2 𝛽𝛽1 + 𝛼𝛼2 𝛽𝛽2 )
𝑖𝑖,𝑗𝑗=1 𝑖𝑖=1 𝑗𝑗=1
2
𝑏𝑏 − 𝑎𝑎 𝑑𝑑 − 𝑐𝑐 2−1 3−4
𝐼𝐼2 = � �� � � 𝜆𝜆𝑖𝑖 𝜆𝜆𝑗𝑗 𝑓𝑓𝑖𝑖𝑖𝑖 = � �� � [(𝜆𝜆1 𝜆𝜆1 𝑓𝑓11 + 𝜆𝜆1 𝜆𝜆2 𝑓𝑓12 ) + (𝜆𝜆2 𝜆𝜆1 𝑓𝑓21 + 𝜆𝜆2 𝜆𝜆2 𝑓𝑓22 )]
2 2 2 2
𝑖𝑖,𝑗𝑗=1

𝑏𝑏 − 𝑎𝑎 𝑏𝑏 + 𝑎𝑎 𝑓𝑓𝑖𝑖𝑖𝑖 = 𝑓𝑓�𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑗𝑗 � = 𝑒𝑒 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥𝑖𝑖 𝑐𝑐𝑐𝑐𝑐𝑐 𝑦𝑦𝑗𝑗

𝑥𝑥𝑖𝑖 = � � 𝑡𝑡𝑖𝑖 + � � = 0.5𝑡𝑡𝑖𝑖 + 1.5 , 𝑖𝑖 = 1, 2
2 2
𝑑𝑑 − 𝑐𝑐 𝑑𝑑 + 𝑐𝑐 𝑓𝑓12 = 𝑓𝑓 (𝑥𝑥1 , 𝑦𝑦2 ) = 𝑒𝑒 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥1 𝑐𝑐𝑐𝑐𝑐𝑐 𝑦𝑦2
𝑦𝑦𝑗𝑗 = � � 𝑡𝑡𝑗𝑗 + � � = 0.5𝑡𝑡𝑗𝑗 + 3.5 , 𝑗𝑗 = 1, 2
2 2
1 1
𝜆𝜆1 = 1 𝜆𝜆2 = 1 𝑡𝑡1 = − 𝑡𝑡2 =
√3 √3

𝑥𝑥1 = 1.211 324 865 𝑦𝑦1 = 3.211 324 865

𝑥𝑥2 = 1.788 675 135 𝑦𝑦2 = 3.788 675 135

(1)(1)(2.544 178 774) + (1)(1)(2.110 342 025) +

𝐼𝐼2 = 0.25 � �
+(1)(1)(2.648 478 467) + (1)(1)(2.179 255 930)

𝐼𝐼2 = 2.370 563 799 ⇽ 𝑆𝑆𝑆𝑆𝑆𝑆.

47
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

4.2.2.-Derivación Numérica mediante las fórmulas de aproximación por Diferencias

4.2.2.1.- Derivada de una Función de una Variable

Dado un conjunto discreto de puntos igualmente espaciados (𝑥𝑥𝑘𝑘 , 𝑦𝑦𝑘𝑘 ) donde 𝑘𝑘 = 0, ±1, ±2, ⋯ , 𝑒𝑒𝑒𝑒𝑒𝑒 ; tal que 𝑓𝑓𝑘𝑘 =
𝑓𝑓 (𝑥𝑥𝑘𝑘 ) = 𝑓𝑓(𝑥𝑥 + 𝑘𝑘ℎ). Los desarrollos de Taylor para esta función alrededor de una abscisa arbitraria𝑥𝑥 para (±ℎ)y
(±2ℎ)están dados por

1 1
𝑓𝑓(𝑥𝑥 + ℎ) = 𝑓𝑓 (𝑥𝑥 ) + ℎ𝑓𝑓 ′ (𝑥𝑥 ) + 2ℎ2 𝑓𝑓 ′′ (𝑥𝑥 ) + 6ℎ3 𝑓𝑓 ′′′ (𝑥𝑥) + 𝑂𝑂(ℎ4 )
1 1
𝑓𝑓(𝑥𝑥 − ℎ) = 𝑓𝑓 (𝑥𝑥 ) − ℎ𝑓𝑓 ′ (𝑥𝑥 ) + 2ℎ2 𝑓𝑓 ′′ (𝑥𝑥 ) − 6ℎ3 𝑓𝑓 ′′′ (𝑥𝑥) + 𝑂𝑂(ℎ4 )
4
(4.2.1)
𝑓𝑓(𝑥𝑥 + 2ℎ) = 𝑓𝑓(𝑥𝑥 ) + 2ℎ𝑓𝑓 ′ (𝑥𝑥 ) + 2ℎ2 𝑓𝑓 ′′ (𝑥𝑥) + 3ℎ3 𝑓𝑓 ′′′ (𝑥𝑥 ) + 𝑂𝑂(ℎ4 )
4
𝑓𝑓 (𝑥𝑥 + 2ℎ) = 𝑓𝑓 (𝑥𝑥 ) − 2ℎ𝑓𝑓 ′ (𝑥𝑥 ) + 2ℎ2 𝑓𝑓 ′′ (𝑥𝑥) − 3ℎ3 𝑓𝑓 ′′′ (𝑥𝑥 ) + 𝑂𝑂(ℎ4 )

Realizando algunas combinaciones con las anteriores expresiones se obtienen las fórmulas:

(a) Para una precisión de primer orden por diferencias progresivas

𝑓𝑓 (𝑥𝑥 + ℎ) − 𝑓𝑓 (𝑥𝑥 )
𝑓𝑓′(𝑥𝑥 ) = + 𝑂𝑂(ℎ)
ℎ (4.2.2)
𝑓𝑓 (𝑥𝑥 + 2ℎ) − 2𝑓𝑓(𝑥𝑥 + ℎ) + 𝑓𝑓 (𝑥𝑥 )
𝑓𝑓′′(𝑥𝑥 ) = + 𝑂𝑂(ℎ)
ℎ

(b) Para una precisión de segundo orden por diferencias centrales

𝑓𝑓 (𝑥𝑥 + ℎ) − 𝑓𝑓 (𝑥𝑥 − ℎ)
𝑓𝑓′(𝑥𝑥 ) = + 𝑂𝑂(ℎ2 )
2ℎ (4.2.3)
𝑓𝑓 (𝑥𝑥 + ℎ) − 2𝑓𝑓 (𝑥𝑥 ) + 𝑓𝑓(𝑥𝑥 − ℎ)
𝑓𝑓′′(𝑥𝑥 ) = + 𝑂𝑂(ℎ2 )
ℎ2

Ejemplo

Dada la función 𝑓𝑓 (𝑥𝑥 ) = 𝑥𝑥 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 , obtenga 𝑓𝑓 ′ (0.4) tome ℎ = 0.000 01

Resolución. - Se tomará 2 puntos y se aproximará la derivada por diferencias centrales, entonces

𝑥𝑥 − ℎ = 0.399 99 ⇒ 𝑓𝑓 (𝑥𝑥 − ℎ) = 0.155 759 758

𝑥𝑥 + ℎ = 0.400 01 ⇒ 𝑓𝑓 (𝑥𝑥 + ℎ) = 0.155 774 915

𝑓𝑓 (𝑥𝑥 + ℎ) − 𝑓𝑓 (𝑥𝑥 − ℎ)
𝑓𝑓′(𝑥𝑥 ) = + 𝑂𝑂(ℎ2 )
2ℎ

0.155 774 915 − 0.155 759 758

𝑓𝑓 ′ (0.4) ≅ = 0.757 871 750 ⇽ 𝑆𝑆𝑆𝑆𝑆𝑆.
2(0.000 01)

48
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

4.2.2.2.- Derivada Parciales de una Función de dos Variables

Sea la función 𝑢𝑢 = 𝑢𝑢(𝑥𝑥, 𝑦𝑦); los desarrollos de Taylor para las longitudes de paso ℎ en la dirección de 𝑥𝑥 ; 𝑘𝑘 en la
dirección de 𝑦𝑦 están dados por

1 2
𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦 + 𝑘𝑘) = 𝑢𝑢 + ℎ𝑢𝑢𝑥𝑥 + 𝑘𝑘𝑢𝑢𝑦𝑦 + �ℎ 𝑢𝑢𝑥𝑥𝑥𝑥 + 2ℎ𝑘𝑘𝑢𝑢𝑥𝑥𝑥𝑥 + 𝑘𝑘 2 𝑢𝑢𝑦𝑦𝑦𝑦 � + 𝑂𝑂(𝐻𝐻3 )
2!
1
𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦) = 𝑢𝑢 + ℎ𝑢𝑢𝑥𝑥 + ℎ2 𝑢𝑢𝑥𝑥𝑥𝑥 + 𝑂𝑂 (ℎ3 )
2!
1 (4.2.4)
𝑢𝑢(𝑥𝑥 − ℎ, 𝑦𝑦) = 𝑢𝑢 − ℎ𝑢𝑢𝑥𝑥 + ℎ2 𝑢𝑢𝑥𝑥𝑥𝑥 + 𝑂𝑂 (ℎ3 )
2!
𝑢𝑢(𝑥𝑥 + 2ℎ, 𝑦𝑦) = 𝑢𝑢 + 2ℎ𝑢𝑢𝑥𝑥 + 2ℎ2 𝑢𝑢𝑥𝑥𝑥𝑥 + 𝑂𝑂(ℎ3 )
1
𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦 − 𝑘𝑘) = 𝑢𝑢 + ℎ𝑢𝑢𝑥𝑥 − 𝑘𝑘𝑢𝑢𝑦𝑦 + �ℎ2 𝑢𝑢𝑥𝑥𝑥𝑥 − 2ℎ𝑘𝑘𝑢𝑢𝑥𝑥𝑥𝑥 + 𝑘𝑘 2 𝑢𝑢𝑦𝑦𝑦𝑦 � + 𝑂𝑂(𝐻𝐻3 )
2!

Tomando combinaciones lineales de las anteriores expresiones se obtienen las fórmulas:

(a) Para una precisión de primer orden por diferencias progresivas

𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦) − 𝑢𝑢(𝑥𝑥, 𝑦𝑦)

𝑢𝑢𝑥𝑥 (𝑥𝑥, 𝑦𝑦) = + 𝑂𝑂(ℎ)
ℎ
𝑢𝑢(𝑥𝑥 + 2ℎ, 𝑦𝑦) − 2𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦) + 𝑢𝑢(𝑥𝑥, 𝑦𝑦)
𝑢𝑢𝑥𝑥𝑥𝑥 (𝑥𝑥, 𝑦𝑦) = + 𝑂𝑂(ℎ) (4.2.5)
ℎ2
𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦 + 𝑘𝑘) − 𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦) − 𝑢𝑢(𝑥𝑥, 𝑦𝑦 + 𝑘𝑘) + 𝑢𝑢(𝑥𝑥, 𝑦𝑦)
𝑢𝑢𝑥𝑥𝑥𝑥 (𝑥𝑥, 𝑦𝑦) = + 𝑂𝑂(𝐻𝐻)
ℎ𝑘𝑘

(b) Para una precisión de segundo orden por diferencias Centrales

𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦) − 𝑢𝑢(𝑥𝑥, 𝑦𝑦)

𝑢𝑢𝑥𝑥 (𝑥𝑥, 𝑦𝑦) = + 𝑂𝑂(ℎ)
ℎ
𝑢𝑢(𝑥𝑥 + 2ℎ, 𝑦𝑦) − 2𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦) + 𝑢𝑢(𝑥𝑥, 𝑦𝑦)
𝑢𝑢𝑥𝑥𝑥𝑥 (𝑥𝑥, 𝑦𝑦) = + 𝑂𝑂(ℎ) (4.2.6)
ℎ2
𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦 + 𝑘𝑘) − 𝑢𝑢(𝑥𝑥 + ℎ, 𝑦𝑦) − 𝑢𝑢(𝑥𝑥, 𝑦𝑦 + 𝑘𝑘) + 𝑢𝑢(𝑥𝑥, 𝑦𝑦)
𝑢𝑢𝑥𝑥𝑥𝑥 (𝑥𝑥, 𝑦𝑦) = + 𝑂𝑂(𝐻𝐻)
ℎ𝑘𝑘

49
ING. EDGAR H. MEDINA TAPIA METODOS NUMERICOS I

Ejemplo
𝜕𝜕𝜕𝜕
Dada la función 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) = 𝑥𝑥 2 𝑦𝑦𝑒𝑒 −𝑥𝑥𝑥𝑥 , obtenga para 𝑥𝑥 = 0.5 , 𝑦𝑦 = 0.2 tome ℎ = 𝑘𝑘 = 0.000 01
𝜕𝜕𝜕𝜕

Resolución. - Se aproximará la derivada parcial empleando diferencias centrales, con 2 puntos,

𝑥𝑥 = 0.5 𝑦𝑦 − 𝑘𝑘 = 0.199 99 𝑓𝑓 (𝑥𝑥, 𝑦𝑦 − ℎ) = 0.045 239 835

𝑥𝑥 = 0.5 𝑦𝑦 = 0.200 00 𝑓𝑓 (𝑥𝑥, 𝑦𝑦) = 0.045 241 870
𝑥𝑥 = 0.5 𝑦𝑦 − 𝑘𝑘 = 0.200 01 ( )
𝑓𝑓 𝑥𝑥, 𝑦𝑦 + ℎ = 0.045 243 906

𝜕𝜕𝜕𝜕 𝑓𝑓 (𝑥𝑥, 𝑦𝑦 + 𝑘𝑘) − 𝑓𝑓 (𝑥𝑥, 𝑦𝑦 − 𝑘𝑘)

(𝑥𝑥, 𝑦𝑦) = + 𝑂𝑂(𝑘𝑘 2 )
𝜕𝜕𝜕𝜕 2𝑘𝑘

𝜕𝜕𝜕𝜕 0.045 243 906 − 0.045 239 835

(0.5, 0.2) ≅ = 0.203 588 77 ⇽ 𝑆𝑆𝑆𝑆𝑆𝑆.
𝜕𝜕𝜕𝜕 2(0.000 01)

50