Integral numérica

Introducción
En los cursos de Cálculo integral, nos enseñan como calcular una integral definida de una función
continua mediante la aplicación del Teorema Fundamental del Cálculo.

Teorema Fundamental del Cálculo

Sea 𝑓(𝑥) una función continua y definida en el intervalo [𝑎, 𝑏] y sea 𝐹(𝑥) una función primitiva de
𝑏
𝑓(𝑥), entonces: 𝐼 = ∫𝑎 𝑓(𝑥)𝑑𝑥 = 𝐹(𝑏) − 𝐹(𝑎)

El problema en la práctica se presenta cuando nos vemos imposibilitados de encontrar la función

1 2
primitiva requerida, aún para integrales aparentemente sencillas como: 𝐼 = ∫0 𝑒 𝑥 𝑑𝑥 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 precisas a integrales como la mencionada anteriormente.

Regla del Punto Medio

Usar la fórmula del punto medio para integrar una función en un intervalo, equivale a sustituir,
dentro de la función integral, la función a integrar por el polinomio de interpolación de grado cero,
es decir, una recta horizontal que pasa por el punto medio de la función (la altura del área del
rectángulo que estamos considerando). Sustituimos la función por una recta e integramos.
Estamos entonces calculando el área de un rectángulo.

Sea f una función continua en [a,b]. La regla del punto medio para aproximar su integral viene
𝑏 𝑏−𝑎 𝑛 𝑏−𝑎
dada por: ∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑𝑖=1 𝑓(𝑥̂𝑖 ) = [𝑓(𝑥̂1 ) + ⋯ + 𝑓(𝑥
̂ 𝑛 )] donde 𝑥𝑖 es el punto medio
𝑎 𝑛 𝑛
1
del i-ésimo subintervalo [𝑥𝑖−1 , 𝑥𝑖 ], es decir, 𝑥𝑖 = 2 (𝑥𝑖−1 , 𝑥𝑖 )

Aquí en el dibujo

MDM. Carmen Sosa Garza

Cálculo Integral
 xi 1  xi
ci  xi 1 cj  xi cij  es decir punto medio de xi y xj
2

x x 
Por lo que se busca es el área del rectángulo Ai de base xi  xi 1  x y altura f  i i 1 
 2 
Es decir
n n
x x 
f ( x)dx   A i   f  i i 1  ( xi  xi 1 )
b
a
i 1 i 1  2 

  f  x*i  x
n
donde x*i punto medio de xi 1 y xi
i 1

ba
  f  x*i 
n

i 1 n
ba n
  f  x*i 
n i 1
ba
  f ( x*1 )  f ( x*2 )  ...  f ( x*n ) 
n

Ejemplo 1
21
Utilizar la regla del punto medio para aproximar ∫1 𝑥 𝑑𝑥 con 𝑛 = 5
𝑏−𝑎
Donde ℎ = 𝑛
, y esto es ℎ = 0.2 (en algunos libros también se le conoce como ∆𝑥)
1 1 1 1
Luego 𝑓(1.1) = 1.1 =0.9 , 𝑓(1.3) = 1.3 = 0.76, 𝑓(1.5) = 1.5 = 0.6, 𝑓(1.7) = 1.7 = 0.58,
1
𝑓(1.9) = 1.9 = 0.52

Cada 0.2 se va marcando una línea es decir 1.2, 1.4,1.6,1.8 y 2 pero por la regla del punto medio
cuando hagamos un rectángulo, no tomaremos la altura de 1.2 ni la de 1.4, 1.6,1.8 ni mucho
MDM. Carmen Sosa Garza
Cálculo Integral
menos 2, la altura será la mitad del segmento FA, AB, BC, CD y DE, que es 1.1, 1.3, 1.5, 1.7, 1.9
respectivamente. (son los segmentos GQ, HM, IN, JO y KP, respectivamente)
Ahora, lo que falta hacer es sacar el área de cada rectángulo y sumarlas:

Sabemos que la base es 0.2, evaluamos los puntos medios porque es donde llega la altura
1 1 1 1
𝑓(1.1) = 1.1 =0.9 , 𝑓(1.3) = 1.3 = 0.76, 𝑓(1.5) = 1.5 = 0.66, 𝑓(1.7) = 1.7 = 0.58, 𝑓(1.9) =
1
1.9
= 0.52,
Luego A≈ 0.2(0.9 + 0.76 + 0.66 + 0.58 + 0.52), A≈ 0.6919
La integral original nos da el A≈ 0.6931.

