UNIVERSIDAD DEL MAGDALENA

FACULTAD DE INGENIERÍA
ANÁLISIS NUMÉRICO
NOCIONES DE INTEGRACIÓN NUMÉRICA
Mg. Alvaro Espinosa Pérez

INTRODUCCIÓN
La regla de Barrow permite integrar de modo exacto las funciones cuya primitiva existe y puede calcularse
(aunque sea a trozos). Muchas veces esto no es posible, bien porque la función solo se conoce en una serie
de puntos, o porque no hay modo de expresar la primitiva de una forma elemental. En estos casos, que son
muy frecuentes, por ejemplo:
3
√𝑥 √1 − 𝑥, √1 − 𝑥 3 , 𝑠𝑒𝑛𝑥 2
Los métodos de integración numérica que vamos a estudiar son de tipo interpolatorio. Esto quiere decir que la
integral de una función 𝑓 en un intervalo [𝑎, 𝑏] se aproximará de la siguiente forma:
𝑏 𝑏
∫ 𝑓(𝑥 )𝑑𝑥 ≅ ∫ 𝑆(𝑥 )𝑑𝑥
𝑎 𝑎
donde 𝑆(𝑥) es una función que interpola a 𝑓(𝑥) en determinados puntos de [𝑎, 𝑏] y que es fácil de integrar en
este intervalo. En estos casos cuando son útiles los métodos aproximados de integración, veamos algunos:

Las fórmulas de NEWTON-COTES son los tipos de integración numérica más comunes. Se basan en la
estrategia de reemplazar una función complicada o datos tabulados por un polinomio de aproximación que es
fácil de integrar. Veamos algunas:

MÉTODO DEL PUNTO MEDIO

𝑏
Sea 𝑓 una funcion continua en [𝑎, 𝑏]. La regla del punto medio para aproximar ∫𝑎 𝑓(𝑥)𝑑𝑥 viene dada por:
𝑏 𝑛
𝑏−𝑎 𝑏−𝑎
∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑓(𝑥̅𝑖 ) = [𝑓(̅̅̅
𝑥1 ) + 𝑓 (̅̅̅
𝑥2 ) + ⋯ + 𝑓(̅̅̅
𝑥𝑛 )]
𝑎 𝑛 𝑛
𝑖=1
1
Con 𝑥̅𝑖 = (𝑥𝑖−1 + 𝑥𝑖 )
2

Ejemplo: Consideremos la siguiente integral definida:

1
1
∫ cos(2𝑥) 𝑑𝑥 = 𝑠𝑒𝑛(2) = 0,4546 …
0 2
1
Si uitilizamos la fórmula de puintos medios con 5 sub-intervalos de igual longitud, esto es: ℎ = = 0,2 y
5
{𝑥𝑖 }6𝑖=1 = {0, 0.2, 0.4, 0.6, 0.8, 1}, tenemos:
1 5
𝑥𝑖 + 𝑥𝑖+1
∫ cos(2𝑥 ) 𝑑𝑥 = ℎ ∑ 𝑐𝑜𝑠 (2 )
0 2
1
1
∫ cos(2𝑥) 𝑑𝑥 = 0,2(cos(0.2) + cos(0.6) + cos(1) + cos(1.4) + cos(1.8)) = 0.4577 …
0

METODO DEL TRAPECIOS

Una manera de aproximar una integral definida es emplear n trapecios. Este método supone que f es una
𝑏
función continua y positiva en [𝑎, 𝑏]. De manera que la integral definida ∫𝑎 𝑓(𝑥)𝑑𝑥 representa el área de la
región limitad por la gráfica de f y las rectas 𝑥 = 𝑎 y 𝑥 = 𝑏.
 (𝑏−𝑎)
Haremos una partición n del intervalo [𝑎, 𝑏] en 𝒏 subintervalos con ancho ∆𝑥 = de manera que:
2
𝑎 = 𝑥0 < 𝑥1 < ⋯ < 𝑥𝑛 = 𝑏

El área de i-ésimo trapecio sería entonces:

𝑓 (𝑥𝑖−1 ) + 𝑓(𝑥𝑖 ) 𝑏 − 𝑎
𝐴𝑖 = [ ]( )
2 𝑛

Luego la suma de las áreas de 𝑛 trapecios es:

𝑏 − 𝑎 𝑓(𝑥0 ) + 𝑓(𝑥1 ) 𝑓 (𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )
𝐴=( )[ +⋯+ ]
𝑛 2 2
𝑏−𝑎
=( ) [𝑓(𝑥0 ) + 𝑓(𝑥1 ) + 𝑓(𝑥1 ) + 𝑓(𝑥2 ) + ⋯ + 𝑓(𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )]
2𝑛
𝑏−𝑎
=( ) [𝑓 (𝑥0 ) + 2𝑓 (𝑥1 ) + 2𝑓 (𝑥2 ) + ⋯ + 2𝑓 (𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )]
2𝑛
Si tomamos el límite cuando 𝑛 → ∞ se tiene que:

𝑏−𝑎
lim ( ) [𝑓(𝑥0 ) + 2𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + ⋯ + 2𝑓(𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )]
𝑛→∞ 2𝑛
(𝑓(𝑥𝑎 ) + 𝑓(𝑥𝑏 ))∆𝑥 𝑛
= lim [ + ∑ 𝑓(𝑥𝑖 )∆𝑥 ]
𝑛→∞ 2 1

(𝑓(𝑥𝑎 ) + 𝑓(𝑥𝑏 ))(𝑏 − 𝑎) 𝑛

= lim [ + ∑ 𝑓(𝑥𝑖 )∆𝑥 ]
𝑛→∞ 2𝑛 1

(𝑓(𝑥𝑎 ) + 𝑓(𝑥𝑏 ))(𝑏 − 𝑎) 𝑛

= lim [ ] + lim ∑ 𝑓(𝑥𝑖 )∆𝑥
𝑛→∞ 2𝑛 𝑛→∞ 1
𝑏
= 0 + ∫ 𝑓(𝑥)𝑑𝑥
𝑎

El resultado se resume en el siguiente teorema:

𝑏
TEOREMA: La regla del trapecio para aproximar ∫𝑎 𝑓(𝑥)𝑑𝑥 viene dada por:
𝑏
𝑏−𝑎
∫ 𝑓(𝑥)𝑑𝑥 = ( ) [𝑓(𝑥0 ) + 2𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + ⋯ + 2𝑓(𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )]
𝑎 2𝑛
𝑏
A medida que 𝑛 → ∞, el lado derecho tiende a aproximarse a ∫𝑎 𝑓(𝑥)𝑑𝑥

Ejemplo: Aproximar con cifras decimales la integra definida utilizando la regla del trapecio con n = 6:
3
𝑑𝑥
∫
0 (16 − 𝑥 2 )
Solución: Como [𝑎, 𝑏] = [0, 3] y 𝑛 = 6, tenemos que:

𝑏−𝑎 3−0 3 1
∆𝑥 = = = =
𝑛 6 6 2
 𝑏−𝑎 3−0 3 1
= = =
2𝑛 2(6) 12 4
De ahí que:
3
𝑑𝑥 1
∫ 2 = [𝑓(𝑥0 ) + 2𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + ⋯ + 2𝑓(𝑥5 ) + 𝑓(𝑥6 )]
0 (16 − 𝑥 ) 4
1
Como 𝑓(𝑥) = 2 , nos ayudaremos con una tabla para hallar la sumatoria de
(16−𝑥 )
los términos dentro del corchete:

Entonces:
3
𝑑𝑥 1
∫ ≈ [0,643] ≈ 0,161
0 (16 − 𝑥 2 ) 4

REGLA DE SIMPSON
La regla de los trapecios se obtiene de aproximar la curva por una poligonal, con lo que la región bajo la curva
se aproxima por una unión de trapecios. En consecuencia, el área de dicha región se aproxima por la suma de
las áreas de los trapecios.

