Contenidos

Integral simple
Integral doble
Integral triple
Cambio de variables
Aplicaciones
Ejercicios resueltos
Ejercicios propuestos

2 Integración múltiple
El cálculo de integrales dobles y triples es uno de los pilares básicos para el estudio
de la Ingeniería Industrial. El cálculo de áreas y volúmenes, los momentos de inercia,
centros de masa, etc., son problemas a los que se enfrentan con frecuencia los ingenie-
ros. Es por ello, que dedicamos este segundo capítulo al estudio de dichas integrales,
haciendo especial hincapié en el estudio de las regiones elementales del plano y del
espacio, así como en los cambios de variables para integrales dobles y triples.

2.1 Integral simple

R
En el cálculo de funciones de una variable se estudia la integral definida ab f (x) dx
de una función f : [a, b] ⊂ R → R acotada sobre el intervalo [a, b]. El procedimiento
para definir la integral de una función de una variable en un intervalo [a, b], también
conocido como integral de Riemann, es el siguiente:
1. Tomamos una partición P del intervalo [a, b] en n subintervalos [xi−1 , xi ], con
i = 1, . . . , n, es decir,

P := {a = x0 < x1 < · · · < xn = b}.

La longitud de cada subintervalo [xi−1 , xi ] viene dada por ∆xi = xi − xi−1 .

2. Se define la norma de la partición P y se nota por kPk al valor kPk := máx ∆xi .
1≤i≤n

3. Se elige un punto xi∗ arbitrario en cada subintervalo [xi−1 , xi ] y se generan las

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

sumas parciales asociadas a la partición P, dadas por

n
S(P) = ∑ f (xi∗ )∆xi .
i=1

Intuitivamente, esta suma parcial es la suma de las áreas de los rectángulos que
rellenan aproximadamente la región delimitada por la gráfica de la función y el
eje de abscisas, tal y como puede verse en la figura 2.1.

y=f (x) P
n
S(P ) = f (x∗i )∆xi
i=1

f (x∗
i)

