INTEGRACIÓN NUMÉRICA

En los cursos de Cálculo Integral, nos enseñan como calcular una integral definida de
una función contínua mediante una aplicación del Teorema Fundamental del Cálculo: 

Teorema Fundamental del Cálculo

Sea una función contínua en el intervalo y sea una antiderivada de
. Entonces: 

El problema en la práctica, se presenta cuando nos vemos imposibilitados de encontrar

la antiderivada requerida, aún para integrales aparentemente sencillas como: 

la cual simplemente es imposible de resolver con el Teorema Fundamental del Cálculo.

En este capítulo estudiaremos diversos métodos numéricos que nos permitirán obtener
aproximaciones bastante exactas a integrales como la mencionada anteriormente.
Esencialmente, veremos dos tipos de integración numérica: las fórmulas de Newton-
Cotes y el algoritmo de Romberg. 

Las fórmulas de Newton-Cotes están conformadas por las bien conocidas reglas del
trapecio y de Simpson (regla de un tercio y de tres octavos). El algoritmo de Romberg
forma parte de un método conocido como método de extrapolación de Richardson. 

Haciendo uso de algunos programas computacionales (por ejemplo, en Mathematica)

es posible discernir sobre las cualidades y defectos de cada uno de los métodos
mencionados arriba.

FORMULAS DE INTEGRACION DE NEWTON-COTES

Estas fórmulas se basan en la idea de integrar una función polinomial en vez de

: 

donde es un polinomio de interpolación de grado para

ciertos datos de que se escogen apropiadamente. 

Es importante observar que estas fórmulas se pueden aplicar inclusive a una tabla de
datos, ya que lo que se usa es un polinomio de interpolación, el cual puede ser
calculado con la tabla. 

Dentro de las fórmulas de Newton-Cotes, existen las formas cerradas y abiertas. En

las formas cerradas se conocen los valores de y ; en caso contrario, se
llaman formas abiertas. 

Nosotros nos remitiremos a estudiar únicamente las formas cerradas, y por lo tanto,
siempre suponemos que conocemos los valores y .

REGLA DEL TRAPECIO 

Corresponde al caso donde , es decir :  

donde es un polinomio de interpolación (obviamente de grado 1) para los

datos:  

Del capítulo anterior, sabemos que este polinomio de interpolación es:

Integrando este polinomio, tenemos que:

Por lo tanto, tenemos que: 

Que es la conocida Regla del Trapecio. Este nombre se debe a la interpretación

geométrica que le podemos dar a la fórmula. El polinomio de interpolación para una
tabla que contiene dos datos, es una línea recta. La integral, corresponde al área bajo
la línea recta en el intervalo , que es precisamente el área del trapecio que se
forma.

Ejemplo 1:
Utilizar la regla del trapecio para aproximar la integral: 

Solución.
Usamos la fórmula directamente con los siguientes datos:  

  Por lo tanto tenemos que: 

Ejemplo 2.  
Usar la regla del trapecio para aproximar la integral:  

 
Solución.  
Igual que en el ejemplo anterior, sustituímos los datos de manera directa en la fórmula
del trapecio. En este caso, tenemos los datos:  

Por lo tanto, tenemos que: 

La regla del trapecio se puede ampliar si subdividimos el intervalo en

subintervalos, todos de la misma longitud . 

Sea la partición que se forma al hacer dicha subdivisión. Usando

propiedades de la integral tenemos que: 

Aplicando la regla del trapecio en cada una de las integrales, obtenemos:  

Ahora bien, ya que todos los subintervalos tienen la misma longitud h, tenemos que:  

Sustituyendo el valor de h y usando la notación sigma, tenemos finalmente: 

 
Esta es la regla del trapecio para n subintervalos. Obviamente, esperamos que entre
más subintervalos usemos, mejor sea la aproximación a la integral. 

Ejemplo 1:  
Aplicar la regla del trapecio para aproximar la integral  

si subdividimos en 5 intervalos.

Solución.  
En este caso, identificamos , y la partición generada es:  

Así, aplicando la fórmula tenemos que:

= 1.48065  

Cabe mencionar que el valor verdadero de esta integral es de 1.4626… 

