Capı́tulo 4

Integración numérica.

En la ciencia y la ingenierı́a existen muchos problemas que se resuelven utilizando aplicaciones

de la integral. Entre estas aplicaciones se pueden mencionar, el cálculo del área de una superficie
con dimensiones irregulares, la cantidad de calor que fluye a través de un material, el movimiento
de una partı́cula a través de un fluido, entre otras.
Los métodos de integración numérica nos permiten integrar funciones que están definidas
analı́ticamante o que se presentan en forma tabular. El principio básico de los métodos de inte-
gración numérica consiste en ajustar un polinomio a puntos de datos funcionales y a continuación
integrarlo. Ası́, podemos derivar muchos métodos de integración distintos modificando la distribu-
ción de las abscisas de los puntos de datos.
Si el grado del polinomio es demasiado elevado, los errores debido al redondeo y a las irregula-
ridades locales pueden ocasionar problemas, por tanto, se usan polinomios de grado bajos.
En las siguientes secciones se abordarán las técnicas de integración numérica que utilizan poli-
nomios de primer (regla del trapecio), segundo (regla de simpson 13 ) y tercer (regla de simpson 83 )
grado.

4.1. Regla del trapecio

La regla del trapecio, también llamada trapezoidal, es la primera de las fórmulas de integración
cerrada de Newton-Cotes. Se basa en la estrategia de reemplazar una función complicada o datos
tabulados con una función que sea fácil de integrar.
La regla trapezoidal corresponde al caso donde la función original se sustituye por un polinomio
de primer orden (de grado 1).
Una aproximación de este tipo es usar una lı́nea recta la cual está dada por

f (b) − f (a)
f (x) ≈ f (a) + (x − a)
b−a

de donde

b Z b 
f (b) − f (a)
Z
f (x)dx ≈ f (a) + (x − a) dx
a a b−a
1
2 CAPÍTULO 4. INTEGRACIÓN NUMÉRICA

desarrollando el lado derecho:

Z b  Z b Z b
f (b) − f (a) f (b) − f (a)
f (a) + (x − a) dx = f (a)dx + (x − a)dx
a b−a a a b−a

b !
f (b) − f (a) x2 
= f (a)(b − a) + − a(b − a)
b−a 2 a

b 2 − a2
 
f (b) − f (a)
= f (a)(b − a) + − a(b − a)
b−a 2

b2 − a2 − 2a(b − a)
 
f (b) − f (a)
= f (a)(b − a) +
b−a 2

(b − a)2
 
f (b) − f (a)
= f (a)(b − a) +
b−a 2

 
2f (a) f (b) − f (a)
= (b − a) +
2 2

 
f (b) + f (a)
= (b − a)
2

El valor de la integral es

Z b  
f (b) + f (a)
f (x)dx ≈ (b − a)
a 2

Geométricamente, la regla trapezoidal es equivalente a aproximar el área del trapezoide bajo la

lı́nea recta que conecta a f (a) y f (b) como se muestra en la Fig. 4.1
4.1. REGLA DEL TRAPECIO 3

Figura 4.1: Aproximación de la integral con la regla del trapecio

Ejemplo 1 Integre numéricamente el siguiente polinomio en el intervalo de a = 0 a b = 0.8

f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5

Solución. Usando la regla del trapecio tenemos que:

Z 0.8  
f (0.8) + f (0)
f (x)dx ≈ (0.8)
0 2
 
0.232 + 0.2
= (0.8)
2
= 0.1728

Ya que la f (x) dada es un polinomio, podemos integrarla y obtener el valor real de la integral
Z
25 2 200 3 675 4 200 6
f (x)dx = 0.2x + x − x + x − 180x5 + x
2 3 4 3

de donde
Z 0.8
0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5 dx = 1.64053333
 
0

por lo cual en la estimación se cometió un error de |0.1728 − 1.64053333| = 1.46773333. La gráfica

de esta aproximación se muestra en la Fig. 4.2.
4 CAPÍTULO 4. INTEGRACIÓN NUMÉRICA

Figura 4.2: Aproximación de la integral de f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5 con
la regla del trapecio