∆x k
z}|{
a=x0 x1 x2 xi−1 xi xn =b
x∗
i

Figura 2.1: Sumas parciales asociadas a la partición P.

Entonces, se define
Z b n

a
f (x) dx = lı́m ∑ f (xi∗ )∆xi
kPk→0 i=1
(2.1)

si existe el límite anterior y su valor es independiente de la elección de los puntos

xi∗ , en cuyo caso se dice que f es integrable en [a, b]. En otras palabras, la función es
integrable cuando el límite de las sumas parciales (cuando la anchura de los intervalos
tiende a cero) existe y no depende de las elecciones de los valores. Dicho límite
coincidirá, por definición, con el área con signo de la región delimitada por la gráfica
y el eje de abscisas. Aquí, cuando hablamos de área con signo, estamos haciendo
referencia a que el área se cuenta como positiva cuando se halla por encima del eje
OX y negativa cuando se halla por debajo.

118 Capítulo 2. Integración múltiple

 Apuntes de Ampliación de Matemáticas

Teorema 2.1.1 (Teorema fundamental del cálculo ) Si f : [a, b] → R es conti-

nua, entonces f es integrable en [a, b], y la función F : [a, b] → R, definida por
Z x
F(x) = f (t) dt,
a

es derivable en [a, b], siendo F 0 (x) = f (x).

Observaciones 2.1.2
R
1. Si f (x) ≥ 0 en [a, b], la integral ab f (x) dx es el área de la región limitada por
la gráfica y = f (x), el eje OX y las rectas x = a y x = b (ver la figura 2.2).
2. Para que f sea integrable en [a, b] no es necesario que f sea continua en
[a, b]. Es suficiente que f esté acotada en [a, b] y sea continua a trozos, es
decir, tenga un número finito de discontinuidades de salto en [a, b].

y=f (x)

Z b
D Área(D) = f (x) dx
a

a b

Figura 2.2: Interpretación geométrica de la integral definida.

2.2 Integral doble

2.2.1 La integral de una función de dos variables sobre un rectángulo

El concepto de función integrable puede extenderse fácilmente a funciones acotadas

f : Q ⊂ R2 → R, donde Q es el rectángulo
Q = [a, b] × [c, d] = {(x, y) ∈ R2 : a ≤ x ≤ b, c ≤ y ≤ d}.

2.2 Integral doble 119

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

1. Tomamos una partición P del rectángulo Q en m × n subrectángulos tomando

los puntos

a = x0 < x1 < · · · < xm = b, c = y0 < y1 < · · · < yn = b,

no necesariamente equidistantes. Es decir, P está formado por los subrectángulos

Qi j = [xi−1 , xi ] × [y j−1 , y j ] para 1 ≤ i ≤ m, 1 ≤ j ≤ n.
2. El área de cada subrectángulo Qi j será ∆Ai j = ∆xi ∆y j .
3. Se define la norma de la partición P, y se nota por kPk, al valor máximo de las
diagonales de los rectángulos Qi j , es decir,
q
kPk := máx ∆xi2 + ∆y2j .
1≤i≤m, 1≤ j≤n

4. Se elige un punto (xi∗ , y∗j ) arbitrario en cada rectángulo Qi j y se generan las

sumas parciales asociadas a la partición P, dadas por
m n m n
V (P) = ∑ ∑ f (xi∗ , y∗j )∆Ai j = ∑ ∑ f (xi∗ , y∗j )∆xi ∆y j .
i=1 j=1 i=1 j=1

Cada sumando es el volumen de un prisma con base Qi j y altura f (xi∗ , x∗j ) (ver
figura 2.3). Es interesante observar la analogía con la integral simple.
Entonces, se define
ZZ m n

Q
f (x, y) dx dy := lı́m
kPk→0
∑ ∑ f (xi∗ , y∗j )∆xi ∆y j , (2.2)
i=1 j=1

si el límite anterior existe y su valor es independiente de la elección de los puntos

(xi∗ , y∗j ), en cuyo caso se dice que f es integrable en el rectángulo Q.

Teorema 2.2.1 (Teorema de Fubini) Sea f : Q = [a, b] × [c, d] ⊂ R2 → R una

función continua. Entonces f es integrable en Q y el valor de la integral puede
obtenerse como
ZZ Z b Z d  Z d Z b 
f (x, y) dx dy = f (x, y) dy dx = f (x, y) dx dy.
Q a c c a

Problema 2.2.2 Calcular la integral de la función f (x, y) = 2xy + y2 sobre el rectán-

gulo Q = [0, 1] × [−1, 2].

120 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

∆Ai,j

f (x∗ ∗
i ,yj )

c yj−1 yj d

a=x0

xi−1
∆xi
xi
∆yj

b=xn

Figura 2.3: Sumas parciales de la función f (x, y) asociadas a la partición P.

S OLUCIÓN :
ZZ Z 2 Z 1  Z 2  1 
2 2
(2xy + y ) dx dy = (2xy + y ) dx dy = x 2 y + y2 x 0
dy
Q −1 0 −1
Z 2 2

2 y2 y3
= (y + y ) dy = +
−1 2 3 −1
   
4 8 1 1 9
= + − − = .
2 3 2 3 2
También podríamos haber intercambiado el orden de integración
ZZ Z 1 Z 2  Z 1  2 !
2 2 2 y3
(2xy + y ) dx dy = (2xy + y ) dy dx = xy + dx
Q 0 −1 0 3 −1
Z 1  2 1
3x 3 9
= (3x + 3) dx = + 3x = + 3 = . 
0 2 0 2 2

+ En lo que sigue, para simplificar las expresiones en las integrales iteradas,

utilizaremos la notación:
Z b Z d  Z b Z d
not
f (x, y) dy dx = dx f (x, y) dy
a c a c

2.2 Integral doble 121

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

y, análogamente,
Z d Z b  Z d Z b
not
f (x, y) dx dy = dy f (x, y) dx.
c a c a
Así, el orden de integración viene establecido de derecha a izquierda.

2.2.2 Regiones elementales del plano

Aunque frecuentemente se presenta el problema de integrar una función de dos

variables en un rectángulo, podemos extender el concepto de integral doble a funciones
definidas sobre regiones más generales de R2 .

Definición 2.2.3 (Región de Tipo I) Sea D ⊂ R2 , decimos que D es una región

elemental de Tipo I en el plano XY (ver figura 2.4 (a)) si:

D = {(x, y) ∈ R2 : a ≤ x ≤ b, φ1 (x) ≤ y ≤ φ2 (x)}

con a, b ∈ R y φ1 , φ2 : [a, b] → R funciones continuas tales que

φ1 (x) ≤ φ2 (x), ∀x ∈ [a, b].

La integral doble sobre una región del plano de Tipo I se calcula como
ZZ Z b Z φ2 (x)  Z b Z φ2 (x)
f (x, y) dx dy = f (x, y) dy dx = dx f (x, y) dy.
D a φ1 (x) a φ1 (x)

Definición 2.2.4 (Región de Tipo II) Sea D ⊂ R2 , decimos que D es una región
elemental de Tipo II en el plano XY (ver figura 2.4(b)) si:

D = {(x, y) ∈ R2 : c ≤ y ≤ d, Φ1 (y) ≤ x ≤ Φ2 (y)}

con c, d ∈ R y Φ1 , Φ2 : [c, d] → R funciones continuas tales que

Φ1 (y) ≤ Φ2 (y), ∀y ∈ [c, d].

La integral doble sobre una región del plano de Tipo II se calcula como
ZZ Z d Z Φ2 (y)  Z d Z Φ2 (y)
f (x, y) dx dy = f (x, y) dx dy = dy f (x, y) dx.
D c Φ1 (y) c Φ1 (y)

122 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

y=φ2 (x)

D
x=Φ1 (y) D x=Φ2 (y)

y=φ1 (x)
c

a b
(a) Región Tipo I. (b) Región Tipo II.

Figura 2.4: Regiones elementales del plano.

Definición 2.2.5 (Región de Tipo III) Sea D ⊂ R2 , decimos que D es una región
elemental de Tipo III en el plano XY si es simultáneamente de Tipo I y de Tipo II.

La integral doble sobre una región del plano de Tipo III se puede calcular como (ver
figura 2.5)
ZZ Z b Z φ2 (x) Z d Z Φ2 (y)
f (x, y) dx dy = dx f (x, y) dy = dy f (x, y) dx.
D a φ1 (x) c Φ1 (y)

Lamentablemente, las regiones del plano que nos encontraremos no son siempre de
uno de estos tres tipos. Por ejemplo, basta pensar en la corona circular

D = {(x, y) ∈ R2 : 1 ≤ x2 + y2 ≤ 4}.

Si queremos integrar una función f sobre D, no podemos aplicar ninguna de las

definiciones anteriores, pero no es difícil darse cuenta de que la corona puede descom-
ponerse en regiones del Tipo I ó del Tipo II. El siguiente resultado nos dice que esto
es suficiente para calcular la integral sobre tales dominios.

2.2 Integral doble 123

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

x=Φ1 (y)
y=φ2 (x)

x=Φ2 (y)
y=φ1 (x)
c

a b
Figura 2.5: Región elemental del plano de Tipo III.

Resultado 2.2.6 (Propiedad de aditividad) Sea D ⊆ R2 una región que puede

descomponerse como D = D1 ∪ D2 ∪ . . . ∪ Dk , siendo cada Di una región de Tipo I,
de Tipo II o de Tipo III. Supongamos además que

int(D1 ) ∩ int(D2 ) ∩ . . . ∩ int(Dk ) = 0.

/

Decimos que una función f es integrable en D cuando sea integrable en todas las
subregiones Di , en cuyo caso definimos
ZZ k ZZ
f (x, y) dx dy = ∑ f (x, y) dx dy.
D i=1 Di

En tales condiciones, la integrabilidad de la función y la definición de integral no

dependen de la descomposición {D1 , D2 , . . . , Dk } elegida.

+ Es importante observar que la condición de que los interiores de las subregiones

tienen intersección vacía no quiere decir que dichas subregiones no se solapen,
pero si se solapan, entonces han de hacerlo en la frontera.
Informalmente, lo que estamos diciendo es que las fronteras de los conjuntos
de Tipos I, II o III no importan a la hora de hacer integrales pues son conjuntos
que podríamos llamar “de área cero”.

124 Capítulo 2. Integración múltiple

 Apuntes de Ampliación de Matemáticas

Figura 2.6: Volumen del sólido por encima de la región D limitado por la superficie z = f (x, y).

En lo sucesivo, aunque no lo indiquemos explícitamente, siempre trabajaremos con

dominios que pueden descomponerse en el sentido anterior.

2.2.3 Aplicaciones de la integral doble

Volúmenes: Si f : D ⊂ R2 → R es una función continua y positiva, entonces la

integral doble de la función f sobre la región D nos proporciona el volúmen del
sólido comprendido entre la región D y la superficie z = f (x, y) (ver figura 2.6):
ZZ
V= f (x, y) dx dy.
D

Áreas:
RR
En particular, si f (x, y) = 1 para todo (x, y) ∈ D, entonces la integral
D f (x, y) dx dy nos proporciona el área de la región D:
ZZ
Área(D) = dx dy.
D

2.2 Integral doble 125

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

Problemas 2.2.7
1. Calcular el área del triángulo T de vértices {(0, 0), (2, 0), (0, 1)}, mediante una
integral doble.
S OLUCIÓN : El área del triángulo D vendrá dada por
ZZ
dx dy,
D
donde D es la región de la figura 2.7(a). La ecuación de la recta que une los
vértices (2, 0) y (0, 1) viene dada por
y2 − y1 2−x
y = y1 + (x − x1 ) = .
x2 − x1 2

Por tanto, T = (x, y) ∈ R2 : 0 ≤ x ≤ 2, 0 ≤ y ≤ 2−x 2 puede considerarse una
región de Tipo I (ver figura 2.7(a)) y el área se puede calcular como sigue:
Z 2 Z 2−x
2 1 2
Área(T ) = dx dy = − (2 − x)2 0 = 1.
0 0 4
Observemos que la región D también puede expresarse como
D = {(x, y) ∈ R2 : 0 ≤ y ≤ 1, 0 ≤ x ≤ 2 − 2y},
es decir, es tabién de Tipo II. La integral se calcula como
Z 1 Z 2−2y Z 2  1
Área(T ) = dy dx = (2 − 2y) dy = 2y − y2 0 = 1. 
0 0 0

2. Calcular el área del círculo de centro (0, 0) y radio r.

S OLUCIÓN : La región de integración será ahora (ver figura 2.7(b))
n p p o
D = (x, y) ∈ R2 : −r ≤ x ≤ r, − r2 − x2 ≤ y ≤ r2 − x2 .

Por tanto, la integral se calcula como

ZZ Z r Z √r2 −x2 Z rp
Área(D) = dx dy = dx √ dy = 4 r2 − x2 dx
D −r − r2 −x2 0
" # Z π
x = r sent 2 2
= = 4r cos2 t dt
dx = r cost dt 0
Z π   π
2 2 1 + cos 2t 2 t sen 2t 2
= 4r dt = 4r + = πr2 .
0 2 2 4 0

126 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

√
y= r 2 −x2

−r 0 r
2−x
y = 2

0 y=0 2 √
y=− r 2 −x2
(a) (b)

x
√ 2a
y =
2
√ a2 −x
y=

b
√
b y= a2 − x2
a
√
0 ( 2−1)a
D
−a a

√
−b y = − ab a2 − x2

(c) (d)
y= √
5 25−
x2
4
/3
4x
y=

0 3 5
(e)

Figura 2.7: Representación de las regiones de los problemas 2.2.7.

2.2 Integral doble 127
 Ana M. Lerma, José M. Quesada y Rafael Sánchez

En este caso, por simetría, la región D también se puede considerar como una
región de Tipo II,
n p p o
D = (x, y) ∈ R2 : −r ≤ y ≤ r, − r2 − y2 ≤ x ≤ r2 − y2 . 

2
x2
3. Calcular el área limitada por la elipse a2
+ by2 = 1.

S OLUCIÓN : Considerando la región de Tipo I dada por (ver figura 2.7(c)):

 
bp 2 bp 2
D = (x, y) ∈ R2 : −a ≤ x ≤ a, − a − x2 ≤ y ≤ a − x2 ,
a a

se calcula el área de la elipse como sigue

ZZ Z a Z b √a2 −x2 Z a Z b √a2 −x2
a a
dx dy = dx √ dy = 4 dx dy
D −a − ab a2 −x2 0 0
Z ap  
4b x = a sent
= a2 − x2 dx =
a 0 dx = a cost dt
Z π
2
= 4ab sen2 t dt = πab. 
0

4. Calcular
ZZ
I= y dx dy,
D

donde D es el recinto limitado por la circunferencia x2 + y2 = a2 , la parábola

y2 = 2ax, el eje OY y el semiplano y ≥ 0 (ver figura 2.7(d)).

S OLUCIÓN : Calculamos los puntos de intersección de la circunferencia con la

parábola, es decir, resolvemos el sitema
 2
x + y2 = a2 ,
y2 = 2ax,

sustituyendo la segunda ecuación en la primera, se tiene

x2 + 2ax − a2 = 0.

128 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

Si ahora completamos cuadrados, se tiene que

(x + a)2 − 2a2 = 0
√
y, por tanto, x = a( 2 − 1).
Como se observa en la figura 2.7(d), el dominio D dado por
n √ √ p o
D = (x, y) ∈ R2 : 0 ≤ x ≤ a( 2 − 1), 2ax ≤ y ≤ a2 − x2 .
es una región de Tipo I. Así pues, se tiene
√
Z a(√2−1) Z √a2 −x2 Z a(√2−1)  2  a2 −x2
y
I= dx √ y dy = √
dx
0 2ax 0 2 2ax
Z a(√2−1)
1 a3 √
= (a2 − x2 − 2ax) dx = (4 2 − 5). 
2 0 6
5. Cambiar el orden de integración en la integral iterada
Z 3 Z √25−x2
I= dx f (x, y) dy.
0 4x/3

S OLUCIÓN : De acuerdo con los límites de integración, el recinto de integración

D, es una región de Tipo I dada por
 p 
2 4x 2
D = (x, y) ∈ R : 0 ≤ x ≤ 3, ≤ y ≤ 25 − x .
3
La intersección de la recta y = 43 x con la circunferencia x2 + y2 = 25 para
x ∈ [0, 3] nos da el punto de coordenadas (3, 4). Luego (ver figura 2.7(e)) para
invertir el orden de integración dividiremos la región D de Tipo I en dos regiones
D1 y D2 de Tipo II:
 
2 3y
D1 = (x, y) ∈ R : 0 ≤ y ≤ 4, 0 ≤ x ≤ ,
4
n p o
D2 = (x, y) ∈ R2 : 4 ≤ y ≤ 5, 0 ≤ x ≤ 25 − y2 .

Por tanto,
ZZ ZZ ZZ
I= f (x, y) dx dy = f (x, y) dx dy + f (x, y) dx dy
D D1 D2
Z 4 Z 3y/4 Z 5 Z √ 25−y2
= dy f (x, y) dx + dy f (x, y) dx. 
0 0 4 0

2.2 Integral doble 129

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

2.2.4 Teorema de Guldin

Teorema 2.2.8 (Teorema de Guldin) El volumen engendrado por una región

elemental Ω del plano XY que gira alrededor del eje Y sin que éste la atraviese,
viene dado por
ZZ
V = 2π x dx dy.
Ω

Problema 2.2.9 Calcular el volumen del toro de radios r y R.

S OLUCIÓN : El toro es la superficie de revolución obtenida por la circunferencia

(x − R)2 + y2 = r2 ,

con R > r, al girar alrededor del eje OY (ver figura 2.8).

Aplicando el teorema de Guldin, se tiene que
ZZ
V = 2π x dx dy,
Ω

siendo Ω el círculo de centro (R, 0) y radio r. Dicha región puede considerarse como
una región de Tipo I:
n q q o
Ω = (x, y) ∈ R2 : R − r ≤ x ≤ R + r, − r2 − (x − R)2 ≤ y ≤ r2 − (x − R)2 .

Y
(x−R)2 +y 2 =r 2

r
R X

Figura 2.8: Superficie de revolución alrededor del eje OY .

130 Capítulo 2. Integración múltiple

 Apuntes de Ampliación de Matemáticas

Por tanto,
Z R+r Z √r2 −(x−R)2 Z R+r q
V = 2π x dx √ dy = 4π x r2 − (x − R)2 dx
R−r − r2 −(x−R)2 R−r
Z R+r q Z R+r q
= 4π 2 2
(x − R) r − (x − R) dx + 4πR r2 − (x − R)2 dx
R−r R−r
 R+r Z π
1 3 2
= −4π (r2 − (x − R)2 ) 2 + 8πr2 R cost 2 dt = 2π 2 r2 R. 
3 R−r 0

2.3 Integral triple

De manera análoga a lo visto para la integral doble, comenzaremos definiendo la

integral triple en funciones definidas sobre paralelepípedos de lados paralelos a los
ejes, para posteriormente extender la definición a dominios más generales.
Sea f : Q ⊂ R3 → R donde Q es el paralelepípedo
Q = [a, b] × [c, d] × [e, f ] = {(x, y, z) ∈ R3 : a ≤ x ≤ b, c ≤ y ≤ d, e ≤ z ≤ f }.
Para definir la integral de f sobre Q procedemos de la siguiente manera:
1. Tomamos una partición P del paralelepípedo Q en m × n × p subparalelepípedos
tomando los puntos
a = x0 < x1 < · · · < xm = b, c = y0 < y1 < . . . yn = d,
e = z0 < z1 < · · · < z p = f
no necesariamente equidistantes. El volumen de cada uno de los subparalele-
pípedos Qi jk = [xi−1 , xi ] × [y j−1 , y j ] × [zk−1 , zk ] será ∆Vi jk = ∆xi ∆y j ∆zk .
2. Se define la norma de la partición P, y se nota por kPk, al valor máximo de las
diagonales de los paralelepípedos Qi jk , es decir,
q
kPk := máx ∆xi2 + ∆y2j + ∆z2k .
1≤i, j,k≤n

3. Se elige un punto (xi∗ , y∗j , z∗k ) arbitrario en cada paralelepípedo Qi jk y se generan

las sumas parciales asociadas a la partición P, dadas por
m n p m n p
F(P) = ∑ ∑ ∑ f (xi∗ , y∗j , z∗k )∆Vi jk = ∑ ∑ ∑ f (xi∗ , y∗j , z∗k )∆xi ∆y j ∆zk .
i=1 j=1 k=1 i=1 j=1 k=1

2.3 Integral triple 131

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

Entonces, se define
ZZZ ZZZ m n p

Q
f dV =
Q
f (x, y, z) dx dy dz := lı́m ∑ ∑ ∑ f (xi∗ , y∗j , z∗k )∆xi ∆y j ∆zk ,
kPk→0 i=1 j=1 k=1

si existe el límite anterior y su valor es independiente de la elección de los puntos

(xi∗ , y∗j , z∗k ), en cuyo caso, se dice que f es integrable en el rectángulo Q.

Teorema 2.3.1 (Teorema de Fubini) Sea f : Q = [a, b] × [c, d] × [e, f ] → R una

función continua. Entonces f es integrable en Q y el valor de la integral puede
obtenerse por cualquiera de las seis posibles integrales sucesivas
ZZZ Z b Z d Z f
f (x, y, z) dx dy dz = dx dy f (x, y, z) dz
Q a c e
Z b Z f Z d
= dx dz f (x, y, z) dy = . . .
a e c

2.3.1 Regiones elementales del espacio

Para poder extender el concepto de integral triple a regiones más generales del espacio
que los paralelepípedos rectangulares definiremos tres tipos especiales de regiones de
R3 , que al igual que en R2 llamaremos de Tipo I, II y III.

Definición 2.3.2 (Región de Tipo I) Una región W ⊂ R3 se dice de Tipo I, si

W = {(x, y, z) ∈ R3 : (x, y) ∈ D, γ1 (x, y) ≤ z ≤ γ2 (x, y)}

con D ⊂ R2 una región elemental de R2 (Tipo I, II o III) y γi : D → R, i = 1, 2, son

funciones continuas sobre D.

Si en la definición anterior invertimos los papeles de x con z (o de y con z) estaremos

ante una región de Tipo II (o de Tipo III).
Sea W una región de Tipo I de R3 y sea f : W → R una función continua sobre W ,
ahora se tiene
ZZZ ZZ Z γ2 (x,y) 
f (x, y, z) dx dy dz = f (x, y, z) dz dx dy.
W D γ1 (x,y)

132 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

En el caso de que D sea una región elemental del plano de Tipo I o de Tipo II, podemos
a su vez desarrollar el miembro de la derecha, obteniendo
Z Z φ2 (x) Z γ2 (x,y)
 b


ZZZ  a dx φ (x) dy γ (x,y) f (x, y, z) dz (D Tipo I),
1 1
f (x, y, z) dx dy dz = Z Z Φ2 (y) Z γ2 (x,y)
W 
 d

 dy dx f (x, y, z) dz (D Tipo II).
c Φ1 (y) γ1 (x,y)

Si D es una región que no es de Tipo I, II o III, pero puede descomponerse como

vimos en la propiedad de aditividad (resultado 2.2.6), la integral puede calcularse sin
más que descomponerla en suma de integrales sobre dominios elementales.

Observación 2.3.3 De la definición de integral triple se deduce que si f (x, y, z) =

1, para todo (x, y, z) ∈ W , entonces
ZZZ
dx dy dz = Vol(W ).
W

Problemas 2.3.4
1. Encontrar el volumen de la esfera x2 + y2 + z2 = r2 .

S OLUCIÓN : Consideramos la región W de la esfera situada en el primer octante

(ver figura 2.9), es decir,

W = {(x, y, z) ∈ R3 : x2 + y2 + z2 ≤ r2 , x ≥ 0, y ≥ 0, z ≥ 0}.

La región W puede considerarse como una región de cualquiera de los tres tipos.
Si, por ejemplo, la consideramos como una región de Tipo I tenemos:
n p
W = (x, y, z) ∈ R3 : 0 ≤ x ≤ r, 0 ≤ y ≤ r2 − x2 ,
p o
0 ≤ z ≤ r2 − x2 − y2 .

Entonces, por simetría:

ZZZ Z r Z √r2 −x2 Z √r2 −x2 −y2
V =8 dx dy dz = 8 dx dy dz
W 0 0 0
Z r Z √
r2 −x2 p
=8 dx r2 − x2 − y2 dy.
0 0

2.3 Integral triple 133

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

Figura 2.9: Primer octante de la esfera x2 + y2 + z2 = r2 .

√
r2 − x2 , resulta
Ahora, si llamamos a =
Z r Z ap  
2 2
y = a sent
V =8 dx a − y dy =
0 0 dy = a cost dt
Z r Z π/2 Z r Z π/2
2 2 2 1 + cos 2t
=8 dx a cos t dt = 8 a dx dt
0 0 0 0 2
Z r
4
= 2π (r2 − x2 ) dx = πr3 . 
0 3

2. Sea W la región acotada por los planos x = 0, y = 0, z = 2 y la superficie

z = x2 + y2 , con x ≥ 0, y ≥ 0. Calcular
ZZZ
x dx dy dz.
W

S OLUCIÓN : Esta región es de Tipo I, II y III simultáneamente. Para compro-

barlo (observar la figura 2.10) vamos a plantear los seis diferentes caminos por

134 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

Figura 2.10: Región acotada por los planos x = 0, y = 0, z = 2 y la superficie z = x2 + y2 ,

con x ≥ 0, y ≥ 0.

los que puede integrarse. Finalmente, la resolveremos por dos de ellos que no
sean simétricos.
a) Como región de Tipo I. La proyección de la superficie sobre el plano XY
nos da la región D1 limitada por el primer cuadrante de la circunferencia
x2 + y2 = 1 (ver figura 2.11 (a)). Por tanto, D1 puede escribirse como
n √ p o
D1 = (x, y) ∈ R2 : 0 ≤ x ≤ 2, 0 ≤ y ≤ 2 − x2
n √ p o
= (x, y) ∈ R2 : 0 ≤ y ≤ 2, 0 ≤ x ≤ 2 − y2 .

