Integración Numérica I

Coordinación de Métodos Numéricos, DMCC

Reglas de integración numérica.

DMCC, Facultad de Ciencia, USACH

Para aproximar una integral de la forma
Z b
f (x) dx,
a

puede aproximarse el integrando f por el polinomio de interpolación Pn (x) que interpola a f en (n + 1) nodos
x0 , . . . , xn e integrar:
Z b Z b
f (x) dx ≈ Pn (x) dx =: In (f ).
a a

El error que se comete al aproximar la integral por In (f ) es

Z b Z b Z b
Rn (f ) := f (x) dx − Pn (x) dx = E (x) dx,
a a a

donde E (x) := f (x) − Pn (x) es el error de interpolación

(x − x0 ) · · · (x − xn ) (n+1)
E (x) := f (x) − Pn (x) = f (ξx ),
(n + 1)!

para algún ξx ∈ (a, b).

2/19
Si se usa la fórmula de Lagrange:
n
X
Pn (x) = f (xi )`i (x),
i=0

para calcular explı́citamente In (f ), se tiene

Z b n Z b
X  n
X
In (f ) = Pn (x) dx = `i (x) dx f (xi ) = Ai f (xi ),
a i=0 a i=0

con Z b
Ai = `i (x) dx, i = 0, . . . , n.
a

Llamaremos:
n
X
Ai f (xi ): regla de integración numérica o regla de cuadratura.
i=0

xi : nodos de la regla de integración,

Ai : coeficientes o pesos de la regla de integración,

3/19
Fórmulas de Newton Cotes
a+b
1. Regla del punto medio (elemental): n = 0, con x0 = .
2
En este caso, A0 = b − a y se obtiene:

a+b
I0 (f ) = (b − a) f (b
x ), con xb := .
2

Si f ∈ C 2 ([a, b]), el error de integración está dado por:

Z b Z b 00
f (ξx ) 2
R0 (f ) := f (x) dx − I0 (f ) = (x − xb) dx.
a a 2
00
Llamando M2 := max |f (x)| y evaluando la integral restante, se obtiene:
x∈[a,b]

M2 3
|R0 (f )| ≤ (b − a) .
24

Regla del punto medio (compuesta). El intervalo [a, b] se divide en n subintervalos iguales:

b−a
xi := a + ih, i = 0, 1, . . . , n, con h := .
n
La regla del punto medio compuesta se obtiene aplicando la regla del punto medio elemental en cada subintervalo
[xi−1 , xi ]:

Z b n Z x
X i
f (x) dx = f (x) dx
a i=1 xi−1
n
X xi−1 + xi
≈ h f (b
xi ), con xbi := .
i=1
2 4/19
La regla del punto medio (compuesta) consiste en aproximar la integral por:

n
X xi−1 + xi
IM (f ) := h f (b
xi ), con xbi := .
i=1
2

Se puede demostrar que el error viene dado por

M2 2 00
|RM (f )| ≤ (b − a)h , con M2 := max |f (x)|.
24 x∈[a,b]

En consecuencia, esta regla es exacta para integrar polinomios de grado menor o igual que uno .

Definición: El grado de exactitud de una fórmula de cuadratura es el máximo entero k ≥ 0 para el cual la integral
aproximada de cualquier polinomioP de grado k es igual a la integral exacta.
n
La fórmula de cuadratura In (f ) = i=0 Ai f (xi ) tiene grado de exactitud k, si

Z b
Pk (x) dx = In (Pk )
a

para todo Pk y no para Pk+1 .

5/19
Fórmulas de Newton Cotes
2. Regla del trapecio (elemental): n = 1, con x0 = a y x1 = b.
b−a
En este caso, A0 = A1 = y se obtiene:
2
b−a
I1 (f ) = [f (a) + f (b)] .
2

Si f ∈ C 2 ([a, b]), para el error de integración se tiene:

f 00 (ξx ) (b − a)3 00
Z b Z b
R1 (f ) := f (x) dx − I1 (f ) = (x − x0 )(x − x1 ) dx = − f (ξ),
a a 2 12

para algún ξ ∈ (a, b).

Regla de los trapecios (compuesta).

El intervalo [a, b] se divide en n subintervalos iguales:

b−a
xi := a + ih, i = 0, 1, . . . , n, con h := .
n
La regla de los trapecios compuesta se obtiene aplicando la regla del trapecio elemental en cada subintervalo
[xi−1 , xi ]:

Z b n Z x n
X i X h  
f (x) dx = f (x) dx ≈ f (xi−1 ) + f (xi )
a i=1 xi−1 i=1
2
 
n−1
f (a) + f (b) X
= h + f (xi ) .
2 i=1 6/19
La regla de los trapecios (compuesta) consiste en aproximar la integral por:

 
n−1
f (a) + f (b) X
IT (f ) := h  + f (xi ) .
2 i=1

Se puede demostrar que el error es del orden h2

M2 2 00
|RT (f )| ≤ (b − a)h , con M2 := max |f (x)|.
12 x∈[a,b]

Definición: Se dice que una regla es de orden hp , y se escribe O(hp ), cuando el error satisface |R| ≤ Chp , para
alguna constante C > 0 independiente de h.

Esta regla es exacta para integrar polinomios de grado menor o igual que uno.

7/19
Fórmulas de Newton Cotes

a+b
3. Regla de Simpson (elemental): n = 2, con x0 = a, x1 = y x2 = b.
2
b−a b−a
En este caso, A0 = A2 = y A1 = 4 , y se obtiene:
6 6

b−a a+b
I2 (f ) := [f (a) + 4f (b
x ) + f (b)] , con xb := .
6 2

Si f ∈ C 4 ([a, b]), para el error de integración se tiene:

5
b−a
Z b
1

(4)
R2 (f ) := f (x) dx − I2 (f ) = − f (ξ),
a 90 2

para algún ξ ∈ (a, b).

8/19
Fórmulas de Newton Cotes

Regla de Simpson (compuesta)

El intervalo [a, b] se divide en 2n subintervalos iguales:

b−a
xi := a + ih, i = 0, 1, . . . , 2n, con h := .
2n

La regla de de Simpson compuesta se obtiene aplicando la regla de Simpson elemental en cada subintervalo
[x2i−2 , x2i ]:

Z b n Z x
X 2i
f (x) dx = f (x) dx,
a i=1 x2i−2
n
X 2h  
≈ f (x2i−2 ) + 4f (x2i−1 ) + f (x2i )
i=1
6
 
n n−1
h X X
= f (a) + f (b) + 4 f (x2i−1 ) + 2 f (x2i ) .
3 i=1 i=1

9/19
La regla de Simpson (compuesta) consiste en aproximar la integral por:

 
n n−1
h X X
IS (f ) := f (a) + f (b) + 4 f (x2i−1 ) + 2 f (x2i ) .
3 i=1 i=1

El error es de orden h4

M2 4 (4)
|RS (f )| ≤= (b − a)h con M2 := max |f (x)|.
180 x∈[a,b]

Esta regla es exacta para integrar polinomios de grado menor o igual que tres.

10/19
Observaciones.

1 Es costumbre referirse a las reglas compuestas IM (f ), IT (f ) e IS (f ) simplemente como regla del punto
medio, regla de los trapecios y regla de Simpson.

2 Existen versiones de las reglas para nodos no equiespaciados.

3 Las reglas del punto medio y de los trapecios son O(h2 ), mientras que la regla de Simpson es O(h4 ).

11/19
Otras Fórmulas de Newton-Cotes (cerradas) en [b − a]

n
b−a X
I (f ) := σi f (xi ).
ns i=0

n σi ns Error Nombre
1 1 1 2 h3 12
1 f (2) (ξ) Regla del Trapecio
2 1 4 1 6 h5 90
1 f (4) (ξ) Regla de Simpson
3 1 3 3 1 8 h5 80
3 f (4) (ξ) Regla de los 3/8
4 7 32 12 32 7 90 h7 945
8 f (6) (ξ) Regla de Milne
5 19 75 50 75 19 288 h7 12096
275 f (6) (ξ) -

12/19
 Z 1
Ejercicio 1: Sea f (x) = 1 + e −x sin(4x). Aproxime el valor de la integral f (x) dx, usando
0

1 la regla compuesta del punto medio

2 la regla compuesta del trapecio,

3 la regla compuesta de Simpson simple.

Compare los resultados obtenidos en cada aproximación con el resultado exacto de la integral

21e − 4 cos(4) − sin(4)

Z 1
f (x) dx = ≈ 1.308250604642669 . . . .
0 17e

Solución: 1. Regla del punto medio con h = 0.5, a = 0, b = 1 y n = 2:

2
xi−1 + xi
X  
IM (f ) = h f , xi = ih
i=1
2
−0.25 −0.75
= 0.5(4 + e sin(4 · 0.25) + e sin(4 · 0.75))
= 1.360999316815140.

13/19
2. Regla del trapecio con h = 0.25, a = 0, b = 1 y n = 4:

 
3
f (0) + f (1) X
IT (f ) = h + f (xi ) , xi = ih
2 i=1

(1) + (1 + e −1 sin(4))
"
−0.25 −0.5
= 0.25 + (1 + e sin(4 · 0.25)) + (1 + e sin(4 · 0.5))
2
i
−0.75
+(1 + e sin(4 · 0.75))
= 1.283577340568086.

3. Regla de Simpson con h = b−a

2n
= 0.25, a = 0, b = 1 y n = 2:

 
2 1
h X X
IS (f ) = f (0) + f (1) + 4 f (x2i−1 ) + 2 f (x2i ) , xi = ih
3 i=1 i=1
h
= [f (0) + f (1) + 4f (x1 ) + 4f (x3 ) + 2f (x2 )]
3
= 1.309384665983771

14/19
Calculemos los errores:

Regla del punto medio

EM = |1.308250604642669 − 1.360999316815140| = 0.102095240321637.

Regla del trapecio

ET = |1.308250604642669 − 1.283577340568086| = 0.012441898556515.

Regla de Simpson

ES = |1.308250604642669 − 1.309384665983771| = 0.001134061341102.

Se puede concluir que con la regla de Simpson se obtiene una mejor aproximación.

15/19
Ejercicio 2: Obtenga el grado de precisión de la fórmula de la cuadratura

Z 1 √ ! √ !
3 3
f (x)dx = f − +f .
−1 3 3

Solución: Consideremos la base {1, x, x 2 , . . . }, calcularemos el valor exacto y el valor usando la fórmula de
cuadratura I (f ).

Para f (x) = 1:
Z 1 Z 1
dx = 2, I (1) = 2 =⇒ dx = I (1).
−1 −1

Para f (x) = x:

√ √
x 2 1
Z 1 Z 1
3 3
x dx =  = 0, I (x) = − + =0 =⇒ x dx = I (x).
−1 2 −1 3 3 −1

Para f (x) = x 2 :

√ !2 √ !2
x 3 1
Z 1 Z 1
2 2 2 3 3 2 2 2
x dx =  = , I (x ) = − + = =⇒ x dx = I (x ).
−1 3 −1 3 3 3 3 −1

Para f (x) = x 3 :

√ !3 √ !3
x 4 1
Z 1 Z 1
3 3 3 3 3 3
x dx =  = 0, I (x ) = − + =0 =⇒ x dx = I (x ).
−1 4 −1 3 3 −1

16/19
 Para f (x) = x 4 :

√ !4 √ !4
x 5 1
Z 1 Z 1
4 2 4 3 3 2 4 4
x dx =  = , I (x ) = − + = =⇒ x dx 6= I (x ).
−1 5 −1 5 3 3 9 −1

La fórmula de cuadratura √ ! √ !
Z 1
3 3
f (x)dx = f − +f .
−1 3 3

es exacta para polinomios de grado menor o igual a 3.

17/19
Ejercicio 3: Considere la fórmula de cuadratura
Z b
b
   
f (x) dx = b A1 f (0) + A2 f + A3 f (b) .
0 3

1 Determine los valores de las constantes A1 , A2 y A3 de modo que la fórmula sea exacta para polinomios de
mayor grado posible. Indique dicho grado.

2 Use la fórmula para aproximar Z 0.5

2 −t 2
I = √ e dt.
π 0

Solución: Para calcular las constantes necesitamos tres ecuaciones, para eso evaluamos en la base {1, x, x 2 , . . . },

Para f (x) = 1:
Z b
dx = b = b(A1 + A2 + A3 ) =⇒ A1 + A2 + A3 = 1,
0

Para f (x) = x:

b2
Z b
b A2 1
x dx = = b(A1 · 0 + A2 + A3 b) =⇒ + A3 = ,
0 2 3 3 2

Para f (x) = x 2 :

b3 b2
Z b
2 2 A2 1
x dx = = b(A1 · 0 + A2 + A3 b ) =⇒ + A3 = ,
0 3 9 9 3

Resolviendo el sistema de ecuaciones se obtiene

3 1
A1 = 0, A2 = A3 = .
4 4 18/19
La fórmula de cuadratura propuesta es

Z b
3 b 1
   
f (x) dx = b f + f (b) .
0 4 3 4

Esta fórmula es exacta para polinómios de grados menor o igual a 2.

2. Usemos la fórmula para aproximar Z 0.5

2 −t 2
I = √ e dt,
π 0

−t 2
donde b = 0.5 y f (t) = e .

Z 0.5
2 2 3 0.5 1
 
−t 2
e dt = 0.5 f( )+ f (0.5)
π 0 π 4 3 4
1 3 −( 0.5 )2 1 −(0.5)2
 
= e 3 + e = 0.521397808
π 4 4

Ası́ el valor de la integral I = 0.521397808.

19/19