Integración numérica, Método del Trapecio
En matemática la regla del trapecio es un método de integración numérica, es decir, un
método para calcular aproximadamente el valor de la integral definida
b
La regla se basa en aproximar el valor de la integral ∫ f ( x ) dx de f(x) por el de la función
a
lineal que pasa a través de los puntos (a, f(a)) y (b, f(b)). La integral de ésta es igual al área
del trapecio bajo la gráfica de la función lineal. Se sigue que
Y el error es:
Siendo ε un número entre a y b
Regla del trapecio compuesta
La regla del trapecio compuesta o regla de los trapecios es una forma de aproximar una
integral definida utilizando n trapecios. En la formulación de este método se supone que f
es continua y positiva en el intervalo [a,b].
b
De tal modo la integral definida ∫ f ( x ) dx representa el área de la región delimitada por la
a
gráfica de f y el eje x, desde x=a hasta x=b. Primero se divide el intervalo [a,b] en n
b−a
subintervalos, cada uno de ancho ∆ x= .
n
Después de realizar todo el proceso matemático se llega a la siguiente fórmula:
b−a
Donde h= y n es el número de divisiones.
n
La expresión anterior también se puede escribir como:
Estudiemos esta regla al utilizar trapecios para aproximar el área bajo la gráfica de la
función f(x) = 3 *ln(x)f(x) en el intervalo [2, 8].
Así es como se ve la regla en un diagrama, donde llamamos al primer trapecio T_1 al
segundo T_2 y al tercero T_3.
Recuerda que el área de un trapecio está dada por h= ( b 1+b2 2 ), donde h es la altura y b1
y b2, son las longitudes de las bases.
Encontrar el área de T1
Debemos pensar que el trapecio está de costado.
La altura h es el 2 en la parte baja de T1, que genera el intervalo comprendido entre x =2 y
x=4.
La primera base, b1, es el valor de 3 *ln(x) en x=2, es decir, 3ln(2)
La segunda base, b2 es el valor de 3 *ln(x) en x=4, es decir, 3ln(4)
Así es como se ve todo esto:
Encontramos el área total de la aproximación al sumar el área de cada uno de los trapecios:
Area total = T1+T2+T3
La respuesta fianl ya simplificada:
Area total= 3(ln2+2ln4+2ln6+ln8) = 7.035
Método de Simpson
Hemos utilizado diferentes métodos numéricos para aproximar el valor de las integrales
definidas. Los métodos que hemos observado hasta el momento utilizar la segmentación
rectangular o trapezoidal del área bajo la curva de la función, pero no son tan exactos para
funciones donde los segmentos de líneas rectas no son la mejor aproximación para la
función curva. Es razonables pensar que otros métodos de aproximación curvas puede ser
más útiles para algunas funciones. La. regla de Simpson es un método que utiliza parábolas
para aproximar la curva en vez trazar segmentos lineales.
Este método consiste en la aproximación del cálculo del ´área plana bajo una curva
utilizando trapecios curvilíneos a partir una interpolación con una función cuadrática:
En este procedimiento, se toma el intervalo de anchura 2h, comprendido entre xi y xi+2, y se
sustituye la función f(x) por la parábola que pasa por tres puntos (xi, yi), (xi+1, yi+1), y (xi+2,
yi+2).
Calculamos la contribución a la integral del primer intervalo (x0, x0+2h) y generalizaremos
al resto de los intervalos.
La ecuación de la parábola y=ax2+bx+c que pasa por los puntos (x0, y0), (x0+h, y1), (x0+2h,
y2) es
Este sistema de tres ecuaciones con tres incógnitas, se reduce a
Despejamos el coeficiente a, y 2ax0+b
Sustituimos la curva por la porción de parábola en el intervalo (x0, x0+2h). La integral
aproximada vale
En general, el valor del área aproximada, en el intervalo (xi, xi+2h) sombreada en la figura,
es
El primer paréntesis, contiene la suma de los extremos, el segundo, la suma de los términos
de índice impar, y el tercero la suma de los términos de índice par. En el método de
Simpson, el número de divisiones n debe de ser par. En el caso de que el usuario introduzca
un número impar el programa lo convierte en el número par siguiente.