Capı́tulo 6: La integral definida para funciones de una variable

Índice
1 La integral definida 1

2 Aplicaciones del cálculo integral 3

3 Métodos numéricos de integración 6

3.0.1 Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

1 La integral definida

0.3
fHxi L

0.2

+ 20 divisiones de @a,bD

0.1

a a b
b
-0.1
- 1.0 1.5 2.0 2.5 xi
Dxi

La integral de una función real entre dos extremos de un intervalo es el área, con signo,
de la zona limitada por la gráfica de la función, el eje de las x y los dos segmentos verticales
que van desde los extremos del intervalo a la gráfica de la función. El signo del área se toma
positivo en las zonas con coordenada y > 0 y negativo en las zonas con coordenada y < 0.
Gráficamente:
Z b
f (x)dx = área con signo (zona rayada)
a
Esta área de la zona rayada se calcula de modo aproximado dividiendo el intervalo [a, b]
en subintervalos [xi , xi+1 ] ( a = x1 < x2 < ... < xN < b) de longitud ∆xi = xi+1 − xi muy
pequeña, calculando el área de los rectángulos de base ∆xi y altura f (xi ) y sumando las
áreas de estos rectángulos. La integral o área rayada se define como el lı́mite de estas sumas
(si existe) cuando ∆xi → 0, esto es:

Z b N  
X Con frecuencia, escribiremos esto con la notación
f (x)dx = lim f (xi )∆xi , Rb P
a ∆xi →0 más condensada a f (x)dx = lim∆x→0 f (x)∆x.
i=1

1
C Integración en una variable 2

o, de un modo más formal, dando a todos los ∆xi el mismo valor,

b N
b−a
Z X
f (x)dx = lim f (xi ) .
a N →∞ N
i=1
Rb
Este lı́mite (y, por tanto, a f (x)dx) existe cuando a, b ∈ R y f es contı́nua, pero también
existe en otros casos en que estas condiciones no se cumplen, por ejemplo:
Z 4
1 3 2
1dx = (1 − 2 3 )
1 (x − 3)3 2

LaRintegral definida verifica las propiedades:

b
i) a λdx = λ(b − a),
Rb Rb Rb
ii) a (λf + µg)dx = λ a f (x)dx + µ a g(x)dx,
Rc Rb Rc
iii) a f (x)dx = a f (x)dx + b f (x)dx, siendo a ≤ b ≤ c.
El teorema fundamental del Cálculo o regla de Barrow relaciona la integral que acabamos
de definir con la primitiva o integral indefinida que definimos en la sección anterior, y dice
lo siguiente:
Teorema fundamental del Cálculo o regla de Barrow. Si f es una función contı́nua
e integrable sobre el intervalo [a, b], entonces
Z b
f (x)dx = F (b) − F (a),
a

siendo F una función diferenciable que verifica F 0 (x) = f (x), es decir, siendo F una primi-
tiva o integral impropia de f .
Idea de la demostración Consideremos la función
Z x
F (x) = f (x)dx,
a

su derivada verificará, teniendo en cuenta la propiedad iii),

R x+∆x Rx R x+∆x
0 F (x + ∆x) − F (x) a f (x)dx − a f (x)dx f (x)dx
F (x) = lim = lim = lim x ,
∆x→0 ∆x ∆x→0 ∆x ∆x→0 ∆x
C Integración en una variable 3

y, como para ∆x pequeño, f varı́a poco en el intervalo [x, x + ∆x], usando la propiedad i)
R x+∆x
se tiene que, aproximadamente, x f (x)dx = f (x)∆x, luego
R x+∆x
f (x)dx f (x)∆x
F 0 (x) = lim x = lim = f (x),
∆x→0 ∆x ∆x→0 ∆x
y, de la propia definición de F (x), se tiene F (a) = 0 y, por tanto,
Z b
F (b) − F (a) = F (b) = f (x)dx,
a
lo que acaba la idea de la demostración del teorema.

2 Aplicaciones del cálculo integral