Fijado (x, y) ∈ D1 , la coordenada z de un punto situado en la región W

satisface x2 + y2 ≤ z ≤ 2. Por tanto, la integral puede calcularse como
ZZZ ZZ Z 2 
I1 = x dx dy dz = dz x dx dy
W D1 x2 +y2
Z √ Z √2−x2 Z 2
2



 0 x dx 0 dy dz =: I11 ,
x2 +y2
= Z√
 Z √2−y2 Z 2


2
 dy x dx dz =: I .
12
0 0 x2 +y2

b) Como región de Tipo II. La proyección de la superficie sobre el plano Y Z

nos da la región D2 limitada por la parábola z = y2 , la recta z = 2 y los

2.3 Integral triple 135

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

(a) Proyección de la superficie z = x2 + y2 sobre el plano XY .

(b) Proyección de la superficie z = x2 + y2 sobre el plano Y Z.

(c) Proyección de la superficie z = x2 + y2 sobre el plano XZ.

Figura 2.11: Proyecciones en los planos coordenados.

136 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

ejes coordenados z = 0 e y = 0 (ver figura 2.11 (b)). Por tanto, D2 puede

escribirse como
 √
D2 = (y, z) ∈ R2 : 0 ≤ y ≤ 2, y2 ≤ z ≤ 2
 √
