Mna1 T05-1

Tema 5
Integración numérica
Dra. Neus Garrido
Máster en Ingenierı́a Matemática y Computación

Escuela Superior en Ingenierı́a y Tecnologı́a
Métodos Numéricos Aplicados I PER 3362 1 / 63

Contenido
1 Introducción
2 Fórmulas cerradas de Newton-Cotes

Método de Trapecios
Método de Simpson
3 Fórmulas abiertas de Newton-Cotes

Método de Punto Medio
4 Cuadratura de Gauss
Cuadratura de Gauss-Legendre
Cuadratura de Gauss-Chebyshev
Cuadratura de Gauss-Laguerre
Cuadratura de Gauss-Hermite
5 Integración múltiple
Nodos equiespaciados
Nodos no equiespaciados

1
Introducción

Introducción
Aproximación de la integral
Longitud de arco Z 1 p
L= 1 + cos2 xdx
0
Función de Bessel Z π
1
Jn (x) = cos (x sin θ − nθ)dθ
x 0
Función error Z x
2 2
erf (x) = √ e−t dt
π 0
Función de distribución normal en un proceso de fabricación

Z x
2 2
φ(x) = √ e−t /2 dt
2π −∞

Introducción
Ejemplo 1. Fuerza total ejercida por el mástil de un velero

La fuerza total ejercida por el mástil de un velero se
define como Z 30
F = f (z)dz,
0
donde z es la distancia vertical a la cubierta.

Se utiliza un modelo a escala en un túnel de viento
para medir la fuerza ejercida por el mástil en
diferentes puntos del mismo. En la siguiente tabla se
observan dichas mediciones en función de la distancia
respecto a la cubierta:
z 0 3 6 9 12 15 18 21 24 27 30
f (z) 190 141 104 77.5 54.7 42.5 31.5 23.3 17.3 12.8 9.5
El cálculo de la fuerza total, es decir, de una integral definida cuyo integrando no posee
una expresión analı́tica, es imprescindible para un correcto diseño del mástil.
¿Qué podemos hacer si nuestra función no tiene una expresión explı́cita y solo conocemos
la función en algunos puntos?
Introducción
INTEGRACIÓN NUMÉRICA
N O D O S E Q U I E S PA C I A D O S N O D O S N O E Q U I E S PA C I A D O S
Newton-Cotes Cuadratura de Gauss
Métodos cerrados Métodos abiertos

Gauss-Legendre Gauss-Hermite
Trapecios Punto medio
Simpson
Gauss-Chebyshev Gauss-Laguerre
I N T E G R A C I Ó N S I M P L E Y M Ú LT I P L E
Objetivos
○ Conocer la expresión general de las técnicas de cuadratura
○ Comprender e implementar las fórmulas derivadas de las expresiones de
Newton-Cotes, cerradas y abiertas
○ Comprender e implementar los casos particulares de la cuadratura de Gauss
○ Aplicar la integración numérica sobre integración múltiple
Introducción
Cuadratura numérica
La técnica de cuadratura numérica consiste en calcular una integral a partir de unos
valores discretos, obtenidos del valor del integrando en los diferentes nodos o en valores
tabulados.
En general, la cuadratura numérica consiste en obtener
Z b Xn
I= f (x) dx ≈ ai f (xi ).
a i=0

Introducción
tabulados.
Z b Xn
a i=0
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 0.5 1 1.5 2 0 0.5 1 1.5 2
x x
Introducción
tabulados.
Z b Xn
a i=0
Polinomio de interpolación de Lagrange de grado n

n n
X Y (x − xj )
ln (x) = Li (x)f (xi ), Li (x) =
(xi − xj )
i=0 j=0
j6=i
f (n+1) (ξ(x)) Qn
Error cometido en f (x) ≈ ln (x): (x) = i=0 (x − xi ).
(n + 1)!

Introducción
Cuadratura numérica y polinomio de Lagrange
Z b Z b
I = f (x) dx = (ln (x) + (x)) dx
a a

Introducción
Z b Z b
I = f (x) dx = (ln (x) + (x)) dx
a a
(
Z b X n n
)
f (n+1) (ξ(x)) Y
= Li (x)f (xi ) + (x − xi ) dx
a i=0
(n + 1)! i=0

Introducción
Z b Z b
I = f (x) dx = (ln (x) + (x)) dx
a a
(
Z b X n n
)
f (n+1) (ξ(x)) Y
= Li (x)f (xi ) + (x − xi ) dx
a i=0
(n + 1)! i=0
n Z b Z b n
X 1 Y
= f (xi ) Li (x) dx + f (n+1) (ξ(x)) (x − xi ) dx.
i=0 a (n + 1)! a i=0

