Cálculo IV

Pedro Villalba

Integración numérica

1 Introducción
La integral definida de una función continua f : [a; b] → R en el intervalo
[a; b] ⊂ R,
Z b
I(f ) = f (x)dx
a
es, geométricamente, el área de la región del plano delimitada por la
gráfica de la función, el eje de abscisas y las rectas verticales x = a, x = b.
Como se puede observar en la figura de abajo.

Este concepto matemático tiene, además del cálculo de áreas, numerosas

aplicaciones, como por ejemplo:

• La longitud del arco de la curva y = f (x) comprendido entre los puntos

(a, f (a)) y (b, f (b)) que viene dada por

1
 Z bp
L= 1 + [f ′ (x)]2 dx
a

• La distancia recorrida por un objeto que se mueve con velocidad va-

riable v = v(t) desde el instante t0 hasta el instante T , que viene dada
por:
Z T
S= v(t)dt
t0

• El centro de gravedad C del arco de la curva y = f (x), comprendido

entre los puntos (a, f (a)) y (b, f (b)) que tiene las coordenadas (L es la
longitud de la curva):
Z b Z b
1 p 1 p
xC = x 1 + (y ′ )2 dx; yC = x 1 + (y ′ )2 dx
L a L a

Si se conoce una primitiva F de la función f , es bien sabido que el valor

de la integral definida se puede calcular mediante la Regla de Barrow:
Z b
f (x)dx = F (b) − F (a)
a
En la mayorı́a de los casos, sin embargo, no se puede utilizar esta fórmula,
ya que no se conoce dicha primitiva. Es posible, por ejemplo, que no se
conozca la expresión matemática de la función f , sino sólo sus valores en
determinados puntos. Pero también hay funciones (de apariencia sencilla)
para las que se puede demostrar que no tienen ninguna primitiva que pueda
escribirse en términos de funciones elementales, por ejemplo
Z 1
2
e−x dx
0

La integración numérica es una herramienta de las matemáticas que

proporciona fórmulas y técnicas para calcular aproximaciones de integrales
definidas. Gracias a ella se pueden calcular, aunque sea de forma aproximada,
valores de integrales definidas que no pueden calcularse analı́ticamente y, so-
bre todo, se puede realizar ese cálculo en un ordenador.

Otra forma de encarar el cálculo de la integral definida

2
 Z b
I(f ) = f (x)dx
a
es mediante los conocidos métodos de Newton-Côtes. Este método se basa
en integrar, en lugar de la función f (x), un polinomio de interpolación que
aproxime a f (x) en [a, b]. Se trata por tanto de toda una familia general de
métodos, según el polinomio de interpolación que se considere, puede elegirse
diferentes grados, diferentes puntos para interpolar, etc. Para el caso de las
interpolaciones lineal y cuadrática, estos métodos se denominan Método de
los Trapecios y Método de Simpson, respectivamente,

2 Fórmulas de cuadratura
Las fórmulas que proporcionan una aproximación del valor de una integral
definida se conocen con el nombre de fórmulas de cuadratura.
En sus versiones más sencillas, estas fórmulas aproximan el área bajo la
curva por el área, parecida, de un paralelogramo. Esto sólo proporciona una
buena aproximación si la base del paralelogramo es pequeña. Por ello, las
fórmulas verdaderamente útiles aproximan la integral definida mediante una
suma infinita de áreas de paralelogramos de base pequeña.
En general, las fórmulas de cuadratura se pueden escribir mediante la
fórmula
n
X
∗
I (f ) = αk f (xk )
k=1

variando una de otras en la forma de elegir los puntos {x1 < x2 < · · · <
xn } en el intervalo [a, b] y los coeficientes, también llamados pesos, αk para
k = 1, 2, · · · , n.
Para determinar el grado de exactitud de una fórmula de cuadratura, es
decir el error que se comete al sustituir la integral definida por la suma finita,
Z b n
X
∗
E(f ) = I(f ) − I (f ) = f (x)dx − αk f (xk )
a k=1