Otra posibilidad puede ser agrupar los puntos de la partición de tres en tres y aproximar la
función por segmentos parabólicos (en lugar de rectilíneos). Entonces la integral de la
función se puede aproximar por la integral de los segmentos parabólicos.

Antes de presentar la regla de Simpson daremos a conocer un teorema necesario para evaluar integrales de
polinomio de grado 2 o menor.

TEOREMA: Si 𝑝(𝑥 ) = 𝐴𝑥 2 + 𝐵𝑥 + 𝐶 , entonces:

𝑏
𝑏−𝑎 𝑎+𝑏
∫ 𝑝(𝑥)𝑑𝑥 = ( ) [𝑝(𝑎) + 4𝑝 ( ) + 𝑝(𝑏)]
𝑎 6 2
Veamos:
𝑏 𝑏 𝑏
𝐴𝑥 3 𝐵𝑥 2
∫ 𝑝(𝑥)𝑑𝑥 = ∫ (𝐴𝑥 2 + 𝐵𝑥 + 𝐶)𝑑𝑥 = [ + + 𝐶𝑥]
𝑎 𝑎 3 2 𝑎

𝐴 (𝑏 3 − 𝑎 3 ) 𝐵 (𝑏 2 − 𝑎 2 )
= + + 𝐶 (𝑏 − 𝑎 )
3 2
𝑏−𝑎
=( ) [2𝐴(𝑎2 + 𝑎𝑏 + 𝑏2 ) + 3𝐵(𝑏 + 𝑎) + 6𝐶 ]
6
Si además los términos dentro del corchete los expresamos:
2
𝑏+𝑎 2 𝑏+𝑎
(𝐴𝑎 + 𝐵𝑎 + 𝐶) + 4 [𝐴 ( ) +𝐵( ) + 𝐶] + (𝐴𝑏 2 + 𝐵𝑏 + 𝐶)
2 2
Si además:
𝑎+𝑏 𝑏+𝑎 2 𝑏+𝑎
𝑃(𝑎) = (𝐴𝑎2 + 𝐵𝑎 + 𝐶 ), 4𝑝 ( ) = 4 [𝐴 ( ) +𝐵( ) + 𝐶] 𝑦 𝑝(𝑏) = (𝐴𝑏2 + 𝐵𝑏 + 𝐶 )
2 2 2
Tendríamos entonces:
𝑏
𝑏−𝑎 𝑎+𝑏
∫ 𝑃(𝑥)𝑑𝑥 = ( ) [𝑝(𝑎) + 4𝑝 ( ) + 𝑝(𝑏)]
𝑎 6 2

Para desarrollar la regla de Simpson considere una función f continua en un intervalo [𝑎, 𝑏] y una partición
𝑃2𝑛 = {𝑥0 , 𝑥1 , . . . , 𝑥2𝑛 } equidistante de [𝑎, 𝑏] con un número par de subintervalos. Ahora, en cada subintervalo
doble [𝑥𝑖−2 , 𝑥𝑖 ] puede aproximarse 𝑓 mediante un polinomio 𝑝 de grado menor o igual a 2; esto es: en e
subintervalo [𝑥0 , 𝑥2 ], la fórmula del Simpson de f para en ese subintervalo es:
𝑥2 𝑥2
𝑥2 − 𝑥0 𝑥0 + 𝑥2
∫ 𝑓 (𝑥 )𝑑𝑥 ≈ ∫ 𝑝(𝑥 )𝑑𝑥 = ( ) [𝑝(𝑥0 ) + 4𝑝 ( ) + 𝑝(𝑥2 )]
𝑥0 𝑥0 6 2
𝑏−𝑎
2[ ]
=( 𝑛 ) [𝑝(𝑥 ) + 4𝑝(𝑥 ) + 𝑝(𝑥 )]
0 1 2
6

𝑏−𝑎
=( ) [𝑝(𝑥0 ) + 4𝑝(𝑥1 ) + 𝑝(𝑥2 )]
3𝑛

Si se da el procedimiento para todo el intervalo [a, b], tenemos:

