3er corte 20%
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
x f(x)
0 0,1
0,2 0,20816215
0,4 0,23470216
0,6 0,28666585
0,8 0,3792961
1 0,2718281
sumatoria 1,48065437
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.
N=10
1-0=1/10=0.1
x f(x)
0 0,05
0,1 0,10100502
0,2 0,10408108
0,3 0,10941743
0,4 0,11735108
0,5 0,12840253
0,6 0,14333293
0,7 0,1632316
0,8 0,18964805
0,9 0,22479074
1 0,13591405
sumatoria 1,4671745
Puede elaborar una tabla en Excel .pero para efectos del examen recuerden dejar
todos los cálculos plamados
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:
(0+0.2)/2=0.1
(0.2+0.4)/2=0.3
(0.4+0.6)/2=0.5
.
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:
PARA CALCULAR LA PARTICION RECUERDEN SE RESTAN LOS LIMITES DE
INTEGRACION (4-2)/N=2/4=0.5. ENTONCES EMPIEZAN A SUMAR 0.5 INICIANDO
EN 2
ES DECIR; 2+0.5=2.5
2.5+0.5=3
3+0.5=3.5
Y ASI SUCESIVAMENTE
Además, los puntos medios de cada subintervalo son:
RECUERDEN PARA SACAR LOS PUNTOS MEDIOS
PM1= (2+2.5)/2=2.25
PM2=(2.5+3)/2=2.75
PM3=(3+3.5)/2=3.25
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:
PARTICION RECUERDEN
(4-1)/N=3/3=1
RNTONCE VAN DE 1 EN 1 FORMANDO LAS PARTICIONES
1+1=2
2+1=3
3+1=4
Al considerar los puntos que dividen en tres partes iguales a cada subintervalo,
tenemos los siguientes datos:
PARA CALCULAR ESTOS SUBINTERVALOS. RESTAN (2-1)/3=0.333333
LUEGO LE SUMAN EL 1
1+0.333333=1.3333333=4/3
1.3333333=4/3+0.33333333=1.6666666=5/3
5/3 Ó 1.6666666+0.33333333=1.99999999= SE APROXIMA A 2
LUEGO SEGUIMOS SUMANDO 0.3333333 Y NOS VAN DANDO LOS OTRO VALORES
POR COMODIDAD ES MEJOR TRABAJAR EN FRACCIONES
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. De la
misma longitud
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 regla del trapecio, en el
intervalo la regla de Simpson de 3/8 y en el intervalo la regla de
Simpson de 1/3. Así, tenemos las siguientes integrales:
0.342
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.
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