Así, vemos que con 5 intervalos, la aproximación no es tan mala. Para hacer cálculos
con más subintervalos, es conveniente elaborar un programa que aplique la fórmula
con el número de subintervalos que uno desee. El lector debería hacer su propio
programa y checar con 50, 500, 1000, 10000 y 20000 subintervalos, para observar el
comportamiento de la aproximación.
REGLA DE SIMPSON DE UN TERCIO
Suponemos que tenemos los datos: 

donde es el punto medio entre y .

En este caso se tiene que:

donde es el polinomio de interpolación para los datos en la tabla anterior.

Usaremos el polinomio de Lagrange.

Así, tenemos que:

Si denotamos , entonces: 

Simplificando términos: 

Vemos que cada uno de los términos anteriores, es esencialmente de la misma forma,
es decir, una constante por

Así, calculamos la siguiente integral por partes: 

Sea:
por lo tanto, 

Usamos esta fórmula para calcular la integral de cada uno de los tres términos de
.
Debido al factor se le conoce como la regla de Simpson de un tercio.

En la práctica, sustituímos el valor de para obtener nuestra fórmula final: 

Ejemplo 1.
Usar la regla de Simpson de 1/3 para aproximar la siguiente integral: 

Solución.
Aplicamos la fórmula directamente, con los siguientes datos: 

Por lo tanto, tenemos que: 

  

Ejemplo 2.
Usar la regla de Simpson de 1/3, para aproximar la siguiente integral: 

Solución.
Igual que en el ejercicio anterior, sustituímos datos adecuadamente:

 
Al igual que con la regla del trapecio, podemos extender la regla de Simpson de 1/3,
si subdividimos el intervalo en subintervalos de la misma longitud

. 

Sea la partición que se forma al hacer la subdivisión, y denotemos

por el punto medio en cada subintervalo. 

Aplicamos primero propiedades básicas de la integral definida: 

Ahora, aplicamos la regla de Simpson de 1/3, en cada una de las integrales de arriba: 

Sustituímos y usamos la notación sigma:

Ejemplo 1.

Aproximar la siguiente integral, aplicando la regla de Simpson de y subdividiendo

en 5 intervalos.  

Solución.
En este caso, tenemos que , y la partición que se genera es: 

Además, los puntos medios de cada subintervalo son: 

 
Por lo tanto, sustituímos los datos en la fórmula para obtener: 

Nótese que esta aproximación ya es exacta hasta el cuarto decimal! 

Ejemplo 2.

Aproximar la siguiente integral, utilizando la regla de Simpson de y subdividiendo

en 4 intervalos. 

Solución.
En este caso, tenemos que , y la partición que se genera es: 

Además, los puntos medios de cada subintervalo son: 

Sustituyendo todos estos datos en la fórmula obtenemos la siguiente aproximación: 

REGLA DE SIMPSON DE TRES OCTAVOS 
Este caso corresponde a , es decir, 

donde es un polinomio de interpolación para los siguientes datos: 

Y donde , y , son los puntos que dividen en tres partes iguales al

intervalo .

Igual que en el caso anterior, se usa el polinomio de interpolación de Lagrange, y

usando el método de integración por partes se llega a la siguiente fórmula: 

donde . Debido al factor es que se le dió el nombre de Regla de

Simpson de 3/8. En la práctica, se sustituye el valor de h para obtener: 

Ejemplo 1.

Aproximar la siguiente integral, usando la regla de Simpson de : 

Solución.
En este caso, tenemos los siguientes datos: 
  Los cuales sustituímos en la fórmula, para obtener: 

Al igual que en los dos casos anteriores, la regla de Simpson de 3/8, se puede
extender si subdividimos el intervalo en intervalos de la misma longitud

. 

Sea la partición determinada de esta forma. Cada subintervalo lo

dividimos en tres partes iguales, y sean y los puntos determinados así: 

Aplicando la regla de en cada uno de los intervalos tenemos: 

 
    

Esta última, es la regla de Simpson de 3/8 para n subintervalos todos de la misma

longitud. 

Ejemplo 2.
Aproximar la siguiente integral: 

aplicando la regla de Simpson de 3/8, y subdiviendo en 3 intervalos. 

Solución.
Identificamos y la partición correspondiente:

Al considerar los puntos que dividen en tres partes iguales a cada subintervalo,
tenemos los siguientes datos: 

Sustituyendo todos los datos en la fórmula, obtenemos: 

