Magister Lord Barrera:

Coordinador del área

de Matemática
Integrales triples

Escuelas Profesionales
de
Ingeniería Industrial
Ingeniería de Sistemas
Ingeniería Empresarial

II
 Lord Barrera

1. Integrales Triples

La Ley de Gravitación Universal es una ley fundamental de la física que in-

volucra masas puntuales, pero cuando queremos calcular la fuerza de atracción
entre la tierra y la luna no podemos aplicarla directamente porque la tierra y
la luna no pueden ser consideradas como puntos. Este problema, abordado por
Newton fué resuelto mediante integración. Un procedimiento natural para resol-
ver este problema consiste en dividir la tierra y la luna en pequeñas partes, ahora
tomamos la fuerza gravitacional entre estas partes pequeñas y añadiendo la con-
tribución de estas partes pequeñas conseguimos la fuerza total. A continuación
ilustraremos este proceso de integración con el siguiente ejemplo:

Con el fín de definir y explicar el concepto de integral, consideremos nueva-

mente el campo gravitacional en * . En lugar de una colección de masas puntuales,
calculemos el campo gravitacional de un punto P de una distribución continua
de masas que son distribuidas en el volumen de la tierra. El punto P es llamado
punto campo. Dividimos la masa grande en N piezas y denotamos por mi a la
masa de la i-ésima pieza en el punto Pi . Para poder escribir la ecuación del campo
para la i-ésima pieza de la masa, debemos asegurarnos que el tamaño mi de la
i-ésima pieza sea suficientemente pequeño. Así escribimos

G  mi
gi ≈ ( r − ri ) .
 r − r i 3

Pi Pi
r - ri r - ri
P P

Cuanto más pequeño sea el tamaño, la expresión se aproxima mejor al campo.

Utilizando el principio de superposición tenemos
n n
G  mi n
G  m ( ri )
g(r) ≈ ∑ g i = − ∑  r − r i 3 ( r − r i ) = − ∑  r − ri 3
( r − ri ) (1.1)
i =1 i =1 i =1
 www.lordbarrera.com.pe

donde en la última igualdad reemplazamos mi por m(ri ). Aparte de un

cambio de notación, esta sustitución hace hincapié en la dependencia de masa
pequeña y su ubicación. En cualquier corte del objeto gravitante, como la Tierra,
cada pieza tiene todavía un tamaño distinto de cero. Esto hace que sea imposi-
ble definir la distancia entre mi y el punto P. Podríamos definir la distancia del
çentro"de mi a P, pero la dificultad se encuentra en la definición del centro de la
pieza. Afortunadamente, resulta que el tamaño de todos los mi son tan peque-
ños en cualquier punto Pi como se muestra en la figura. Esto nos lleva a tomar el
límite de la ecuación (1.1) cuando el tamaño de todas las piezas tiende a cero. Si
este límite existe, lo llamamos la integral del campo gravitatorio y denotamos de
la siguiente manera:

G  m ( ri )
n
G dm(r )
g(r) = lı́m ∑ ( r − ri ) = − (r − r ) (1.2)
 m → 0 i =1  r − r i 3 Ω  r − r  3
N→ ∞

y un procedimiento general lleva a la fórmula del campo electrostático y potencial

 
k e dq(r ) k e dq(r )
E= ( r − r  ), Φ= . (1.3)
Ω  r − r  3 Ω r − r 

Las ecuaciones (1.2) y (1.3) se utilizan con frecuencia en diferentes situaciones

físicas. Tenga en cuenta que las ecuaciones (1.2) y (1.3) son independientes de
cualquier sistema de coordenadas como son todas las leyes físicas.

En las expresiones anteriores, Ω es la región de integración, por ejemplo, el

volumen de la Tierra en el que reside la distribución de madsas y dm(r) es lla-
mado elemento de masa situado en el punto P cuyo vector posición es r . P se
llama punto origen, ya que es la ubicación fuente del campo gravitatorio, es decir,
el elemento de masa en ese punto. También se llama punto de integración. Las
coordenadas de r sobre la cual el elemento de masa depende se puede expresar
se llaman las variables de integración, y lo que multiplica los productos de los
diferenciales de estas variables se llama el integrando.