Introducción
Z b Z b
I = f (x) dx = (ln (x) + (x)) dx
a a
(
Z b X n n
)
f (n+1) (ξ(x)) Y
= Li (x)f (xi ) + (x − xi ) dx
a i=0
(n + 1)! i=0
n Z b Z b n
X 1 Y
i=0 a (n + 1)! a i=0
Por tanto: n
X Z b
I≈ f (xi ) Li (x) dx,
i=0 a
denominándose error de cuadratura al término

Z b n
1 Y
f (n+1) (ξ(x)) (x − xi ) dx.
(n + 1)! a i=0

Introducción
Z b Z b
I = f (x) dx = (ln (x) + (x)) dx
a a
(
Z b X n n
)
f (n+1) (ξ(x)) Y
= Li (x)f (xi ) + (x − xi ) dx
a i=0
(n + 1)! i=0
n Z b Z b n
X 1 Y
i=0 a (n + 1)! a i=0
Por tanto: n
X Z b
I≈ f (xi ) Li (x) dx,
i=0 a
denominándose error de cuadratura al término

Z b n
1 Y
f (n+1) (ξ(x)) (x − xi ) dx.
(n + 1)! a i=0
ú nodos equiespaciados fórmulas de Newton-Cotes

ú nodos NO equiespaciados cuadratura de Gauss
2
Fórmulas cerradas de Newton-Cotes

Xn Z b
𝑥𝑖
I ≈𝑥 f (xi ) Li (x) dx,
𝑖+1 a
i=0
𝑎 𝑥 b−a
Nodos equiespaciados: x ∈ [a, b], h = n
, {x0 = a, x1 , x2 , . . . , xn = b}
𝑎 𝑎+ℎ 𝑎 + 2ℎ 𝑏−ℎ 𝑏
Teorema 1 (Error de las fórmulas

𝑥𝑖−1 𝑥𝑖 cerradas
𝑥𝑖+1 de Newton-Cotes)
𝑎 𝑎+ℎ 𝑎 + 2ℎ 𝑏−ℎ
Pn
Sea i=1 ai f (xi ) la fórmula de Newton-Cotes cerrada de n + 1 puntos donde x0 = a,
xn = b y h = b−a
n
. Entonces existe un ξ ∈ (a, b) tal que
para valores pares de n y si f ∈ C n+2 [a, b] se cumple que
b n n
hn+3 (n+2)
Z X Z
f (x) dx = ai f (xi ) + f (ξ) τ 2 (τ − 1) · · · (τ − n) dτ.
a i=0
(n + 2)! 0
para valores impares de n y si f ∈ C n+1 [a, b] se cumple que

b n n
hn+2 (n+1)
Z X Z
f (x) dx = ai f (xi ) + f (ξ) τ (τ − 1) · · · (τ − n) dτ.
a i=0
(n + 1)! 0


Rb
Aproximación de I = a
f (x) dx:
Trapecios:
h h3 00
I≈ (f (a) + f (b)) , =− f (ξ)
2 12
Simpson:
h5 (iv)

h a+b
I≈ f (a) + 4f + f (b) , =− f (ξ)
3 2 90
Simpson 3/8:
3h5 (iv)

3h 2a + b a + 2b
I≈ f (a) + 3f + 3f + f (b) , =− f (ξ)
8 3 3 80
Milne:

2h 3a + b a+b a + 3b
I≈ 7f (a) + 32f + 12f + 32f + 7f (b) ,
45 4 2 4
8h7 (vi)
=− f (ξ)
945
Contenidos
1 Introducción

Método de Simpson

Desarrollo
Puntos (n = 1): x0 = a, x1 = b
Z b Z b
I ≈ f (a) L0 (x) dx + f (b) L1 (x) dx
a a
Z b Z b
x−b x−a
= f (a) dx + f (b) dx
a a−b b−a a
f (a) b f (b) b
Z Z
= − (x − b) dx + (x − a) dx
h a h a
h
= (f (a) + f (b)) .
2
h
I≈ (f (a) + f (b))
2
3
Error: = − h12 f 00 (ξ(x))
Reducción del error

1. Dividir el intervalo [a, b] en subintervalos
2. Aplicar el método de trapecios sobre cada subintervalo
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 1 2 0 1 2
x x
Método de trapecios compuesto
n−1
hX h
I≈ (f (xi ) + f (xi+1 )) = (f (x0 ) + 2f (x1 ) + 2f (x2 ) + · · · + 2f (xn−1 ) + f (xn ))
2 i=0 2
2
Error: = − h12 (b − a)f 00 (ξ(x))
Trapecios.m
function I=Trapecios(f,a,b,n)
% fórmula de Trapecios compuesta
h=(b-a)/n;
x=a:h:b;
pesos=[1 2*ones(1,n-1) 1];
I=h/2*sum(pesos.*f(x));
end