De acuerdo a los ejemplos vistos, resulta evidente que la regla de Simpson de 3/8, es
más exacta que la de 1/3 y a su vez, ésta es más exacta que la regla del trapecio. En
realidad, pueden establecerse cotas para los errores que se cometen en cada uno de
estos métodos.

Puesto que no es nuestra intención justificar formalmente cada uno de los teoremas,
los siguientes resultados se mencionan para completar la información, pero omitimos
las demostraciones correspondientes.

REGLA FORMULA ERROR DONDE...

Trapecio

Simpson

Simpson

INTEGRACIÓN EN INTERVALOS DESIGUALES

  Cuando la longitud de los subintervalos no es igual, se usa una combinación de
la regla Trapezoidal y las reglas de Simpson, procurando seguir el siguiente orden
jerárquico: 

1 .- Simpson

Esta se aplica, si contamos con 4 puntos igualmente espaciados. 

2 .- Simpson

Esta se aplica si falla (1) y contamos con 3 puntos igualmente

espaciados. 

3 .- Regla Trapezoidal

Solo se aplica si no se cumple y  

Ejemplo 1.

Evaluar , usando la siguiente tabla :

Solución.
Vemos que en el intervalo podemos aplicar la regal del trapecio, en el intervalo
 la regal de Simpson de 3/8 y en el intervalo la regal de Simpson de
1/3. Así, tenemos las siguientes integrales:

Finalmente, la integral buscada es la suma de las tres integrales anteriores: 

Ejemplo 2.

Calcula la integral , usando la siguiente tabla de datos:

Solución.
En este caso, vemos que podemos aplicar la regla de Simpson de 1/3 en el intervalo
, la regla del trapecio en el intervalo y la regla de Simpson de 3/8 en el
intervalo . Así, tenemos las siguientes integrales:

Por lo tanto, la integral buscada es la suma de las tres integrales anteriores: 

 
Vale la pena comentar que no siempre tiene que suceder que se apliquen exactamente
las tres reglas. En realidad, esto depende de cómo se encuentran espaciados los
intervalos de la tabla de datos.

MÉTODO DE INTEGRACIÓN DE ROMBERG 

Sea el valor de la integral que aproxima a , mediante una

partición de subintervalos de longitud y usando la regla del trapecio.

Entonces, 

donde es el error de truncamiento que se comete al aplicar la regla. 

El método de extrapolación de Richardson combina dos aproximaciones de

integración numérica, para obtener un tercer valor más exacto.

El algoritmo más eficiente dentro de éste método, se llama Integración de

Romberg , la cual es una fórmula recursiva.
Supongamos que tenemos dos aproximaciomnes : e  

Se puede demostrar que el error que se comete con la regla del trapecio para n
subintervalos está dado por las siguientes fórmulas: 

donde es un promedio de la doble derivada entre ciertos valores que pertenecen a

cada uno de los subintervalos.

Ahora bien, si suponemos que el valor de es constante, entonces : 

  

Sustituyendo esto último en nuestra primera igualdad, tenemos que: 

De aquí podemos despejar : 

En el caso especial cuando (que es el algoritmo de Romberg), tenemos : 

Esta fórmula es solo una parte del algoritmo de Romberg. Para entender el método, es
conveniente pensar que se trabaja en niveles de aproximación. En un primer nivel, es
cuando aplicamos la regla del Trapecio, y para poder usar la fórmula anterior,
debemos de duplicar cada vez el número de subintervalos: así, podemos comenzar con
un subintervalo, luego con dos, cuatro, ocho, etc, hasta donde se desee. 

Posteriormente, pasamos al segundo nivel de aproximación, que es donde se usa la

fórmula anterior, tomando las parejas contiguas de aproximación del nivel anterior, y

que corresponden cuando .

Después pasamos al nivel tres de aproximación, pero aquí cambia la fórmula de

Romberg, y así sucesivamente hasta el último nivel, que se alcanza cuando solo
contamos con una pareja del nivel anterior.  

Desde luego, el número de niveles de aproximación que se alcanzan, depende de las

aproximaciones que se hicieron en el nivel 1. En general, si en el primer nivel,
iniciamos con n aproximaciones, entonces alcanzaremos a llegar hasta el nivel de
aproximación n. 

Hacemos un diagrama para explicar un poco más lo anterior. 

Ejemplo 1.
Usar el algoritmo de Romberg, para aproximar la integral 

usando segmentos de longitud .

