Capı́tulo 6

Integrales Múltiples

“God does not care about our mathematical difficulties. He integrates empirically.”
- Albert Einstein 6.1

Queremos seguir reconstruyendo el cálculo que conocemos pero en más variables. Ahora que el lector
domina las técnicas de diferenciación en Rn , el paso siguiente es la integración.
Como el lector supone, cuando n = 1 la integral no es realmente “múltiple”, sino que es la ya conocida
integral de Riemann sobre R. Bien, como buenos matemáticos, haremos un poco de inducción antes de
saltar al caso general para un natural cualquiera. Primero estudiaremos el caso n = 2, luego para n = 3
y después generalizaremos la integral de Riemann a Rn .
El lector podrá notar que los procesos de construcción de la integral de Riemann son siempre iguales y por
lo tanto, es útil partir con una situación más simple para comprender ı́ntegramente lo que está ocurriendo
(además del hecho que visualizar las cosas puede ser un poco más difı́cil cuando n > 3).
Queremos hacer hincapié en que existe más de un sentido de integración, pero en este apunte sólo nos
limitaremos a la integral de Riemann. Lamentablemente, esta integral es la más picante de todas y se
invita al lector a tomar el curso de Teorı́a de la Integración, donde se construye una integral muchı́simo
más poderosa que la que estudiaremos nosotros —hablamos de la integral de Lebesgue. No es que
no le tengamos aprecio a Riemann; simplemente, es más entretenido matar hormigas con dinamita...
estocástica.
Sin embargo, para efectos de calcular edificios, máquinas y cosas por el estilo (en general, que no sean
muy indignas, en el sentido de la continuidad, o que no sean muy aleatorias) nos basta con la integral
de Riemann.
Finalmente, advertimos que para no insultar la inteligencia del lector, en algunas ocasiones omitiremos
el cálculo de integrales simples. Además, si el lector está aquı́, deberı́a ser capaz de hacer las cuentas
sin mayores problemas. ¡Manos a la obra!

6.1
Albert Einstein (1879-1955) fue un fı́sico de origen alemán, nacionalizado posteriormente suizo y estadounidense. Es
considerado como el cientı́fico más importante del siglo XX, además de ser el más conocido. Realizó importantes avances
dentro de la fı́sica moderna, entre ellos están la Teorı́a de la Relatividad (que engloba a la Teorı́a general de la relatividad
y a la Teorı́a especial de la relatividad ), el Movimiento browniano y el Efecto fotoeléctrico. Aunque es considerado el
“padre de la bomba atómica”, abogó en sus escritos por el pacifismo, el socialismo y el sionismo. Fue proclamado como el
“personaje del siglo XX” y como el más preeminente cientı́fico por la célebre revista Time.
120 Integrales Dobles

6.1. Integrales Dobles

Empezamos por buscar nuestros dx y dy. Dados dos intervalos I1 = [a, b], I2 = [c, d] ⊆ R, consideremos
el rectángulo de R2 dado por

Q := I1 × I2 = [a, b] × [c, d] = (x, y) ∈ R2 : a ≤ x ≤ b, c ≤ y ≤ d .



Como I1 y I2 son intervalos de R, podemos realizar una partición de cada uno, digamos

P1 = {a = x0 < x1 < . . . < xn = b} , P2 = {c = y0 < y1 < . . . < yn = d} .

Ahora, trazando las rectas x = xi e y = yj , dividimos el rectángulo Q en un montón de subrectángulos