Ejemplo 2

1 2
Utilizar la regla del punto medio para aproximar ∫0 𝑒 𝑥 𝑑𝑥 con 𝑛 = 4

Solución
1
2 1−0 0 + 0.25 0.25 + 0.5 0.5 + 0.75 0.75 + 1
∫ 𝑒 𝑥 𝑑𝑥 ≈ [𝑓 ( )+𝑓( )+𝑓( )+𝑓( )]
4 2 2 2 2
0

1
= [𝑓(0.125) + 𝑓(0.375) + 𝑓(0.625) + 𝑓(0.875)]
4
1 2 2 2 2
= [𝑒 0.125 + 𝑒 0.375 + 𝑒 0.625 + 𝑒 0.875 ]
4
1
= [1.01574 + 1.15099 + 1.4799 + 2.1503]
4
1
= [5.79498] ≈ 1.448745
4

REGLA DEL TRAPECIO

𝑏 𝑏
Corresponde al caso donde n=1, es decir: 𝐼 = ∫𝑎 𝑓(𝑥)𝑑𝑥 ≈ ∫𝑎 𝑃1 (𝑥)𝑑𝑥 donde 𝑃1 (𝑥) es un
polinomio de grado 1.

En el gráfico trazamos la recta que une los puntos: (𝑎, 𝑓(𝑎)) 𝑦 (𝑏, 𝑓(𝑏)) obteniendo un trapecio
cuya superficie será, aproximadamente, el valor de la integral 𝐼.

MDM. Carmen Sosa Garza

Cálculo Integral
 𝑏 𝑏 𝑏 𝑓(𝑎)+𝑓(𝑏)
Así tendremos : 𝐼 = ∫𝑎 𝑓(𝑥)𝑑𝑥 ≈ ∫𝑎 𝑃1 (𝑥)𝑑𝑥 = ∫𝑎 (𝑎0 + 𝑎1 𝑥)𝑑𝑥 = (𝑏 − 𝑎) 2
conocida
como la Regla del Trapecio.

El error que se llega a cometer con esta forma de aplicación puede ser significativo. Una mayor
aproximación se obtiene dividiendo el intervalo de integración en subintervalos y aplicando en cada
uno de ellos la regla trapecial. A este procedimiento se le conoce como Regla trapecial compuesta.

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.

MDM. Carmen Sosa Garza

Cálculo Integral
ci  xi 1 cj  xi Ci  f (ci) y Cj  f (cj )

Por lo que se busca es el área del trapecio Ti de base

xi  xi 1  x y base mayor y menor f  xi 1  y f ( xi )

Es decir
n n
 f ( xi 1 )  f ( xi ) 
f ( x)dx   T i   
b
a
i 1 i 1  2  ( xi  xi 1 )

n
 f ( xi 1 )  f ( xi ) 
   x donde x*i punto medio de xi 1 y xi
i 1  2 
n
 f ( xi 1 )  f ( xi )  b  a
   n
i 1  2 
  f ( x1 )  f ( x2 )   b  a   f ( x2 )  f ( x3 )   b  a   f ( xn1 )  f ( xn )   b  a  
      ....   
 2  n  2  n  2  n 
ba
  f ( x1 )  f ( x2 )    f ( x2 )  f ( x3 )   ....   f ( xn 1 )  f ( xn )  
2n 
ba
  f ( x1 )  2 f ( x2 )  2 f ( x2 )  ...  2 f ( xn1 )  f ( xn )
2n

Ejemplo 3
1 2
Utilizar la regla del trapecio para aproximar la integral: ∫0 𝑒 𝑥 𝑑𝑥

Solución: Usamos la fórmula directamente con los siguientes datos:

2
𝑎 = 0, 𝑏 = 1 𝑦 𝑓(𝑥) = 𝑒 𝑥

1 2 𝑓(0)+𝑓(1) 1+𝑒
Por lo tanto, tenemos que: ∫0 𝑒 𝑥 𝑑𝑥 ≈ (1 − 0) [ ] = ≈ 1.85914
2 2

Ejemplo 4
4 𝑒𝑥
Usar la regla del trapecio para aproximar la integral: ∫2 𝑥
𝑑𝑥

Solución: Igual que en el ejemplo anterior, sustituimos los datos:

MDM. Carmen Sosa Garza

Cálculo Integral
 𝑒𝑥
