Integracion numerica

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 [pic] una función contínua en el intervalo [pic] y sea [pic] una antiderivada de [pic]. Entonces: 
[pic] 
El problema en la práctica, se presenta cuando nos vemosimposibilitados de encontrar la antiderivada requerida, aún para integrales aparentemente sencillas como: 
[pic] 
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 deintegració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, enMathematica) 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 [pic]: 
[pic] 
donde [pic] es un polinomio de interpolación de grado [pic] para ciertos datos de [pic] que se escogen apropiadamente. 
Es importante observarque 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 [pic] y [pic]; en caso contrario, se llaman formas abiertas. 
Nosotros nos remitiremos aestudiar únicamente las formas cerradas, y por lo tanto, siempre suponemos que conocemos los valores [pic] y [pic].
REGLA DEL TRAPECIO 
Corresponde al caso donde [pic] , es decir :  
[pic] 
donde [pic] es un polinomio de interpolación (obviamente de grado 1) para los datos:  
[pic]
Del capítulo anterior, sabemos que este polinomio de interpolación es:
[pic] 
Integrando este polinomio,tenemos que:
| |[pic] |
| |[pic] |
| |[pic]|
| |[pic] |
| |[pic] |


Por lo tanto, tenemos que: 
[pic] 
Que es la conocida Regla del Trapecio. Este nombre se debe a la interpretacióngeomé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 [pic], que es precisamente el área del trapecio que se forma.
[pic]
Ejemplo 1:
Utilizar la regla del trapecio para aproximar la integral: 
[pic] 
Solución.
Usamos la fórmula directamente con lossiguientes datos:  
[pic]
  Por lo tanto tenemos que: 
[pic] 
Ejemplo 2.  
Usar la regla del trapecio para aproximar la integral:  
[pic]  
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:  
[pic]
Por lo tanto, tenemos que: 
[pic]
La regla del trapecio se puede...