= (y, z) ∈ R2 : 0 ≤ z ≤ 2, 0 ≤ y ≤ z .

Fijado un punto (y, z) ∈ D2 , la coordenada x de un punto situado en la

región W satisface que
p
0 ≤ x ≤ z − y2 .
Por tanto, la integral puede calcularse como
ZZZ ZZ Z √
!
z−y2
I2 = x dx dy dz = x dx dy dz
W D2 0
Z √ Z 2 Z √z−y2


2

 dy dz x dx =: I21 ,
0 y2 0
=

Z 2 Z √
z Z √z−y2

 dz dy x dx =: I22 .
0 0 0

c) Como región de Tipo III

La proyección de la superficie sobre el plano XZ nos da la región D3
limitada por la parábola z = x2 , la recta z = 2 y los ejes coordenados z = 0
y x = 0 (ver figura 2.11 (c)). Por tanto, D3 puede escribirse como
 √
D3 = (x, z) ∈ R2 : 0 ≤ x ≤ 2, x2 ≤ z ≤ 2
 √
= (x, z) ∈ R2 : 0 ≤ z ≤ 2, 0 ≤ x ≤ z .

Fijado un punto (x, z) ∈ D3 , la coordenada

√ y de un punto situado en la
región W satisface que 0 ≤ y ≤ z − x2 . Por tanto, la integral puede
calcularse como
ZZZ ZZ Z √
!
z−x2
I3 = x dx dy dz = dy x dx dz
W D3 0
Z √ Z 2 Z √z−x2

 2

 x dx dz dy =: I31 ,