Contenidos
1 Introducción

Método de Simpson

Método de Simpson
Desarrollo
a+b
Puntos (n = 2): x0 = a, x1 = 2
, x2 = b
Z b Z b Z b
a+b
I ≈ f (a) L0 (x) dx + f L1 (x) dx + f (b) L2 (x) dx
a 2 a a

h a+b
= f (a) + 4f + f (b)
3 2
Método de Simpson

h a+b
I≈ f (a) + 4f + f (b)
3 2
h5 (iv)
Error: = − f (ξ(x))
90

Método de Simpson
Método de Simpson compuesto

Aplicando la regla de Simpson a cada subintervalo [xj−1 , xj+1 ], j = 1, 3, 5, . . . , n − 1,
con n par:
h
I ≈ (f (x0 ) + 4f (x1 ) + 2f (x2 ) + 4f (x3 ) · · · + 2f (xn−2 ) + 4f (xn−1 ) + f (xn ))
3
n −1 n

2 2
h X X
= f (x0 ) + 2 f (x2j ) + 4 f (x2j−1 )
3 j=1 j=1
h4
Error: = − (b − a)f (iv) (ξ(x))
180
Simpson.m
function I=Simpson(f,a,b,n)
% fórmula de Simpson compuesta
h=(b-a)/n; x=a:h:b;
pesos=ones(1,n+1); pesos(2:2:n)=4; pesos(3:2:n-1)=2;
end
Método de Simpson
Simpson.m
function I=Simpson(f,a,b,n)
% fórmula de Simpson compuesta
h=(b-a)/n;
x=a:h:b;
pesos=ones(1,n+1);
pesos(2:2:n)=4; pesos(3:2:n-1)=2;
end

Método de Simpson
Rπ
Ejemplo 2. Calcula I = 02 sin(x)e−x dx, utilizando los métodos de Trapecios y
Simpson tomando 4 y 8 subintervalos. Para cada caso, calcula el error cometido
−π
sabiendo que el resultado analı́tico es 1−e2 2 .
>> f=@(x) sin(x).*exp(-x);

>> a=0; b=pi/2; n1=4; n2=8;
Trapecios Simpson T S
n=4 0.380591 0.395839 0.015469 0.000221
n=8 0.392183 0.396047 0.003877 0.000013
Para un mismo método, el error se reduce cuando disminuye el tamaño de h

Las fórmulas de Simpson tienen menor error que las de Trapecios

3
Fórmulas abiertas de Newton-Cotes

n
X Z b
I≈ f (xi ) Li (x) dx
i=0 a
Nodos equiespaciados: x ∈ [a, b]

El primer y último nodo no se toman: {a = x−1 , x0 , x1 , x2 , . . . , xn , xn+1 = b}
ℎ 𝑎 + 2ℎ 𝑏−ℎ 𝑏 b−a
h=
n+2
𝑥𝑖+1
𝑥−1 𝑥0 𝑥1 𝑥𝑛 𝑥𝑛+1
𝑎 𝑎+ℎ 𝑎 + 2ℎ 𝑏−ℎ 𝑏

Teorema 2 (Error de las fórmulas abiertas de Newton-Cotes)

Sea n
P
i=1 ai f (xi ) la fórmula de Newton-Cotes abierta de n + 1 puntos donde x−1 = a,
b−a
xn+1 = b y h = n+2 . Entonces existe un ξ ∈ (a, b) tal que
para valores pares de n y si f ∈ C n+2 [a, b] se cumple que
b n n+1
hn+3 (n+2)
Z X Z
f (x) dx = ai f (xi ) + f (ξ) τ 2 (τ − 1) · · · (τ − n) dτ.
a i=0
(n + 2)! −1
para valores impares de n y si f ∈ C n+1 [a, b] se cumple que

b n n+1
hn+2 (n+1)
Z X Z
f (x) dx = ai f (xi ) + f (ξ) τ (τ − 1) · · · (τ − n) dτ.
a i=0
(n + 1)! −1


Rb
Aproximación de a
f (x) dx:
Punto medio, n = 0:
h3 00

a+b
I ≈ (b − a)f , = f (ξ)
2 3
n = 1:
3h3 00

b−a 2a + b a + 2b
I≈ f +f , = f (ξ)
2 3 3 4
n = 2:
14h5 (iv)

b−a 3a + b a+b a + 3b
I≈ 2f +f + 2f , = f (ξ)
3 4 2 4 45

Contenidos
1 Introducción


Desarrollo
a+b
Puntos (n = 0): x−1 = a, x0 = 2
, x1 = b
Z b
a+b
I ≈ f (x0 ) dx = (b − a)f (x0 ) = (b − a)f
a 2