No es difícil abstraer el concepto de integral a partir del ejemplo concreto de

la gravedad. En lugar de la forma específica de la integral en las ecuaciones (1.2)
y (1.3), usamos f (r, r ), y en lugar del elemento de masa, utilizamos el elemento
de la otra cantidad genérica como dQ(r ). De este modo, escribimos
n 
h(r) = lı́m ∑ f (r, ri )  Q(ri ) = f (r, r ) dQ(r )
 m → 0 i =1 Ω
N→ ∞

donde h(r) es el resultado de la integral la función r indica la posición del vector

P cuyas coordenadas son llamadas parámetros de integración.

2
 Lord Barrera

2. Integrales Triples sobre Cubos

Consideremos el cubo

E = {( x, y, z) : a ≤ x ≤ b, c ≤ y ≤ d, r ≤ z ≤ s}
o también denotado por

E := [ a, b] × [c, d] × [r, s]

Geométricamente:
z

r
a y
c d

El volumen de estre cubo es el número

vol( E) = (b − a)(d − c)(s − r ) .

Una partición de E := [ a, b] × [c, d] × [e, f ] es un conjunto finito

 
P := ( xi , y j , zk ) ∈ E : 0 ≤ i ≤ m, 0 ≤ j ≤ n, y 0 ≤ k ≤ p (2.4)

de puntos en E tal que

a = x0 < x1 < . . . < x m = b

c = y0 < y1 < . . . < y n = d

r = z0 < z1 < . . . < z p = s

3
 www.lordbarrera.com.pe

Los puntos ( xi , y j , zk ), donde 0 ≤ i ≤ m, 0 ≤ j ≤ n y 0 ≤ k ≤ p, son

llamados puntos red. El cubo

Eij k := [ xi−1 , xi ] × [y j−1 , y j ] × [zk−1 , zk ]

es llamado ij k-ésimo subcubo inducido por P . Esto se hace dividiendo el inter-

valo [ a, b] en m subintervalos [ xi−1 , xi ] de igual longitud, dividiendo [c, d] en n
subintervalos [y j−1 , y j ], y el intervalo [r, s] en p subintervalos, todos de igual lon-
gitud, respectivamente. Su volumen se denota por

ΔEij k := Δxi Δy j Δzk

donde

Δxi := xi − xi−1 , Δy j := y j − y j−1 , y Δzk := zk − zk−1 .

Definición 2.1. (Riemann). Sea f : E ⊆ R3 → R una función definida en el cubo

E. Decimos que f es Riemann integrable sobre E si existe el número
 m n p

E
f ( x, y, z) dV := lı́m
m,n,p→ +∞
∑∑∑ f (ri , s j , tk )ΔEij k
i =1 j =1 k =1

para cualquier elección de los puntos (ri , s j , tk ) ∈ Eij k .

Definición 2.2. Sea f : E ⊆ R3 → R una función definida sobre el cubo E :=

[ a, b] × [c, d] × [r, s]. Las integrales iteradas de f sobre E son
 b d s  b  
f ( x, y, z) dz dy dx := f ( x, y, z) dA dx
a c r a [c,d]×[r,s]
 d s b  d  
f ( x, y, z) dx dz dy := f ( x, y, z) dA dy
c r a c [r,s]×[ a,b]

y
 s b d  s  
f ( x, y, z) dy dx dz := f ( x, y, z) dA dz
r a c r [ a,b]×[c,d]

4
 Lord Barrera

Teorema 2.1. (Fubini). Sea f : E ⊆ R3 → R una función continua sobre el cubo

R := [ a, b ] × [ c, d ] × [ r, s ]. Entonces
  b d s
f ( x, y, z) dV = f ( x, y, z) dz dy dx
E a c r
 d s b
= f ( x, y, z) dx dz dy
c r a
 s b d
= f ( x, y, z) dy dx dz.
r a c

Ejemplo 2.2. Calcular la integral triple


xyz2 dV
E

donde
E = {( x, y, z) : 0 ≤ x ≤ 1, −1 ≤ y ≤ 2, 0 ≤ z ≤ 3} .

Solución.
  3 2  1
xyz2 dV = xyz2 dx dy dz
E 0 −1 0
 3 2  1