se suele utilizar el concepto de orden.

Se dice que una fórmula de cuadratura es de orden m o bien que es exacta
para polinomios de grado m, si el error de dicha fórmula verifica que

E(xk ) = 0; k = 0, 1, · · · , m
E(xm+1 ) ̸= 0;

3
3 Fórmulas de cuadraturas simples
3.1 Fórmula de los rectángulos
La fórmula de cuadratura más sencilla son las fórmulas de los rectángulos.
Z b
f (x)dx = I1 (f ) = (b − a)f (a)
a

En el primer caso se aproxima la integral por el área del rectángulo de base

[a; b] y altura f (a), y
Z b
f (x)dx = I2 (f ) = (b − a)f (b)
a

en el segundo por el de altura f (b).

Es claro que ambas son de orden cero, es decir, exactas para polinomios
constantes, ya que:

4
 Z b
E(1) = dx − I3 (1) = (b − a) − (b − a) = 0
a
Z b h1 i h i
E(x) = xdx − I3 (x) = (b2 − a2 ) − (b − a)a
a 2
2 2
b a
= − − ab + a2
2 2
a2 b2
= − ab +
2 2
1
= (a − b)2 ̸= 0
2

3.2 Fórmula del punto medio

Es similar a las anteriores pero tomando como altura del rectángulo el valor
de f en el punto medio del intervalo.
Z b
a+b
f (x)dx = I3 (f ) = (b − a)f ( )
a 2

Esta fórmula es de orden 1, ya que

Z b
E(1) = dx − I3 (1) = (b − a) − (b − a) = 0
a
Z b h1
i h a + bi
E(x) = (b2 − a2 ) − (b − a)
xdx − I3 (x) = =0
a 2 2
Z b h1 i h1 i
2 2 2 3 3 2 2
E(x ) = x dx − I3 (x ) = (b − a ) − (b − a )(a + b) ̸= 0
a 3 4

5
3.3 Fórmula del trapecio
La fórmula del trapecio se aproxima la integral por el área del trapecio, como
se ve en la figura de abajo. Esta fórmula también es de orden 1 (verificar!).
Z b
b − a 
f (x)dx = I4 (f ) = f (a) + f (b)
a 2

3.4 Fórmula de Simpson

La última de las fórmulas elementales que veremos aquı́ es la fórmula de
Simpson. En ésta, se aproxima la integral de f por el área encerrada bajo un
arco de parábola que coincide con f en tres puntos; los extremos del intervalo
[a, b] y su punto medio.
Z b
b − ah a + b i
f (x)dx = I5 (f ) = f (a) + 4f + f (b)
a 6 2

6
 Razonando como antes, podemos comprobar que la fórmula de cuadratura
de Simpson de orden 3, es decir, es exacta para polinomios de grado 3.

Ejemplo 3.1.

Consideremos la función f (x) = cos 2x, sabemos que

Z 1
1
I(f ) = cos 2xdx = sen(2) = 0, 4546
0 2

• Fórmula del rectángulo:

IR1 (f ) = (1 − 0) cos 2(0) = cos(0) = 1

IR2 (f ) = (1 − 0) cos 2(1) = cos(2) = −0, 4161

• Fórmula del punto medio:

0 + 1
IP M (f ) = (1 − 0) cos 2 = cos(1) = 0, 5403
2

• Fórmula del trapecio:

1 − 0  1 
IT (f ) = cos 2(0) + cos 2(1) = cos(0) + cos(2)
2 2
= 0, 2919

• Fórmula de Simpson:
1 − 0 0 + 1 
IS (f ) = cos 2(0) + 4 cos 2 + cos 2(1)
6 2
1 
= cos(0) + 4 cos(1) + cos(2)
6
= 0, 4575

Es posible proponer otras fórmulas de cuadratura del estilo de las anteri-

ores, por ejemplo utilizando el valor de la función en más puntos, o también,
eligiendo los puntos de manera óptima para conseguir que la fórmula de
cuadratura asociada sea del mayor orden posible.