a+b
I ≈ (b − a)f
2
h3 00
Error: = − f (ξ(x))
3

Reducción del error

1. Dividir el intervalo [a, b] en subintervalos
2. Aplicar el método de punto medio a cada subintervalo
El método de punto medio abarca los puntos anterior y posterior:
x−1 ←- x0 ,→ x1
ú Tomaremos solo los puntos pares para no solapar

ú El número de subintervalos n debe ser par
Método de Punto Medio compuesto
n/2
X
I ≈ (x2i+1 − x2i−1 )f (x2i )
i=0
= 2h (f (x0 ) + f (x2 ) + · · · + f (xn−2 ) + f (xn ))
b−a 2 00
Error: = 6
h f (ξ)
PuntoMedio.m
function I=PuntoMedio(f,a,b,n)
% fórmula de punto medio compuesta
h=(b-a)/(n+2);
x=a+h:h:b-h;
I=2*h*sum(f(x(1:2:end)));
end

Ejemplo 3. El trabajo W necesario para trasladar un cuerpo de una posición a otra

Rb
se puede calcular como W = a F (x) cos(α(x)) dx. Con los datos de la tabla,
calcula el trabajo realizado utilizando el método de Punto Medio.
𝐹
x 0.00 1.52 3.04 4.56 6.08 7.60 9.12
F (x) 0.00 40.04 57.83 62.28 46.71 53.38 22.24 𝛼
α(x) 0.50 1.40 0.75 0.90 1.30 1.48 1.50


Rb
𝐹
x 0.00 1.52 3.04 4.56 6.08 7.60 9.12
F (x) 0.00 40.04 57.83 62.28 46.71 53.38 22.24 𝛼
α(x) 0.50 1.40 0.75 0.90 1.30 1.48 1.50
40
Como h = 1.52, para tomar todos los puntos definimos: 30
20
1.52
xi = i, i = −1, 0, . . . , 13, 10
2 0
-1.52 0 1.52 3.04 4.56 6.08 7.6 9.12 10.64
x


Rb
𝐹
x 0.00 1.52 3.04 4.56 6.08 7.60 9.12
F (x) 0.00 40.04 57.83 62.28 46.71 53.38 22.24 𝛼
α(x) 0.50 1.40 0.75 0.90 1.30 1.48 1.50
40
20
1.52
xi = i, i = −1, 0, . . . , 13, 10
2 0
-1.52 0 1.52 3.04 4.56 6.08 7.6 9.12 10.64
x
1.52 1.52 1.52 + 3.04

x−1 = − , x0 = 0, x1 = , x2 = 1.52, x3 = , x4 = 3.04,
2 2 2
3.04 + 4.56 4.56 + 6.08 6.08 + 7.60
x5 = − , x6 = 4.56, x7 = , x8 = 6.08, x9 = ,
2 2 2
7.60 + 9.12 1.52
x10 = 7.60, x11 = − , x12 = 9.12, x13 = 9.12 + .
2 2

Rb
𝐹
x 0.00 1.52 3.04 4.56 6.08 7.60 9.12
F (x) 0.00 40.04 57.83 62.28 46.71 53.38 22.24 𝛼
α(x) 0.50 1.40 0.75 0.90 1.30 1.48 1.50
40
20
1.52
xi = i, i = −1, 0, . . . , 13, 10
2 0
-1.52 0 1.52 3.04 4.56 6.08 7.6 9.12 10.64
x
Aplicando el método de punto medio compuesto:
I = 2h (f (x0 ) + f (x2 ) + · · · + f (x10 ) + f (x12 )) = 162.24638.

4
Cuadratura de Gauss

Cuadratura de Gauss
Objetivo
Obtención de los nodos x1 , x2 , · · · , xn y de los coeficientes c1 , c2 , . . . , cn que minimizan
el error obtenido en la aproximación:
Z b Xn
w(x)f (x) dx ≈ ci f (xi ),
a i=1
w(x): función peso, y cumple w(x) > 0, x ∈ [a, b]

Para nodos no equiespaciados
Obtención de los coeficientes ci

1. Determinar unos polinomios p1 (x), . . . , pn (x) tales que
hpi , pj i = 0, i 6= j.
2. Los nodos x1 , . . . , xn son las raı́ces del polinomio pn (x).

A partir de estas dos condiciones, se pueden obtener los coeficientes c1 , . . . , cn y, por
tanto, la fórmula de cuadratura.

Cuadratura de Gauss
Fórmulas de cuadratura
w(x) = 1, [a, b] = [−1, 1] ú Cuadratura de Gauss-Legendre
w(x) = √ 1 2 , [a, b] = [−1, 1] ú Cuadratura de Gauss-Chebyshev
1−x
w(x) = e−x , [a, b] = [0, +∞) ú Cuadratura de Gauss-Laguerre

