Integración aproximada: Regla del trapecio

Algunas veces no es nada sencillo calcular la antiderivada de una función dada. En esos casos es mejor hacer una aproximación al valor del área debajo de la curva utilizando métodos numéricos ampliamente conocidos. La regla del trapecio consiste utilizar trapeciod en lugar de rectángulos al hacer la aproximación del área bajo la curva.

Ya hicimos la primera aproximación del valor del área bajo la parábola $y = x^2$ desde $x=0$ hasta $x=1$ en alguna lección previa. Ahora vamos a hacer la misma aproximación usando trapecios.

Ahora, la comparación de las aproximaciones usando rectángulos por un lado y trapecios por otro, se muestra enseguida:

$\begin{minipage}{0.5\linewidth} \begin{center} \textbf{Usando rect\'angulos:}\\ \begin{tabular}{lll}\hline \textbf{i} & $A_{inf}$ & $A_{sup}$\\\hline 1 & 0.0 & 0.0080 \\ 2 & 0.0080 & 0.032 \\ 3 & 0.032 & 0.072 \\ 4 & 0.072 & 0.128 \\ 5 & 0.128 & 0.2 \\ \hline \textbf{Totales:} & 0.24 & 0.44 \\ \hline \end{tabular} \end{center} \end{minipage} % + + + + + + + + + + + + + + + + + + + \begin{minipage}{0.5\linewidth} \begin{center} \textbf{Usando trapecios:}\\ \begin{tabular}{ll}\hline \textbf{i} & $dA$ \\\hline 1 & 0.0040 \\ 2 & 0.02 \\ 3 & 0.052 \\ 4 & 0.1 \\ 5 & 0.164 \\ \hline \textbf{Totales:} & 0.34 \\ \hline \end{tabular} \end{center} \end{minipage}$

Solo para ver la diferencia, recuerda que el área bajo la curva es exactamente de $1/3$ . Si comparas esta gráfica con el uso de rectángulos en lugar de trapecios, verás por qué esta nueva aproximación es mucho mejor: los trapecios se acercan mucho mejor a la curva que los rectángulos. Esa es la idea que está detrás de la regla del trapecio.

Contenido

1 Ejemplo
2 Ejemplo
3 Ejemplo
4 Ejemplo
5 Ejemplo
6 Ejemplo

Ejemplo

Calcula el área bajo la curva de la función:

$f(x) = \frac{1}{\sqrt{2\,\pi}}\,e^{-x^2/2}$

desde $x=0$ hasta $x = 1$ .

Geométricamente, tenemos que calcular la siguiente área:

Si intentamos calcular la antiderivada para aplicar el teorema fundamental del cálculo tendremos serias dificultades. En términos de integral definida, el problema se puede pronunciar como:

Ejemplo

Calcula:

$\begin{equation*} \frac{1}{\sqrt{2\,\pi}}\int\limits_{0}^{1}\!e^{-x^2/2}\,dx \end{equation*}$

Calcular la antiderivada de esta función es posible a través de métodos de cálculo avanzado. Así que lo más recomendable por ahora es hacer una aproximación. Podemos hacer una primera aproximación ocupando rectángulos como hicimos en la lección Área bajo la curva. Pero podemos mejorar la solución si en lugar de dibujar rectángulos en su lugar utilizamos trapecios como se muestra enseguida:

Usando dos trapecios en el intervalo obtenemos los siguientes resultados:

$\begin{tabular}{ll}\hline \textbf{i} & $dA$ \\\hline 1 & 0.18775 \\ 2 & 0.14851 \\ \midrule \enc{Total:} & 0.33626 \\ \hline \end{tabular} \end{center}$

Igual, podemos aumentar el número de particiones y comparar los resultados. Vamos a utilizar $n=5$ para hacer la comparación con la aproximación anterior ( $n=2$ ).

Y las áreas de cada una de las particiones es:

$\begin{tabular}{ll}\hline \textbf{i} & $dA$ \\\hline 1 & 0.079 \\ 2 & 0.07593 \\ 3 & 0.07015 \\ 4 & 0.06229 \\ 5 & 0.05317 \\ \hline \enc{Total:} & 0.34054 \\ \hline \end{tabular}$

Si ahora hacemos $n=25$ obtenemos una muy buena aproximación al valor real:

En la siguiente tabla se muestran las áreas de cada trapecio dibujado en la gráfica anterior, y la aproximación del área bajo la curva en este caso es: 0.34131 unidades de área.

Ejemplo

Utiliza la regla del trapecio para aproximar la integral definida:

$\begin{equation*} \int\limits_{-1}^{1}\!e^{x^2}\,dx \end{equation*}$

El cálculo de esta integral de manera analíica es imposible por los métodos que hemos estudiado. Así que empezamos definiendo $n=10$ :

Y haciendo $n = 20$ , obtenemos:

Y el valor de esta aproximación es: 2.93435 unidades de área. El valor del área buscada (correcto a 8 decimales) es de 2.925303492 unidades de área.

Observa que para calcular el área $\Delta A$ de un trapecio utilizamos la fórmula:

$\begin{equation*} \Delta A = \left(\frac{f(x_i) + f(x_i + \Delta x)}{2}\right)\cdot \Delta x \end{equation*}$

y para calcular la aproximación del área debajo de la curva usando trapecios en lugar de rectángulos ocupamos la sumatoria de todas las áreas de los $n$ trapecios que hemos dibujado debajo de la curva:

$\begin{equation*} \int\limits_{a}^{b}\!f(x)\,dx \approx \sum\limits_{i=1}^{n}\left(\frac{f(x_i) + f(x_i + \Delta x)}{2}\right)\cdot \Delta x \end{equation*}$

donde $\Delta x = (b - a)/n$ .

Si definimos: $x_{i+1} = x_{i} + \Delta x$ , podemos reescribir la expresión anterior de la siguiente forma:

$\begin{eqnarray*} \int\limits_{a}^{b}\!f(x)\,dx &\approx& \sum\limits_{i=1}^{n}\left(\frac{f(x_i) + f(x_i + \Delta x)}{2}\right)\cdot \Delta x\\ &=& \sum\limits_{i=1}^{n}\left(\frac{f(x_i) + f(x_{i+1})}{2}\right)\cdot \Delta x\\ &=& \frac{\Delta x}{2}\cdot\left([f(x_1) + f(x_2)] + [f(x_2) + f(x_3)] + \cdots + [f(x_{n-1}) + f(x_n)]\right)\\ &=& \frac{\Delta x}{2}\cdot\left[f(x_1) + 2\,f(x_2) + 2\,f(x_3) + \cdots + 2\,f(x_{n-1}) + f(x_n)\right]\\ &=& (\Delta x)\cdot\left(f(x_2) + f(x_3) + \cdot + f(x_{n-1})\right) + (\Delta x)\cdot\left(\frac{f(x_1) + f(x_n)}{2}\right) \end{eqnarray*}$

Mucho software computacional que se utiliza para calcular integrales utiliza algún método como el que se acaba de explicar. La computadora es programada para realizar los cálculos. El usuario solamente debe indicar el número de intervalos que desea utilizar.

Ejemplo

Utiliza la regla trapezoidal para aproximar:

$\begin{equation*} \int\limits_{1}^{2}\!\frac{dx}{x} \end{equation*}$

haciendo $n = 10$ .

Nosotros haremos $n=10$ , luego:

$\begin{equation*} \Delta x = \frac{2-1}{10} = 0.1 \end{equation*}$

De acuerdo a la regla trapezoidal, tenemos:

$\begin{eqnarray*} \int\limits_{a}^{b}\!f(x)\,dx &\approx& (\Delta x)\cdot\left(f(x_2) + f(x_3) + \cdot + f(x_{n-1})\right) + (\Delta x)\cdot\left(\frac{f(x_1) + f(x_n)}{2}\right)\\ &=& (0.1)\left[f(1.1) + f(1.2) + \cdots + f(1.9)\right] + (0.1)\left(\frac{f(1) + f(2)}{2}\right) \end{eqnarray*}$

En este caso particular, $f(x) = 1/x$ , así que:

$\begin{eqnarray*} \int\limits_{1}^{2}\!\frac{dx}{x} &\approx& (0.1)\left[f(1.1) + f(1.2) + \cdots + f(1.9)\right] + (0.1)\left(\frac{f(1) + f(2)}{2}\right)\\ &=& (0.1)\left(\frac{1}{1.1} + \frac{1}{1.2} + \cdots + \frac{1}{1.9}\right) + (0.1)\left(\frac{1 + 0.5}{2}\right)\\ &=& (0.1)(6.187714032) + (0.1)(0.75)\\ &=& 0.6187714032 + 0.075 = 0.6937714032 \end{eqnarray*}$

Nosotros ya sabemos que:

$\begin{equation*} \int\limits_{1}^{2}\frac{dx}{x} = \left.\textcolor{white}{\frac{}{}}\!\!\!\ln x\right\vert_{1}^{2} = \ln(2) - \ln(1) \approx \ln 2 = 0.6931471806 \end{equation*}$

Así que hemos una aproximación correcta hasta 3 decimales. Representa geométricamente la aproximación que hemos calculado usando la regla del trapecio en papel milimétrico.

Ejemplo

Aplicando la regla del trapecio aproxima:

$\begin{equation*} \int\limits_{0}^{\pi^2}\!\sin\left(\sqrt{x}\right)\,dx \end{equation*}$

Geométricamente, tenemos que calcular el área sombreada en la siguiente gráfica:

Observa que el integrando se hace cero exactamente en $\pi^2$ :

$\begin{equation*} \sin\left(\sqrt{\pi^2}\right) = \sin\pi = 0 \end{equation*}$

Empezamos la primera aproximación haciendo $n = 10$ . Entonces, $\Delta x = 0.1\,\pi^2$ , y sustituyendo en la regla del trapecio obtenemos:

$\begin{eqnarray*} \int\limits_{0}^{\pi^2}\!\sin\left(\sqrt{x}\right)\,dx &\approx& \left(f(x_2) + f(x_3) + \cdot + f(x_{n-1})\right) + \left(\frac{f(x_1) + f(x_n)}{2}\right)\\ &=& \left(\sin\sqrt{0.2\,\pi^2} + \cdots + \sin\sqrt{0.9\,\pi^2}\right) + \frac{\sin(0.1\,\pi^2) + \sin(\pi^2)}{2}\\ &=& 6.0706 \mbox{ unidades de \'area.} \end{eqnarray*}$

Si usamos $n = 20$ obtenemos una mejor aproximación: 6.20862 unidades de área. Por otra parte, para $n = 50$ obtenemos: 6.26451 unidades de área. El área buscada es exactamente: $2\,\pi \approx 6.283185307$ unidades cuadradas.

Ejemplo

Utiliza la regla del trapecio para aproximar:

$\begin{equation*} \int\limits_{0}^{\pi}\!\frac{\sin x}{x}\,dx \end{equation*}$

Geométricamente, tenemos:

Observa que el integrando no está definido para $x = 0$ . Pero eso no es problema, porque podemos acercarnos tanto como queramos a $x=0$ , pero sin llegar a cero. Podemos, por ejemplo, empezar desde $x=0.0001$ y el error cometido al hacer la aproximación es muy pequeño, además de que el integrando sí está definido en ese punto.

Empezamos haciendo $n = 10$ , con lo que obtenemos la primera aproximación: 1.84921 unidades de árera. Para $n = 20$ obtenemos: 1.85118 unidades de área. Para $n = 50$ obtenemos: 1.85173 unidades de área. El valor aproximado de esta integral definida es: 1.851837052 unidades de área.

VER TODO

Add a note

TÚ

Añadir tu comentario

Integración aproximada: Regla del trapecio

Aprenderás a calcular el valor de una Integral Definida de manera aproximada por medio de la regla del trapecio.

Ejemplo

Ejemplo

Ejemplo

Ejemplo

Ejemplo

Ejemplo