𝑏 𝑛−1 𝑛
𝑏−𝑎
∫ 𝑓(𝑥)𝑑𝑥 ≈ ( ) [𝑓(𝑎) + 2 ∑ 𝑓(𝑥2𝑘 ) + 4 ∑ 𝑓(𝑥2𝑘−1 ) + 𝑓(𝑏)]
𝑎 3𝑛
𝑘=1 𝑘=1

𝑏
TEOREMA: La regla De Simpson para aproximar ∫𝑎 𝑓(𝑥)𝑑𝑥 viene dada por:
𝑏
𝑏−𝑎
∫ 𝑓(𝑥)𝑑𝑥 = ( ) [𝑓(𝑥0 ) + 4𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + 4𝑓(𝑥3 ) + ⋯ + 4𝑓(𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )]
𝑎 3𝑛
𝑏
A medida que 𝑛 → ∞, el lado derecho tiende a aproximarse a ∫𝑎 𝑓(𝑥)𝑑𝑥

Ejemplo: Aproximar con cifras decimales la integra definida utilizando la regla de Simpson con n = 4:
1
𝑑𝑥
∫
0 𝑥+4
Solución: Como [𝑎, 𝑏] = [0, 1] y 𝑛 = 4, tenemos que:

𝑏−𝑎 1−0 1
∆𝑥 = = =
𝑛 4 4
𝑏−𝑎 1−0
=
3𝑛 3(4)
𝑏−𝑎 1
=
3𝑛 12
De ahí que:
1
𝑑𝑥 1
∫ = [𝑓(𝑥0 ) + 4𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + 4𝑓(𝑥3 ) + 𝑓(𝑥4 )]
0 𝑥 + 4 12
1
Como 𝑓(𝑥) = , nos ayudaremos con una tabla para hallar la sumatoria de los
𝑥+4
términos dentro del corchete:

Entonces:
3
𝑑𝑥 1
∫ ≈ [2,678] ≈ 0,223
0 𝑥 + 4 12
REGLA DE SIMPSON 1/3
La regla de Simpson 1/3 resulta cuando un polinomio de interpolación de segundo grado se sustituye en la
ecuación:
𝑏 𝑏
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑃2 (𝑥)𝑑𝑥
𝑎 𝑎
Si se designan 𝑎 y 𝑏 como 𝑥0 y 𝑥2 , y 𝑃2 (𝑥) se representa por un polinomio
de Lagrange de segundo grado, la integral se transforma en:

Después de la integración y de las manipulaciones algebraicas, se obtiene la siguiente fórmula:

𝑏 𝑏 ℎ
∫𝑎 𝑓 (𝑥 )𝑑𝑥 = ∫𝑎 𝑃2 (𝑥 )𝑑𝑥= 3 [𝑓 (𝑥0 ) + 4𝑓 (𝑥1 ) + 𝑓 (𝑥2 )]

𝑏–𝑎
Donde, en este caso, ℎ = . Esta ecuación se conoce como regla de Simpson 1/3, y es la segunda fórmula
2
de integración cerrada de Newton-Cotes. La especificación “1/3” se origina del hecho de que h está dividida
entre 3 en la ecuación.

Se puede demostrar que la aplicación a un solo segmento de la regla de Simpson 1/3 tiene un error de
truncamiento de:
1
𝐸𝑟 = − ℎ5 𝑓 (4) (𝜉)
90
Donde ξ está en algún lugar en el intervalo de [𝑎, 𝑏].

LA REGLA DE SIMPSON 1/3 DE APLICACIÓN MÚLTIPLE

Así como en la regla del trapecio, la regla de Simpson se mejora al dividir el intervalo de integración en varios
𝑏−𝑎
segmentos de un mismo tamaño ℎ = . La integral total se puede representar como:
𝑛

Al sustituir la regla de Simpson 1/3 en cada integral se obtiene:

Para definirse la variante “compuesta” del método para el cual se utilizan N puntos que corresponden a 𝑛 =
𝑁 − 1 subintervalos. Este caso requiere que la cantidad de subintervalos sean pares (el caso simple utiliza
dos, por lo tanto ´este debe ser un múltiplo). Por lo tanto, se define un valor h que corresponde al ancho del
𝑏−𝑎
subintervalo o el paso que hay entre puntos. Se calcula como ℎ = . Finalmente, la aproximación se puede
𝑛
calcular como:
Los coeficientes “4” y “2” en la ecuación a primera vista parecerían peculiares. No obstante, siguen en forma
natural la regla de Simpson 1/3. Los puntos impares representan el término medio en cada aplicación y, por lo
tanto, llevan el peso de 4 de la ecuación. Los puntos pares son comunes a aplicaciones adyacentes y, por lo
tanto, se cuentan dos veces.

Un error estimado en la regla de Simpson de aplicación múltiple se obtiene de la misma forma que en la regla
del trapecio: sumando los errores individuales de los segmentos y sacando el promedio de la derivada para
llegar a lo siguiente:
(𝑏 − 𝑎)5
𝐸𝑟 = max |𝑓 (4) (𝜉)|
180𝑛4 𝑎≤𝜉≤𝑏

REGLA DE SIMPSON 3/8

De manera similar a la obtención de la regla del trapecio y Simpson 1/3, es posible ajustar un polinomio de
Lagrange de tercer grado a cuatro puntos e integrarlo:
𝑏 𝑏
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑃3 (𝑥)𝑑𝑥
𝑎 𝑎
Para obtener:
3ℎ
= [𝑓(𝑥0 ) + 3𝑓 (𝑥1 ) + 3𝑓 (𝑥2 ) + 𝑓 (𝑥3 )]
8

𝑏–𝑎
Donde ℎ = . Esta ecuación se llama regla de Simpson 3/8 debido a que h se multiplica por 3/8. Ésta es la
3
tercera fórmula de integración cerrada de Newton-Cotes. La regla 3/8 se expresa también en la forma de la
ecuación:
𝑏 𝑏 3ℎ
∫𝑎 𝑓 (𝑥 )𝑑𝑥 = ∫𝑎 𝑃3 (𝑥 )𝑑𝑥 = 8
[𝑓(𝑥0 ) + 3𝑓 (𝑥1 ) + 3𝑓(𝑥2 ) + 𝑓 (𝑥3 )]

LA REGLA DE SIMPSON 3/8 DE APLICACIÓN MÚLTIPLE

Es más exacta que la regla de Simpson 3/8 simple, ya que divide el intervalo de integración en más
subintervalos. Se expresa de la siguiente forma:

𝑏−𝑎
Tomando ℎ = donde n es el número de subintervalos, con la condición de que n sea múltiplo de 3 y que
𝑛
en cada sumatorio se tomen los valores de 𝑖 = 𝑖 + 3:
𝑛 𝑛
𝑏 3ℎ −1 −2
∫𝑎 𝑓 (𝑥 )𝑑𝑥= 8
[𝑓(𝑥0 ) + 3 ∑𝑘=0
3
𝑓(𝑥3𝑘+2 ) + 2 ∑𝑘=0
3
𝑓(𝑥3𝑘+3 ) + 𝑓(𝑥4 )]

Para el cálculo del error, se obtiene la cuarta derivada de la función y tomando en cuenta que 𝜉 debe pertenecer
al intervalo de integración, se aplica la siguiente fórmula
𝑛ℎ5 (4)
𝐸𝑟 = 𝑓 (𝜉)
80
Algunos datos importantes:
1. Los métodos de Simpson 1/3 y 3/8 solo se pueden aplicar a segmentos equidistantes.
2. Es posible aplicar Simpson 1/3 o 3/8 a grupos de segmentos equidistantes y el aporte del resto se
calcula con la regla del Trapecio
3. Considerar que el método de Simpson 3/8 requiere de al menos 4 puntos equidistantes (3 intervalos)
CUADRATURA DE GAUSS
En análisis numérico un método de cuadratura es una aproximación de una integral definida de una función.

Una cuadratura de Gauss n, es una cuadratura que selecciona los puntos de la evaluación de manera óptima
y no en una forma igualmente espaciada, construida para dar el resultado de un polinomio de
grado 2𝑛 − 1 o menos, elegibles para los puntos 𝑥𝑖 y los coeficientes 𝑤𝑖 para 𝑖 = 1, . . . , 𝑛.