−x2
w(x) = e , [a, b] = (−∞, +∞) ú Cuadratura de Gauss-Hermite
Error
b
f (2n) (ξ(x))
Z
= p2n (x)w(x) dx, ξ ∈ (a, b).
(2n)! a

Contenidos
1 Introducción

Z 1 n
X
f (x) dx ≈ ci f (xi )
−1 i=1
ú w(x) = 1
ú [a, b] = [−1, 1]
ú Polinomios pi (x):
p0 (x) = 1,
p1 (x) = x,
1
pk+1 (x) = k+1
[(2k + 1)xpk (x) − kpk−1 (x)] , k = 1, 2, . . . , n − 1.
ú Raı́ces de pn (x):
1 1 4k − 1
xi = 1 − 2 + 3 cos π
8n 8n 4n + 2
ú Coeficientes ci :
2
ci =
(1 − x2i )(p0n (xi ))2

i 1 2 3 4 5
xi -0.577350 0.577350 - - -
n=2
ci 1.000000 1.000000 - - -
xi 0.000000 -0.774597 0.774597 - -
n=3
ci 0.888889 0.555556 0.555556 - -
xi -0.339981 -0.861136 0.339981 0.861136 -
n=4
ci 0.652145 0.347855 0.652145 0.347855 -
xi 0.000000 -0.538469 -0.906180 0.538469 0.906180
n=5
ci 0.568889 0.478629 0.236927 0.478629 0.236927
Tabla: Nodos y coeficientes de la cuadratura de Gauss-Legendre
Ejemplo 4. Cuadratura de Gauss-Legendre para n = 4

Z 1
f (x) dx ≈ c1 f (x1 ) + c2 f (x2 ) + c3 f (x3 ) + c4 f (x4 )
−1
= 0.652145f (−0.339981) + 0.347855f (−0.861136)
+0.652145f (0.339981) + 0.347855f (0.861136)

Intervalo genérico [a, b]

[a, b] [−1, 1]
Utilizamos el cambio de variable
b−a b+a b−a
x= y+ , dx = dy,
2 2 2
por lo que podremos aproximar
Z b n
b−a X b−a b+a
I= f (x) dx ≈ ci f xi + .
a 2 i=1 2 2

R 1.5 2
Ejemplo 5. Obtén el valor de I = 1 e−x dx con la cuadratura de
Gauss-Legendre utilizando n = 2 y n = 3
Transformamos el intervalo [1, 1.5] en el [−1, 1]:
1.5 − 1 1.5 + 1 y+5 dy
x= y+ = , dx =
2 2 4 4
de modo que Z 1.5 Z 1 (y+5)2
2 1
I= e−x dx = e− 16 dy.
1 4 −1
n = 2:
Z 1 (y+5)2

(−0.577350+5)2 (0.577350+5)2

1 1
e− 16 dy ≈ e− 16 + e− 16 = 0.109400.
4 −1 4
n = 3:
Z 1
(y+5)2
1
4
e− 16 dy ≈
−1
(0+5)2 (−0.774597+5)2 (0.774597+5)2
≈ 14 0.888889e− 16 + 0.555556e− 16 + 0.555556e− 16
= 0.109364.
Contenidos
1 Introducción

1
ú w(x) = √
1−x2
ú [a, b] = [−1, 1]
p0 (x) = 1,
p1 (x) = x,
pk+1 (x) = 2xpk−1 (x) − pk−2 , k = 2, 3, . . . , n.
2k − 1
xi = cos π
2n
π
ci =
n
Aproximamos la integral:
Z 1 n n
f (x) X πX
√ dx ≈ ci f (xi ) = f (xi )
−1 1 − x2 i=1
n i=1

i 1 2 3 4 5
n=2 xi -0.707107 0.707107 - - -
n=3 xi -0.866025 0.000000 0.866025 - -
n=4 xi -0.923880 -0.382683 0.382683 0.923880 -
n=5 xi -0.951057 -0.587785 0.000000 0.587785 0.951057
Tabla: Nodos de la cuadratura de Gauss-Chebyshev
Ejemplo 6. Cuadratura de Gauss-Chebyshev para n = 4

Z 1
f (x) π
√ dx ≈ (f (−0.923880) + f (−0.382683) + f (0.382683) + f (0.923880))
−1 1 − x2 4

R1 x
Ejemplo 7. Calcula la integral I = √e dx con seis decimales exactos
−1 1−x2
Tenemos que determinar el valor de n tal que el error satisface
2π
= eξ < 10−6 , ξ ∈ (−1, 1).
22n (2n)!
Acotemos el error :

2π 2π
|eξ | ≤ e ξ

