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 3:
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 4.
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 5.
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 6.
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 7.
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 8.
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 9.
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:

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

Ejemplo10.

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 11.

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.

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