7
4 Fórmulas de cuadratura compuestas
Cuando el número de puntos aumenta, es decir, n grande, las fórmulas de
cuadratura simples consideradas en la sección anterior, en general no pro-
porcionan aproximaciones muy fiables de la integral. En la práctica, se usan
las fórmulas de cuadratura compuestas, cuya idea de base es descomponer la
integral definida en una suma de integrales sobre subintervalos pequeños y
aplicar las fórmulas anteriores sobre cada uno de ellos:

• Sean {a = x1 < x2 < · · · < xn = b} un conjunto de n puntos del

intervalo [a, b].
• Por la propiedad de las integrales se tiene que:
Z b Z x2 Z x3 Z xn
f (x)dx = f (x)dx + f (x)dx + · · · + f (x)dx
a x1 x2 xn−1
n−1 Z
X xi+1
= f (x)dx
i=1 xi

• Se utiliza alguna de lasRfórmulas elementales, digamos Ik , en cada subin-

x
tervalo, para calcular xii+1 f (x)dx.

Las fórmulas de la sección anterior, dan lugar ası́, a las siguientes fórmulas
compuestas.

4.1 Fórmula de los rectángulos compuestas

Z b n−1 
X 
f (x)dx = xi+1 − xi f (xi ), o bien
a i=1
n−1 
X 
= xi+1 − xi f (xi+1 )
i=1

En el caso particular en que todos los subintervalos tienen la misma lon-

b−a
gitud h, es decir h = , las fórmulas anteriores toman la forma:
n
Z b n−1 
X  n−1
X n−1
X
f (x)dx = xi+1 − xi f (xi ) = hf (xi ) = h f (xi )
a i=1 i=1 i=1

8
 O bien

Z b n−1 
X  n−1
X n−1
X
f (x)dx = xi+1 − xi f (xi+1 ) = hf (xi+1 ) = h f (xi+1 )
a i=1 i=1 i=1

4.2 Fórmula del punto medio compuesta

Z b n−1 
X  x + x 
i i+1
f (x)dx = xi+1 − xi f
a i=1
2

De vuelta, en el caso en que todos los subintervalos son de la misma

b−a
longitud h = = xi+1 − xi , obtenemos:
n
9
 Z b n−1 x + x  n−1 
X i i+1
X xi + xi+1 
f (x)dx = hf =h f
a i=1
2 i=1
2

4.3 Fórmula de los trapecios compuesta

Z b n−1 
X  f (x ) + f (x )
i i+1
f (x)dx = xi+1 − xi
a i=1
2

y en el caso de subintervalos de igual longitud;

Z b n−1 n−1
X f (xi ) + f (xi+1 ) h X 
f (x)dx = h = f (xi ) + f (xi+1 )
a i=1
2 2 i=1
n−1
!
h X
= f (x1 ) + 2 f (xi ) + f (xn )
2 i=2

10
4.4 Fórmula de Simpson compuesta

n−1
!
b
(xi+1 − xi )
Z X x + x 
i i+1
f (x)dx = f (xi ) + 4f + f (xi+1 )
a i=1
6 2

que en el caso de subintervalos de la misma longitud;

n−1
!
Z b X h x + x 
i i+1
f (x)dx = f (xi ) + 4f + f (xi+1 )
a i=1
6 2
n−1
!
hX x + x 
i i+1
= f (xi ) + 4f + f (xi+1 )
6 i=1 2
" n−1 n
#
h X X  xi + xi+1 
= f (x1 ) + 2 f (xi ) + 4 f + f (xn )
6 i=2 i=1
2

Ejemplo 4.1.

Del ejemplo anterior, recordemos que

Z 1
I(f ) = cos(2x)dx = 0, 4546
0

Utilicemos las fórmulas del trapecio y de Simpson, con 5 subintervalos.

• Trapecio:
Recordemos que
n−1
!
h X
IT (f ) = f (x1 ) + 2 f (xi ) + f (xn )
2 i=1