⇔ || = 2n
e ≤
2n
e .
2 (2n)! 2 (2n)!
Calculamos para diferentes valores de n, obteniendo:
n = 1 : || < 2.134933555
n = 2 : || < 0.044477782
n = 3 : || < 3.70648 · 10−4
n = 4 : || < 1.65468 · 10−6
n = 5 : || < 4.59633 · 10−9 ⇒ para n = 5 se garantiza que el error es menor a 10−6
π −0.951057
+ e−0.587785 + e0 + e0.587785 + e0.951057 = 3.977463

I≈ e
5
Contenidos
1 Introducción

Z +∞ n
X
e−x f (x) dx ≈ ci f (xi )
0 i=1
ú w(x) = e−x
ú [a, b] = [0, +∞)
p0 (x) = 1,
p1 (x) = 1 − x,
pk+2 (x) = (2k + 3 − x)pk+1 (x) − (k + 1)2 pk (x), k = 0, 1, . . . , n − 2.
2 2

j0i −2 + j0k
xi = 1+ ,
4hn 48h2n
siendo hn = n + 21 y j0i la raı́z i-ésima de la función J0 (x), es decir, la función de
Bessel de primera especie y orden.
(n!)2 xi
ci = 2 .
pn+1 (xi )
R +∞
Ejemplo 8. Calcula la integral I = 0
e−10x sin(x) dx utilizando la cuadratura
de Gauss-Laguerre para n = 3
R +∞
El dominio de la integral está en [0, +∞), pero no es de la forma I = 0
e−x f (x) dx
ú Cambio de variable:
y = 10x ⇒ dy = 10 dx
entonces Z +∞ Z +∞
1 y
I= e−10x sin(x) dx = e−y sin dy.
0 10 0 10
Polinomios de Laguerre hasta n = 4:
p0 (x) = 1,
p1 (x) = 1 − x,
p2 (x) = x2 − 4x + 2,
p3 (x) = −x3 + 9x2 − 18x + 6,
p4 (x) = x4 − 16x3 + 72x2 − 96x + 24,

R +∞
Ejemplo 8. Calcula la integral I = 0
e−10x sin(x) dx utilizando la cuadratura
de Gauss-Laguerre para n = 3
Calculamos las raı́ces de p3 (x): Obtenemos los coeficientes:

>> syms x >> p4=x.∧4-16*x.∧3+72*x.∧2-96*x+24;
>> p3=-x.∧3+9*x.∧2-18*x+6; >> ci=(factorial(3))∧2*xi./...
>> xi=double(solve(L3==0)) (double(subs(p4,x,xi))).∧2
xi = ci =
0.415774556783479 0.711093009929173
2.294280360279042 0.278517733569241
6.289945082937479 0.010389256501586
Por tanto, el valor de la integral para n = 3 será:

Z +∞ 3
1 y 1 X
I= e−y sin dy ≈ ci f (xi )
10 0 10 10 i=1
>> f=@(y) sin(y/10);

>> I=1/10*sum(f(xi).*ci)
I=
0.009900991829812
Contenidos
1 Introducción

Z +∞ n
2 X
e−x f (x) dx ≈ ci f (xi )
−∞ i=1
2
ú w(x) = e−x
ú [a, b] = (−∞, +∞)
p0 (x) = 1,
p1 (x) = 2x,
pk+2 (x) = 2xpk+1 (x) − 2(k + 1)pk (x), k = 0, 1, . . . , n − 2.
ú Coeficientes ci : √
2n−1 n! π
ci = 2
.
n2 Hn−1 (xi )

R +∞ 2
Ejemplo 9. Calcula la integral I = −∞
e−4x |x| dx utilizando la cuadratura de
Gauss-Hermite para n = 4
Cambio de variable:
y = 2x ⇒ dy = 2 dx
entonces Z +∞ Z +∞
2 1 2
I= e−4x |x| dx = e−y |y| dy.
−∞ 4 −∞
Calculamos los polinomios de Hermite hasta n = 4:
p0 (x) = 1,
p1 (x) = 2x,
p2 (x) = 4x2 − 2,
p3 (x) = 8x3 − 12x,
p4 (x) = 16x4 − 48x2 + 12.

R +∞ 2
Ejemplo 9. Calcula la integral I = −∞
e−4x |x| dx utilizando la cuadratura de
Gauss-Hermite para n = 4
Calculamos las raı́ces de p4 (x):
>> xi=roots([16 0 -48 0 12])
xi =
-1.650680123885785
1.650680123885786
-0.524647623275290
0.524647623275290
Obtenemos los coeficientes:
>> syms x
>> p3=8*x.∧3-12*x;
>> ci=2∧(4-1)*factorial(4)*sqrt(pi)/4∧2./(double(subs(p3,x,xi)).∧2);
El valor de la integral para n = 4 será:
f=@(y) abs(y);
>> I=1/4*sum(ci.*f(xi))
I =
0.278258872775874