De la propia definición se deduce que las integrales sirven para calcular áreas. Ası́, el
área de la región limitada por las gráficas de las funciones y = sen 5x y y = − sen 5x entre
x = 0 y x = π/5 (zona rayada de la figura adjunta, es
Z π/5 Z π/5
Area = sen 5xdx − (− sen 5x)dx
0 0
π/5 π/5
cos 5x cos 5x 2 2 4
=− − = − (− ) = .
5 0 5 0 5 5 5

El área de la región limitada por las gráficas de las funciones cos x y sen x desde el primer
x positivo en que se cortan hasta el segundo (ver figura) es
Z 5π/4 Z 5π/4
Area = sin xdx − cos xdx
π/4 π/4
Z5π/4
= (sin x − cos x)dx
π/4
5π/4 √
= − cos x − sin x|π/4 = 2 2
El área de la región R de la figura siguiente, limitada por la elipse
x2 y 2
+ 2 =1
a2 b
la calculamos dividiendo la elipse en dos zonas R+ y R− de modo que cada una de ellas sea
la gráfica de una función, lo que es fácil de conseguir despejando y de la ecuación anterior
de la elipse r
x2
y = ±b 1 − 2 ,
a
C Integración en una variable 4

región R limitada por la elipse región R+ limitada por región R− limitada por
p p
y = b 1 − (x2 /a2 ) y = −b 1 − (x2 /a2 )

Area(R) = Area(R+ ) + Area(R− )

Z a r Z a r ! r
a
x2 x2 x2
Z
= b 1 − 2 dx + −b 1 − 2 dx = 2 b 1− dx,
−a a −a a −a a2

y resolvemos esta última integral haciendo el cambio de variable

r
x x2
= sen u, lo que da 1 − 2 = cos u y dx = a cos udu,
a a
de donde resulta que si x varı́a entre −a y a, entonces u varı́a entre arcsen(−1) = −π/2 y
arcsen 1 = π/2, y
Z ar Z π/2
x2 a π/2 a
1 − 2 dx = a cos2 (u)du = (sen u cos u + u)|−π/2 = π,
−a a −π/2 2 2

que, sustituido en la expresión de Area(R) da

a
Area(R) = 2b π = πab.
2

La integral definida también permite el cálculo del volumen de un sólido de rev-

olución S, que es volumen de R3 limitado por la superficie obtenida por la rotación de la
gráfica de una función f alrededor del eje X, y dos planos de ecuaciones x = a y x = b.

La idea para aplicar el cálculo integral es que, dividiendo el intervalo [a, b] en que está
definida la función f en subintervalos de longitud pequeña ∆x, el volumen que se quiere
calcular vendrá bien aproximado por la unión de cilindros de radio f (x) y altura ∆x, de
modo que
X Z b
2
V olume(S) = lim πf (x) ∆x = πf (x)2 dx.
∆x→0 a
Como ejemplo, podemos calcular el volumen del elipsoide de revolución obtenido por
rotación de la elipse del ejercicio anterior. Por lo que ya vimos, este elipsoide E se puede
C Integración en una variable 5

p
obtener por revolución de la gráfica de la función y = b 1 − (x2 /a2 ). Por lo tanto, aplicando
la fórmula anterior
Z a a
x2 x3
 
4
V olume(E) = πb 1 − 2 dx = πb2 x − 2
2
= πb2 a.
−a a a 3 −a 3
.

También se puede usar la integración para

calcular la longitud de una curva que es la
gráfica Γ de una función f del siguiente modo:
Dividiendo de nuevo el intervalo [a, b] donde ∆s
∆y
esté definida la función en subintervalos de lon-
(x, f (x))
gitud pequeña ∆x, la longitud del segmento de ∆x
gráfica entre (x, f (x)) y (x + ∆x, f (x + ∆x))
viene aproximada por el segmento de recta
∆s entre esos mismos puntos, y los segmen-
tos ∆x, ∆y = f (x + ∆x) − f (x)) y ∆s for-
man un triángulo rectángulo, luego (∆s)2 =
(∆x)2 + ∆y)2 , y
X Xp
Longitud(Γ) = lim ∆s = lim (∆x)2 + ∆y)2
∆x→0 ∆x→0
s s Z bp
X (∆x)2 (∆y)2 X (∆y)2
= lim + ∆x = lim 1+ ∆x = 1 + (f 0 (x))2 dx.
∆x→0 (∆x)2 (∆x)2 ∆x→0 (∆x)2 a