El dominio de tal cuadratura por regla es de [−1, 1] dada por:

Tal cuadratura dará resultados precisos solo si 𝑓(𝑥) es aproximado por un polinomio dentro del rango
[−1, 1]. Si la función puede ser escrita como 𝑓(𝑥) = 𝑊(𝑥)𝑔(𝑥), donde 𝑔(𝑥) es un polinomio aproximado
y 𝑊(𝑥) es conocido.

Podemos calcular los coeficientes están dados (método de Gauss-

Legendre) de la siguiente manera:

Donde 𝑃(𝑛) son polinomios de Legendre en el intervalo [−1, 1].

Los cambios de intervalos van de [−1, 1] después de aplicar la

cuadratura de Gauss:

Después de aplicar la cuadratura la aproximación es:

Lista de coeficientes de 𝑤𝑖 y puntos 𝑥𝑖 para 𝑛 = 1, . . . . ,5

INTEGRACIÓN DE ROMBERG
La integración de Romberg es una técnica diseñada para obtener integrales numéricas de funciones
de manera eficiente. Se basa en aplicaciones sucesivas de la regla del trapecio. Sin embargo, a través
de las manipulaciones matemáticas, se alcanzan mejores resultados con menos trabajo.

El método de Romberg genera una matriz triangular cuyos elementos son estimaciones numéricas de
la integral definida:
𝑏
∫ 𝑓 (𝑥 )𝑑𝑥
𝑎
Usando la extrapolación de Richardson de forma reiterada en la regla del trapecio.

Si se aplica la regla de los Trapecios sucesivamente para tamaños de intervalo ℎ𝑘 y requiere conocer
la función o disponer de 𝑛 = 𝑚𝑘 + 1 = 2𝑘−1 + 1 puntos:
Las Fórmulas recursiva de Romberg para trapecios son:

1. Regla compuesta de trapecios:

a. Un trapecio:

b. Dos trapecios:

c. Cuatro trapecios:

En general, 2𝑘−1 trapecios:

Las extrapolaciones que se hacen son, a partir de Richardson,

Con una fórmula de recurrencia para cada coeficiente 𝑗𝑘, 𝑘 > 1:

Al final se llega a una tabla general:

La primera columna está dada por las aproximaciones de la integral dada por la regla recursiva de trapecio. El
resto de las columnas se obtienen por la fórmula de extrapolación de Richardson

Ejemplo: Use la técnica de integración de Romberg para aproximar con 𝑚 = 3 la siguiente integral:
1
2
∫ 𝑒 𝑥 𝑑𝑥
0
Solución. Primero se aplica la regla del trapecio para aproximar la integral con tamaños de paso (segmentos)
1 2
de longitud ℎ1 = 1, ℎ2 = 𝑦 ℎ3 = 1/4. Se tiene que 𝑓(𝑥) = 𝑒 𝑥 , 𝑎 = 0, 𝑏 = 1.
2

ℎ1
𝑅1,1 = [𝑓(𝑎) + 𝑓(𝑏)] = 1.859140914,
2
1
Para ℎ2 = se tiene que 𝑥1 = 𝑎 + ℎ2 = 0.5,
2

ℎ2
𝑅2,1 = [𝑓(𝑎) + 𝑓(𝑏) + 2𝑓(𝑥1 )] = 1.571583165,
2

Para ℎ3 = 1/4, 𝑥1 = 𝑎 + ℎ3 = 0.25, 𝑥2 = 𝑥1 + ℎ3 = 0.50 𝑦 𝑥3 = 𝑥2 + ℎ3 = 0.75,

ℎ3
𝑅3,1 = [𝑓(𝑎) + 𝑓(𝑏) + 2𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 )] = 1.490678862,
2
Ahora tenemos que:

Donde tenemos que 𝑅3,3 = 1.46290944 es la mejor estimación para la integral:

1
2
∫ 𝑒 𝑥 𝑑𝑥
0