5
Integración múltiple

Integración sobre más de una variable: f : Rn → R

Consideraremos n = 2:
1
0.9
0.8
0.5
2
f : R → R, (x, y) ∈ [a, b] × [c, d] 0.7
0
2
0.5 0.6
0
0
y -2 -0.5
x
Objetivo
Aproximación de integrales dobles utilizando:
Nodos equiespaciados (método de Trapecios)
Nodos no equiespaciados (fórmulas de cuadratura de Gauss)

Contenidos
1 Introducción

Podemos escribir la integral doble como

x Z b Z d
I= f (x, y) dA = f (x, y) dy dx.
R a c
Rd
Si nombramos g(x) = c
f (x, y) dy, entonces
Z b
I= g(x) dx.
a
d−c
b−a
Con la notación h = , k=
, aproximamos la integral como
n m
Z b n−1 n−1 Z d Z d
hX hX
I= g(x) dx ≈ g(xi ) + g(xi+1 ) = f (xi , y) dy + f (xi+1 , y) dy
a 2 i=0 2 i=0 c c
n−1
( m−1 m−1
)
hX k X k X
≈ f (xi , yj ) + f (xi , yj+1 ) + f (xi , yj ) + f (xi , yj+1 ) .
2 i=0 2 j=0 2 j=0

Aproximación de la integral doble con la regla de Trapecios

Z b Z d
I= f (x, y) dy dx
a c
n−1
(m−1 )
hk X X
≈ f (xi , yj ) + f (xi , yj+1 ) + f (xi+1 , yj ) + f (xi+1 , yj+1 ) ,
4 i=0 j=0
b−a d−c
siendo h = n
, k= m
.
Ejemplo 10. Regla de Trapecios para n = 2 y m = 2

𝑦2 1 2 1
hk
I≈ [f (x0 , y0 ) + 2f (x0 , y1 ) + f (x0 , y2 )+ 𝑦1 2 4 2
4
+ 2f (x1 , y0 ) + 4f (x1 , y1 ) + 2f (x1 , y2 )+
𝑦0 1 2 1
+ f (x2 , y0 ) + 2f (x2 , y1 ) + f (x2 , y2 )]
𝑥0 𝑥1 𝑥2

𝑦𝑚 1 2 2 2 1
𝑦𝑚−1 2 4 4 4 2
𝑦2 2 4 4 4 2
2 4 4 4 2
𝑦1
1 2 2 2 1
𝑦0
𝑥0 𝑥1 𝑥2 𝑥𝑛−1 𝑥𝑛
Figura: Distribución de los nodos para cualesquiera n y m

Contenidos
1 Introducción

ú Nodos no equiespaciados ⇒ cuadratura de Gauss

ú Transformar el recinto de integración [a, b] × [c, d] en el que corresponda a la
cuadratura

Enunciado
Calcula la integral
s
x ∂f
2
∂f
2
I= (x, y) + (x, y) dA,
∂x ∂y
R
donde f (x, y) es la superficie de la semiesfera
x2 + y 2 + f 2 (x, y) = 9,
y el recinto de integración es
R = {(x, y) ∈ R2 : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1}.
Utiliza los métodos de Simpson con n = m = 8 y de Gauss-Legendre con n = 4.

Solución
p
La función f que aparece en el integrando es f (x, y) = 9 − x2 − y 2 , y sus derivadas
parciales
∂f x ∂f y
(x, y) = − p , (x, y) = − p ,
∂x 9−x −y 2 2 ∂y 9 − x2 − y 2
de modo que el integrando es
s 2 2 s
∂f ∂f x2 + y 2
(x, y) + (x, y) = .
∂x ∂y 9 − x2 − y 2
El recinto de integración es un cuadrado [0, 1] × [0, 1], por lo que la integral a resolver es
Z 1Z 1s
x2 + y 2
I= dy dx.
0 0 9 − x2 − y 2

Solución (Método de Simpson con n = m = 8)

La distribución de nodos del método de Simpson es la expuesta en la siguiente tabla para
n = m = 8:
1 4 2 4 2 4 2 4 1
4 16 8 16 8 16 8 16 4
2 8 4 8 4 8 4 8 2
4 16 8 16 8 16 8 16 4
2 8 4 8 4 8 4 8 2
4 16 8 16 8 16 8 16 4
2 8 4 8 4 8 4 8 2
4 16 8 16 8 16 8 16 4
1 4 2 4 2 4 2 4 1
El resultado es
I ≈ 0.267814255559730.

Solución (Método de Gauss-Legendre con n = 4)