También se puede calcular el área de una superficie de revolución obtenida por

rotación de la gráfica de una función f entre los valores x = a y x = b. Un razonamiento
similar al empleado para el cálculo de volúmenes de revolución, teniendo en cuenta el modo
en que hemos deducido la expresión anterior para la longitud de una curva permite ver que
una tal área viene dada por
Z b p
Area = 2π f (x) 1 + (f 0 (x))2 dx.
a
C Integración en una variable 6

3 Métodos numéricos de integración

A pesar de los muchos métodos de integración que existen, para la mayor parte de las
funciones no es posible encontrar una primitiva. En esos casos, para encontrar un valor de
Z b
f (x) dx no queda más remedio que acudir a métodos numéricos que dan valores aproxi-
a
mados de la misma. Estos métodos numéricos son los que emplean los distintos programas
que usan los ordenadores. Para hacerse una idea de como trabajan esos programas, va-
mos a describir los dos métodos más sencillos, que se basan muy directamente en la misma
definición de integral.
Método de los trapecios. Para calcular la integral entre a y b de una función f , se
divide el intervalo [a, b] en N partes iguales, de longitud h = b−a N ; es decir, se hace la partición
a = x0 < x1 < x2 < ... < xN −1 < xN = b del intervalo [a, b] con xi − xi−1 = h = b−a N , para
1 ≤ i ≤ N , y se toma como valor aproximado de la integral la suma de las áreas Ai de los
trapecios de vértices (xi−1 , 0), (xi , 0), (xi−1 , f (xi−1 )), (xi−1 , f (xi )). El área de cada uno de
estos trapecios es (ver figura)
   
f (xi ) − f (xi−1 ) f (xi−1 ) + f (xi )
Ai = (xi − xi−1 ) f (xi−1 ) + =h ,
2 2
por lo tanto, el valor aproximado de la integral IT (f ) dado por el método de los trapecios es
N  
X f (xi−1 ) + f (xi )
IT (f ) = h
2
i=1

f (x0 ) + f (x1 ) f (x1 ) + f (x2 ) f (x2 ) + f (x3 )
=h + + + ...
2 2 2

f (xN −2 ) + f (xN −1 ) f (xN −1 ) + f (xN )
... + +
2 2
 
f (x0 ) f (xN )
=h + f (x1 ) + f (x2 ) + f (x3 ) + ... + f (xN −2 ) + f (xN −1 ) +
2 2
f (xi )

f (xi )−f (xi−1 )

2
f (xi−1 )

a = x0 x1 x2 x3 x4 x5 x6 = b
xi−1 xi
C Integración en una variable 7

Se puede demostrar usando el desarrollo de Taylor que el error cometido por este pro-
cedimiento al calcular una integral es
1
(b − a) h2 f 00 ,


error aproximado =
12
donde (f 00 ) representa un valor medio de f 00 en [a, b].
Rb
Como ejemplo, vamos a escribir como se harı́a el cálculo de la integral a f (x) dx usando
el método de los trapecios con Mathematica. La regla, para una subdivisión de [a, b] en n
intervalos serı́a:
n−1
X b − a b − a f [a] + f [b] b − a
trapeciorule[n− ] := f [a + i ] + .
n n 2 n
i=1

Si la aplicamos al ejemplo f[x− ]:= Log[Sin[x]], tomando como lı́mites de integración a =