0 x2 0
= Z Z √z Z √z−x2

 2

 dz x dx dy =: I32 .
0 0 0

2.3 Integral triple 137

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

Resolvamos ahora la integral por I11 y por I22 :

Z √2 Z √2−x2 Z 2 Z √2 Z √2−x2
I11 = x dx dy dz = x dx (2 − x2 − y2 ) dy
0 0 x2 +y2 0 0
Z √2 √
2x 2 32 8 2
= (2 − x ) dx = ,
0 3 15
Z 2 Z √ Z √ Z 2 Z √z
z z−y2 1
I22 = dz dy x dx = dz (z − y2 ) dy
0 0 0 2 0 0
Z 2 √
1 3/2 8 2
= z dz = . 
3 0 15

2.4 Cambio de variables

En la práctica, el cálculo de algunas integrales dobles y triples puede resultar muy

tedioso y complicado debido a la complejidad de la región de integración o a la
dificultad que conlleva en sí misma la propia integral. Para intentar paliar estos
problemas haremos uso de una técnica llamada cambio de variables.

2.4.1 Cambio de variables en las integrales dobles

Supongamos que D y D∗ son dos regiones del plano y T : D∗ → D dada por

T (u, v) = (x(u, v), y(u, v)),
una función biyectiva y de clase C1 .

T
D∗ D

(u, v) (x, y)

Sea f : D → R una función integrable en D y tomemos la función compuesta f ◦ T :

T f
D∗ −→ D −→ R
(u, v) 7−→ (x(u, v), y(u, v)) 7−→ f (x(u, v), y(u, v))

138 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

Teorema 2.4.1 (Cambio de variables para integrales dobles) Sea T : D∗ → D

una función biyectiva de clase C1 , tal que a un punto (u, v) ∈ D∗ le asocia un punto
(x, y) ∈ D, dado por

x = x(u, v), y = y(u, v).

Entonces para cualquier función f : D → R integrable sobre D se tiene que

ZZ ZZ
∂ (x, y)
f (x, y) dx dy = f (x(u, v), y(u, v)) du dv.
D D∗ ∂ (u, v)

