Integracion Numerica
Integracion Numerica
Integracion Numerica
INTEGRACIN NUMRICA
Sea de
y sea
una antiderivada
El problema en la prctica, se presenta cuando nos vemos imposibilitados de encontrar la antiderivada requerida, an para integrales aparentemente sencillas como:
La cual simplemente es imposible de resolver con el Teorema Fundamental del Clculo. Las frmulas de Newton-Cotes estn 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 mtodo conocido como mtodo de extrapolacin de Richardson. FORMULAS DE INTEGRACION DE NEWTON-COTES Estas frmulas se basan en la idea de integrar una funcin polinomial en vez :
de
Dentro de las frmulas 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 .
Donde datos:
Que es la conocida Regla del Trapecio. Este nombre se debe a la interpretacin geomtrica que le podemos dar a la frmula. El polinomio de interpolacin para una tabla que contiene dos datos, es una lnea recta. La integral, corresponde al rea bajo la lnea recta en el intervalo , que es precisamente el rea del trapecio que se forma.
Solucin. Igual que en el ejemplo anterior, sustituimos los datos de manera directa en la frmula del trapecio. En este caso, tenemos los datos:
La regla del trapecio se puede ampliar si subdividimos el intervalo subintervalos, todos de la misma longitud .
en
Sea la particin que se forma al hacer dicha subdivisin. Usando propiedades de la integral tenemos que:
Ahora bien, ya que todos los subintervalos tienen la misma longitud que:
h, tenemos
Esta es la regla del trapecio para n subintervalos. Obviamente, esperamos que entre ms subintervalos usemos, mejor sea la aproximacin a la integral. Ejemplo 1: Aplicar la regla del trapecio para aproximar la integral
= 1.48065 Cabe mencionar que el valor verdadero de esta integral es de 1.4626 As, vemos que con 5 intervalos, la aproximacin no es tan mala. Para hacer clculos con ms subintervalos, es conveniente elaborar un programa que aplique la frmula con el nmero de subintervalos que uno desee. El lector debera hacer su propio programa y checar con 50, 500, 1000, 10000 y 20000 subintervalos, para observar el comportamiento de la aproximacin.
Donde
Al igual que con la regla del trapecio, podemos extender la regla de Simpson, si subdividimos el intervalo en subintervalos de la misma longitud .
Sustituimos
Ejemplo 1. Aproximar la siguiente integral, aplicando la regla de Simpson compuesta y subdividiendo en 5 intervalos.
Ntese que esta aproximacin ya es exacta hasta el cuarto decimal! Ejemplo 2. Aproximar la siguiente integral, utilizando la regla de Simpson compuesta y subdividiendo en 4 intervalos.