x2
= yz 2
dy dz
0 −1 2 0
 3 2
yz2
= dy dz
0 2−1
 3  2 2
y
= z2 dz
0 4 −1
 3 
z2
= z2 − dz
0 4
 3 2  3
3z 3 z2 27
= dz = = .
0 4 4 3 0 4

Proposición 2.3. Sean φ : [ a, b] → R, ψ : [c, d] → R y δ : [r, s] → R funciones

continuas y definimos f : [ a, b] × [c, d] × [r, s] → R por f ( x, y, z) = φ( x )ψ(y)δ(z)
para ( x, y, z) ∈ [ a, b] × [c, d] × [r, s]. Entonces f es continua en [ a, b] × [c, d] × [r, s] y
  b   d   s 
f ( x, y, z) dV = φ( x ) dx ψ(y) dy δ(z) dz .
[ a,b]×[c,d]×[r,s] a c r

5
 www.lordbarrera.com.pe

Ejemplo 2.4. Evaluar la integral


xz sen y dV
E

donde E = {( x, y, z) ∈ R3 : 0 ≤ x ≤ 1, 0 ≤ y ≤ π/2, 0 ≤ z ≤ 1}.

Solución. De acuerdo a la proposición anterior tenemos

  1   π/2   1

xz sen y dV = x dx sen y dy z dz
E 0 0 0
 1
π/2  z2 1
x2
= − cos y
2 0 0 2 0
1   1
= × 0 − (−1) ×
2 2
1
= .
4

3. Integrales Triples sobre Regiones Generales

Sea D ⊆ R3 una región acotada y f : D → R una región acotada. Encerramos

a D en un cubo E, o sea, D ⊆ E y definimos

f ( x, y, z) si ( x, y, z) ∈ D
f ( x, y, z) =
0 si ( x, y, z)
∈ D

Decimos que f es integrable sobre D si f es integrable sobre E, en este caso defi-

nimos  
f ( x, y, z) dV = f ( x, y, z) dV .
D E

Así como en  variables, la integrabilidad de la función f y el valor de la inte-

 dos
gral triple f ( x, y, z) dV son independientes de la elección del cubo E conte-
D
niendo a D.

Cuando f ( x, y, z) = 1, entonces el número 1 dV es precisamente el vo-
D
lumen del sólido D.

Si D es un subconjunto acotado de R3 y f : D → R es integrable, también

denotamos la integral triple
 
f ( x, y, z) dV simplemente por f.
D D

6
 Lord Barrera

Proposición 3.1. Sea D un subconjunto acotado de R3 y f , g : D → R funciones

integrables. Sea también a ∈ R y n ∈ N. Entonces
  
(a) f + g es integrable y ( f + g) = f+ g.
D D D
 
(b) a f es integrable y (a f ) = a f.

(c) f g es integrable.

4. Integrando sobre Regiones Elementales

4.1. Regiones elementales de tipo I

El teorema de Fubini es válido para regiones que son más generales que las
regiones rectangulares. Más precisamente, también es válido para regiones que
se describen a continuación:

Definición 4.1. (Región elemental de tipo I). Una región plana es de tipo I si está
entre dos funciones dependiendo de x; es decir
 
D = ( x, y) ∈ R2 : a ≤ x ≤ b, g1 ( x ) ≤ y ≤ g2 ( x )

donde g1 y g2 son funciones continuas en [ a, b].

Algunos ejemplos de regiones de tipo I se muestran en la figura de abajo:

y y = g2(x ( y y
y = g2(x (
y = g2(x (

D D
D

y = g1(x ( y = g1(x ( y = g1(x (

a b x a b x a b x

7
 www.lordbarrera.com.pe

Ejemplo 4.1. Expresar las regiones de las figuras de abajo como regiones del
tipo I.

(i) y (ii) y (iii) y

y = 4 - x2
4 4 4
y = x2

D2 D3
D1 1 1

2 x -2 2 x 1 5 x

Solución. (i) En este caso, la curva que limita superiormente a la región, es

la recta y = 2x, mientras que la curva que limita inferiormente a la región, es la
parábola y = x2 . De esta manera se tiene

D1 : 0≤x≤2 y x2 ≤ y ≤ 2x

(ii) La región D2 se expresa mediante

D2 : −2 ≤ x ≤ 2 y 1 ≤ y ≤ 4 − x2

(iii) Para este caso debemos conocer la ecuación de la recta que pasa por los
puntos (1, 1) y (5, 4). La pendiente de esta recta es

4−1 3
m= =
5−1 4
o sea que la ecuación de dicha recta es

3
y= x+b
4
Desde que (1, 1) es un punto de la recta, tenemos

3 1
1= (1) + b que implica b=
4 4
Por tanto,
3 1
x+y=
4 4
Finalmente, la región D3 se expresa mediante

3 1
D3 : 1≤x≤5 y 1≤y≤ x+
4 4

8
 Lord Barrera

El siguiente resultado nos permite calcular integrales dobles sobre regiones

de tipo I.

Teorema 4.2. (Teorema de Fubini para regiones de tipo I). Si f es una función
continua en una región D de tipo I tal que

D = {( x, y) ∈ R2 : a ≤ x ≤ b, g1 ( x ) ≤ y ≤ g2 ( x )}

entonces   b  g2 ( x )

f ( x, y) dA = f ( x, y) dy dx
D a g1 ( x )

Ejemplo 4.3. Evaluar las siguiente integral


( x + 2y) dA
R

donde R es la región R : 0≤x≤1 e 0≤y≤x

Solución. Claramente la región R es de tipo I. Entonces

  1 x
( x + 2y) dA = ( x + 2y) dy dx
R 0 0
 1
 x
2 
= xy + y  dx
0 0
 1    

= ( x )( x ) + ( x )2 − ( x )(0) + (0)2 dx
0 1
= 2x2 dx
0
1
2x3  2
=  =
3 0 3

Ejemplo 4.4. Evaluar las siguiente integral


( x3 + 2y) dA
R

donde R es la región

R: 0≤x≤2 e x2 ≤ y ≤ 2x

9
 www.lordbarrera.com.pe

Solución. Claramente la región R es de tipo I. Entonces

  2 2x  2
2x
( x + 2y) dA =
3
( x + 2y) dy dx =
3
x y + y  dx
3 2
R 0 x2 0 2 x
 2    

= x3 (2x ) + (2x )2 − x3 ( x2 ) + ( x2 )2 dx
0
  2 
2
x5 4x3 x6 
= ( x + 4x − x ) dx =
4 2
+ − 5
0 5 3 6 0
   
25 4(2)3 26 05 4(0)3 06 32
= + − − + − =
5 3 6 5 3 6 5

4.2. Regiones elementales de tipo II

A continuación estudiaremos regiones de tipo II e integraremos sobre estas

regiones. El tratamiento es similar al caso visto en el párrafo anterior. Para este
caso, la variable y varía entre constante, mientras que la variable x varía entre
dos funciones.

Definición 4.2. (Región elemental de tipo II). Una región plana es de tipo II si
está entre dos funciones dependiendo de y; es decir
 
D = ( x, y) ∈ R2 : c ≤ y ≤ d, h1 (y) ≤ x ≤ h2 (y)

donde h1 y h2 son funciones continuas en [c, d].

Algunos ejemplos de regiones de tipo II se muestran en la figura de abajo:

y y x = h ( y( y
1
d d
d
x = h1( y( x = h1( y(
D D
D
c c
x = h2( y( x = h2( y( c x = h2( y(
x x x

10
 Lord Barrera

Ejemplo 4.5. Expresar las regiones de las figuras de abajo como regiones del
tipo II.

(i) y (ii) y (iii) y

y = 4 - x2
4 4 4
y = x2

D2 D3
D1 1 1

2 x -2 2 x 1 5 x

Solución.

El siguiente resultado nos permite calcular integrales dobles sobre regiones

de tipo I.

Teorema 4.6. (Teorema de Fubini para regiones de tipo II). Si f es una función
continua en una región D de tipo II tal que

D = {( x, y) ∈ R2 : c ≤ y ≤ d, h1 (y) ≤ x ≤ h2 (y)}

entonces    
d h2 ( y )
f ( x, y) dA = f ( x, y) dx dy
D c h1 ( y )

Ejemplo 4.7. Evaluar las siguiente integral


( x + 2y) dA
R

donde R es la región R : 0≤x≤1 e 0≤y≤x

Solución.

11
 www.lordbarrera.com.pe

Ejemplo 4.8. Evaluar las siguiente integral


( x3 + 2y) dA
R

donde R es la región

R: 0≤x≤2 e x2 ≤ y ≤ 2x

Solución. Claramente la región R es de tipo I. Entonces

  2 2x
( x + 2y) dA =
3
( x3 + 2y) dy dx
R 0 x2
 2
2x
= x y + y  dx
3 2
0 2 x
 2    

= x3 (2x ) + (2x )2 − x3 ( x2 ) + ( x2 )2 dx
0 2
= ( x4 + 4x2 − x5 ) dx
0
  2
x5 4x3 x6 
= + −
5 3 6 0
   
25 4(2)3 26 05 4(0)3 06 32
= + − − + − =
5 3 6 5 3 6 5

Ejemplo 4.9. Evaluar las siguiente integral


1
dA
R y4 +1
√
donde R es la región R : 0≤x≤8 e 3
x ≤ y ≤ 2.

Solución. La región R es claramente una

región de tipo I. Su gráfica se muestra en la
figura de la derecha. Por otro lado, la inte-
gral no puede ser evaluada directamente,
por lo que cambiaremos el orden de inte-
gración.

12
 www.lordbarrera.com.pe

Cambiando el orden de integración se tiene

 8  2   2   y3 
dy dx
√
dx = dy
0 3 x y4 + 1 0 0 y +1
4

 2  y3 
x 
= dy
0 y4 + 1 0
 2
y3
= dy
0 y4 + 1

ln(y4 + 1) 2
=  = ln 17
4  4
0

Ejemplo 4.10. Calcular la siguiente integral

  
1 1
2 xy
x e dx dy
0 y

Solución. La región de integración es el

conjunto

D = {( x, y) ∈ R2 : y ≤ x ≤ 1, 0 ≤ y ≤ 1}

como se muestra en la figura:

Cambiando el orden de integración se tiene

  1 1
   1 x 
2 xy
x e dx dy = 2 xy
x e dy dx
0 y 0 0
 1  x 

= xe xy 
0
0 1  
2
= xe x − x dx
0 2  1
ex x2 
= −
2 2 0
e−2
=
2

13
 www.lordbarrera.com.pe

Ejemplo 4.11. Una pieza sólida de rompecabezas tiene la forma de la figura

abajo, la proyección sobre el plano xy tambien se indica a la derecha. Más preci-

z y

R4 R4

y R3
1

R2 R1 x
x -1 1

samente, el sólido es el conjunto

E = {( x, y, z) : ( x, y) ∈ B, y 0 ≤ z ≤ 2 − x2 }

Calcular el volumen del sólido.

Solución. Desde que el sólido tiene como techo a la superficie gráfico de la

función f ( x, y) = 2 − x2 , su volumen viene dado por

volumen = (2 − x2 ) dx dy
R

Notemos que R es la unión de las regiones R1 , R2 , R3 , R4 , R5 . Además

1
B1 : −1 ≤ x ≤ 1, 0≤y≤ ( x + 1)
2
tiene volumen
 1 ( x +1)/2
V1 = (2 − x2 ) dy dx
−1 0
 1
1
= (2 − x2 ) ( x + 1) dx
−1 2
 1
1
= (− x3 − x2 + 2x + 2) dx
2 −1
5
=
3
Por otro lado,
1 1
B2 : −1 ≤ x ≤ 0, ( x + 1) ≤ y ≤
2 2

14
tiene volumen
 0 1/2
V2 = (2 − x2 ) dy dx
−1 ( x +1)/2
 0  
1 1
= (2 − x )
2
− ( x + 1) dx
−1 2 2
3
=
8
El volumen de R3 viene de
 0 5/2
10
V3 = (2 − x2 ) dy dx =
−1 1/2 3

Por simetría, el volumen V5 es igual al volumen V1 y el volumen V2 es igual

al volumen V4 . Luego

5 3 10 3 5 89
V = V1 + V2 + V3 + V4 + V5 = + + + + =
3 8 3 8 3 12

15