1−0
Como n = 5, a = 0 y b = 1, entonces h = = 0, 2. Denotemos por
5
P la patición de puntos {x1 , x2 , · · · , x6 } de [0, 1], ası́:

P = {0; 0, 2; 0, 4; 0, 6; 0, 8; 1}
| {z }
×2

De esta manera,

11
 Z 1
IT (f ) = cos(2x)dx
0
(
0, 2 h
= f (0) + 2 f (0, 2) + f (0, 4) + f (0, 6)
2
)
i
+f (0, 8) + f (1)

0, 2 
= cos 2(0) + 2 cos 2(0, 2) + cos 2(0, 4) + cos 2(0, 6)
2
!

+ cos 2(0, 8) + cos 2(1)
   
= 0, 1 cos(0) + 2 cos(0, 4) + cos(0, 8) + cos(1, 2) + cos(1, 6) + cos(2)
= 0, 4486

• Simpson: Recordemos que

" n−1 n
#
h X X x + x 
i i+1
IS (f ) = f (x1 ) + 2 f (xi ) + 4 f + f (xn )
6 i=2 i=1
2

Para aplicar la fórmula de Simpson, es bueno considerar dos particiones,

uno
P = {0; 0, 2; 0, 4; 0, 6; 0, 8; 1}
| {z }
×2

y otra partición de los puntos medios xmi de cada subintrvalo [xi , xi+1 ]

PM = {0, 1; 0, 3; 0, 5; 0, 7; 0, 9}
| {z }
×4

Ahora, reemplacemos en la fórmula

12
 Z 1
IS (f ) = cos(2x)dx
0
0, 2 n h i
= f (0) + 2 f (0, 2) + f (0, 4) + f (0, 6) + f (0, 8)
6h i
+4 f (0, 1) + f (0, 3) + f (0, 5) + f (0, 7) + f (0, 9)
o
+f (1)
0, 2 n h i
= cos 2(0) + 2 cos 2(0, 2) + cos 2(0, 4) + cos 2(0, 6) + cos 2(0, 8)
6h i
+4 cos 2(0, 1) + cos 2(0, 3) + cos 2(0, 5) + cos 2(0, 7) + cos 2(0, 9)
o
+ cos 2(1)
1n h i
= cos(0) + 2 cos(0, 4) + cos(0, 8) + cos(1, 2) + cos(1, 6)
30 h i o
+4 cos(0, 2) + cos(0, 6) + cos(1) + cos(1, 4) + cos(1, 8) + cos(2)
= 0, 4546

5 Métodos de Newton-Côtes
Como hemos mencionado, los métodos de Newton-Côtes se basan en la inte-
gración medainte la interpolación. Por ejemplo, el Método de los Trapecios
en un Método de Newton-Côtes basado en la interpolación lineal.
La idea esencial por tanto, de cara a integrar f (x) desde el punto (a, f (a))
hasta (b, f (b)), es aproximar f (x) por su polinomio de interpolación lineal en
[a, b]. En ese sentido, sabemos que, por Lagrance:
x−a x−b
f (x) ≈ P1 (x) = f (a) + f (b)
b−a a−b
de esta manera
b b
b − a
Z Z 
I(f ) = f (x)dx ≈ P1 (x)dx =
f (a) + f (b)
a a 2
En definitiva se trata de aproximar el valor de la integral I por el área
del trapecio que determinan las rectas x = a, x = b, el eje de abscisas y la
recta que une los puntos (a, f (a)) y (b, f (b)), como ya habı́amos visto antes.

13
 Si el intervalo en el que se realiza la integración es grande, el Método
de los Trapecios simple suele ser muy impreciso. Para mejorar la exactitud,
debemos subdividir el intervalo en otros más pequeños y aplicar en cada uno
de ellos el método simple, ası́ como lo hemos hecho para las fórmulas de
cuadratura.
Si dividimos el intervalo [a, b] en n subintervalos obtenemos una partición
b−a
P = {x0 = a, x1 , x2 , · · · , xn = b} de manera que xi+1 − xi = h y h = ,
n
entonces:

Z b Z x1 Z x2 Z xn
f (x)dx = f (x)dx + f (x)dx + · · · + f (x)dx
a x0 x1 xn−1

En cada integral de cada sumando aplicamos el método de los trapecios

simple y obtenemos que:
Z b n−1
!
h X
f (x)dx ≈ f (a) + 2 f (xi ) + f (b)
a 2 i=1

El método de Simpson se puede ver también como un método de Newton-

Côtes de segundo orden, es decir, basado en integrar un polinomio de inter-
polación de segundo grado, de la forma siguiente:
Dada la función f (x) en el intervalo [a, b], tomamos como tercer punto
a+b
para la interpolación el punto medio del intervalo, es decir, xm = , y
2
b−a
denominamos h = a la mitad de la longitud del intervalo [a, b]. De
2
esta forma el polinomio de interpolación de grado 2 que pasa por los untos
(a, f (a)), (xm , f (xm )) y (b, f (b)) es:

f (xm ) − f (a) f (a) + f (b) − 2f (xm )

P2 (x) = f (a) + (x − a) + (x − a)(x − xm )
h 2h2
Al integrar P2 (x) entre a y b, obtenemos:

Z b Z b
f (x)dx ≈ P2 (x)dx
a a
h 
= f (a) + 4f (xm ) + f (b)
3
Ejemplo 5.1.

14
 Calcular el valor aproximado de la integral
Z 1
xdx
0 (x + 1)(x + 2)

1. Utilizando la regla de los trapecios compuesta con n = 8 subintervalos.

2. Utilizando la regla de Simpson compuesta con n = 8.

3. Evaluar exactamente el valor de la integral y compare con los valores

aproximados obtenidos.

Ejercicios 1.

Aproxima el valor de las siguientes integrales definidas por los métodos

del rectángulo, del punto medio, del trapecio y de Simpson, tomando para
todos los casos el mismo valor de h = 0, 1. Calcula el error que se comete en
cada caso en relación con el valor exacto que se proporciona.
Z 3
1. 3x2 dx = 8
0
Resp.: 7,41; 7,995; 8,01; 8
Z 1
2. ex dx = 1, 71828
0
Resp.: 1,6337994; 1,717566086; 1,6337994; 1,661006721
Z 1
3. (x + 2x2 − x3 + 5x4 )dx = 3, 3333
−1
Resp.: 3,3733 ; 3,3133625; 3,3733; 3,333466667
Z 3
1
4. √ dx = 0, 828427
2 x−1
Resp.: 0,843340902; 0,828292655; 0,828696241; 0,828428056
Z 2
2x + 1
5. 2
dx = 1, 09861
1 x +x
Resp.: 1,132685544; 1,098242635; 1,09935221; 1,098615505
Z 1
1
6. 2
dx = 0, 785398
0 1+x
Resp.: 0,809981497; 0,785606496; 0,784981497; 0,785398153

15
 Ejemplo:
Un agricultor desea conocer la superficie aproximada de un prado limitado
por una carretera, dos caminos perpendiculares a ella y la ribera de un rı́o, de
manera que si colocamos unos ejes cartesianos sobre la carretera (eje OX) y
uno de los caminos (eje OY, abscisa x = 0), el segundo camino serı́a la recta
vertical x = 2 (unidades en cientos de metros). Se toman varias medidas
desde la carretera hasta la ribera, obteniéndose las siguientes coordenadas
para los puntos de la ribera: (0;1.5), (0.5;1.8), (1;2.1), (1.5;1.75), (2;1.3).
Calcular aproximadamente el área de dicho terreno utilizando las reglas de
los trapecios y de Simpson. Determinar el área si extendemos el terreno hasta
la abscisa x = 2.5 sabiendo que el rı́o en tal caso pasa por el punto (2.5;1.1).

16