Problemas 2.4.2
1. Sea P el paralelogramo acotado por las rectas y = x, y = x + 1, y = 2x − 2 e
y = 2x. Realizar
RR
un cambio de variables que convierta P en un rectángulo y
calcular P xy dx dy.
S OLUCIÓN : Observemos que el dominio P (ver figura 2.12) se expresa como

P = (x, y) ∈ R2 : x ≤ y ≤ x + 1, 2x − 2 ≤ y ≤ 2x

= (x, y) ∈ R2 : 0 ≤ y − x ≤ 1, −2 ≤ y − 2x ≤ 0 .

Esta última expresión de P nos sugiere el cambio de variables

u = y − x, v = y − 2x,

que transforma el recinto P en el recinto P∗ dado por


P∗ = (u, v) ∈ R2 : 0 ≤ u ≤ 1, −2 ≤ v ≤ 0 .

Despejando x e y en función de u, v, se tiene

x = u − v, y = 2u − v,

y podemos definir la aplicación biyectiva T : P∗ → P, dada por

T (u, v) = (x(u, v), y(u, v)) = (u − v, 2u − v).

El jacobiano será
∂ (x, y) 1 −1
= = 1.
∂ (u, v) 2 −1

2.4 Cambio de variables 139

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

2x
v

y=

x
y=
0 1
u
T
1 P P∗
x+
y=
−2
2x

−2
y=

Figura 2.12: Paralelogramo acotado por las rectas y = x, y = x + 1, y = 2x − 2 e y = 2x

Por tanto, aplicando el teorema del cambio de variables se tiene que

ZZ ZZ
∂ (x, y)
xy dx dy = (u − v)(2u − v) du dv
P P∗ ∂ (u, v)
Z 0 Z 1
= dv (u − v)(2u − v) du = 7. 
−2 0

2. Calcular
ZZ
log(x2 + y2 ) dx dy
D
donde D es la región situada en el primer cuadrante entre los arcos de las
circunferencias x2 + y2 = a2 y x2 + y2 = b2 , con 0 < a < b.
S OLUCIÓN : Un punto P ∈ D puede expresarse de forma unívoca en coordena-
das polares mediante
 h πi
x = ρ cos θ ,
θ ∈ 0, , a ≤ ρ ≤ b.
y = ρ sen θ , 2
Esto nos permite definir una transformación biyectiva T : D∗ → D, tal que
T (ρ, θ ) = (x, y) = (ρ cos θ , ρ sen θ ),

140 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas


donde D∗ = (ρ, θ ) ∈ R2 : a ≤ ρ ≤ b, 0 ≤ θ ≤ π
2 , (ver figura 2.13).
y

b x2
+y 2
=b 2 ρ

b
a T
D D∗
x
2
+
y
2
=

a
a
2

θ
a b 0 π θ
2

Figura 2.13: Cambio a coordenadas polares.

El jacobiano vendrá dado por:

∂ (x, y) cos θ −ρ sen θ
= = ρ > 0.
∂ (ρ, θ ) sen θ ρ cos θ
Luego la integral quedará como sigue:
ZZ ZZ
2 2 ∂ (x, y)
log(x + y ) dx dy = log(ρ 2 ) dρ dθ
D D∗ ∂ (ρ, θ )
ZZ Z π/2 Z b
= ρ log(ρ 2 ) dρ dθ = dθ
2ρ log(ρ) dρ
D∗ 0 a
 
π 1
= b2 log b − a2 log a − (b2 − a2 ) . 
2 2
RR
3. Calcular la integral R (x2 + y2 ) dx dy, donde R es la región del plano XY
limitada por las curvas x2 − y2 = 1, x2 − y2 = 9, xy = 2 y xy = 4, en el primer
cuadrante.
S OLUCIÓN : La región R puede expresarse en la forma (ver figura 2.14)

R = (x, y) ∈ R2 : 1 ≤ x2 − y2 ≤ 9, 2 ≤ xy ≤ 4 .

2.4 Cambio de variables 141

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

1
=
T

2
y
−
x2
v

=9
−y 2
4

x2 R∗
R
x y=4 2
x y=1

1 3 x 0 1 9 u

Figura 2.14: Región situada en el primer cuadrante limitada por las curvas x2 − y2 = 1,
x2 − y2 = 9, xy = 2 y xy = 4.

Lo que nos sugiere el cambio de variables

u = x2 − y2 , v = xy,

que transforma la región R en la región R∗ dada por


R∗ = (u, v) ∈ R2 : 1 ≤ u ≤ 9, 2 ≤ v ≤ 4 .

Dado que las variables x e y no están despejadas en función de u y v, podemos

recurrir al teorema de la función inversa para calcular el jacobiano:

∂ (u, v) 2x −2y ∂ (x, y) 1

= = 2(x2 + y2 ) ⇒ = .
∂ (x, y) y x ∂ (u, v) 2(x + y2 )
2

Sustituyendo este resultado en la integral se tiene finalmente:

ZZ ZZ
∂ (x, y)
(x2 + y2 ) dx dy = (x2 + y2 ) du dv
R R∗ ∂ (u, v)
ZZ
1
= (x2 + y2 ) 2 du dv
R ∗ 2(x + y2 )
Z Z Z 4
1 1 9
= du dv = du dv = 8.
2 R∗ 2 1 2

Observemos que en este caso no ha sido necesario expresar el integrando x2 + y2

en función de u y de v ya que se cancela con el valor del jacobiano. 

142 Capítulo 2. Integración múltiple

 Apuntes de Ampliación de Matemáticas

2.4.2 Cambio de variables en las integrales triples

Teorema 2.4.3 (Cambio de variables para integrales triples) Sea T : B∗ → B

una función biyectiva de clase C1 , tal que a un punto (u, v, w) ∈ B∗ le asocia un
punto (x, y, z) ∈ B, dado por

x = x(u, v, w), y = y(u, v, w), z = z(u, v, w).

Entonces, para cualquier función f : B → R integrable sobre B, se tiene que

ZZZ
f (x, y, z) dx dy dz
B
ZZZ
∂ (x, y, z)
= f (x(u, v, w), y(u, v, w), z(u, v, w)) du dv dw.
B∗ ∂ (u, v, w)

Cambio a coordenadas esféricas

Dado un punto P(x, y, z) en R3 , las coordenadas esféricas (ρ, θ , φ ) se obtienen de la

siguiente forma (ver figura 2.15):
1. Tomamos ρ = d(O, P), la distancia del punto P al origen de coordenadas,
p
ρ = d(O, P) = x2 + y2 + z2 .

2. Proyectamos el punto P(x, y, z) sobre el plano XY y obtenemos el punto P0 (x, y, 0).