de la forma
Qi,j = [xi−1 , xi ] × [yj−1 , yj ], i = 1, . . . , n, j = 1, . . . , m.
Es decir, se tiene que
n [
[ m
Q= Qi,j .
i=1 j=1

Además, el área del subrectángulo Qi,j está dada por

A (Qi,j ) = (xi − xi−1 ) (yj − yj−1 ) .

A partir de lo anterior, definimos la partición P de Q como

P := P1 × P2 = (xi , yj ) ∈ R2 : x0 = a < x1 < . . . < xn = b, y0 = c < y1 < . . . < ym = d .



Ahora, sean P, P 0 dos particiones de Q. Diremos que P 0 es un refinamiento de P si P 0 ⊆ P, es decir,

si todos los subrectángulos generados por P 0 están completamente contenidos por un subrectángulo
generado por P.
6.1 Definición (Suma superior e inferior). Sea Q ⊆ R2 un rectángulo y P una partición de éste. Sea
f : Q → R acotada y consideremos

Mi,j (f ) := sup {f (x, y)} , mi,j (f ) := ı́nf {f (x, y)} .

(x,y)∈Qi,j (x,y)∈Qi,j

Definimos la suma superior y la suma inferior de f con respecto a P como

n X
X m
S (f, P) = Mi,j (f ) (xi − xi−1 ) (yj − yi−1 ) ,
i=1 j=1

n X
X m
s (f, P) = mi,j (f ) (xi − xi−1 ) (yj − yi−1 ) ,
i=1 j=1

respectivamente.
6.2 Observación. Es claro que
s (f, P) ≤ S (f, P) .
Ahora, consideremos P, P 0 dos particiones de Q tales que P 0 ⊆ P, entonces si S es un subrectángulo
generado por P , existen subrectángulos Sk generados por P 0 tales que
[
S= Sk ,
k
Integrales Múltiples 121

y además, X
A(S) = A(Sk ).
k
También, por propiedades de ı́nfimos y supremos, se tiene que

mS (f ) ≤ mSk (f ) , MS (f ) ≥ MSk (f ), ∀k.

Luego,
s (f, P) ≤ s f, P 0 ≤ S f, P 0 ≤ S (f, P) .



(6.1)

Con esto, estamos en condiciones de probar el siguiente Lema.

6.3 Lema. Sean P, Q dos particiones cualquiera del rectángulo Q ⊆ R2 . Sea f : Q → R acotada,
entonces
s (f, P) ≤ S (f, Q) .

Demostración. Consideremos la partición T = P ∪ Q, entonces T es un refinamiento tanto de P como

de Q y por (6.1), sigue que

s (f, P) ≤ s (f, T ) ≤ S (f, T ) ≤ S (f, Q) .

Seguimos nuestro camino hacia la integral. El siguiente paso es hacer un buen refinamiento y para
formalizar qué queremos decir con esto, definimos la norma de una partición P, como

kPk := máx {∆xi = xi − xi−1 , ∆yj = yj − yj−1 } ,

i,j

es decir, la máxima altura o ancho de un subrectángulo. Luego, si kPk → 0, estamos diciendo que cada
una de las alturas y anchos de los subrectángulos están tendiendo a cero y por lo tanto, cada vez las
particiones son más finas.
6.4 Definición (Integrabilidad). Sea f : Q ⊆ R2 → R acotada, con Q rectángulo. Diremos que f es
integrable en Q si existe I ∈ R tal que

lı́m S (f, P) = lı́m s (f, P) = I.

kPk→0 kPk→0

En tal caso, escribiremos f ∈ R (Q) y

¨ ¨
I= f (x, y)dA = f (x, y)dxdy.
Q Q

6.5 Proposición. Sea f : Q ⊆ R2 → R acotada, con Q rectángulo. Entonces, f es integrable en Q si y

sólo si para cada  > 0 existe una partición P de Q tal que

S (f, P) − s (f, P) < .

6.6 Ejemplo. Sea f (x, y) = c constante en Q = [a, b] × [c, d]. Entonces para una partición P de Q, se
tiene que
Mi,j (f ) = mi,j (f ),
y por lo tanto,
S (f, P) − s (f, P) = 0 < , ∀ > 0.
122 Integrales Dobles

Luego, f es integrable y se tiene que

¨ Xn X m n X
X m
f dA = S (f, P) = Mi,j (f )∆xi ∆yj = c ∆xi ∆yj = c(b − a)(d − c).
Q i=1 j=1 i=1 j=1

2
6.7 Observación.
  Sea f : Q ⊆ R → R acotada, con Q rectángulo. Sea P una partición de Q, entonces
para todo x∗i , yj∗ ∈ Qi,j se tiene que

mi,j (f ) ≤ f x∗i , yj∗ ≤ Mi,j (f ).

Ası́,
n X
X m
f x∗i , yj∗ ∆xi ∆yj ≤ S (f, P) ,


s (f, P) ≤
i=1 j=1

Por lo tanto, para la suma de Riemann, se tiene que

n X
X m ¨
x∗i , yj∗


lı́m f ∆xi ∆yj = f (x, y)dxdy.
kPk→0 Q
i=1 j=1

6.8 Teorema. Sea f : Q ⊆ R2 → R, con Q rectángulo. Si f es continua, entonces es integrable.

6.9 Ejercicio. Sea f (x, y) = g(x)h(y) con g : [a, b] → R, h : [c, d] → R continuas. Muestre que f es
integrable y
¨ ˆ b  ˆ d 
f (x, y)dxdy = g(x)dx h(y)dy ,
Q a c

con Q = [a, b] × [c, d].

6.10 Proposición. Sean f, g : Q ⊆ R2 → R integrables y c ∈ R. Entonces

(a) f + g es integrable y
¨ ¨ ¨
(f + g)(x, y)dA = f (x, y)dA + g(x, y)dA.
Q Q Q

(b) cf es integrable y ¨ ¨
cf (x, y)dA = c f (x, y)dA.
Q Q

(c) Si f (x, y) ≥ g(x, y), ∀x, y ∈ R, entonces

¨ ¨
f (x, y)dA ≥ g(x, y)dA.
Q Q

(d) |f | es integrable y ¨ ¨
f (x, y)dA ≤ |f (x, y)| dA.
Q Q

(e) Si M = sup f (x, y), m = ı́nf f (x, y), entonces

(x,y)∈Q (x,y)∈Q
¨
m A (Q) ≤ f (x, y)dA ≤ M A (Q) ,
Q

con A (Q) = (b − a)(d − c).

Integrales Múltiples 123

Demostración. Análogas a las propiedades mostradas en Cálculo I. Se dejan como ejercicio al lector.

6.11 Ejemplo. Calculemos ¨

I= 16 − x2 − y 2 dxdy,
Q

con Q = [0, 1] × [0, 2]. Usando la Proposición 6.10, tenemos que

¨ ¨ ¨ ¨ ¨
2 2 2
I= 16dxdy − x dxdy − y dxdy = 16 A (Q) − x dxdy − y 2 dxdy.
Q Q Q Q Q

Luego, si el lector hizo el Ejercicio 6.9, podemos usarlo y sigue que

ˆ 1  ˆ 2  ˆ 1  ˆ 2 
2 2 56
I = 32 − x dx dy − dx y dy = .
0 0 0 0 3

Con esto hemos construido gran parte de lo que es la integral de Riemann en R2 . Pero ¿qué significa
todo esto? El lector recordará que la integral en R representa el área bajo la curva si es que el integrando
es positivo, en caso contrario es la suma de áreas signadas.
Pues bien, si le sumamos una dimensión al asunto, ya no integramos sobre una lı́nea (el eje real) sino que
sobre un rectángulo y por lo tanto, si cada dA está siendo ponderado por una altura f (x, y)... ¡adivine
qué! Si f (x, y) ≥ 0 sobre Q entonces ¨
f (x, y)dA
Q

representa el volumen del sólido encerrado, es decir, el volumen de

Ω = (x, y, z) ∈ R3 : a ≤ x ≤ b, c ≤ y ≤ d, 0 ≤ z ≤ f (x, y) ,


que denotaremos por V (Ω). Análogamente, si f cambia de signo, entonces la integral será la suma
alternada de volumenes signados. Es importante notar que en este último caso,
¨
|f (x, y)|dA
Q

equivale al volumen total del sólido encerrado.

Con esa inquietud filosófica resuelta, empezamos a desarrollar la maquinaria que nos permitirá calcular
integrales más generales.

6.12 Teorema (Fubini). Sea f : Q ⊆ R2 → R integrable, con Q rectángulo. Entonces

¨ ˆ b ˆ d  ˆ d ˆ b 
f (x, y)dxdy = f (x, y)dy dx = f (x, y)dx dy.
Q a c c a

6.13 Ejemplo. Calculemos ¨

I= yexy dxdy,
Q

con Q = [1, 2] × [2, 4]. Tenemos que

ˆ 4ˆ 2 ˆ 4 ˆ 4 4
e2y e8 3e4
I= yexy dxdy = exy |21 dy = e2y − ey dy = − ey = − + e2 .
2 1 2 2 2 2 2 2
124 Integrales Dobles

6.14 Observación. Muchas veces ocurrirá que es más cómodo evaluar la integral en un orden que en
otro. De hecho, alguna vez ni siquiera se podrá evaluar en uno de los órdenes. Si hubiesemos calculado
la integral del Ejemplo 6.13, en el orden dydx, tendriamos que
ˆ 2ˆ 4
I= yexy dydx.
1 2
Pero ˆ ˆ
4 4
y xy 4 exy 4e4x 2e2x e4x e2x
yexy dy = e − dy = − − 2 + 2,
2 x 2 2 x x x x x
es decir, ˆ 2
4e4x 2e2x e4x e2x
−I= − 2 + 2 dx.
1 x x x x
Y hasta ahı́ llegamos pues el integrando no tiene primitiva.

Pero integrar sobre rectángulos es muy fome. Queremos algo más general. Para ello, primero es necesario
definir un par de cosas.
6.15 Definición (Despreciable). Sea A ⊆ R2 . Diremos que A es despreciable o de medida nula si
existe una sucesión de rectángulos (Ui )i∈N tales que
∞
[
(a) A ⊆ Ui .
i=1
∞
X
(b) A (Ui ) < .
i=1
6.16 Observación. ∅ es trivialmente despreciable.
6.17 Ejemplo. El conjunto A = {a1 , . . . , am } con ak ∈ R2 es despreciable pues basta con tomar
( 
 

ai − 4m , ai + 4m , si i ≤ m
Ui = ,
∅ , si i > m
ya que claramente contienen a A y además,
∞ m
X X  
A (Ui ) = = < .
2m 2
i=1 i=1
S∞
6.18 Teorema. Si A = j=1 Aj , con Aj despreciable, entonces A es despreciable.

Demostración. Sea  > 0, entonces para cada Aj , existe (Uj,i )i∈N tales que
∞
X 
A (Uj,i ) < .
2j
i=1

Luego,
∞ [
[ ∞
A⊆ Uj,i ,
j=1 i=1
y
∞ X
∞ ∞
X X 
A (Uj,i ) < = .
2j
j=1 i=1 j=1
Integrales Múltiples 125

6.19 Corolario. Sea A = {a1 , . . . , ak , . . .} ⊆ R2 . Entonces A es despreciable.

6.20 Ejercicio. Traduzca el concepto de medida nula a R y sabiendo que los números racionales pueden
escribirse como una sucesión, pruebe que Q es despreciable.

6.21 Proposición. Sea ϕ : [a, b] → R continua. Entonces Γ(ϕ) es despreciable.

6.22 Teorema (Du Bois-Reymond). Sea f : Q ⊆ R2 → R acotada. Entonces, f es integrable en Q si y

sólo si el conjunto
Nf := {x ∈ Q : f es discontinua en x} ,

es despreciable.

6.23 Definición (Función caracterı́stica). Sea A ⊆ R2 . Definimos su función caracterı́stica como 1A :

R2 → {0, 1} dada por
(
1 , si x ∈ A
1A (x) := .
0 , si x 6∈ A

6.24 Ejemplo. Sea A = Q2 ∩ [0, 1] × [0, 1]. Entonces, 1A no es integrable pues N1A = [0, 1] × [0, 1], el
cual no es de medida nula.

6.25 Observación. Sea f : Q ⊆ R2 → R integrable. Consideremos A ⊆ Q, entonces podemos definir

la función f 1A : Q → R dada por
(
f (x) , si x ∈ A
f 1A (x) := .
0 , si x 6∈ A

Ahora sı́, finalmente, somos capaces de dar una definición general para la integral doble.

6.26 Definición (Integral Doble). Sea f : Q ⊆ R2 → R integrable y A ⊆ Q tal que 1A es integrable.

Se define la integral doble de f sobre A en el sentido de Riemann como
¨ ¨
f (x, y)dA := f 1A (x, y)dA.
A Q

Esta definición nos presenta un nuevo dilema: ¿cuándo 1A es integrable? Las siguientes proposiciones
nos ayudaran a determinarlo con mayor facilidad.

6.27 Proposición. Sea A ⊆ R2 acotado. Entonces, 1A es integrable si y sólo si ∂A es despreciable.

Demostración. N1A = ∂A. Luego, por el Teorema 6.38, se concluye lo pedido.

6.28 Proposición. Sea f : Q ⊆ R2 → R integrable. Sea D ⊆ Q y tal que ∂D es una unión de curvas
suaves finitas. Entonces f es integrable en D.

Demostración. Por la Proposición 6.21, se tiene que cada curva suave finitas es despreciable. Luego,
usando el Teorema 6.18 se tiene que ∂D es despreciable y por la Proposición 6.27, sigue que 1D es
integrable.

6.29 Definición. A continuación definiremos las regiones tipo I y II.

126 Integrales Dobles

(a) Sean ϕ1 , ϕ2 : [a, b] → R continuas y tales que

ϕ2 (t) ≤ ϕ1 (t), ∀t ∈ [a, b].

Llamaremos región tipo I a

D = (x, y) ∈ R2 : a ≤ x ≤ b, ϕ2 (x) ≤ y ≤ ϕ1 (x) .



(b) Sean ψ1 , ψ2 : [c, d] → R continuas y tales que

ψ2 (t) ≤ ψ1 (t), ∀t ∈ [c, d].

Llamaremos región tipo II a

D = (x, y) ∈ R2 : c ≤ y ≤ d, ψ2 (y) ≤ x ≤ ψ1 (y) .



6.30 Observación. Todas las regiones a estudiar se podrán descomponer en un número finito de cortes
en regiones del tipo I y II.

6.31 Observación. Sea D una región del tipo I. Entonces,

ˆ b ˆ bˆ ϕ1 (x)
A (D) = ϕ1 (x) − ϕ2 (x)dx = dydx.
a a ϕ2 (x)

6.32 Ejercicio. Escriba el análogo para una región del tipo II.

6.33 Ejemplo. Calculemos ¨ p

1 − y 2 dA,
D
n √ o
donde D = (x, y) ∈ R2 : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1 − x2 . Notemos que D corresponde al cuarto de cir-
cunferencia unitaria en el primer cuadrante. Luego, por Fubini,
¨ p ˆ √
1ˆ 1−y 2 ˆ 1 1
2
p
2 2 y3 2
1 − y dA = 1 − y dxdy = 1 − y dy = y − = .
D 0 0 0 3 0 3

Nuevamente, notamos que era mucho más cómodo integrar en el orden dxdy que en dydx.

6.34 Observación. Al igual que antes ¨

|f (x)|dA
D

entrega el volumen del sólido encerrado.

6.35 Ejercicio. Calculemos el volumen del sólido S bajo z = xy sobre la región del primer cuadrante
limitada por y 2 = x, y 2 = 8 − x, y = 0. Tenemos que la intersección de las dos primeras curvas ocurre
en (4, 2), luego
ˆ 2 ˆ 8−y2 ˆ
1 2
2 
V (S) = xydxdy = 8 − y 2 − y 4 ydy = 32.
0 y2 2 0
Integrales Múltiples 127

6.2. Integrales Triples

El procedimiento que permite construir la integral triple de Riemann es análogo al de la integral doble.
La única diferencia es que ahora en vez de un rectángulo, consideraremos un caja

Ω = [a, b] × [c, d] × [e, f ],

entonces ahora las sumas superiores e inferiores están dadas por

n X
X p
m X
S (f, P) = Mi,j,k (f ) (xi − xi−1 ) (yj − yi−1 ) (zk − zk−1 ) ,
i=1 j=1 k=1

n X
X p
m X
s (f, P) = mi,j,k (f ) (xi − xi−1 ) (yj − yi−1 ) (zk − zk−1 ) ,
i=1 j=1 k=1

respectivamente.
6.36 Definición (Integrabilidad). Sea f : Ω ⊆ R3 → R acotada, con Ω caja. Diremos que f es integrable
en Ω si existe I ∈ R tal que
lı́m S (f, P) = lı́m s (f, P) = I.
kPk→0 kPk→0

En tal caso, escribiremos f ∈ R (Ω) y

˚ ˚
I= f (x, y, z)dV = f (x, y, z)dxdydz.
Ω Ω

Y el resto, sigue igual. Rescatamos lo siguiente:

6.37 Observación. Sea f (x, y, z) = f1 (x)f2 (y)f3 (z) con f1 : [a, b] → R, f2 : [c, d] → R, f3 : [e, f ] → R
integrables. Entonces
¨ ˆ b  ˆ d  ˆ f 
f (x, y)dxdy = f1 (x)dx f2 (y)dy f3 (z)dz ,
Ω a c e

con Ω = [a, b] × [c, d] × [e, f ].

6.38 Teorema (Du Bois-Reymond). Sea f : Ω ⊆ R3 → R acotada. Entonces, f es integrable en Ω si y
sólo si el conjunto
Nf := {x ∈ Ω : f es discontinua en x} ,
es despreciable.
6.39 Definición (Integral Triple). Sea f : Ω ⊆ R3 → R integrable y A ⊆ Ω tal que 1A es integrable.
Se define la integral triple de f sobre A en el sentido de Riemann como
˚ ˚
f (x, y, z)dV := f 1A (x, y, z)dV.
A Ω

6.40 Ejemplo. Calculemos ˚

xydzdydx,
S
donde S es el sólido limitado por los planos coordenados y el plano x+2y+3z = 12. Primero, despejamos
12 − x − 2y
z= .
3
128 Integrales n-ésimas

Luego, si z = 0, se tiene que

12 − x
y=
2
y si y = 0, entonces x = 12. Ası́, la integral queda dada por
˚ ˆ 12 ˆ 12−x ˆ 12−x−2y
2 3
xydzdydx = xydzdydx
S 0 0 0
ˆ 12 ˆ 6− x  
2 x 2y
= xy 4 − − dydx
0 0 3 3
ˆ 12 
2
3
x 2 x2 6 − x2 2x 6 − x2 864
= 2x 6 − − − dx = .
0 2 6 9 5

Toda la complejidad de las integrales triples está en cambiar el orden de integración lo cual puede llegar
a ser bastante complicado.

6.41 Ejemplo. Dada la integral

ˆ 2 ˆ 4 ˆ 2− y2
f (x, y, z)dzdydx,
−2 x2 0

cambiemos el orden de integración a dxdydz. Notemos que el dominio de integración es

n yo
Ω = (x, y, z) ∈ R3 : −2 ≤ x ≤ 2, x2 ≤ y ≤ 4, 0 ≤ z ≤ 2 − .
2

Notamos que sobre el plano z = 0, la curva la proyección está dada por la región encerrada entre y = x2
e y = 4. Por otro lado, podemos despejar del plano

y = 4 − 2z.

Luego,
˚ ˆ 2 ˆ 4−2z ˆ √
y
f (x, y, z)dV = √
f (x, y, z)dxdydz.
Ω 0 0 − y

6.3. Integrales n-ésimas

El lector debiera ser capaz de reescribir todo lo hecho para R2 (y lo poco en R3 ), en Rn , teniendo en
consideración que ahora serán rectángulos n-dimensionales de la forma

Ψ = [a1 , b1 ] × [a2 , b2 ] × . . . × [an , bn ].

Lo proponemos como un excelente ejercicio que permitirá repasar toda la materia y entenderla a caba-
lidad. Lo más importante de la integral sobre Rn , es que siguen siendo válidos todos los teoremas antes
mencionados. De ahora en adelante, todas las integrales, a menos que queramos enfatizar su dimensión,
las escribiremos simplemente como ˆ
f dn x,
Ψ

donde dn x denota al diferencial de volumen de n dimensiones.

Integrales Múltiples 129

6.4. Teorema del cambio de variable

Recordemos que el teorema del cambio de variable para R dice
ˆ g(b) ˆ b
f (u)du = (f ◦ g)(x)g 0 (t)dt,
g(a) a

con f, g 0 continuas y g 0 (t) 6= 0 para todo t ∈]a, b[.

Notemos que si g 0 (t) 6= 0, al ser continua, se tiene que g 0 (t) > 0 o g 0 (t) < 0, para todo t ∈]a, b[.
Por lo tanto, g(t) es estrictamente creciente o estrictamente decreciente. Luego,

g (]a, b[) = (g(a), g(b)) ,

si g 0 (t) > 0, ó bien,

g (]a, b[) = (g(b), g(a)) ,
si g 0 (t) < 0.

Sea I = g (]a, b[), entonces podemos reescribir el teorema como

ˆ ˆ
f (u)du = (f ◦ g)(x)g 0 (t)dt,
I g −1 (I)

ó bien, ˆ ˆ
f (u)du = − (f ◦ g)(x)g 0 (t)dt.
I g −1 (I)

Es decir, en cualquiera de los dos casos,

ˆ ˆ
f (u)du = (f ◦ g)(x) g 0 (t) dt.
I g −1 (I)

Formalizaremos esta idea y la llevaremos a más variables.

6.42 Definición (Difeomorfismo). Sea ϕ : U1 → U2 , con U1 , U2 ⊆ Rn abiertos. Diremos que ϕ es un

difeomorfismo si es una biyección y tanto ϕ como ϕ−1 son funciones diferenciables.

Si además, ϕ, ϕ−1 ∈ C1 , entonces diremos que es ϕ es un C1 -difeomorfismo.

6.43 Observación. Si ϕ ∈ C1 y Jϕ no se anula en A, entonces por el Teorema de la Función Inversa,

se tiene que ϕ es un C1 -difeomorfismo.

6.44 Teorema (Cambio de variable). Sea ϕ : U1 → U2 , con U1 , U2 ⊆ Rn abiertos, un C1 -difeomorfismo.

Sea f : U2 → R continua y acotada, entonces
ˆ ˆ
f dn x = (f ◦ ϕ)(x) |Jϕ | dn x,
U2 U1

donde |Jϕ |, en este caso, es el módulo del jacobiano.

6.45 Observación. Alguna de las condiciones del Teorema 6.44 se puede violar sobre conjuntos des-
preciable y no nos importará.
130 Teorema del cambio de variable

6.46 Observación. Una buena nemotecnia es que si ϕ(u, v, w) = (x, y, z), entonces

∂(x, y, z)
dxdydz = |Jϕ |dudvdw = dudvdw.
∂(u, v, w)

6.47 Observación. Es bueno recordar que para

 
a b
T =
c d

invertible, entonces  
−1 1 d −b
T = .
|A| −c a

Con esto terminamos de entregar todas las técnicas de integración en varias variables que se ven en este
curso. Sugerimos al lector repasar todas los cambios de coordendas vistos en la Sección 4.6 pues seran
usados constantemente.

6.48 Ejemplo. Calculemos ¨

xydA,
D
donde D es el cuadrilátero acotado por y = 2x, y = x, y = x + 1, y = 2x − 2. Sean u, v dados por
    
u 2 −1 x
= .
v −1 1 y

Entonces,
∂(u, v) 2 −1
= = 1.
∂(x, y) −1 1
Por lo tanto,
∂(x, y)
= 1.
∂(u, v)
Para obtenemos los lı́mites de integración, notemos que

2x = y ⇐⇒ u = 2x − y = 0,

e
y = x ⇐⇒ v = y − x = 0.
Haciendo lo mismo para las otras dos rectas, se obtiene que
¨ ˆ 2ˆ 1 ˆ 2
2 3
xydA = (u + v)(2v + u)dvdu = + u + u2 du = 7.
D 0 0 0 3 2

6.49 Ejemplo. Calculemos ¨

log(x2 + y 2 )dxdy,
D

donde D es el anillo entre los cı́rculos x2 + y 2 = a2 y x2 + y 2 = b2 , con 0 ≤ a ≤ b. Hacemos el cambio a

coordenadas polares, y sigue que
¨ ˆ 2π ˆ b ˆ 2π  ˆ b  ˆ b
2 2 2
log(x + y )dydx = log r rdrdθ = dθ 2r log rdr = 2π 2r log rdr.
D 0 a 0 a a
Integrales Múltiples 131

Pero
ˆ b ˆ b
b2 a2
   
b 1 1
2r log rdr = r2 log r a
− rdr = b2 log b − a2 log a − + = b2 log b − − a2 log a − .
a a 2 2 2 2

Ası́, ¨    
2 2 1 2 1 2
log(x + y )dydx = 2π b log b − − a log a − .
D 2 2
2
6.50 Ejemplo (Integral de Gauss). Desde Cálculo I, siempre que vemos la integral de e−x , sabemos
que debemos tratar de hacerle el quite pues no esa función no posee primitiva. Pues bien, ahora podemos
calcular el valor de ˆ ∞
2
I= e−x dx.
0
Notemos que por el Teorema de Fubini,
ˆ ∞  ˆ ∞  ˆ ∞ˆ ∞
−x2 −y 2 2 −y 2
2
I = e dx e dy = e−x dxdy.
0 0 0 0

Luego, cambiando a polares,

ˆ π ˆ ∞
2 2 π −r2 ∞ π
I = 2
e−r rdrdθ = e = .
0 0 4 0 4
Es decir, ˆ √
∞
−x2 π
e dx = . (6.2)
0 2
6.51 Ejemplo. Calculemos el volumen del sólido Ξ dentro de x2 + y 2 = 1 sobre z = 1 − x2 − y 2 y bajo
z = 4 + x2 + y 2 . Usando coordenadas cilı́ndricas, tenemos que
ˆ 2π ˆ 1 ˆ 4+r2 ˆ 1
r 3 + 2r2 dr = 4π.


V (Ξ) = rdzdrdθ = 2π
0 0 1−r2 0

6.52 Ejemplo. Calculemos ˚

x2 + y 2 dV,
Σ

donde Σ es la región limitada por x2 + y 2 + z 2 ≤ 1 con x, y, z ≥ 0. Usamos coordenadas esféricas y sigue

que
˚ ˆ πˆ πˆ 1
2 2
2 2
r2 cos2 θ sin2 ϕ + r2 sin2 θ sin2 ϕ r2 sin ϕdrdϕdθ


x + y dV =
Σ 0 0 0
ˆ 1  ˆ π !
π 2 π
= r4 dr sin3 ϕdϕ = .
2 0 0 15

6.5. Aplicaciones
Considere un sólido con forma de región Ω ⊆ R3 y cuya densidad está dada por la función integrable
ρ(x, y, z). Entonces, tenemos que su masa total está dada por
˚
M (Ω) = ρ(x, y, z)dV. (6.3)
Ω
132 Problemas Resueltos

Luego, su densidad media está dada por

˚
ρ(x, y, z)dV
Ω˚ M (Ω)
ρ(x, y, z) = = . (6.4)
V (Ω)
dV
Ω

6.53 Definición (Centro de gravedad). Llamaremos centro de gravedad del sólido Ω al punto r =
(x, y, z), donde
˚ ˚ ˚
1 1 1
x= xρ(x, y, z)dV , y = yρ(x, y, z)dV , z = zρ(x, y, z)dV.
M (Ω) Ω M (Ω) Ω M (Ω) Ω

Además, si ρ(x, y, z) ≡ c, r será el centroide.

Notar que el centro de gravedad depende de la posición del origen y los ejes coordenados.
6.54 Observación. Las definiciones son análogas para R2 (y Rn ).
6.55 Ejemplo. Calculemos el centroide del semidisco unitario con ρ(x, y) ≡ 1. Tomaremos y ≥ 0. Es
claro que x = 0 por la simetrı́a del disco. Por otro lado,
¨ ˆ πˆ 1 ˆ π  ˆ 1 
1 2 2 4
y= ydA = 1 r sin θrdrdθ = sin θdθ r dr = .
A (D) D 2π 0 0 π 0 0 3π

6.56 Definición (Momento de inercia). El momento de inercia en torno a un eje conocido puede
calcularse como ˆ
I= ρ(x, y, x)d2 (x, y, z)dV
Ω
donde la función d representa la distancia del punto al eje en cuestión.
6.57 Observación. Cuando calcule momentos de inercia y no se diga otra cosa, asuma que los ejes
pasan por el centro de gravedad del sólido.

6.6. Problemas Resueltos

6.1 Problema. Sean 0 < a < b. Considere f : Q → R continua, donde Q = [a, b] × [a, b].

(a) Justifique que f es integrable y muestre que

ˆ bˆ x ˆ bˆ b
I= f (x, y)dydx = f (x, y)dxdy.
a a a y

(b) Deduzca que si f (x, y) = f (y, x) sobre Q, entonces

¨
1
I= f (x, y)dxdy.
2 Q

Solución:

(a) Dado que f es continua sobre Q tenemos que f es integrable allı́. Luego, por Fubini, se tiene que
ˆ bˆ x ˆ bˆ b
I= f (x, y)dydx = f (x, y)dxdy.
a a a y
Integrales Múltiples 133

(b) Como f (x, y) = f (y, x), tenemos que

ˆ bˆ b ˆ bˆ b
I= f (x, y)dxdy = f (y, x)dxdy.
a y a y

Luego, como las variables son mudas, podemos cambiar x con y para obtener
ˆ bˆ b
I= f (x, y)dydx.
a x

Finalmente,
ˆ bˆ x ˆ bˆ b ˆ b ˆ x ˆ b 
2I = f (x, y)dydx + f (x, y)dydx = f (x, y)dy + f (x, y)dy dx,
a a a x a a x

es decir, ¨ ¨
1 1
I= f (x, y)dydx = f (x, y)dxdy.
2 Q 2 Q

6.2 Problema. Considere la región D = (x, y) ∈ R2 : 0 ≤ y ≤ 1, 0 ≤ x ≤ y . Calcule



¨ p
x y 2 − x2 dA.
D

Solución: Tenemos que

¨ p ˆ 1ˆ y ˆ ˆ
p 1 1 2 3 y 1 1
1
x y 2 − x2 dA = x y 2 − x2 dxdy = − (y − x2 ) 2 dy = y 3 dy = .
D 0 0 3 0 0 3 0 12

6.3 Problema. Calcule ˆ 1ˆ 1

dydx
√
p .
0 x 1 + y3
Solución: Usamos Fubini para cambiar el orden de integración y sigue que
ˆ 1ˆ 1 ˆ 1 ˆ y2 ˆ 1 √

dxdy dydx y2 2 2−1
√
p = p = p dy = .
0 x 1 + y3 0 0 1 + y3 0 1 + y3 3

6.4 Problema. Calcule el volumen del sólido limitado por la superficie z = x2 − y 2 , el plano z = 0 y
los planos x = 1, x = 3.
Solución: Dado que el sólido es simétrico respecto a y = 0, el volumen de éste es
ˆ 3 ˆ x ˆ x2 −y2 ˆ 3ˆ x ˆ 3
2 2 4 3 3 1 80
V =2 dzdydx = 2 x − y dydx = x dx = x4 = .
1 0 0 1 0 3 1 3 1 3

6.5 Problema. Calcule ˆ 1ˆ 1

x3 sin(y 3 )dydx.
0 x2

Solución: Usamos Fubini y sigue que

ˆ 1ˆ 1 ˆ 1ˆ √
y ˆ 1
3 3 3 1 3 1
x sin(y )dydx = x sin(y )dxdy = y 2 sin(y 3 )dy = (1 − cos 1) .
0 x2 0 0 4 0 12
134 Problemas Resueltos

6.6 Problema. Usando que

ˆ b
xb − xa
xy dy = , (6.5)
a log x
calcule ˆ 1
xb − xa
dx.
0 log x
Solución: Usando la ayuda y Fubini, tenemos que
ˆ 1 b ˆ 1ˆ b ˆ bˆ 1 ˆ b
x − xa
 
y y dy b+1
dx = x dydx = x dxdy = = log .
0 log x 0 a a 0 a y+1 a+1

6.7 Problema. Calcule ˆ ∞

cos x − 1
dx.
0 x2
Hint: Recuerde que ˆ ∞
sin x π
dx = .
0 x 2

Solución: Notemos que

ˆ ∞ ˆ ∞ˆ 1 ˆ 1ˆ ∞
cos x − 1 sin(yx) sin(yx)
2
dx = − dydx = − dxdy.
0 x 0 0 x 0 0 x
Sea t = xy, entonces dt = ydx y sigue que
ˆ ∞ ˆ 1ˆ ∞ ˆ ∞
cos x − 1 sin t dt sin t π
2
dx = − t dy = − dt = − .
0 x 0 0 y
y 0 t 2

6.8 Problema. Calcule ˆ ∞

arctan(πx) − arctan x
dx.
0 x

Solución: Se tiene que

ˆ ∞ ˆ ∞ˆ π ˆ πˆ ∞ ˆ
arctan(πx) − arctan x dydx dxdy π π dy π
dx = 2
= 2
= = log π.
0 x 0 1 1 + (xy) 1 0 1 + (xy) 2 1 y 2
p
6.9 Problema. Calcule el volumen del sólido encerrado entre z = x2 + y 2 , z = −1 y x2 + y 2 = 2ay.
Solución: Notemos que
x2 + y 2 = 2ay ⇐⇒ x2 + (y − a)2 = a2 ,
es decir, esta superficie corresponde a un cilindro de centro (0, a) y radio a. Escribamos nuestro problema
en coordenadas polares, tenemos que
p
z = x2 + y 2 ⇐⇒ z = r,

y además,
x2 + y 2 = 2ay ⇐⇒ r2 = 2ar sin θ ⇐⇒ r = 2a sin θ.
Integrales Múltiples 135

Luego, el volumen pedido es

ˆ π ˆ 2a sin θ π ˆ π ˆ 2a sin θ ˆ
(2a sin θ)3 (2a sin θ)2 2
V = (r − (−1)) rdrdθ =
r + rdrdθ = + drdθ
0 0 0 0 0 3 2
ˆ ˆ π π
8a3 π 8a3 cos3 θ

2 2
= sin θ − cos θ sin θdθ + a 1 − cos 2θdθ = − cos θ + + a2 π
3 0 0 3 3 0
32a3
= + a2 π.
9

6.10 Problema. Considere la región D = (x, y, z) ∈ R3 : 0 ≤ z ≤ 1 − y, 0 ≤ x ≤ y 2 , y ≥ 0 . Calcule



˚
(x + y + z)dV.
D

en el orden dydxdz.
Solución: Tenemos que
˚ ˆ 1 ˆ (1−z)2 ˆ 1−z
1
(x + y + z)dV = √
x + y + zdydxdz = .
D 0 0 x 12

6.11 Problema. Calcule ˆ 3 ˆ x+1

(x + y)ey−x dydx.
0 x−2

Solución: Sean u, v dados por     

u 1 1 x
= .
v −1 1 y
Entonces,     
x 1 1 −1 u
= .
y 2 1 1 v
Por lo tanto,
∂(x, y) 1
= .
∂(u, v) 2
Para obtener los lı́mites de integración, notemos que

y = x − 2 ⇐⇒ v = y − x = −2 , y = x + 1 ⇐⇒ v = y − x = 1,

x = 0 ⇐⇒ u = v , x = 3 ⇐⇒ u = v + 6.
Ası́, la integral pedida es igual a
ˆ 1 ˆ v+6 ˆ 1 ˆ 1
(v + 6)2 − v 2
 
v dudv 1 v 1
ue = e dv = ev (3v + 9) dv = 9e.
−2 v 2 4 −2 2 2 −2

6.12 Problema. Sea f, g : R2 → R continuas. Muestre que

ˆ a ˆ a−y ˆ u
−(x+y)
e f (x)g(y)dxdy = e−u f (u − v)g(v)dvdu.
0 0 a
136 Problemas Resueltos

Solución: Sean u, v dados por     

u 1 1 x
= .
v 0 1 y
Entonces,     
x 1 −1 u
= .
y 0 1 v
Por lo tanto,
∂(x, y)
= 1.
∂(u, v)
Para obtener los lı́mites de integración, notemos que

x = a − y ⇐⇒ u = x + y = a , x = 0 ⇐⇒ u = v,

y = a ⇐⇒ v = a , y = 0 ⇐⇒ v = 0.
Ası́,
ˆ a ˆ a−y ˆ aˆ v ˆ aˆ u
−(x+y) −u
e f (x)g(y)dxdy = e f (u − v)g(v)dudv = e−u f (u − v)g(v)dvdu.
0 0 0 a 0 0

6.13 Problema. Calcule ¨

yexy dxdy,
D
donde D es la región determinada por las hipérbolas xy = 1, xy = 3 y las rectas y = 1, y = 3.
Solución: Sean u = y, v = xy, entonces
∂(u, v) 0 1
= = y.
∂(x, y) y x

Por lo tanto,
∂(x, y) 1
= .
∂(u, v) u
Ası́, ¨ ˆ 3ˆ 3 ˆ 3
dudv
yexy dxdy = uev ev dv = 2e e2 − 1 .


=2
D 1 1 u 1

6.14 Problema. Sea ∆ =]0, 1[×]0, 1[×] − π, π[ y sea φ : R3 → R3 la aplicación dada por

φ(u, v, w) = (u, uv cos w, v sin w).

Calcule el volumen de φ(∆).

Solución: Tenemos por un lado el sólido ∆, en el espacio de coordenadas u, v, w. La imagen de este
conjunto bajo la función φ corresponde a un sólido en R3 con coordenadas x, y, z. El volumen de éste
puede calcularse integrando ˚
V (φ(∆)) = dxdydz.
φ(∆)

En analogı́a a los casos anteriores, hacemos un cambio de variables, determinado por la función (x, y, z) =
φ(u, v, w). Ası́, ˚ ˚
∂(x, y, z)
V (φ(∆)) = dxdydz = dudvdw.
φ(∆) ∆ ∂(u, v, w)
Integrales Múltiples 137

Tenemos además que

1 0 0
∂(x, y, z)
abs = v cos w u cos w −uv sin w = uv.
∂(u, v, w)
0 sin w v cos w
Por lo tanto, ˆ 1ˆ 1ˆ π
π
V (φ(∆)) = uvdwdvdu = .
0 0 −π 2

6.15 Problema. Sea

ϕ : [1, 2] × [1, 2] →  R2  .
2 −v 2
(u, v) 7→ uv, u 2

Calcule el área de la imagen de [1, 2] × [1, 2] a través de ϕ.

Solución: Sea  = [1, 2] × [1, 2], entonces
¨ ¨
∂ (x, y)
A (ϕ()) = dxdy = dudv.
ϕ()  ∂ (u, v)
Pero
∂ (x, y) v u
= = −(u2 + v 2 ).
∂ (u, v) u −v
Ası́, ¨ ˆ 2ˆ 2
2 2 14
A (ϕ()) = u + v dudv = u2 + v 2 dudv = .
 1 1 3

6.16 Problema. Considere un disco D de radio a > 0 y un punto A ubicado en su frontera. La densidad
de masa por unidad de área en cualquier punto P del disco es igual a la distancia de P con A. Encontrar

(a) La masa del disco

(b) El centro de masa del disco
(c) El momento de inercia del disco con respecto a un eje perpendicular que pasa por el punto A

Solución: Sin pérdida de generalidad, podemos centrar el disco en (0, a), es decir, éste queda dado por
la ecuación
x2 + (y − a)2 = a2 ⇐⇒ x2 + y 2 = 2ay.
Si lo escribimos en coordenadas polares, tenemos que la ecuación es r = 2a sin θ para θ ∈ [0, π].

(a) La densidad queda dada por ρ(r, θ) = r. Ası́,

ˆ π ˆ 2a sin θ ˆ π ˆ 2a sin θ ˆ
2 8a3 π 3
M (D) = ρ(r, θ)rdrdθ = r drdθ = sin θdθ
0 0 0 0 3 0
ˆ ˆ
8a3 π 3 8a3 π 32a3
= sin θdθ = sin θ − cos2 θ sin θdθ = .
3 0 3 0 9

(b) Por simetrı́a, es claro que x̄ = 0. Por otro lado,

ˆ π ˆ 2a sin θ ˆ π ˆ 2a sin θ
1 9
ȳ = yρ(r, θ)rdrdθ = r3 sin θdrdθ
M (D) 0 0 32a3 0 0
ˆ
9a π 5 6a
= sin θdrdθ = .
8 0 5
138 Problemas Resueltos

(c) Tenemos que

ˆ π ˆ 2a sin θ ˆ π ˆ 2a sin θ
2
I= r ρ(r, θ)rdrdθ = r4 drdθ
0 0 0 0
ˆ π
32a5 512a5 48a2
= sin5 θdθ = = M (D) .
5 0 75 25

6.17 Problema. Calcule el área de la región D, la cual está encerrada por

2
x2 + y 2 = 2a2 x2 − y 2 ,


a > 0.

Solución: Usamos coordenadas polares y sigue que

√
r4 = 2a2 r2 cos2 θ − sin2 θ ⇒ r = a 2 cos 2θ,

con la restricción cos 2θ > 0. Como para θ ∈ [0, 2π[ hay 4 intervalos de largo π4 donde cos 2θ > 0, sigue
que
ˆ π ˆ a√2 cos 2θ ˆ π
4 4
A (D) = 4 rdrdθ = 4 a2 cos 2θdθ = 2a2 .
0 0 0

6.18 Problema. Cambie el orden de integración de

ˆ 1 ˆ 1−2x ˆ 2−4x−2y
2
f (x, y, z)dzdydx
0 0 0

a dxdzdy.
Solución: Tenemos que la región Ω de integración está dada por
 
1
Ω = (x, y, z) ∈ R3 : 0 ≤ x ≤ , 0 ≤ y ≤ 1 − 2x, 0 ≤ z ≤ 2 − 4x − 2y ,
2

es decir, la región del primer octante que queda bajo el plano 4x + 2y + z = 2.

Al dibujar la región (¡hágalo!), se ve que y varı́a en [0, 1], es decir,

0 ≤ y ≤ 1.

Luego, para obtener cómo varı́a z, nos fijamos en la proyección de la región sobre el plano Y Z, es decir,
imponemos x = 0. Notamos que z varı́a entre 0 y la intersección de nuestra “tapa” y el plano Y Z es
decir,
0 ≤ z ≤ 2 − 2y.
Finalmente, x debe variar de tal forma que sigamos estando bajo el plano 4x + 2y + z = 2, es decir,

2 − 2y − z
0≤x≤ .
4
Ası́,
ˆ 1 ˆ 1−2x ˆ 2−4x−2y ˆ 1 ˆ 2−2y ˆ 2−2y−z
2 4
f (x, y, z)dzdydx = f (x, y, z)dxdzdy.
0 0 0 0 0 0
Integrales Múltiples 139

xy
6.19 Problema. Calcule el volumen del sólido bajo z = 1+x2 y 2
, sobre z = 0 y limitada por xy =
1, xy = 5, x = 1, x = 5.
Solución: Sea Θ el sólido en cuestión. Notemos que
˚ ¨
xy
V (Θ) = dzdxdy = dxdy,
Θ D 1 + x2 y 2
donde D es la región limitada por xy = 1, xy = 5, x = 1, x = 5. Sean u = x, v = xy, entonces

∂(u, v) 1 0
= = x.
∂(x, y) y x

Por lo tanto,
∂(x, y) 1
= .
∂(u, v) u
Ası́, ¨ ˆ 5ˆ 5 ˆ 5  ˆ 5 
xy v dudv du v
V (Θ) = dxdy = = dv ,
D 1 + x2 y 2 1 1 1 + v2 u 1 u 1 1 + v2
es decir,
1
V (Θ) = log 5 (3 log 3 − log 2) .
2
√ √
6.20 Problema. Considere la región Ω del espacio limitada por x, y, z ≥ 0, x+ y + 3 z ≤ 1. Determine
V (Ω).
Solución: Sean x = u, v 2 = y, w3 = z, entonces la región puede ser descrita como u, v, w ≥ 0, u+v+w ≤
1. El jacobiano es
1 0 0
∂(x, y)
= 0 2v 0 = 6vw2
∂(u, v)
0 0 3w2
y se anula si y sólo si vw = 0. Por lo tanto,
˚
V (Ω) = 6vw2 dudvdw.
Ω

Ası́ ˆ 1 ˆ 1−u ˆ 1−u−v ˆ 1 ˆ 1−u

2
V (Ω) = 6vw dwdvdu = 2v(1 − u − v)3 dvdu,
0 0 0 0 0
pero
ˆ 1−u 1−u ˆ 1−u 1−u
v(1 − u − v)4 1 (1 − u − v)5
2v(1 − u − v)3 dv = − + (1 − u − v)4 dv = − .
0 2 0 2 0 10 0

Luego, ˆ 1
1 1
V (Ω) = (1 − u)5 du = .
10 0 60

6.21 Problema. Calcule ¨ p

x2 + y 2 dA,
Q

con Q = [0, 1] × [0, 1].

140 Problemas Resueltos

Solución: Notamos que la integral es simétrica respecto a y = x, luego

¨ p ¨ p
x2 + y 2 dA = 2 x2 + y 2 dA,
Q D

donde D es la región limitada por 0 ≤ x ≤ 1 y 0 ≤ y ≤ x. Haciendo el cambio a coordenadas polares,

sigue que x ≤ 1 implica que
r cos θ ≤ 1 ⇒ r ≤ sec θ.
Luego,
¨ p ˆ π ˆ sec θ ˆ π π

2 2
4
2 1 4
3 1 4
x + y dA = r drdθ = sec θdθ = (log (sec θ + tan θ) + sec θ tan θ) .
D 0 0 3 0 3 0
Ası́, ¨ p
2  √  √ 
x2 + y 2 dA = log 2+1 + 2 .
Q 3

6.22 Problema. Sea V = B (0, 1) r C, con C = (x, y, z) ∈ R3 : x2 + y 2 ≤ z 2 . Calcule el volumen de



V.
Solución: Calcularemos el volumen encerrado por Ω = C ∩ B (0, 1) y luego se lo restaremos al vo-
lumen de B (0, 1). Dada la simetrı́a del problema analizaremos sólo la parte positiva de Ω y luego la
multiplicaremos por 2. Buscamos la intersección despejando z y obtenemos que
1
x2 + y 2 = ,
2
es decir, corresponde a una circunferencia a en el plano
1
z=√ .
2
Luego,
ˆ 2π ˆ √1 ˆ √
1−r2 ˆ √1
+

2 2
p  π2
V Ω = rdzdrdθ = 2π 1 − r2 − r rdr = .
0 0 r 0 4
Ası́,

4π π 2
V (V ) = V (B (0, 1)) − 2 V Ω+ = − .
3 2
6.23 Problema. Calcule la masa de B(0, 1), si su densidad está dada por
1
ρ(x, y, z) = p .
1 + x2 + y 2 + z 2
Solución: La masa de la bola está dada por
˚
M (B(0, 1)) = ρ(x, y, z)dV.
B(0,1)

Usando coordenadas esféricas y Fubini, tenemos que

ˆ 2π ˆ π ˆ 1 2 ˆ 2π  ˆ π  ˆ 1 2 
r sin ϕ r
M= drdϕdθ = dθ sin ϕdϕ dr
0 0 0 1+r 0 0 0 1+r
ˆ 1 2 ˆ 1 2 ˆ 1
r r −1 1 1
= 4π dr = 4π + dr = 4π r−1+ dr
0 1+r 0 1+r 1+r 0 1+r
 
1
= 4π log 2 − .
2
Integrales Múltiples 141

6.24 Problema. Sea 2

p la densidad de carga ρ(x, y, z) = z. Calcule la carga total encerrada entre x +
y 2 + z 2 = 2z y z = x2 + y 2 .
Solución: Sea Ω la región encerrada. Entonces la carga total está dada por
˚
Q= zdV.
Ω

Notemos que
x2 + y 2 + z 2 = 2z ⇐⇒ x2 + y 2 + (z − 1)2 = 1,
es decir, la superficie corresponde a una esfera trasladada en k̂. Luego, la semiesfera por sobre el cono
se puede escribir en coordenadas cilı́ndricas como
p
z = 1 + 1 − r2 .

Ası́,
ˆ 2π ˆ 1 ˆ 1+ 1−r2
√ ˆ 1 p 
Q= zrdzdrdθ = π 1 − r2 + 1 − r2 2rdr.
0 0 r 0
En esta última integral hacemos t = 1 − r2
⇒ dt = −2rdr para concluir que
ˆ 1 √ 7π
Q=π t + tdt = .
0 6

6.25 Problema. Considere Ω = (x, y, z) ∈ R3 : x2 + y 2 + z 2 ≤ 4, z ≥ 0 .



(a) Calcule la masa del cuerpo definido por Ω si su densidad está dada por
s
x2 + y 2
ρ(x, y, z) = .
x2 + y 2 + z 2

(b) Calcule el centro de gravedad del sólido.

Solución:

(a) Usando coordenadas esféricas, tenemos que

˚ ˆ 2π ˆ π ˆ 2
r
2 r2 sin2 ϕ 2
M (Ω) = ρ(x, y, z)dV = r sin ϕdrdϕdθ
Ω r2
ˆ π ! 0ˆ 0 0
2 
2 4
= 2π sin2 ϕdϕ r2 dr = π2.
0 0 3

(b) Por simetrı́a, es claro que x = y = 0. Además,

˚ ˆ 2π ˆ π ˆ 2
r
2 r2 sin2 ϕ 2
z M (Ω) = zρ(x, y, z)dV = r cos ϕ r sin ϕdrdϕdθ
Ω 0 0 0 r2
ˆ π
! ˆ
2 
2
2 3 8
= 2π sin ϕ cos ϕdϕ r dr = π.
0 0 3

Luego,    
1 8 2
r= 0, 0, π = 0, 0, .
M (Ω) 3 π
142 Problemas Resueltos

6.26 Problema. Calcule los momentos de inercia Ix , Iy , Iz de un cilindro grueso y homogéneo de

altura h y masa m, como el que se describe en la figura.
Nota: Ix corresponde al momento de inercia respecto al eje x. Asuma que el origen está en el centro de
gravedad del cuerpo.

m
Solución: Llamemos C al cilindro. Como es homogéneo, su densidad es simplemente ρ = V (C) .
Para un punto cualquiera P = (x, y, z) del cilindro, la distancia al eje z es
p
d(P, z) = x2 + y 2 .

Por ende, ˆ p
Iz = ( x2 + y 2 )2 ρdV.
C
Usando coordenadas cilı́ndricas, tenemos que
ˆ h ˆ 2π ˆ r2 r2
2 r4
Iz = ρ r2 rdrdθdz = 2πρh .
−h 0 r1 4 r1
2

Por geometrı́a elemental, el volumen del sólido es V (C) = πh(r22 − r12 ). Entonces,
1
Iz = m(r22 + r12 ).
2
Notemos que la figura tiene simetrı́a respecto al eje z. Por lo tanto, Ix = Iy . Si usamos nuevamente
coordenadas cilı́ndricas, obtenemos
ˆ h ˆ 2π ˆ r2
2 m
Ix = ρ (z 2 + r2 sin2 (θ))rdrdθdz = ... = (3(r22 + r12 ) + h2 ).
−h 0 r1 12
2

Haciendo r1 = 0, obtenemos además los momentos de inercia de un cilindro no hueco.

6.27 Problema. Calcule ˚
xyz
dV,
Σ x4 − y 4
donde Σ es la región determinada por 1 ≤ x2 − y 2 ≤ 4, 1 ≤ xy ≤ 3 y x2 + y 2 ≤ z ≤ 2(x2 + y 2 ).
Solución: Sean
u = x2 − y 2 ,
v = xy,
z
w= .
x2 + y2
Integrales Múltiples 143

Entonces,

2x −2y 0
∂(u, v, w) 1 2x −2y 2(x2 + y 2 )
= y x 0 = = = 2.
∂(x, y, z) − (x22xz − (x22yz 1 x2 + y 2 y x x2 + y 2
+y 2 )2 +y 2 )2 x2 +y 2

Luego, por Fubini,

˚ ˆ 4ˆ 3ˆ 2 ˆ 4  ˆ 3  ˆ 2 
xyz vw 1 1 du
4 4
dV = dwdvdu = vdv wdw = 6 log 2.
Σ x −y 1 1 1 u 2 2 1 u 1 1

p
6.28 Problema. Sea Ω la región de B (0, 2) al interior z = x2 + y 2 . Calcule
˚
x2 + y 2 dV.
Ω

Solución: Utilizamos coordenadas esféricas. La intersección entre la esfera y el cono se escribe

π
q
r cos ϕ = r2 cos2 θ sin2 ϕ + r2 sin2 θ sin2 ϕ = r sin ϕ ⇒ tan ϕ = 1 ⇒ ϕ = .
4
Luego,
˚ ˆ 2π ˆ π ˆ 2 ˆ 2  ˆ π
! √

2 2
4
4 3 4
4
3 16 8 − 5 2 π
x + y dV = r sin ϕdrdϕdθ = 2π r dr sin ϕdϕ = .
Ω 0 0 0 0 0 15

6.29 Problema. Calcule ¨  

y−x
exp dxdy,
A y+x
donde A es el triángulo entre la recta x + y = 2 y los ejes cartesianos.
Solución: Sean u, v dados por     
u −1 1 x
= := T (x, y).
v 1 1 y
Como T es una transformación lineal, para obtener el conjunto T (A), basta con mirar las imágenes de
los vértices del triángulo. Tenemos que

T (0, 0) = (0, 0) , T (2, 0) = (−2, 2) , T (0, 2) = (2, 2).

Por lo tanto, T (A) corresponde a un triángulo de vértices (0, 0), (−2, 2), (2, 2). Luego, como

∂ (u, v)
= |T | = −2,
∂ (x, y)

sigue que
¨ ¨ ˆ ˆ
1 2 v
 
y−x  u  ∂ (x, y)
exp dxdy = exp dudv = exp (uv) dudv
A y+x T (A) v ∂ (u, v) 2 0 −v
ˆ
e − e−1 2
= vdv = e − e−1 .
2 0
144 Problemas Resueltos

6.30 Problema. Calcule ˆ 1ˆ 1

dxdy
.
0 0 (log x + log y) (1 + x2 y 2 )

Solución: Sean u = y, v = xy, entonces

∂(u, v) 0 1
= = y.
∂(x, y) y x

Por lo tanto,
∂(x, y) 1
= .
∂(u, v) u
Luego, como v = xu, entonces
ˆ 1ˆ u ˆ 1ˆ 1 ˆ 1
dvdu dudv − log v π
= = dv = − .
0 0 log v(1 + v 2 )u 0 v log v(1 + v 2 )u 0
2
log v(1 + v ) 4

6.31 Problema. Sean 0 < a < b y D = {(x, y) ∈ R2 : x, y ≥ 0, y ≥ x, a ≤ xy ≤ b, y 2 − x2 ≤ 1}. Calcule

¨
(y 2 − x2 )xy (x2 + y 2 )dA.
D

Solución: Sean u = y 2 − x2 , v = xy, entonces

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


=
∂(x, y) y x

Luego,
¨ ¨ ˆ bˆ 1
(x2 + y 2 )
 
1 1 b+1
(y 2 − x2 )xy (x2 + y 2 )dxdy = uv dudv = uv dudv = log .
D D 2(x2 + y 2 ) 2 a 0 2 a+1

6.32 Problema. Calcule ¨

y2
dA,
R x4
x2
donde R es la región determinada por 1 ≤ y ≤ 2, 0 ≤ xy ≤ 1.
x2
Solución: Sean u = y , v = xy, entonces

2
∂(u, v) 2x
− xy2 3x2
= y = .
∂(x, y) y x y

Por lo tanto,
∂(x, y) 1
= .
∂(u, v) 3u
Luego,
¨ ˆ 1ˆ 2 ˆ 2
y2 dudv 1 du 1
dA = = = .
R x4 0 1 3u3 3 1 u3 8
Integrales Múltiples 145

6.33 Problema. Calcule ˚

2
1 + x2 + y 2 dV,
Ψ
p
donde Ψ es la región encerrada entre z = x2 + y 2 y B (0, 1).
Solución: La intersección de las superficies ocurre cuando

1
x2 + y 2 = .
2
Usamos coordenadas cilı́ndricas, donde el cono se escribe z = r, y sigue que
˚ ˆ 2π ˆ 1
√ ˆ √
1−r2
2 2 2 2
1 + r4 rdzdrdθ



1+ x +y dV =
Ψ 0 0 r
ˆ √1
√
2
p 42
 344 − 180 2
= 2π 1+r 1 − r − r rdr = π.
0 420

6.34 Problema. Calcule ˚

2
e(x+y) dV,
Ω

donde Ω es la región del espacio limitada por los planos z = x + y, z = x − y; lateralmente por los planos
x − y = −1, x − y = 1, x + y = 1, x + y = 2.
Solución: Sean u = x + y, v = x − y, w = z, entonces

1 1 0
∂(u, v, w)
= 1 −1 0 = 2.
∂(x, y, z)
0 0 1

Por lo tanto,
˚ ˆ 2ˆ 1 ˆ u ˆ 2ˆ 1 ˆ 2ˆ 1
2 2 dwdvdu 1 2 1 2
e(x+y) dV = eu = ueu dvdu − veu dvdu
Ω 1 −1 v 2 2 1 −1 2 1 −1
ˆ 2 u2
2
2 e 1 4
ueu du =


= = e −e .
1 2 2
1

6.35 Problema. Sea Ω la región interior al sólido limitado por

x2 y 2 z 2
+ 2 + 2 = 1.
a2 b c
Calcule
˚   2  23
x y2 z2
1− + 2 + 2 dV.
Ω a2 b c

Solución: Sean
x = ar cos θ sin ϕ,
y = br sin θ sin ϕ, (6.6)
z = cr cos ϕ.
146 Problemas Resueltos

Entonces,
a cos θ sin ϕ −ar sin θ sin ϕ ar cos θ cos ϕ
∂(x, y, z)
= b sin θ sin ϕ br cos θ sin ϕ br sin θ cos ϕ
∂(r, θ, ϕ)
c cos ϕ 0 −cr sin ϕ
(6.7)
cos θ sin ϕ −r sin θ sin ϕ r cos θ cos ϕ
= abc sin θ sin ϕ r cos θ sin ϕ r sin θ cos ϕ = abcr2 sin ϕ,
cos ϕ 0 −r sin ϕ
pues el último jacobiano corresponde al de coordenadas esféricas. Ası́,
˚   2  23 ˆ 2π ˆ π ˆ 1
x y2 z2
3
1− 2
+ 2 + 2 dV = 1 − r2 2 r2 sin ϕdrdϕdθ
Ω a b c 0 0 0
ˆ 2π  ˆ π  ˆ 1 

3
2 2 2
= dθ sin ϕdϕ 1−r r dr
0 0 0
ˆ 1

3
= 4π 1 − r2 2 r2 dr.
0

Para esta última integral, hacemos r = sin t, dr = cos tdt y sigue que6.2
ˆ 1 ˆ π 3 1√ 1√
1 Γ 52 Γ 32
 


1 5 3 1 2 2 π2 π π

3
2 2 2
2
4 2
1−r r dr = cos t sin tdt = B , = = = .
0 0 2 2 2 2 Γ(4) 2 3! 32
Ası́,
˚   2  23
x y2 z2 π2
1− + + dV = .
Ω a2 b2 c2 8

6.36 Problema. Consideremos el sólido C limitado por a2 x2 + b2 y 2 = 1, z = 0 y z = 1. Además, sea S

la superficie dada por z 2 = a2 x2 + b2 y 2 . Definimos S1 como el sólido bajo S al interior de C y S2 como
V(S1 )
el sólido sobre S al interior de C. Demuestre que V(S 2)
no depende de a ni de b.
Solución: Notemos que C es un cilindro de base elipsoidal y altura 1. Los puntos de su interior cumplen
a2 x2 + b2 y 2 ≤ 1, luego S está bajo el plano z = 1 y sobre z = 0 en el interior de C.
Hacemos el cambio de variables
r cos θ r sin θ
x= , y= , z = z.
a b
Notar que,
a2 x2 + b2 y 2 = 1 ⇐⇒ r = 1.
De esta forma, se tiene que
  
r cos θ r sin θ
C= , , z ∈ R3 : r ∈ [0, 1], θ ∈ [0, 2π[, z ∈ [0, 1]
a b
  
r cos θ r sin θ 3
S= , , z ∈ R : r ∈ [0, 1], θ ∈ [0, 2π[, z = r
a b
cos θ −r sin θ
a a 0
∂(x, y, z) sin θ r cos θ r
= b b 0 =
∂(r, θ, z) ab
0 0 1
6.2
Es posible resolver la integral sólo los conocimientos de Cálculo I pero es más largo y tedioso.
Integrales Múltiples 147

Entonces, ˆ ˆ
2π 1ˆ r
r 2π
V (S1 ) = dzdrdθ =
0 0 0 ab 3ab
ˆ 2π ˆ 1ˆ 1
r π
V (S2 ) = dzdrdθ =
0 0 r ab 3ab
V(S1 )
En consecuencia, V(S2 ) = 2.

6.37 Problema. Calcule el volumen del sólido Ω limitado por las superficies x2 + y 2 + z 2 = 9 y
z = 7 − (x2 + y 2 ).
Solución: Notamos que las superficies corresponden a: una esfera centrada en el origen de radio 3 y un
paraboloide cuya cúspide es (0, 0, 7). Buscamos las intersecciones

z = 7 − 9 − z 2 ⇐⇒ z 2 − z − 2 = 0 ⇐⇒ z = 2 ∨ z = −1.

Pero si z = 2 ⇒ x2 + y 2 = 5 y si z = −1 ⇒ x2 + y 2 = 8. Por lo tanto, las intersecciones corresponden a

circunferencias.
Para calcular
√ V (Ω) usaremos coordenadas cilı́ndricas para dos intervalos√ radiales √distintos. Cuando
0 ≤ r ≤ 5, el sólido queda atrapado entre las dos semiesferas y cuando 5 ≤ r ≤ 8 el sólido queda
atrapado entre el paraboloide y la semiesfera inferior. En consecuencia,
ˆ 2π ˆ √
5ˆ
√
9−r2 ˆ 2π ˆ √
8 ˆ 7−r2
V (Ω) = √ rdzdrdθ + √ √ rdzdrdθ
0 0 − 9−r2 0 5 − 9−r2
ˆ √5 p ˆ √8  p 
= 4π 9 − r2 rdr + 2π √ 7 − r2 + 9 − r2 rdr.
0 5

Haciendo u = 9 − r2 ⇒ du = −2rdr en cada integral sigue que

ˆ √
5p ˆ 9√ 9
4π 3 76π
4π 9− r2 rdr = 2π udu = u2 = ,
0 4 3 4 3
ˆ √
8 ˆ 4 4
√ u2 2 3
 
2
p 
37π
2π √ 7 − r + 9 − r2 rdr = π u+ u − 2 du = π + u 2 − 2u = .
5 1 2 3 1 6
Finalmente,
76π 37π 63π
V (Ω) = + = .
3 6 2

6.38 Problema. Calcule ¨

xyzdV,
Ω
donde Ω es la región común entre x2 + y 2 + z 2 ≤ 1 y x2 + y 2 + z 2 ≤ 2x.
Solución: Primero pondremos “derecho” el dibujo. Para ello, cambiaremos el eje x con el eje z y ası́,
las regiones son
x2 + y 2 + z 2 ≤ 1 , x2 + y 2 + z 2 ≤ 2z.
Notamos que la última es equivalente a x2 + y 2 + (z − 1)2 ≤ 1. Es decir, Ω corresponde a la región
encerrada entre dos esferas desplazadas. Luego, la intersección entre sus fronteras está dada por
1
x2 + y 2 + z 2 = 1 ∧ x2 + y 2 + z 2 = 2z ⇒ z = .
2
148 Problemas Resueltos

p √
1 3
Pero, si z = 2 ⇒ x2 + y 2 = 2 . Por lo tanto, la integral en coordenadas cilı́ndricas es

¨ ˆ 2π ˆ √
3 ˆ √
1−r2
2
xyzdV = √ (r cos θ) (r sin θ) zrdzdrdθ = 0,
Ω 0 0 1− 1−r2

pues
ˆ 2π
cos θ sin θdθ = 0.
0

6.39 Problema. Calcule el volumen del sólido Υ limitado por

x2 y 2 z 2 x2 y 2 (z − 2)2
+ + =1 , + + = 1.
4 9 4 4 9 4

Solución: Sean u = 2x, v = 3y, w = 2z, entonces

2 0 0
∂(u, v, w)
= 0 3 0 = 12.
∂(x, y, z)
0 0 2

Sea Υ0 el nuevo sólido encerrado por las superficies, que ahora se escriben

u2 + v 2 + w2 = 1 , u2 + v 2 + (w − 1)2 = 1,

respectivamente. Luego, la intersección de éstas se da cuando

1
w2 = (w − 1)2 ⇐⇒ w = .
2
Reemplazando en cualquiera de las dos ecuaciones, concluimos que
3
u2 + v 2 = .
4
Por lo tanto, usando coordenadas cilı́ndricas, se tiene que
¨ ¨ ˆ 2π ˆ √
3 ˆ √
1−r2
2
V (Υ) = dxdydz = 12dudvdw = 12 √ rdwdrdθ
Υ Υ0 0 0 1− 1−r2
ˆ √
3
2
 p 
= 24π r 2 1 − r2 − 1 dr = 5π.
0

6.40 Problema. Se define la función zeta de Riemann ζ como

∞
X 1
ζ(x) = , x > 0. (6.8)
kx
k=1

Muestre que
ˆ 1 ˆ 1ˆ 1
dx1 dx2 . . . dxn
ζ(n) = ... .
0 0 0 1 − x1 x2 . . . xn
Integrales Múltiples 149

Solución: En efecto,
ˆ 1 ˆ 1ˆ 1 ˆ 1 ˆ 1ˆ 1X
∞
dx1 dx2 . . . dxn
... = ... (x1 x2 . . . xn )k dx1 dx2 . . . dxn
0 0 0 1 − x1 x2 . . . xn 0 0 0 k=0
∞ ˆ 1
X ˆ 1ˆ 1
= ... xk1 xk2 . . . xkn dx1 dx2 . . . dxn
k=0 0 0 0
∞
XY n ˆ 1 
= xki dxi
k=0 i=1 0

X∞  ˆ 1 n
k
= x dx
k=0 0
∞
X 1
=
(k + 1)n
k=0
∞
X 1
=
kn
k=1
= ζ(n).