Tema 7

Integrales múltiples

7.1. Definición.

En esta sección estudiamos el cálculo de la integral de una función real de dos variables
denominada integral doble. Se puede utilizar el esquema del tema anterior para la integral
de Riemann.

Definición 7.1.1. Sea f : A = [a, b]×[c, d] ⊆ R2 → R acotada en el rectángulo [a, b]×[c, d]

Si P1 = {x0 , x1 , · · · , xn } es una partición de [a, b] y P2 = {y0 , y1 , · · · , ym } es una

partición de [c, d], se obtiene una partición P = P1 × P2 de [a, b] × [c, d] formada por
los rectángulos de la forma:

Ri = [xi , xi+1 ] × [yj , yj+1 ], i = 1, 2, · · · , n, j = 1, 2, · · · , m.

Sobre cada uno de los rectángulos construimos dos paralelepı́pedos de alturas:

mij = inf{f (x, y) : (x, y) ∈ [xi , xi+1 ] × [yj , yj+1 ],

Mij = sup{f (x, y) : (x, y) ∈ [xi , xi+1 ] × [yj , yj+1 ]

Sumamos los volúmenes de los paralelepı́pedos de alturas mij y Mij :

n X
X m n X
X m
U (f, P ) = Mij (xi+1 −xi )(yj+1 −yj ), L(f, P ) = mij (xi+1 −xi )(yj+1 −yj )
i=1 j=1 i=1 j=1

105
Curso 2019/2020 Matemáticas (Grado en Quı́mica)

Decimos que f es Integrable (Riemann) sobre A y definimos la integral de Rie-

mann de f (x, y) sobre A como:
ZZ
f (x, y)dxdy = inf {U (f, P )} = sup{L(f, P )}.
A

Nota 7.1.2. Si f : A ⊂ R2 → R es una función continua , al igual que ocurrı́a en una

variable, f es [Link] que las funciones que utilizaremos en este tema serán todas
continuas,todas serán también integrables.

7.2. Propiedades

1.- Si f y g son integrables en A, entonces f ± g es integrable en A y

ZZ ZZ ZZ
(f ± g)(x, y)dxdy = f (x, y)dxdy ± g(x, y)dxdy
A A A

2.- Si f es integrable en A y α ∈ R, entonces αf es integrable en A y:

ZZ ZZ
(αf )(x, y)dxdy = α f (x, y)dxdy
A A

3.- Si f y g son integrables en A y f (x, y ≤ g(x, y)∀(x, y) ∈ A,

ZZ ZZ
f (x, y)dxdy ≤ g(x, y)dxdy
A A

ZZ ZZ
4.- Si f es integrable en A, |f | es integrable en A y f (x, y)dxdy ≤ |f (x, y)|dxdy
A A

5.- Sean f, g : A ⊆ R2 → R y A1 , A2 tales que A1 ∪ A2 = A, A1 ∩ A2 = ∅.

f es integrable en A si y sólo si lo es en A1 y A2 . En este caso:

ZZ ZZ ZZ
f (x, y)dxdy = f (x, y)dxdy + f (x, y)dxdy
A A1 A2

7.3. Teorema de Fubini

El Teorema que vamos a enunciar nos proporciona una importante herramienta para el
cálculo de integrales múltiples , ya que permite reducir el cálculo de una integral múltiple
sobre Rn al cálculo de n integrales ordinarias.

106
Grupos A y B Curso 2019/2020

Teorema 7.3.1. (Teorema de Fubini). Sea f : D ⊆ R2 → R integrable en D.

Si D = [a, b] × [c, d], entonces:

!
ZZ Z b Z d
f (x, y)dxdy = f (x, y)dy dx
D a c

Si D = {(x, y) ∈ R2 : a ≤ x ≤ b, g1 (x) ≤ y ≤ g2 (x)} , con g1 , g2 : [a, b] → R

continuas en [a, b] y tales que g1 (x) ≤ g2 (x) ∀x ∈ [a, b], entonces:
!
ZZ Z Z b g2 (x)
f (x, y)dxdy = f (x, y)dy dx
D a g1 (x)

Si D = {(x, y) ∈ R2 : c ≤ y ≤ d, h1 (y) ≤ x ≤ h2 (y)} , con h1 , h2 : [c, d] → R

continuas en [c, d] y tales que h1 (y) ≤ h2 (y) ∀y ∈ [c, d], entonces:
!
ZZ Z Z d h2 (y)
f (x, y)dxdy = f (x, y)dx dy
D c h1 (y)

7.4. Cambio de variable en integrales dobles

Teorema 7.4.1. Sean f : D ⊆ R2 → R integrable en D y h : R2 → R2 de clase C 1 en

D y biyectiva.

x = h1 (u, v)
Si el cambio de variables
y = h2 (u, v)

transforma D0 (región del plano (u, v)) en D (región del plano (x, y)) se tiene:
ZZ ZZ
f (x, y)dxdy = f (h1 (u, v), h2 (u, v))|Jh(u, v)|dudv
D D0
 ∂x ∂x 
 ∂u ∂v 
donde J(u, v) = 
 

∂y ∂y
 
∂u ∂v

Nota 7.4.2. El cambio de variable más común en R2 es el cambio a coordenadas polares:


 x = ρ cos θ
 y = ρ sen θ

107
Curso 2019/2020 Matemáticas (Grado en Quı́mica)

∂x ∂x
∂ρ ∂θ
cos θ −ρ sen θ
El determinante de su Jacobiano es: |Jh(ρ, θ)| = = =ρ
∂y ∂y sen θ ρ cos θ
∂ρ ∂θ
Por tanto: ZZ ZZ
f (x, y)dxdy = ρf (ρ, θ)dρdθ
D D0

7.5. Ejercicios resueltos

ZZ
1. Calcular I = (2x + y)dxdy, siendo J = [1, 3] × [0, 1].
J

SOLUCIÓN: Aplicando el Teorema de Fubini, resulta que

Z 3 Z 1 Z 3 y=1 Z 3
y2
I= dx (2x + y)dy = 2xy + dx = (2x + 1/2)dx =
1 0 1 2 y=0 1
 x x=3
= x2 + = (9 + 3/2) − (1 + 1/2) = 9.
2 x=1
ZZ
2. Calcular I = cos xdxdy, siendo J = [1, 2] × [0, 2].
J

SOLUCIZÓN: Z Aplicando el Teorema deZFubini, resulta

Z 2 que
2 2 Z 2 2
I= dx cos xdy = cosxdx dy = 2 cos xdx =
1 0 1 0 1
x=2
= 2 sen x = 2(sen 2 − sen 1).
x=1
ZZ
y
3. Calcular I = dxdy, siendo
D x + 1
D = {(x, y) ∈ R : 0 ≤ x ≤ 1, 0 ≤ y, x ≤ 2 − y 2 }.
2

SOLUCIÓN: El recinto D es el siguiente:

1.5

1.0 x= 2 - y2

0.5
D

0.5 1.0 1.5 2.0

√
Por lo tanto, si x ∈ [0, 1], entonces 0 ≤ y ≤ 2 − x; y aplicando el Teorema de
Fubini, se tiene que
Z 1 Z √2−x Z 1 Z √2−x
y 1
I= dx dy = dx ydy =
0 0 x+1 0 x+1 0

108
Grupos A y B Curso 2019/2020

1 y=√2−x 1 1
y2
  
2−x
Z Z Z
1 1 1 3
= dx = dx = − 1 dx =
0 x+1 2 y=0 2 0 x+1 2 0 x+1
x=1
1 3 ln 2 − 1
= (3 ln |x + 1| − x) = .
2 x=0 2
ZZ
y
4. Calcular I = dxdy, siendo D = {(x, y) ∈ R2 : x, y ≥ 0, x2 + y 2 ≤ 4}.
D x+2

SOLUCIÓN: El recinto D es el siguiente:

2.0

x2 + y2= 4
1.5

1.0
D
0.5

0.5 1.0 1.5 2.0

√
Por lo tanto, si x ∈ [0, 2], entonces 0 ≤ y ≤ 4 − x2 ; y aplicando el Teorema
de Fubini, se √tiene que
Z 2 Z 4−x2 Z 2 Z √4−x2
y 1
I= dx dy = dx ydy =
0 0 x +√2 0 x+2 0
Z 2  2 y= 4−x2
1 2 4 − x2 1 2
Z Z
1 y
= dx = dx = (2 − x)dx =
0 x+2 2 y=0 2 0 x+2 2 0
x=2
x2
 
1 1
= 2x − = (4 − 2) = 1.
2 2 x=0 2

5. Calcular el área del recinto D = {(x, y) ∈ R2 : x2 ≤ y ≤ x}.

SOLUCIÓN: El recinto D es el siguiente:

1.5

y ‡ x2
y‡x
1.0

0.5
D
0.2 0.4 0.6 0.8 1.0 1.2

ZZ
2
Por lo tanto, si x ∈ [0, 1], entonces x ≤ y ≤ x; y el área será A = dxdy.
D
Aplicando el Teorema de Fubini, se tiene que

109
Curso 2019/2020 Matemáticas (Grado en Quı́mica)

1 x 1 y=x 1 x=1
x2 x3
Z Z Z Z 
1 1
A= dx dy = y dx = (x−x2 )dx = − = − =
0 x2 0 y=x2 0 2 3 x=0 2 3
1
.
6
Z +∞ √
2
6. Demostrar que e−x dx = π.
−∞
ZZ
2
−y 2
SOLUCIÓN: Para ello, vamos a calcular la siguiente integral doble e−x dxdy
R2
por dos métodos diferentes.
Por un lado, por el Teorema de Fubini resulta que
ZZ Z +∞ Z +∞ 
−x2 −y 2 −x2 −y 2
e dxdy = e dy dx =
R2 −∞ −∞

Z +∞  Z +∞  Z +∞ 2
−x2 −y 2 −x2
= e dx e dy = e dx .
−∞ −∞ −∞

Por otro lado, si realizamos el cambio a polares resulta que

ZZ Z
2 2 2
e−x −y dxdy = e−r rdrdα =
R2 (0,+∞)×(−π,π)

π +∞ +∞
−1 −r2
Z Z
2
= dα re−r dr = 2π · e = π.
−π 0 2 0

Por consiguiente, uniendo ambas partes, tenemos que

Z +∞ 2 Z Z
−x2 2
−y 2
e dx = e−x dxdy = π
−∞ R2

Z +∞ √
2
de donde se deduce que e−x dx = π.
−∞

7.6. Ejercicios propuestos

1.- Calcular las siguientes integrales dobles en el rectángulo que se indica:

ZZ ZZ
y
(a) e dxdy, I = [0, 1] × [1, 3] (b) (x + 1) dxdy,I = [0, 3] × [0, 1]
I I
ZZ ZZ
(c) (2x + y) dxdy, I = [1, 3] × [0, 1]. (d) cos x dxdy, I = [1, 2] × [0, 2].
I I

2.- Calcular las siguientes integrales dobles en el recinto que se indica:

110
Grupos A y B Curso 2019/2020

ZZ
(a) x2 y dxdy, D = {(x, y) ∈ R2 : 0 ≤ x ≤ 1 0 ≤ y ≤ x}
D
ZZ
(b) xey dxdy, D = {(x, y) ∈ R2 : 0 ≤ x ≤ 1 0 ≤ y ≤ 2x2 }
D
ZZ
(c) xy dxdy, D = {(x, y) ∈ R2 : x2 + y 2 ≤ 4 ; x2 + y 2 ≥ 1 y ≥ 0 , x ≥ 0}
D
x2 y2
ZZ  
(d) x dxdy, D = (x, y) ∈ R2 : 0 ≤ x ; 0 ≤ y ; + ≤1
D 4 25
ZZ n x o
(e) sen x dxdy, D = (x, y) ∈ R2 : 0 ≤ x ; 0 ≤ y ; x − y ≤ 0 ; y − ≤ 1
D 2
ZZ
1 
(f) dxdy, D = (x, y) ∈ R2 : 1 ≤ x ; 0 ≤ y ; x ≤ 1 − y 2
D 1 + y
ZZ
1 
(g) dxdy, D = (x, y) ∈ R2 : y ≥ 3 ; y ≤ 4 − x2
D x−2
ZZ
y
(h) dxdy, D = {(x, y) ∈ R2 : 0 ≤ x ≤ 1, 0 ≤ y, x ≤ 2 − y 2 }.
D x+1
ZZ
y
(i) dxdy, D = {(x, y) ∈ R2 : x, y ≥ 0, x2 + y 2 ≤ 4}.
D x+2

3.- Aplicando el Teorema de Fubini, calcular las siguientes integrales:

ZZ
(a) |y − sen x|dxdy, A = [0, π] × [0, 1].
A

ZZ x si x > y

(b) f (x, y)dxdy, A = [0, 1] × [0, 1], f (x, y) =
A y 2 si x ≤ y.


4.- Construye el dominio cuyas áreas están expresadas por las siguientes integrales:
Z a Z x Z a Z √a2 −y2 Z e Z ln x
(a) dx dy (b) dy dx (c) dx dy.
0 0 0 a−y 1 0

5.- Calcula las siguientes integrales

ZZ
(a) x dxdy, donde A = {(x, y) : x2 + y 2 ≤ 2x}.
A
ZZ
x
(b) e y dxdy, donde A = {(x, y) : y 3 ≤ x ≤ y 2 , y ≥ 0}.
Z ZA
(c) ex+y dxdy, donde A = {(x, y) : |x| + |y| ≤ 1}.
A

6.- Expresar por integrales dobles y calcular las áreas acotadas por las siguientes curvas:

(a) xy = 4, y = x, x = 4, y = 0.

(b) y = x2 , 4y = x2 , y = 4.

(c) y = log x, x − y = 1, y = −1.

111
 Curso 2019/2020 Matemáticas (Grado en Quı́mica)

7.- Calcular mediante una integral doble el área de los siguientes recintos

(a) A = {(x, y) ∈ R2 : 1 ≤ x ; x − y ≤ 0 ; y + x ≤ 4}.

(b) B = {(x, y) ∈ R2 : y ≤ 0 ; x + y ≤ 3 ; y 2 ≤ 4x}.

(c) C = {(x, y) ∈ R2 : x2 ≤ y ≤ x}.

ZZ
8.- Cambiar a coordenadas polares la integral f (x, y)dxdy en los siguientes casos
D

(a) D es el cı́rculo x2 + y 2 ≤ ax, con a > 0.

(b) D es el recinto del primer cuadrante limitado por las curvas x+y = 1 y x2 +y 2 = 1.

(c) D es el cuadrado unidad [0, 1] × [0, 1].

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

ZZ
9.- Utilizando el cambio a coordenadas polares, calcular f (x, y)dxdy en los siguientes
A
casos:

(a) f (x, y) = x2 + y 2 , A = {(x, y) ∈ R2 : x2 + y 2 ≤ 2y, x2 + y 2 ≤ 1, x ≥ 0}.

x
(b) f (x, y) = 2 , A = {(x, y) ∈ R2 : x2 + y 2 − 2x ≤ 0, x ≥ 1}.
x + y2
p
(c) f (x, y) = a2 − x2 − y 2 , siendo A el conjunto cuya frontera es la curva (x2 +
y 2 )2 = a2 (x2 − y 2 ), a > 0.

(d) f (x, y) = (1 + x” + y 2 )3/2 , siendo A el cı́rculo unidad.

10.- Sea D el recinto del plano D = {(x, y) ∈ R2 : (x − 1)2 + y 2 ≤ 1, 0 ≤ y ≤ 2 − x} .

ZZ
(a) Calcular xy dxdy.
D
(b) Expresar la integral anterior en función de las coordenadas polares.
√
Z 2 Z x 3
11.- Transformar la integral x dydx a coordenadas polares y resolverla.
0 x
Z 4 Z 2
1
12.- Dada la integral doble √
p dydx, se pide:
0 x x + y2
(a) Dibujar el recinto de integración.

(b) Cambiar el orden de integración en dicha integral y finalmente, calcularla.

13.- Calcular mediante una integral doble el área de los siguientes recintos

(a) A = {(x, y) ∈ R2 : 1 ≤ x ; x − y ≤ 0 ; y + x ≤ 4}.

112
 Grupos A y B Curso 2019/2020

(b) B = {(x, y) ∈ R2 : y ≤ 0 ; x + y ≤ 3 ; y 2 ≤ 4x}.

(c) C = {(x, y) ∈ R2 : x2 ≤ y ≤ x}.

ZZ
14.- Cambiar a coordenadas polares la integral f (x, y)dxdy en los siguientes casos
D

(a) D es el cı́rculo x2 + y 2 ≤ ax, con a > 0.

(b) D es el recinto del primer cuadrante limitado por las curvas x+y = 1 y x2 +y 2 = 1.

(c) D es el cuadrado unidad [0, 1] × [0, 1].

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

ZZ
15.- Utilizando el cambio a coordenadas polares, calcular f (x, y)dxdy en los siguientes
A
casos:

(a) f (x, y) = x2 + y 2 , A = {(x, y) ∈ R2 : x2 + y 2 ≤ 2y, x2 + y 2 ≤ 1, x ≥ 0}.

x
(b) f (x, y) = 2 , A = {(x, y) ∈ R2 : x2 + y 2 − 2x ≤ 0, x ≥ 1}.
x + y2
p
(c) f (x, y) = a2 − x2 − y 2 , siendo A el conjunto cuya frontera es la curva (x2 +
y 2 )2 = a2 (x2 − y 2 ), a > 0.

(d) f (x, y) = (1 + x” + y 2 )3/2 , siendo A el cı́rculo unidad.

113
Curso 2019/2020 Matemáticas (Grado en Quı́mica)

114