Llamamos θ al ángulo que forma el segmento OP0 con la dirección positiva del
eje OX.
3. Finalmente, definimos φ como el ángulo que forma el segmento OP con la
dirección positiva del eje OZ.
De la definición se deduce que
ρ ≥ 0, 0 ≤ θ ≤ 2π, 0 ≤ φ ≤ π.
Por otra parte, en la figura 2.15, se observa que
z = OS = ρ cos φ , OP0 = PS = ρ sen φ .
Por consiguiente,
x = OQ = OP0 cos θ = ρ cos θ sen φ ,
y = OR = QP0 = OP0 sen θ = ρ sen θ sen φ .

2.4 Cambio de variables 143

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

P (x,y,z)
φ
ρ
O
y R Y

θ
x

P ′ (x,y,0)
Q
X

Figura 2.15: Cambio a coordenadas esféricas.

Finalmente, la relación entre coordenadas cartesianas y esféricas viene dada por

x = ρ cos θ sen φ , y = ρ sen θ sen φ , z = ρ cos φ ,

de donde se tiene que

sen φ cos θ −ρ sen φ sen θ ρ cos φ cos θ

∂ (x, y, z)
= sen φ sen θ ρ sen φ cos θ ρ cos φ sen θ = −ρ 2 sen φ .
∂ (ρ, θ , φ )
cos φ 0 −ρ sen φ

Como φ varia entre 0 y π, se tiene que sen φ ≥ 0. Por lo tanto,

∂ (x, y, z)
= ρ 2 sen φ .
∂ (ρ, θ , φ )

Problema 2.4.4 Calcular