0.11, b = π − 0.1, obtenemos, usando la instrucción de Mathematica R π−0.1
Table[trapeciorule[n],n,2,100], la siguiente lista de valores aproximados de 0.11 ln[sin[x] dx,
para el número de divisiones n del intervalo variando de 2 a 100:
-3.30799, -2.44855, -2.08898, -1.90219, -1.79206, -1.72143, -1.67335, -1.6391, -1. 61382, -
1.59463, -1.57971, -1.56788, -1.55835, -1.55055, -1.54409, -1.53868, -1.5341, - 1.5302, -1.52684,
-1.52393, -1.52139, -1.51917, -1.51721, -1.51547, -1.51392, -1. 51254, -1.5113, -1.51018, -
1.50916, -1.50825, -1.50741, -1.50665, -1.50595, -1.50531, - 1.50472, -1.50418, -1.50368, -
1.50322, -1.50279, -1.50239, -1.50202, -1.50167, -1. 50135, -1.50105, -1.50076, -1.5005, -
1.50025, -1.50001, -1.49979, -1.49959, -1.49939, - 1.4992, -1.49903, -1.49886, -1.49871, -
1.49856, -1.49842, -1.49828, -1.49816, -1. 49803, -1.49792, -1.49781, -1.4977, -1.4976, -1.49751,
-1.49742, -1.49733, -1.49724, - 1.49716, -1.49709, -1.49701, -1.49694, -1.49687, -1.49681, -
1.49675, -1.49669, -1. 49663, -1.49657, -1.49652, -1.49647, -1.49642, -1.49637, -1.49632, -
1.49628, -1.49623, -1.49619, -1.49615, -1.49611, -1.49607, -1.49604, -1.496, -1.49597, -1.49594,
-1.4959, -1.49587, -1.49584, -1.49581, -1.49578, -1.49576.
Si hacemos el cálculo con Mathematica con NIntegrate, obtenemos
NIntegrate[Log[Sin[x]],x,0.11,Pi-0.1]
-1.4944.
Como se ve, para funciones complicadas es necesario tomar un valor grande de n para
encontrar un valor de la integral próximo al real.
C Integración en una variable 8

Método de Simpson. El método de los trapecios usa interpolación lineal para aprox-
imar los valores de f (x) entre cada dos puntos de los extremos de los subintervalos en que
se divide [a, b] (es decir, aproxima f (x) entre xi−1 y xi por los de una función cuya gráfica
es una recta que pasa por los puntos (xi−1 , f (xi )) y (xi−1 , f (xi−1 ))). El método de Simpson
aproxima f por arcos de parábola que pasan por los mismos puntos. Para ello se divide
el intervalo [a, b] en el que se vaya a hacer la integración en 2N partes iguales, es decir, se
considera la partición a = x0 < x1 < x2 < ... < x2N −1 < x2N = b del intervalo [a, b] con
xi − xi−1 = h = b−a 2N . Para cada i impar se define la función aproximación g de la función f
que se desea integrar de la siguiente manera:
g(xi + y) = f (xi ) + a y + b y 2 , con y ∈ [−h, h] y 1 ≤ i ≤ 2N − 1,
donde a y b son las constantes necesarias para que se verifique
f (xi+1 ) = g(xi + h) = f (xi ) + a h + b h2 y
2
f (xi−1 ) = g(xi − h) = f (xi ) − a h + b h .
Sumando las dos expresiones de (??) se obtiene
1
b h2 = (f (xi+1 ) + f (xi−1 ) − 2f (xi )) . (1)
2
Aproximando cada trozo de la integral de f entre xi−1 y xi+1 por la integral Ai entre los
mismos lı́mites de g, tenemos
Z xi+1 Z h
2
f (xi ) + a y + b y 2 dy = 2 h f (xi ) + b h3 ,


Ai = g(x) dx =
xi−1 −h 3
y, sustituyendo (1) en la última expresión,
2 1 1
Ai = 2 h f (xi ) + h (f (xi+1 ) + f (xi−1 ) − 2f (xi )) = h (4 f (xi ) + f (xi+1 ) + f (xi−1 )) ,
3 2 3
Rb
de donde resulta que el valor aproximado IS de a f (x)dx por el método de Simpson es
N N −1 N −1
!
X h X X
IS = A2i+1 = f (x0 ) + f (x2N ) + 4 f (x2i+1 ) + 2 f (x2i ) .
3
i=0 i=0 i=1
Como en el caso del método de los trapecios, se puede demostrar usando el desarrollo de
Taylor que el error cometido por este procedimiento al calcular una integral es
1  
error aproximado = (b − a) h4 f (4) ,
180
(4)


donde f representa un valor medio de la cuarta derivada de f en [a, b].
La regla, usando Mathematica, para calcular una integral de una función f entre los
extremos a y b por el método de Simpson, usando 2n subdivisiones del intervalo [a, b] serı́a
n−1
1b−a X b−a
simpsonrule[n− ] := f [a] + f [b] + 4 f [a + (2i + 1) ]
3 2n 2n
i=0
n−1
!
X b−a
+2 f [a + (2i) ] .
2n
i=1
C Integración en una variable 9

Si aplicamos el método de Simpson al mismo ejemplo al que aplicamos el método de los

trapecios, f[x− ]:= Log[Sin[x]], tomando como lı́mites de integración a = 0.11, b = π −
0.1, obtenemos, usando la instrucción de Mathematica Table[simpsonorule[n],n,2,100], la
R π−0.1
siguiente lista de valores aproximados de 0.11 ln[sin[x] dx, para el número de divisiones 2n
del intervalo variando de 2 a 100:
-2.20534, -1.68265, -1.57322, -1.53481, -1.5177, -1.50893, -1.50398, -1.501, -1.4991, -
1.49785, -1.49698, -1.49638, -1.49594, -1.49561, -1.49537, -1.49518, -1.49504, -1. 49493, -
1.49484, -1.49477, -1.49471, -1.49466, -1.49463, -1.49459, -1.49457, -1.49455, -1.49453, -
1.49451, -1.4945, -1.49449, -1.49448, -1.49447, -1.49446, -1.49445, -1. 49445, -1.49444, -
1.49444, -1.49443, -1.49443, -1.49443, -1.49443, -1.49442, -1.49442, -1.49442, -1.49442, -
1.49442, -1.49441, -1.49441, -1.49441, -1.49441
donde se ve que, en este ejemplo, el método de Simpson da mejores resultados que el de los
trapecios. Generalmente suele ser ası́.

3.0.1 Ejercicios
x2 y 2
(1) Encuentra el área limitada por el arco de la hipérbola a2
− b2 = 1 y la recta x = m > a.
(2) Encuentra la proporción entre las áreas de las regiones en que la curva y = 3x2 divide
al triángulo formado por el eje x, la recta y = λx y la ordenada correspondiente al punto de
intersección de la recta con la curva.
(3) Encuentra el área limitada por la curva x2 = 64y 3 y la recta x = 8y.
(4) Encuentra el área limitada por las curvas y 2 = 4x y y 2 = x + 3. Solución: 8.
(5) Encuentra el área limitada por la curva y = x2 y la recta x = y.
(6) Encuentra el área limitada por la curva 8y = x3 y la recta 2x = y.
(7) Encuentra el área limitada por la curva y 2 = 9x y la recta 2x = y.
(8) Encuentra el área limitada por la curva y = 4x − x2 y la recta y + 2x − 8 = 0.
(9) Encuentra el área limitada por la curva y = 2x2 y la recta 2x = y − 4.
(10) Encuentra el área limitada por las curvas y = x3 y y 2 = x.
(11) Encuentra el área limitada por la curva y 2 = x + 1 y la recta x = y.
(12) Encuentra el área limitada por las curvas y = x2 y y = 2 − x2 .
(13) Encuentra el área limitada por las curvas y = x3 y y = x4 .
(14) Encuentra el volumen del elipsoide de revolución que se obtiene al hacer girar una
elipse de semiejes a y b. Como consecuencia, encuentra el volumen de una esfera de radio a.
(15) Encuentra el volumen de la figura que se obtiene al girar sobre el eje x el área
limitada por las siguientes curvas: y = x3 , y = 0, x = 1, x = 2.