Para aplicar el método de Gauss-Legendre, aplicamos los cambios de variables:
u = 2x − 1 ⇒ du = 2 dx, v = 2y − 1 ⇒ dv = 2 dy.
Por tanto,
s r
Z 1 Z 1 Z 1 Z 1
x2 + y 2 1 u2 + 2u + v 2 + 2v + 2
I= dy dx = dv du
0 0 9 − x2 − y 2 4 −1 −1 34 − u2 − 2u − v 2 − 2v
Los coeficientes y los nodos de la cuadratura de Gauss-Legendre para n = 4 son
u1 = v1 = −0.339981, c1 = 0.652145;
u2 = v2 = −0.861136, c2 = 0.347855;
u3 = v3 = −u1 , c3 = c1 ;
u4 = v4 = −u2 , c4 = c2 .
De este modo,
s
4 4
1 XX u2i + 2ui + vj2 + 2vj + 2
I≈ ci cj ≈ 0.267770529696778.
4 i=1 j=1 34 − u2i − 2ui − vj2 − 2vj
Para finalizar...
i Lecciones magistrales
L Material complementario: A fondo
Bibliografı́a recomendada
...Y por supuesto:
TEST DE APRENDIZAJE!!


Mna1 T05-1

Cargado por

Copyright:

Formatos disponibles

Mna1 T05-1

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Mna1 T05-1

Cargado por

Copyright:

Formatos disponibles

Tema 5

Dra. Neus Garrido

Máster en Ingenierı́a Matemática y Computación

Métodos Numéricos Aplicados I PER 3362 1 / 63

2 Fórmulas cerradas de Newton-Cotes

3 Fórmulas abiertas de Newton-Cotes

Métodos Numéricos Aplicados I PER 3362 2 / 63

Métodos Numéricos Aplicados I PER 3362 3 / 63

Función de distribución normal en un proceso de fabricación

Métodos Numéricos Aplicados I PER 3362 4 / 63

Ejemplo 1. Fuerza total ejercida por el mástil de un velero

donde z es la distancia vertical a la cubierta.

Métodos cerrados Métodos abiertos

Métodos Numéricos Aplicados I PER 3362 7 / 63

Polinomio de interpolación de Lagrange de grado n

Métodos Numéricos Aplicados I PER 3362 7 / 63

Métodos Numéricos Aplicados I PER 3362 8 / 63

Métodos Numéricos Aplicados I PER 3362 8 / 63

Métodos Numéricos Aplicados I PER 3362 8 / 63

denominándose error de cuadratura al término

Métodos Numéricos Aplicados I PER 3362 8 / 63

denominándose error de cuadratura al término

ú nodos equiespaciados fórmulas de Newton-Cotes

Fórmulas cerradas de Newton-Cotes

Métodos Numéricos Aplicados I PER 3362 9 / 63

Teorema 1 (Error de las fórmulas

para valores impares de n y si f ∈ C n+1 [a, b] se cumple que

Métodos Numéricos Aplicados I PER 3362 10 / 63

Fórmulas cerradas de Newton-Cotes

2 Fórmulas cerradas de Newton-Cotes

3 Fórmulas abiertas de Newton-Cotes

Métodos Numéricos Aplicados I PER 3362 12 / 63

Reducción del error

Método de trapecios compuesto

Métodos Numéricos Aplicados I PER 3362 15 / 63

2 Fórmulas cerradas de Newton-Cotes

3 Fórmulas abiertas de Newton-Cotes

Métodos Numéricos Aplicados I PER 3362 16 / 63

Métodos Numéricos Aplicados I PER 3362 17 / 63

Método de Simpson compuesto

Métodos Numéricos Aplicados I PER 3362 19 / 63

>> f=@(x) sin(x).*exp(-x);

Para un mismo método, el error se reduce cuando disminuye el tamaño de h

Métodos Numéricos Aplicados I PER 3362 20 / 63

Fórmulas abiertas de Newton-Cotes

Métodos Numéricos Aplicados I PER 3362 21 / 63

Nodos equiespaciados: x ∈ [a, b]

Métodos Numéricos Aplicados I PER 3362 22 / 63

Teorema 2 (Error de las fórmulas abiertas de Newton-Cotes)

para valores impares de n y si f ∈ C n+1 [a, b] se cumple que

Métodos Numéricos Aplicados I PER 3362 23 / 63

Fórmulas abiertas de Newton-Cotes

Métodos Numéricos Aplicados I PER 3362 24 / 63

2 Fórmulas cerradas de Newton-Cotes

3 Fórmulas abiertas de Newton-Cotes

Métodos Numéricos Aplicados I PER 3362 25 / 63

Método de Punto Medio

Métodos Numéricos Aplicados I PER 3362 26 / 63