ZZZ
2 +y2 +z2 )3/2
e(x dV,
W

donde W es la bola unitaria de R3 .

144 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

P (x,y,z)
φ
ρ

Figura 2.16: Bola unitaria en R3 en el primer octante

S OLUCIÓN : La región W puede expresarse en coordenadas cartesianas como

W = {(x, y, z) ∈ R3 : x2 + y2 + z2 ≤ 1}.
Realizando el cambio de coordenadas cartesianas (x, y, z) a coordenadas esféricas
(ρ, θ , φ ) (ver figura 2.16), se tiene
x = ρ sen φ cos θ , y = ρ sen φ sen θ , z = ρ cos φ .
Con este cambio de variables la región W se transforma en la región

W ∗ = (ρ, θ , φ ) ∈ R3 : 0 ≤ ρ ≤ 1, 0 ≤ θ ≤ 2π, 0 ≤ φ ≤ π .
Por tanto,
ZZZ ZZZ
2 +y2 +z2 )3/2 3 ∂ (x, y, z)
e(x dV = eρ dρ dθ dφ
W W∗ ∂ (ρ, θ , φ )
ZZZ
3
= eρ ρ 2 sen φ dρ dθ dφ
W∗
Z 2π Z π Z 1
3 4
= dθ sen φ dφ eρ ρ 2 dρ = π(e − 1). 
0 0 0 3

2.4 Cambio de variables 145

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

Cambio a coordenadas cilíndricas

Dado un punto P(x, y, z) en R3 , las coordenadas cilíndricas del punto P vienen dadas
por (ρ, θ , z) donde (ρ, θ ) son las coordenadas polares del punto (x, y) en el plano XY
(ver figura 2.17). Se tiene que

x = ρ cos θ , y = ρ sen θ , z = z,

donde
p
ρ= x2 + y2 , 0 ≤ θ ≤ 2π.

Además,

cos θ −ρ sen θ 0
∂ (x, y, z)
= sen θ ρ cos θ 0 = ρ.
∂ (ρ, θ , z)
0 0 1

P (x,y,z)

O
y Y
ρ
θ
x

P ′ (x,y,0)

Figura 2.17: Cambio a coordenadas cilíndricas

146 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

Problemas 2.4.5
1. Calcular el volumen limitado por el cilindro ρ = 2 cos θ , el paraboloide z = ρ 2
(dados en coordenadas cilíndricas) y el plano z = 0.
S OLUCIÓN : Tal y como nos indican en el enunciado, las ecuaciones de ambas
superficies vienen dadas en coordenadas cilíndricas:
x = ρ cos θ , y = ρ sen θ , z = z.
Observemos que x2 + y2 = ρ 2 . Por tanto,
z = ρ 2 = x2 + y2 .
Por otro lado, multiplicando por ρ a ambos lados de la igualdad ρ = 2 cos θ , se
tiene que ρ 2 = 2ρ cos θ , luego
x2 + y2 = ρ 2 = 2ρ cos θ = 2x.
Esto muestra que la ecuación z = ρ 2 se corresponde con el paraboloide z =
x2 + y2 , y la ecuación ρ = 2 cos θ con el cilindro circular (x − 1)2 + y2 = 1.
Sea W la región limitada por el cilindro y el paraboloide (ver figura 2.18). Se
trata, por tanto, de calcular
ZZZ
Vol(W ) = dx dy dz.
W

De la figura, se observa que los puntos situados sobre el sólido verifican

π π
− ≤θ ≤ , 0 ≤ ρ ≤ 2 cos θ , 0 ≤ z ≤ ρ 2.
2 2
Por tanto, la región W se transforma en la región
n π π o
W ∗ = (ρ, θ , z) ∈ R3 : − ≤ θ ≤ , 0 ≤ ρ ≤ 2 cos θ , 0 ≤ z ≤ ρ 2 .
2 2
Por tanto,
ZZZ ZZZ
∂ (x, y, z)
Vol(W ) = dx dy dz = dρ dθ dz
W W∗ ∂ (ρ, θ , z)
ZZZ Z π/2 Z 2 cos θ Z ρ2
= ρ dρ dθ dz = 2 dθ ρ dρ dz
W∗ 0 0 0
Z π/2 3
=8 cos4 θ dθ = π. 
0 2

2.4 Cambio de variables 147

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

y
P (x,y,0)
sθ
2 co
ρ=

θ
1 x

x2 +y 2 =2x

Figura 2.18: Superficie limitada por el paraboloide z = ρ 2 y el cilindro ρ = 2 cos θ . Región

de integración sobre el plano XY .

2. Hallar el volumen determinado por el paraboloide 4z = x2 + y2 , el cilindro

x2 + y2 = 8y y el plano z = 0.
S OLUCIÓN : Sea W la región limitada por el cilindro y el paraboloide (ver
figura 2.20). Se trata, por tanto, de calcular
ZZZ
Vol(W ) = dx dy dz.
W

El cambio a coordenadas cilíndricas

x = ρ cos θ , y = ρ sen θ , z = z,
actuando sobre el paraboloide nos proporciona
4z = x2 + y2 = ρ 2 .

Por tanto, el paraboloide en coordenadas cilíndricas viene dado por z = 14 ρ 2 .

Actuando de la misma forma sobre el cilindro x2 + y2 = 8y, se tiene
ρ 2 = 8ρ sen θ ,

148 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

con θ ∈ [0, π] (ver figura 2.19) y, simplificando, el cilindro en coordenadas

cilíndricas viene dado por
ρ = 8 sen θ , θ ∈ [0, π] .
Con este cambio de variables, la región W se transforma en la región

x2 +y 2 =8x

θ P (x,y,0)
sen
ρ =8
θ
4 x

Figura 2.19: Región limitada por el paraboloide 4z = x2 + y2 y el cilindro x2 + y2 = 8y.

Región de integración sobre el plano XY .

 
∗ 3 1 2
W = (ρ, θ , z) ∈ R : 0 ≤ θ ≤ π, 0 ≤ ρ ≤ 8 sen θ , 0 ≤ z ≤ ρ .
4

∂ (x,y,z)
Además, ∂ (ρ,θ ,z) = ρ, luego
ZZZ ZZZ
∂ (x, y, z)
Vol(W ) = dx dy dz = dρ dθ dz
W W∗ ∂ (ρ, θ , z)
ZZZ Z π Z 8 sen θ Z 1 ρ2
4
= ρ dρ dθ dz = dθ ρ dρ dz
W∗ 0 0 0
Z π
= 256 sen4 θ dθ = 96π. 
0

2.4 Cambio de variables 149

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

2.5 Aplicaciones de las integrales múltiples

Las integrales dobles y triples tienen su aplicación fundamental, como hemos visto,
en el cálculo de áreas y volúmenes. No obstante, existen otras aplicaciones en la
Matemática y en la Física, como son el cálculo de valores medios, centros de masas y
momentos de inercia.

2.5.1 Valor medio de una función

Si f : D ⊂ R2 → R se define el valor medio de f en D como

ZZ
f (x, y) dx dy
WZZ
fD = .
dx dy
W

Para el caso de funciones de 3 variables, f : W ⊂ R3 → R el valor medio de f en W

viene dado por
ZZZ
f (x, y, z) dx dy dz
W
ZZZ
fW = .
dx dy dz
W

Problema 2.5.1 Sea W el cuerpo de R3 limitado por la superficie (x2 + y2 + z2 )2 =

x2 + y2 . Sea ρ : W → R, con ρ(x, y, z) = (x2 + y2 + z2 )3/2 , la función densidad en W .
Calcular la densidad media de W .

S OLUCIÓN : Realizando el cambio a coordenadas esféricas,

x = ρ sen φ cos θ , y = ρ sen φ sen θ , z = ρ cos φ ,

la ecuación de la superficie se transforma en

ρ = sen φ , 0 ≤ θ ≤ 2π, 0 ≤ φ ≤ π.

El hecho de que el ángulo θ no aparezca en la ecuación de la superficie en coordenadas

esféricas, ρ = sen φ , nos indica que se trata de una superficie de revolución. El corte de
la superficie con el plano Y Z (x = 0) nos proporciona la circunferencia z2 + y2 − y = 0,

150 Capítulo 2. Integración múltiple

Apuntes de Ampliación de Matemáticas

Z
z
P (0,y,z)

φ φ
n
se
ρ=
1 y
2

Figura 2.20: Gráfica de la superficie (x2 + y2 + z2 )2 = x2 + y2 .

de centro (0, 12 , 0) y radio 12 . Dicha circunferencia genera la superficie por revolución

alrededor del eje Z (ver figura 2.20).
La densidad media del cuerpo vendrá dada por
ZZZ
ρ(x, y, z) dx dy dz
W
ZZZ
ρW = .
dx dy dz
W

Calculamos la integral del denominador

ZZZ Z 2π Z π Z sen φ Z π
2 1
dx dy dz = dθ sen φ dφ ρ 2 dρ = π sen4 φ dφ = π 2
W 0 0 0 3 0 4
y la del numerador
ZZZ Z 2π Z π Z sen φ
ρ(x, y, z) dx dy dz = dθ sen φ dφ ρ 5 dρ
W 0 0 0
Z π
2 7 32π
= π sen φ dφ = .
6 0 105
Finalmente,
ZZZ
ρ(x, y, z) dx dy dz 32π
WRRR 105 128
ρW = = π2
= . 
W dx dy dz 105π
4

2.5 Aplicaciones 151

 Ana M. Lerma, José M. Quesada y Rafael Sánchez

2.5.2 Centro de masa

Si D ⊂ R2 representa una palanca bidimensional con densidad de masa ρ(x, y), defini-
mos los momentos de masa de D como
ZZ ZZ
Mx = xρ(x, y) dx dy, My = yρ(x, y) dx dy,
D D

mientras que la masa de D es la integral de la densidad:

ZZ
M= ρ(x, y) dx dy.
D

Se define el centro de masa de la palanca como

 
Mx My
(x, y) = , .
M M

Problema 2.5.2 Hallar el centro de masa del rectángulo D = [0, 1] × [0, 1] de densi-
dad de masa ρ(x, y) = ex+y .

S OLUCIÓN : Calculamos en primer lugar la masa, que viene dada por

ZZ Z 1 Z 1
M= e x+y
dx dy = x
e dx ey dy = (e − 1)2 .
D 0 0

A continuación, calculamos los momentos de masa

ZZ Z 1 Z 1
My = Mx = xex+y dx dy = ey dy xex dx = e − 1.
D 0 0

Por último, usando lo calculado previamente, se obtiene el centro de masa como sigue
 
1 1
(x, y) = , . 
e−1 e−1

Para regiones en el espacio con densidad de masa ρ(x, y, z) estas fórmulas se genera-
lizan fácilmente a tres dimensiones obteniéndose
ZZZ ZZZ
M= ρ(x, y, z) dx dy dz, Mx = xρ(x, y, z) dx dy dz, . . .
W W
 
Mx My Mz
y (x, y, z) = M, M, M .

152 Capítulo 2. Integración múltiple

 Apuntes de Ampliación de Matemáticas

2.5.3 Momentos de inercia

Si el sólido W tiene una función de densidad ρ(x, y, z), el momento de inercia alrededor
del eje X se define por
ZZZ
Ix = (y2 + z2 )ρ(x, y, z) dx dy dz,
W

alrededor del eje Y por

ZZZ
Iy = (x2 + z2 )ρ(x, y, z) dx dy dz
W

y alrededor del eje Z por

ZZZ
Iz = (x2 + y2 )ρ(x, y, z) dx dy dz.
W

Problema 2.5.3 Calcular el momento de inercia con respecto del eje Z del sólido
situado sobre el plano XY acotado por el paraboloide z = x2 + y2 y el cilindro
x2 + y2 = a2 con densidad de masa constante ρ = k.

S OLUCIÓN : Realizando un cambio a coordenadas cilíndricas, resulta

ZZZ Z 2π Z a Z ρ2
kπa6
Iz = k (x2 + y2 ) dx dy dz = k dθ ρ 3 dρ dz = . 
W 0 0 0 3

2.5 Aplicaciones 153