4.1.1. Aplicación múltiple de la regla trapezoidal

Una forma de mejorar la exactitud de la regla trapezoidal es dividir el intervalo de integración
desde a hasta b en un número de segmentos y aplicar el método a cada uno de ellos. Las áreas de
los segmentos individuales se suman para calcular el área total del intervalo. La fórmula resultante
es llamada fórmula de integración de múltiple aplicación o fórmula de integración compuesta.
Supongamos que dividimos el segmento [a, b] en n partes iguales, de modo que tenemos n+1 puntos
base igualmente espaciados (x0 , x1 , . . . , xn ), la anchura de cada segmento de la forma (xi , xi+1 ) es

b−a
h=
n
La integral total se expresa como:
Z x1 Z x2 Z xn
I= f (x)dx + f (x)dx + · · · + f (x)dx
x0 x1 xn−1

Al aplicar la regla trapezoidal en cada uno de estos segmentos tenemos

f (x1 ) + f (x0 ) f (x2 ) + f (x1 ) f (xn ) + f (xn−1 )

I≈h +h + ··· + h
2 2 2
o bien,
h
I≈ [f (x1 ) + f (x0 ) + f (x2 ) + f (x1 ) + · · · + f (xn ) + f (xn−1 )]
2
4.1. REGLA DEL TRAPECIO 5

reordenando
h
I≈ [f (x0 ) + f (x1 ) + f (x1 ) + f (x2 ) + · · · + f (xn ) + f (xn−1 )]
2

usando los signos de sumatoria

" n−1
#
h X
I≈ f (x0 ) + f (xn ) + 2 f (xi )
2 i=1

Ejemplo 2 Integre numéricamente

f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5

desde a = 0 a b = 0.8, divida el intervalo en 5 partes iguales.

Solución. Sabemos que

" n−1
#
Z 0.8
h X
f (x)dx ≈ f (x0 ) + f (xn ) + 2 f (xi )
0 2 i=1

0.8−0
para este caso n = 5, h = 5
= 0.16, x0 = 0 y x5 = 0.8 por lo tanto, la igualdad anterior se
reescribe como:
" 4
#
Z 0.8 X
f (x)dx ≈ 0.08 f (0) + f (0.8) + 2 f (xi )
0 i=1

donde x1 = 0.16, x2 = 0.32, x3 = 0.48, x4 = 0.64

Z 0.8
f (x)dx ≈ 0.08 [0.2 + 0.232 + 2(1.296919 + 1.743393 + 3.186014 + 3.181928)]
0
= 1.5398806

En esta ocasión el error cometido fue de |1.64053333 − 1.5398806| = 0.10065269 como se puede
observar en la Fig. 4.3.
6 CAPÍTULO 4. INTEGRACIÓN NUMÉRICA

Figura 4.3: Aproximación de la integral de f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5 con
la regla compuesta del trapecio (5 intervalos)

4.2. Regla de Simpson

Además de aplicar la regla trapezoidal con segmentación más fina, otra forma de obtener una
estimación más exacta de la integral es con el uso de polinomios de orden superior para conectar
los puntos. Las fórmulas que resultan al usar polinomios de grado dos y tres son conocidas como
Reglas de Simpson.

1
4.2.1. Regla de Simpson de 3
La regla de Simpson de 13 resulta cuando una interpolación polinomial de segundo orden es
usada para aproximar a la función en cuestión. Si a = x0 , x1 = b+a
2
y b = x2 ; usando el polinomio
de lagrange se tiene:
(x − x1 )(x − x2 ) (x − x0 )(x − x2 ) (x − x0 )(x − x1 )
f (x) = f (x0 ) + f (x1 ) + f (x2 )
(x0 − x1 )(x0 − x2 ) (x1 − x0 )(x1 − x2 ) (x2 − x0 )(x2 − x1 )
Después de integrar y realizar algunos arreglos algebraicos, se obtienen la siguiente fórmula:
h
I≈ [f (x0 ) + 4f (x1 ) + f (x2 )]
3
b−a
donde h = 2

Ejemplo 3 Integre numéricamente

f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5
4.2. REGLA DE SIMPSON 7

desde a = 0 a b = 0.8