Solución.
Primero calculamos las integrales del nivel 1, usando la regla del trapecio para las
longitudes de segmentos indicadas: 
Con estos datos, tenemos:

Ahora pasamos al segundo nivel de aproximación donde usaremos la fórmula que se

dedujo anteriormente: 

donde es la integral menos exacta (la que usa menos subintervalos) e es

la más exacta (la que usa el doble de subintervalos).

En un diagrama vemos lo siguiente: 

Para avanzar al siguiente nivel, debemos conocer la fómula correspondiente. De forma

similar a la deducción de la fórmula, 
  

se puede ver que la fórmula para el siguiente nivel de aproximación (nivel 3) queda
como sigue: 

donde: 

es la integral más exacta

es la integral menos exacta 

En el siguiente nivel (nivel 4) se tiene la fórmula  

En el ejemplo anterior, obtenemos la aproximación en el nivel 3 como sigue: 

Así, podemos concluir que el valor de la aproximación, obtenido con el método de

Romberg en el ejemplo 1, es: 

Ejemplo 2.
Usar el algoritmo de Romberg para aproximar la integral: 

Agregando a la tabla anterior donde .

Solución.
Calculamos con la regla del trapecio:
Tenemos entonces la siguiente tabla: 

De donde concluímos que la aproximación buscada es: 

Ejemplo 3.
Aproximar la siguiente integral: 

usando el método de Romberg con segmentos de longitud

, , ,  

Solución.
Igual que arriba, primero usamos la regla del trapecio (con los valores de h
indicados) para llenar el nivel 1. Tenemos entonces que:
A continuación, usamos las fórmulas de Romberg para cada nivel y obtenemos la
siguiente tabla: 

De donde concluímos que la aproximación buscada es: 

Podemos escribir una fórmula general para calcular las aproximaciones en cada uno de
los niveles como sigue: 

ALGORITMO DE INTEGRACIÓN DE ROMBERG 

Los coeficientes en cada una de las fórmulas en el método de Romberg, deben

sumar 1.
Así se tiene la siguiente fórmula recursiva: 

donde: 

es la integral más exacta

es la integral menos exacta 

y el indice k indica el nivel de integración o de aproximación. Por ejemplo, digamos

que , entonces tenemos: 

que es nuestra fórmula del nivel 2 de aproximación. 

Como todo proceso iterativo, éste se detiene cuando se obtiene una aproximación
suficientemente buena. En este caso se pide que: 
  

donde es la cota suficiente. 

Ejemplo 1.
Aplicar el algoritmo de integración de Romberg a la integral: 

tomando

Solución.
En este caso no sabemos exactamente cuantas aproximaciones debemos hacer con la
regla del trapecio. Así que para comenzar hacemos los cálculos correspondientes a
uno, dos, cuatro y ocho subintervalos:

Con estos datos, podemos hacer los cálculos hasta el nivel 4. Tenemos la siguiente
tabla: 

Haciendo los cálculos de los errores, nos damos cuenta que efectivamente la
aproximación se obtiene hasta el nivel 4, donde .
Por lo tanto, concluímos que la aproximación buscada es: 
EJERCICIOS 
1. Usar la regla del trapecio para aproximar, 

i)                   Dividiendo en un solo intervalo.

ii)                 Dividiendo en 6 intervalos. 

Soluciones: i) 3.4115 ii) 0.36907 

2. Usar la regla de Simpson 1/3 para aproximar, 

i)                   Dividiendo en un solo intervalo.

ii)                 Dividiendo en 4 intervalos. 

Soluciones: i) 82.60511 ii) 76.94497 

3. Usar la regla de Simpson 3/8 para aproximar, 

i)                   Dividiendo en un solo intervalo.

ii)                 Dividiendo en 4 intervalos. 

Soluciones: i) 2.76591 ii) 2.76501 

4. Integrar las siguientes tablas de datos: 

i)                    
 ii)                

Soluciones: i) -17.11458 ii) 9.425

5. Usar el algoritmo de integración de Romberg para aproximar, 

i)                   Usando 1, 2 y 4 intervalos.

ii)                 Agregando al inciso anterior, 8 intervalos.

Soluciones: i) 9.156626413 ii) 9.153287278 

6. Aproxime la integral del ejercicio anterior, tomando como cota

suficiente. 

Solución. 9.153112082