𝑎 = 2, 𝑏 = 4 𝑦 𝑓(𝑥) =
𝑥
4 𝑒𝑥 𝑓(2)+𝑓(4) 𝑒2 𝑒4
Por lo tanto, tenemos que: ∫2 𝑥
𝑑𝑥 ≈ (4 − 2) [ 2
] = 2
+ 4
≈ 17.3441

La regla del trapecio se puede ampliar si subdividimos el intervalo [𝑎, 𝑏] en 𝑛 subintervalos todos
𝑏−𝑎
de la misma longitud ℎ = 𝑛
.

Sea 𝑃 = {𝑥0 , 𝑥1 , … , 𝑥𝑛 } la partición que se forma al hacer dicha subdivisión. Usando propiedades
𝑏 𝑥 𝑥 𝑥
de la integral tenemos que: ∫𝑎 𝑓(𝑥)𝑑𝑥 = ∫𝑥 1 𝑓(𝑥)𝑑𝑥 + ∫𝑥 2 𝑓(𝑥)𝑑𝑥 + … + ∫𝑥 𝑛 𝑓(𝑥)𝑑𝑥
0 1 𝑛−1

Aplicando la regla del trapecio en cada una de las integrales obtenemos:

𝑏 𝑓(𝑥0 )+𝑓(𝑥1 ) 𝑓(𝑥𝑛−1 )+𝑓(𝑥𝑛 )
∫𝑎 𝑓(𝑥)𝑑𝑥 ≈(𝑥1 − 𝑥0 ) [ 2
]+ ⋯ + (𝑥𝑛 − 𝑥𝑛−1 ) [
2
]

Ahora bien, ya que todos los subintervalos tienen la misma longitud ℎ, tenemos que:
𝑏
ℎ
∫ 𝑓(𝑥)𝑑𝑥 ≈ [𝑓(𝑥0 ) + 2𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + ⋯ + 2𝑓(𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )]
𝑎 2
(𝑏−𝑎)
Sustituyendo el valor de ℎ = 𝑛
, tenemos finalmente:
𝑏
(𝑏 − 𝑎)
∫ 𝑓(𝑥)𝑑𝑥 ≈ [𝑓(𝑥0 ) + 2𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + ⋯ + 2𝑓(𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )]
2𝑛
𝑎

Esta es la regla del trapecio para 𝑛 subintervalos. Cuantos más subintervalos se usen, mejor será
la aproximación a la integral, hasta que la importancia de los errores por redondeo comience a
tomar relevancia.

Ejemplo 5
2
Aplicar la regla del trapecio para aproximar la integral: ∫0 𝑠𝑒𝑛 𝑥 2 𝑑𝑥 si subdividimos en 4
intervalos.

Solución
𝑏−𝑎 2−0
En este caso, tenemos que 𝑛 = 4 y ℎ = 𝑛
, en este caso b= 2, a=0, luego se tiene que ∆𝑥 = 4
=0.5
2
2−0
∫ 𝑠𝑒𝑛 𝑥 2 𝑑𝑥 ≈ [𝑓(0) + 2𝑓(0.5) + 2𝑓(1) + 2𝑓(1.5) + 𝑓(2)]
2(4)
0

2−0
[𝑓(0) + 2𝑓(0.5) + 2𝑓(1) + 2𝑓(1.5) + 𝑓(2)] ≈
2(4)

MDM. Carmen Sosa Garza

Cálculo Integral
 2
≈ [𝑠𝑒𝑛 02 + 2(𝑠𝑒𝑛 (0.5)2 ) + 2(𝑠𝑒𝑛 (1)2 ) + 2(𝑠𝑒𝑛 (1.5)2 ) + 𝑠𝑒𝑛(2)2 ]
8
2
≈ [0 + 0.494 + 1.68 + 1.55 − 0.756]
8
2
≈ [2.968] ≈ 0.742𝑢2
8

Ejercicios:

Usar la regla del Punto Medio y la del Trapecio: comparar con el valor exacto de la integral.
3
1) 
1
x 2 dx n6
1 2
2)   x  2
0 2
dx n4

2
3)  0
x x2  1dx n4

Usar la regla del Punto Medio y la del Trapecio:

1
4) 
0
x 1  xdx n4

 /4
5) 
0
x tan xdx n4

2 1
6) 
0
1  x3
dx n8

MDM. Carmen Sosa Garza

Cálculo Integral
MDM. Carmen Sosa Garza
Cálculo Integral