1
Solución. Ahora calcularemos el valor de esta integral usando la regla de simpson de 3
como se
muestra a continuación:
Z 0.8
h
f (x)dx = [f (x0 ) + 4f (x1 ) + f (x2 )]
0 3
2
= [0.2 + 4(2.456) + 0.232]
15
= 1.367466

En esta ocasión el error cometido es de |1.64053333 − 1.367466| = 0.273066 que como puede
observarse es mucho menor que el cometido utilizando la regla simple del trapecio (integrando sólo
con un trapecio). La Fig. 4.4 muestra la gráfica de esta aproximación.

Figura 4.4: Aproximación de la integral de f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5 con
la regla de Simpson de 31

1
4.2.2. Regla de Simpson de 3 compuesta
Al igual que la regla del trapecio, la regla de simpson se puede aplicar en forma compuesta.
Aunque de manera general, es posible aplicar la regla se simpson para segmentos de diferente
tamaño, los cálculos son más sencillos si dividimos el intervalo en un número de segmentos de
igual anchura. Sea h = b−an
, la anchura de los intervalos, entonces la integral se puede expresar de
8 CAPÍTULO 4. INTEGRACIÓN NUMÉRICA

la siguiente manera:

Z x2 Z x4 Z xn
I= f (x)dx + f (x)dx + · · · + f (x)dx
x0 x2 xn−2

Al sustituir la regla de simpson en cada integral de la ecuación anterior, se tiene:

f (x0 ) + 4f (x1 ) + f (x2 ) f (x2 ) + 4f (x3 ) + f (x4 ) f (x − 2) + 4f (xn−1 ) + f (xn )

I ≈ 2h + 2h + · · · + 2h
6 6 6

Observe que el extremo final de cada suma es el extremo inicial de la suma siguiente; aprovechando
este patrón y recordando que h = b−a
n
, podemos reescribir la expresión anterior como:

Pn−1 Pn−2
f (x0 ) + f (xn ) + 4 i=1,3,5 f (xi ) + 2 i=2,4,6 f (xi )
I = (b − a)
3n

Este método está limitado a una situación donde hay un número par de subintervalos
(número impar de puntos) y a longitudes de subintervalos igualmente espaciados.

Ejemplo 4 Integre numéricamente

f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5

desde a = 0 a b = 0.8, divida el intervalo en 4 partes iguales.

Solución. Calculamos el valor de h = 0.84

= 0.2 de donde los valores de los puntos xi son x0 = 0,
x1 = 0.2, x2 = 0.4, x3 = 0.6, x4 = 0.8. El valor de la integral está dado por:

0.2 + 0.232 + 4(1.288 + 3.464) + 2(2.456)

I = (0.8)
12
= 1.623466

En esta ocasión hemos cometido un error de |1.64053333 − 1.623466| = 0.017066 lo cual mejora
el cálculo hecho con la regla del trapecio compuesta. Observe que se empleó un número menor de
puntos. En la Fig. 4.5, se puede observar la diferencia entre el valor de la integral real y el valor
de la integral aproximado.
4.2. REGLA DE SIMPSON 9

Figura 4.5: Aproximación de la integral de f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5 con
la regla de Simpson de 31 Compuesta

3
4.2.3. Regla de Simpson de 8
Esta regla difiere de la anterior en que en lugar de usar tres puntos para el polinomio de
Lagrange, se usan cuatro, es decir se usa un polinomio de grado tres para hacer los ajustes entre
los puntos. Repitiendo el proceso expuesto con anterioridad se puede mostrar que la aproximación
para este método está dada por:
3h
I= [f (x0 ) + 3f (x1 ) + 3f (x2 ) + f (x3 )] (4.1)
8
b−a
donde h = 3
.

Ejercicio 5 Integre numéricamente

f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5

desde a = 0 a b = 0.8
0.8
Solución. Calculemos el valor de h, h = 3
= 0.26, por lo cual, x0 = 0, x1 = 0.26, x2 = 0.53, x3 =
0.8 de donde
3(0.26)
I= [0.2 + 3(1.432724) + 3(3.487176) + 0.232]
8
= 1.51917
10 CAPÍTULO 4. INTEGRACIÓN NUMÉRICA

En esta aproximación el error cometido fue de |1.64053333 − 1.51917| = 0.121363. En la Fig.

4.6, se puede observar la diferencia entre el valor de la integral real y el valor de la integral
aproximado.

Figura 4.6: Aproximación de la integral de f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5 con
la regla de Simpson de 31 Compuesta

3
4.2.4. Regla de Simpson de 8 compuesta
De la misma forma en la que se aplicaron las reglas compuestas del trapecio y Simpson 13 , es
necesario dividir el intervalo de integración en n subintervalos. Estos subintervalos pueden tener un
tamaño arbitrario, sin embargo, en este curso analizaremos el caso donde los subintervalos tienen
el mismo tamaño h, donde h = b−a n
.
Debido a que se necesitan cuatro puntos para construir un polinomio de grado 3, el método de
Simpson 38 , se aplica a tres subintervalos adyacentes a la vez (los primeros tres, el cuarto, quinto y
sexto y ası́ sucesivamente). Por lo que, el intervalo completo de integración debe de dividirse en un
número de intervalos, que sea de al menos 6 y divisible entre 3. De la misma forma que se obtuvo
una ecuación general para las reglas compuestas del trapecio y Simpson 13 , se obtiene la ecuación
4.2, la cual nos ayuda a obtener mejores aproximaciones de la integral a evaluar.
" n−2 n−3
#
3h X X
I≈ f (x0 ) + f (xn ) + 3 [f (xi ) + f (xi+1 )] + 2 f (xi ) (4.2)
8 i=1,4,7 i=3,6,9

Ejemplo 6 Integre numéricamente

f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5
desde a = 0 a b = 0.8, divida el intervalo en 6 partes iguales.
4.3. EJERCICIOS 11

Solución. Calculamos el valor de h = 0.8

6
= 0.13 de donde los valores de los puntos xi son x0 = 0,
x1 = 0.13, x2 = 0.26, x3 = 0.4, x4 = 0.53, x5 = 0.66, x6 = 0.8. El valor de la integral está dado
por:

I = 0.05[0.2 + 0.232 + 3(1.3101 + 1.4327 + 3.4871 + 2.8749) + 2(2.456)]

= 1.632948

En esta ocasión hemos cometido un error de |1.64053333 − 1.632948| = 0.00758533 lo cual

mejora los cálculos hechos con los métodos revisados en este capı́tulo. En la Fig. 4.7, se puede
observar la diferencia entre el valor de la integral real y el valor de la integral aproximado.

Figura 4.7: Aproximación de la integral de f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5 con
la regla de Simpson de 83 Compuesta

4.3. Ejercicios
Ejercicio 7 Aproxime el valor de las siguientes integrales, luego compare con el valor real
Z 3
1 − e−x dx


1.
0

Z 4
1 − x − 4x3 + x5 dx


2.
−2

Z π
2
3. (8 + 4 sin (x)) dx
0
12 CAPÍTULO 4. INTEGRACIÓN NUMÉRICA

Ejercicio 8 Evalúe la integral de los siguientes datos tabulados:

x 0 0.1 0.2 0.3 0.4 0.5

f (x) 1 7 4 3 5 2

1. Usando la regla trapezoidal

2. Usando la regla de simpson

Ejercicio 9 Evalue la integral de los siguientes datos tabulados:

x -3 -1 1 3 5 7 9 11
f (x) 1 -4 -9 2 4 2 6 -3

1. Usando la regla trapezoidal

2. Usando la regla de simpson

Ejercicio 10 Evalúe las siguientes integrales por la regla del trapecio extendida con n = 2, 4, 8 y
16 intervalos:
Z π
4
1. tan (x) dx
0
Z 1
2. ex dx
0
Z 1
1
3. dx
0 2+x

Ejercicio 11 Con la tabla de la función que se da más abajo, evalúe

Z 0.8
f (x) dx
0

por la regla trapezoidal extendida con h = 0.4, 0.2 y 0.1

x 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

f (x) 0 2.122 3.0244 3.2568 3.1399 2.8579 2.514 2.1639 1.8358