INTEGRALES MÚLTIPLES

UNIVERSIDAD NACIONAL DE SAN JUAN

Facultad de Ingeniería
Departamento de Matemática

DEPARTAMENTO DE MATEMATICA

CÁLCULO II
ALIMENTOS – QUIMICA – INDUSTRIAL

CURSO: SEGUNDO AÑO

INTEGRALES MULTIPLES

Equipo de cátedra:

Mg. Ing. Patricia Cuadros Dra. Bioing. Lorena Orosco

Ing. Nicolás Sardiña

Año 2021

1
 INTEGRALES MÚLTIPLES

1. Introducción.
Este módulo trata el segundo concepto básico del cálculo, la integral, para campos escalares y
vectoriales, de dos o más variables. Se definen tres tipos de integrales: Múltiples, de Línea y de
Superficie. Estas definiciones se realizan a partir de un concepto geométrico o físico según
corresponda.
También se analizan algunas aplicaciones como: masa, centro de gravedad, momento de inercia, flujo
y circulación.

En Cálculo I se trabaja en el cálculo de áreas con funciones continuas y definidas en [a ,b], intervalo
o subconjunto de  (o sea f :  → ) definiendo el área bajo la curva, gráfica de y = f(x) como
y
y = f(x)
𝑏
𝐴 = lim [∑𝑛𝑖=1 𝑓(𝑥𝑖 )∆𝑥] = ∫𝑎 𝑓(𝑥)𝑑𝑥
||𝑃||→0
A

0 a b x

Se generaliza el concepto de integral para funciones de dos variables basándose en esta idea.
El intervalo unidimensional [a,b] se reemplaza por un conjunto R bidimensional, llamado región de
integración, contenido en 2 .
El integrando es una función o campo escalar f: 2 →  definido y acotado en R   2 .
La integral que resulta se llama integral doble y se representa por:

 f ( x, y)dA ,  f ( x, y) dx dy
R R

Los símbolos diferenciales dx, dy indican las variables sobre las que se integra y se realiza la
transformación de integrales.

2
 INTEGRALES MÚLTIPLES

2. Integral doble

2.1. Definición de Integral doble

Sea z = f(x, y), una función definida y acotada en el rectángulo R : [a, b] x [c, d] , dividiendo R
en n-subrectángulos obtenemos una partición, de R.
La gráfica de z es
una superficie en R3
z
z = f(x, y) z

c y d y
a k
c d xk
a y b
x Ak
b R
x

Primero suponemos f(x,y) ≥ 0. La gráfica de f es una superficie con ecuación z = f(x,y). Sea S el
sólido que se encuentra encima de R y debajo de la gráfica de f.

El área del k-ésimo rectángulo es ΔAk =Δxk Δyk

Multiplicando Ak por el valor inferior de f(x, y) en el k-ésimo subrectángulo, mk , formamos la

n
suma inferior s p =  mk Ak
k =0

z z = f(x, y)

b x Ak
y
c d y
x

Multiplicando Ak por el valor superior de f(x, y) en el k-ésimo subrectángulo, M k , formamos la

n
suma superior S p =  M k Ak
k =0
3
 INTEGRALES MÚLTIPLES

z z = f(x, y)

b x Ak
y
c d y
x

Nota: Se verifica s p = Sp solamente cuando f (x, y) es constante en cada rectángulo

La suma inferior sp nos da un volumen contenido en el cuerpo limitado por la superficie

z = f(x,y), el plano xy y el cilindro cuya base es Ak; cuerpo que se llama cilindroide por analogía
al trapezoide.
La suma superior Sp nos da el volumen continente de dicho cilindroide.

Se trata de aproximar el volumen del cilindroide por defecto y por exceso, mediante la suma de
volúmenes contenidos y continentes correspondientes a cada partición. En cada subrectángulo de la
partición se cumple (ver figuras anteriores):

mk Ak  Vk  M k Ak

Si hacemos esto para cada Ak que pertenece a la partición tenemos aproximadamente el volumen
del cuerpo de base R y altura variable medida sobre z = f(x, y):

n
V   Vk
k =0

Comparando el volumen V con las sumas inferior y superior se observa que s p  V  S p para cada

partición P. Si tenemos una sucesión de particiones P1  P2    Pn   tales que cuando

n →  y la norma de la partición tiende a cero, es decir P → 0 ( P = máxima diagonal de los n-

def
rectángulos), entonces tendremos s p n  V  Sp n  n y como:

lim s p n = lim S Pn =
P →0 P →0 def
 f (x, y) dx dy
R

ya que z = f(x, y) es integrable por ser definida y acotada, tendremos :

4
 INTEGRALES MÚLTIPLES

V =  f ( x, y) dx dy
R

Otra forma de definir integral doble.

Los cinco pasos, desarrollados en la siguiente tabla, conducen a la definición de integral doble. Se
utiliza la suma de Riemann y el concepto de volumen. Se desarrolla haciendo la comparación con
la integral de funciones de una variable, y = f(x)

y = f (x) z = f (x,y) Interpretación gráfica

1. Considérese que f está Sea f definida y acotada en
definida y acotada en un una región cerrada R
intervalo cerrado [a,b]
2. Forme una partición P del Por medio de una cuadricula
intervalo [a,b] en n (o red de rectas verticales y
subintervalos de longitudes horizontales paralelas a los
xk ejes coordenados), forme una
partición P de R en n
subregiones rectangulares Rk
de áreas Ak = xk yk
contenidas totalmente en R.
3. Sea P la norma de la Sea P la norma de la
partición, o sea la longitud del partición, o sea la longitud de
subintervalo mayor. la diagonal mayor de las Rk.de
área Ak .
4. Elíjase un xk* intermedio en Selecciónese un punto
cada subintervalo [xk-1,xk] ( xk *, yk *) arbitrario inter-
medio en cada subregión Rk
de área Ak.
5. Evalúe la suma Evalúe la suma
n n
 f ( x k *) x k  f ( x k *, y k *) Ak
k =1 k =1

De esta manera se tiene la siguiente definición.

Definición: Sea f una función de dos variables definida y acotada en una región cerrada R. Entonces
la integral doble de f en R está dada por


R
f ( x, y ) dA = lim
P →0  f ( x *, y *) A
k =0
k k k

Nota:
5
 INTEGRALES MÚLTIPLES

n
1 -  f ( x k *, y k *) A k se llama suma de Riemann.
k =0

2 - Si f es continua en R entonces f es también integrable.

3 - La integral 
R
f ( x, y ) dx dy existe si f está definida y acotada en todos los puntos de R., también

existe la integral si los puntos de discontinuidades forman un conjunto finito de puntos.

4- 
R
f ( x, y ) dx dy calcula un volumen cuyo techo es z = f (x,y).

5 – Si z = f (x,y) = 1,  dx dy calcula el área plana de la región R.

R

2.2. INTEGRALES ITERADAS

2.2.1. Cálculo de integrales dobles de funciones continuas por medio de integrales reiteradas
sobre un rectángulo.

De manera similar al proceso de derivación parcial podemos definir una integración iterada (o
parcial). Este concepto es la clave para el método práctico de evaluación de la integral múltiple.

Sea z = f (x,y), talque f : 2 →  continua en el rectángulo R : [a,b] x [c,d]

z = f(x, cte.) z

z = f(x, y)

c y = cte d y
a
b
x

Por lo que sabemos de primer año se puede calcular el área de la región sombreada, en la gráfica
anterior

6
 INTEGRALES MÚLTIPLES

z
f(x, cte.) b
A( y) =  f ( x, y) dx
a
A(y)

0 a b x

A cada yi  c, d  la integral anterior le asocia un valor A(yi) los que representamos en el siguiente
sistema cartesiano:
A(y)
(yi, A(yi))

0 c yi d y

d
El área sombreada en la figura anterior se calcula como  A( y) dy .
c
b
Reemplazando A( y) por  f (x, y) dx
a
podemos calcular el volumen definido por z = f (x,y), sobre el

rectángulo R : [a,b] x [c,d]

d d b

V =  A( y) dy =   f ( x, y) dx dy
c c a

De igual manera, pero comenzando el trabajo con x = cte, podremos obtener el mismo volumen:

b b d

V =  A( x) dx =   f ( x, y) dy dx
a a c

Nota: Cuando la región de integración es rectangular los límites de integración son todos
constantes.

2.2.2. Teorema de Fubini.

Si f(x,y) es continua en la región rectangular R, con a ≤ x ≤ b , c ≤ y ≤ d entonces

𝑏 𝑑 𝑑 𝑏

∫ ∫ 𝑓(𝑥, 𝑦) 𝑑𝑦 𝑑𝑥 = ∫ ∫ 𝑓(𝑥, 𝑦)𝑑𝑥 𝑑𝑦

𝑎 𝑐 𝑐 𝑎

7
 INTEGRALES MÚLTIPLES

2.2.3. Integrales dobles extendidas a regiones más generales.

Hasta aquí la integral doble solo se ha definido para regiones de integración rectangulares. No
obstante, no hay dificultad para extender el concepto a regiones más generales. I) Sea T una región
acotada e incluyamos T en un rectángulo R.
z

z = f(x, y)

R -T
x

Sea f una función definida y acotada en T. Definamos una nueva función f en R del siguiente
modo:

f(x,y) si (x,y)  T
f (x,y) = 
0 si (x,y)  R- T

Es decir, extendemos la definición de f a todo el rectángulo R haciendo que la función valga cero
fuera de T.
Como f es integrable en T , por lo tanto:

 f (x, y) dx dy =  f (x, y) dx dy
T R

Para resolver una integral doble sobre un recinto cualquiera, se considera la variación de las
variables, en el recinto, de dos maneras diferentes:

a) Definamos el recinto R, dominio de integración de la función f(x,y), considerando a la variable

x como independiente y a la variable y como dependiente, entonces es

R = { (x, y) / a  x  b  θ 1(x)  y  θ 2 (x) }

Se interpreta en un gráfico el recinto R y su variación para visualizar lo antes definido. Se

tiene en cuenta la definición anterior, luego la integral doble queda planteada de la siguiente manera

8
 INTEGRALES MÚLTIPLES

a
. y

 1(x)
x
b .  2(x)

Se tiene en cuenta la definición anterior, luego la integral doble queda planteada de la siguiente
manera

b 2 ( x)

 f ( x, y) dx dy =   f ( x, y ) dy dx
R a 1 ( x )

Ejemplo 1: Sea R la región definida por: R = {(x,y) : 0 ≤ x ≤ 1, 0 ≤ y ≤ x² }, donde

f(x,y) = 2x. Calcule  f(x, y ) dx dy .
R

Solución:
2 1
 f ( x, y ) dxdy =  2 x dydx = 0  0 2 x dydx =
1 x
R R
2

b) Definamos el recinto R, dominio de integración de la función f(x,y), considerando a la variable

y como independiente y a la variable x como dependiente, entonces es

R = { (x, y) / c  y  d  ψ 1( y)  x  ψ 2 ( y) } .

La interpretación gráfica de R es:

9
 INTEGRALES MÚLTIPLES

c d y
 1(y)
. .
x  2(y)
Con esta definición de R la integral doble queda expresada de la siguiente manera:

d 2 ( y )
 f (x, y) dx dy =   f (x, y) dx dy
R2 c 1 ( y)

Ejemplo 2: Sea R la región definida por:

R = { (x,y) : 0 ≤ y ≤ 1, 0 ≤ x ≤ y }, donde
f(x,y) = 2x. Plantee  f(x, y ) dx dy .
R
Solución:
1 y
 f ( x, y) dxdy =  2x dydx =   2x dxdy
R R 0 0

Generalizando: Dada R una región cerrada y acotada que tiene por frontera a una curva  con la
siguiente propiedad, que toda recta corta a  en a lo sumo dos puntos, como lo indica la figura, y una
función f(x.y) definida y continua en R,

z
z

c d c d
y y
a
a
. .  1(y)
.
 1(x)

x
b .  2(x)
x
b
 2(y)

Se puede plantear:

10
 INTEGRALES MÚLTIPLES

b  2 (x) d  2 ( y)
 f (x, y) dx dy =   f ( x, y) dy dx =   f (x, y) dx dy
R a 1 ( x ) c  1( y)

Nota: Observe que el resultado no depende del orden de integración.

Ejemplo 3: Sea R la región limitada por: y = x, y = x2 donde f(x,y) = 2x. Plantee  f(x, y ) dx dy .
R
Solución:
1 x
 f ( x , y ) dydx =  2 x dydx =   2x dydx
R R 0 x2

2.3. Cálculo de áreas planas mediante integrales dobles.

Para calcular el área de la región plana indicada en la Figura [1], utilizando integrales definidas para
funciones de una variable, debemos trabajar de la siguiente manera:

b b b b f 2 (x)
A =  f 2 ( x ) dx −  f1 ( x ) dx =  (f 2 ( x ) − f1 ( x )) dx =   dy dx
a a a a f1 ( x )

La última expresión se obtiene aplicando el concepto vertido por Barrow para el cálculo de
integrales definidas.

y = f2(x)
Figura [1]

y = f1(x)

a b x

11
 INTEGRALES MÚLTIPLES

Podemos observar en el desarrollo anterior que, cuando f(x,y) = 1, la integral doble representa el área
del recinto de integración

Area =  dx dy
R

Ejemplo 4: Plantear la integral que permite calcular el área de la región limitada por: y = x,
y= x.

 dxdy =  dydx = 0 x dydx

Solución: 1 x
R R

Nota: por integrales dobles podemos calcular áreas planas y volúmenes, en este curso utilizaremos
estas integrales solamente para calcular áreas planas .

V =  f ( x, y) dx dy , Area =  dx dy
R R

Ejemplo 5:
1 x
Determinar la región R e invertir el orden de integración: a)   f ( x, y) dy dx
0 x2

Solución: y = x ; y = x2

1 y
  f ( x, y) dx dy
0 y

b) Ídem para:
12
 INTEGRALES MÚLTIPLES

y
3 25− y 2
  f ( x, y) dx dy
0 4 5
y
3
4
Solución: x = 25 − y 2 x= y
3
R
3 4 5 x
x
4 4 5 25− x 2
  f ( x, y) dy dx +   f ( x, y) dy dx
0 0 4 0

Ejemplo 6:

Calcular el área del dominio limitado por las curvas y = 2 - x ², y = x

Solución:

Primero se determinan los puntos de intersección.

x = 2 - x² => x² + x - 2 = 0

(−1  1 + 8 )
x1,2 = ; x1 = 1, x2 = -2
2

P1 (1,1) y P2 (-2,-2)

1 2− x 2
A=   dy dx
−2 x

1
1 x3 x2 1 1   8   27
A =  (2 − x - x) dx = 2x −
2
− = 2 − − −  − 4 +   − 2 =
−2 3 2
−2
3 2   3  6

22.5. Cambio de coordenadas en integrales.

.
En teoría de integración unidimensional el método de sustitución nos permite calcular integrales
complicadas transformándolas en otras más sencillas o en tipos de integrales que pueden calcularse
más fácilmente, quedando así:

b t2
 f ( x ) dx =  f [g( t )] g ( t ) dt con x = g(t) (1)
a t1

13
 INTEGRALES MÚLTIPLES

donde a = g(t1) y b = g(t2) . Suponiendo que g tiene derivadas continuas en [t1, t2] y que f es continua
en el conjunto de valores que forma g(t) al variar t en [t1, t2].

2.5.1. Cambio de coordenadas en integrales dobles. Coordenadas curvilíneas

En dos dimensiones existe análogo problema para las integrales dobles. Se transforma una integral
doble de la forma

 f ( x , y) dx dy
R

extendida a una región R del plano xy, en otra integral doble  F(u , v) du dv extendida a una nueva
T
región T del plano uv.

Se va a estudiar a continuación la relación entre las áreas definidas por las regiones R y T, también
de los integrandos f(x,y) y F(u,v). El método de sustitución en integrales dobles es más laborioso
que en las simples debido a que existen dos sustituciones formales a efectuar, una respecto a x; y otra
respecto a y. Esto significa que en lugar de una función g que aparece en (1), tenemos ahora dos
funciones X e Y que relacionan x,y con u,v del modo siguiente :

 x = X ( u , v)
Sea  la aplicación  :  (2)
 y = Y ( u , v)

v y
R
(x,y)
T
(u,v)

u x
}

 x = X ( u , v)
Geométricamente, puede considerarse que  definen una "aplicación" que hace
 y = Y ( u , v)
corresponder a un punto (u,v) del plano uv, el punto imagen (x,y) del plano xy. Es decir que un
conjunto T de puntos del plano uv es aplicado a otro conjunto R del plano xy, como se representa en
la figura.
La aplicación también se expresa mediante una función vectorial. En el plano xy se traza el radio
vector que une el origen con un punto genérico (x,y) de R.
r(u, v) = X(u, v) i + Y(u, v) j ecuación vectorial de la aplicación.

14
 INTEGRALES MÚLTIPLES

Como (u,v) son puntos de T, el vértice de r(u, v) = X(u, v) i + Y(u, v) j describe puntos de R.
Algunas veces puede expresarse “u” y “v” en función de “x” e “y” quedando:

u = U(x, y)
 -1: 
v = V(x, y)

Estas ecuaciones definen una aplicación del plano xy en el plano uv, llamada aplicación inversa de la
definida en (2), ya que transforma los puntos de R en los de T.
Esta aplicación es la que se usa en la practica por eso se pide que X(u,v) ; Y(u,v) tengan inversa, para
lo cual  debe ser uno a uno o biunívoca..
Las aplicaciones uno a uno son de especial importancia. Estas transforman puntos distintos de T en
puntos distintos de R; tales aplicaciones establecen una correspondencia biunívoca entre los puntos
de T y los correspondientes de R y permite (por lo menos teóricamente) regresar de R a T por la
aplicación inversa (que naturalmente es uno a uno).
v y

dA
U=cte. T
(x0, y0)
(u0,v0) V=cte.
x
u

Si consideramos un segmento horizontal en el plano uv, sobre dicho segmento v es constante. La

función vectorial r aplica este segmento sobre una curva (llamada curva u) en el plano xy.
     
Si r ( u , v) = X ( u , v) i + Y ( u , v) j luego dr = ru du + rv dv
  
Si v = v 0 = cte entonces d u r = ru du pues rv dv = 0
  
Si u = u 0 = cte entonces d v r = rv dv pues ru du = 0
     
Se observa que d u r = ru du es paralelo a ru y que d v r = rv dv es paralelo a rv

La región rectangular del plano uv se transforma en una porción del plano xy que es casi un
 
paralelogramo cuyos lados son los vectores ru du y rv dv como se observa en la figura anterior. El
área de dicho paralelogramo es el módulo del producto vectorial de ambos vectores.
   
Luego el dA = ru du x rv dv = ru x rv du dv

Donde el módulo del producto vectorial es:

15
 INTEGRALES MÚLTIPLES

  
i j k
  x yu  (x, y)
ru x rv = x u yu 0 = u = = J(u, v)
xv yv  (u, v)
xv yv 0

 (x,y) x yu
Pues por definición el Jacobiano de (u,v) es J(u,v) = = u
 (u,v) x v yv
Recordar y relacionar con matriz jacobiana vista en Unidad I, parágrafo 10.7 pag.41.

Por lo que el área:

xu yu  (x, y)
dA = du dv = du dv
xv yv  (u, v)

Es decir que:  f ( x, y) dx dy =  f [X(u, v), Y(u, v)]  J(u, v) du dv

R T

• Si J(u,v) = 1 para todos los puntos de T, el "paralelogramo" tiene la misma área que el
rectángulo y la aplicación conserva las áreas. Si no es así, para obtener el área del
paralelogramo se debe multiplicar el área del rectángulo por J(u,v). .
Se ve que el jacobiano es un factor de proporcionalidad de las áreas de las regiones.
 
• Si J(u,v) = 0 en un punto (u,v), los vectores ru du y rv dv son paralelos (ya que su producto
vectorial es el vector nulo) y el paralelogramo degenera.
Tales puntos se llaman puntos singulares de la aplicación. Ya sabemos que la forma de
transformación es válida también si existe un número finito de puntos singulares o más general
cuando tales puntos forman un conjunto de medida nula.

De lo estudiado anteriormente tenemos

Teorema : 
Sea z = f(x,y) integrable en una región R. Sea  una transformación biunívoca

 :  2 →  2 dada por x = X (u,v) , y = Y (u,v) tal que J(u,v)  0 ,

para todo (u,v)   −1 ( R ) = T. Entonces se verifica

 f (x, y) dx dy =  f [X(u, v), Y(u, v)] J(u, v) du dv

R T

Observación: la formula anterior vale también si J(u,v) = 0 sobre un número finito de puntos
singulares.

16
 INTEGRALES MÚLTIPLES

2.5.2. Cambio de coordenadas en integrales dobles. Coordenadas polares

x =  cos 
Consideremos:  Y
 y =  sen 
()
x (, ) =  cos  
Es decir que: 
y(, ) =  sen  y= sen 


con  > 0 y 0    0 + 2 x= cos  x

para que la aplicación sea uno a uno.

 x  y
J (, ) =     = cos  sen 
=  (cos 2  + sen 2 ) =   J (, ) =  = 
 x  y −  sen   cos 
   

Quedando dA =   d   d 

La fórmula de transformación es:

 f ( x, y) dx dy =  f ( cos,  sen)  d d
R T

Las curvas  = cte. son rectas por el origen y las curvas  = cte. son círculos concéntricos en el origen
en xy. La imagen del rectángulo en el plano  es un "seudo paralelogramo" en el plano xy limitado
por dos radios y dos arcos de circulo.
y


d d

x


Las coordenadas polares son particularmente importantes cuando la región de integración tiene
fronteras a lo largo de las cuales  o  son constante. En el caso de contornos circulares o elípticos
en el plano xy, se puede transformar en rectas paralelas a los ejes coordenados en el plano ,.

17
 INTEGRALES MÚLTIPLES

Ejemplo 7: Calcular el área de un cuarto de círculo de radio a.

Solución: Se realiza el cálculo utilizando integrales dobles

y coordenadas polares.
y

A =  dx dy
R

x=0
x =  cos 
 J(, ) = 
 y =  sen  y=0 a
x
En el plano xy se toma un cuarto de círculo, en el primer cuadrante y se ve en que se transforma
en el plano (,), transformamos cada una de las fronteras:

Si x = 0 => 0 =  cos  →  = /2

Si y = 0 => 0 =  sen  →  = 0

Si x² + y² = a² → ² = a² →  = ± a,  sólo toma valores positivos   = a

Si x = 0  0 =  cos  
→=0
Si y = 0  0 =  sen 

Luego 0    a El cuarto de circulo se transforma en el rectángulo [0,a] x [0,/2]

 y

 =/2

=a x=0
=0

=0 a y=0 a
 x
Luego el área es:


2 a

A =  d d  =   d d  = 4 a
2

T 0 0

Ejemplo 8:

Calcular el área del círculo de ecuación x 2 + y 2  4x en coordenadas polares.

18
 INTEGRALES MÚLTIPLES

Solución: utilizando MAPLE :

1
 4 cos 

 f (x, y) dx dy =  
R
2
1
− 
2
0
 d d = 4

Ver desarrollo de la solución en Anexo de Integrales Múltiples.

2.5.3. Cambio de coordenadas en integrales dobles. Coordenadas polares generalizadas

x
 a =   cos  0   1
 J(, ) = a  b  
 y =   sen  0    2
 b

! Queda para el lector la demostración del Jacobiano para estas coordenadas y el planteo del dA.

2.6.
2.6.Aplicaciones
Aplicacionesde
delas integrales
integrales dobles.
múltiples
2.6.1. Valor medio de una función.

Sea f(x,y) una función continua en las variables x e y, el valor medio de f(x,y) en una región plana R
está dado por :
∬𝑹 𝒇(𝒙, 𝒚)𝒅𝑨
𝑽𝒂𝒍𝒐𝒓 𝑴𝒆𝒅𝒊𝒐 =
∬𝑹 𝒅𝑨

2.6.1. Masa total de una lámina.

En estas aplicaciones se trabaja con una lámina en 2
Se define utilizando integrales dobles a la masa total de una lámina con densidad  (x,y), como

19
 INTEGRALES MÚLTIPLES

m(R ) =  ( x, y) dx dy
R

  ( x, y) dx dy
masa
el cociente = R
se llama densidad media de la lámina.
area  dx dy
R

Si R no es una lámina, sino una figura geométrica de dos dimensiones, el anterior cociente se
llama promedio o valor medio de la función  sobre la región R.

2.6.2. Centro de gravedad de una lámina.

Se define como centro de gravedad ( x, y ) de una lámina a

 x ( x, y) dx dy  y ( x, y) dx dy
x= R
; y= R
m( R ) m( R )

Cuando la densidad es constante  (x,y) = c , entonces :

 x dx dy  y dx dy
x =c R
; y=c R

 dx dy  dx dy
R R

donde  dx dy es el área de R.
R
En este caso ( x , y ) se denomina centroide de la lámina (o de la región R).

2.6.3. Momento de inercia de una lámina.

L es una recta en el plano de la lámina y d(x,y) la distancia desde el punto (x,y) de R a la recta
L.

IL =  d ( x, y ) ( x, y ) dx dy
2

Se llama momento de inercia de la lámina respecto a L.

Nota: Si  (x,y) = 1, IL se denomina momento de inercia o segundo momento de la región R respecto

de L.
Los momentos de inercia respecto a los ejes x e y se designan por Ix e Iy

20
 INTEGRALES MÚLTIPLES

I x =  y 2  ( x, y) dx dy , I y =  x 2  ( x, y) dx dy
R R

La suma de estos dos se llama momento polar de inercia I0 respecto del origen: I0 = Ix + Iy

Ejemplo 9: Calcular el centroide de la región plana limitada por un arco de sinusoide,

con 0    .

Solución: Como (x,y) no es dato se considera constante.

Por simetría


x= = 1,57
2

 y dx dy
y= R
 dx dy
R

 senx  

  dx dy =  senx dx = − cos x
0 0 0
=2
0

 senx   
2 sen x
1 (sen 2x )  
   
y 1
y dx dy = dx = sen 2 x dx = x − =
2 2 4 2  4
0 0 0 0 0 0

Luego:

 y dx dy 
y= R
= = 0,39
 dx dy 8

21
 INTEGRALES MÚLTIPLES

3. INTEGRALES TRIPLES

3.1. DEFINICIÓN DE INTEGRALES TRIPLES

Siguiendo los cinco pasos utilizados para definir integrales dobles, se define la integral
definida tridimensional o integral triple. La que notamos como

 F(x, y, z) dV
R

1. Considere que F(x ,y ,z) está definida y acotada en una región R cerrada en el espacio.

2. Por medio de una red tridimensional de planos verticales y horizontales paralelos a los planos
coordenados, forme una partición P de R en subregiones (cajas) Rk de volúmenes Vk contenidas
totalmente en R, donde P es la norma de la partición, o sea la longitud de la diagonal mayor de los
Rk .

3. Elija un punto intermedio arbitrario ( x *k , y *k , z *k ) en cada subregión R k.

4. Forme la suma de Riemann:

n
 F ( x k , y k , z k ) Vk
* * *
k =0

5. Tome el límite de la sumatoria anterior cuando la norma de la partición tiende a cero para obtener
la definición buscada.

 F(x, y, z)dV = lim  F(x *k , y*k , z*k )Vk

P →0 k = 0
R

22
 INTEGRALES MÚLTIPLES

Definición: Sea f una función de tres variables definida y acotada en una región cerrada R del espacio.
Entonces la integral triple de F en R esta dada por

n
 F ( x, y, z) dV = lim  F( x k *, y k *, z k *) Vk
R P →0 k = 0

siempre y cuando el límite exista.

Se trabaja con funciones continuas sobre una región R   3 simplemente conexa del tipo:

R = {(x, y, z) : a  x  b ;  1 (x)  y   2(x) ;  1(x, y)  z   2(x, y)}

Supóngase que este dominio especial tridimensional R, limitado por una superficie cerrada S tiene
las siguientes propiedades:

1) Que toda paralela al eje z, trazada por un punto interior del dominio R (es decir por un punto que
no pertenece a la frontera S), corta a la superficie S en dos puntos.

2) Todo el dominio R se proyecta sobre el plano xy en forma de un dominio regular (de dos
dimensiones) D.

z z = 2 (x,y)

z = 1(x,y)

y
a
y = 1(x) y = 2(x)
b
x (x)

El volumen V es la diferencia del volumen que limita z =  2 (x ,y) menos el volumen que limita z
=  1 (x ,y). Luego el volumen total es igual a:

b  2 (x) b  2 (x) b  2 (x)

V=   2 (x, y) dy dx −   1 (x, y) dy dx =   [ 2 (x, y) − 1(x, y)] dy dx
a 1 ( x ) a 1 ( x ) a 1 ( x )

23
 INTEGRALES MÚLTIPLES

b  2 ( x )  2 ( x , y)
V=   dz dy dx
a  1 ( x ) 1 ( x , y)

Las integrales triples no tienen interpretación geométrica cuando el integrando no es constante, es

decir del tipo:

b  2 ( x )  2 ( x , y)
   f (x, y, z) dz dy dx
a 1 ( x )  1 ( x , y)

x y z
Ejemplo 10: Calcular el volumen comprendido entre los tres planos coordenados y + + =1
4 4 8

Solución: Despejando z, z = 8 - 2x - 2y

4 4 − x 8 − 2x − 2 y 4 4− x 4 4− x
2
V=   dz dy dx =   (8 − 2x − 2 y) dy dx =  8y − 2xy − y dx =
0 0 0 0 0 0 0

4
2
 [8 (4 − x ) − 2 x (4 − x ) − (4 − x ) ] dx
0
4 4
x3 64
V =  (16 − 8x + x 2 )dx = 16x − 4x 2 + =
0
3 3
0

24
 INTEGRALES MÚLTIPLES

3.2. Cambio de coordenadas en integrales triples.

3.2.1. Cambio de coordenadas en integrales triples. Coordenadas curvilíneas

Teorema:
Sea f = f(x ,y ,z) integrable en R   3 ,  una aplicación uno a uno tal que,  : 3→ 3 y
x = X(u, v, w )

 =  y = Y(u, v, w )
z = Z(u, v, w )


tal que J(u ,v ,w)  0 para todo (u ,v ,w)   -1 ( R ) = R' entonces se verifica:

 f ( x, y, z) dx dy dz =  f [ X (u, v, w), Y (u, v, w), Z (u, v, w)] J (u, v, w) du dv dw

R R

Observación:

1 - La fórmula anterior vale aún si, sobre un conjunto de 3de medida nula, J = 0

2 – Nótese que la expresión del diferencial volumen, dV, en el nuevo sistema de coordenadas es
  
dV = J (u, v, w) du dv dw = ru  rv  rw du dv dw (1)

Con du > 0, dv > 0, dw > 0 escalares positivos.

Se recuerda que el vector posición de la superficie es

r( u,v,w ) = x( u,v,w )i + y( u,v,w ) j + z( u,v,w )k , por lo que

   
ru = x u i + y u j + z u k

   
rv = x v i + y v j + z v k

   
rw = x w i + y w j + z w k

Reemplazando en (1) resulta:

25
 INTEGRALES MÚLTIPLES

Xu Yu Zu
  
dV = ru  rv  rw du dv dw = X v Yv Z v du dv dw
Xw Yw Zw

Este determinante es el jacobiano o factor de proporcionalidad de volúmenes,

Xu Yu Zu
J(u, v, w ) = X v Yv Zv
Xw Yw Zw

Existen dos sistemas de coordenadas curvilíneas particularmente importantes, las coordenadas

cilíndricas y las esféricas.

3.2.2. Cambio de coordenadas en integrales triples. Coordenadas cilíndricas

Se ubica un punto en el espacio y sus

correspondientes coordenadas , , z. z

x =  cos  P (x, y, z) o P (, , z)


La transformación  :  y =  sen 
z = z z

con 0    2 ,  > 0,  y

(, )
x

recibe el nombre de sistema de coordenadas cilíndricas, se calcula su Jacobiano y se ve donde es

válida dicha transformación.

x  y z cos  sen  0
 ( x, y, z)
J(, , z) = = x y z = −  sen   cos  0 =  ;
 (, , z)
x z yz zz 0 0 1

Se obtiene la expresión del jacobiano y del diferencial de volumen correspondiente

J(,,z) =  y dV =  d d dz

El jacobiano es distinto de cero (J(,,z)  0) si y solo si   0, luego para todo punto distinto del
origen existe  -1 y por lo tanto para estos puntos la terna (,,z) es un sistema de coordenadas.
Al aplicar la transformación tenemos:

26
 INTEGRALES MÚLTIPLES

 f (x, y, z)dxdydz =  f ( cos  ,  sen  , z)  d d dz

R R

Gráficamente:

d
z d
dz
dV=  d d dz

0  d
x d

3.2.3. Cambio de coordenadas en integrales triples. Coordenadas esféricas

x = r cos 

La transformación  :  y = r sen 
z =  cos  z


con r =  sen  (ver gráfica)

P (x, y, z) o (, , )
da lugar a la transformación 

 cos 
x =  sen cos  0    2 y

 y =  sen sen 0  r cos  r
z =  cos 
  0 r =  sen 
r sen 
x
llamada sistema de coordenadas esféricas.

En la siguiente figura se visualiza la transformación de un diferencial de volumen cuando se le aplica

coordenadas esféricas.

27
 INTEGRALES MÚLTIPLES

z
z

d

d

0 
y

 d

Se determina la expresión del jacobiano para las coordenadas esféricas

x y z sen  cos  sen  sen  cos 

J(, , ) = x  y z  = −  sen sen  sen cos  0
x y z  cos  cos   cos  sen  −  sen 

= − 2 sen3 cos2  −  2 sen sen2 cos2  − [  2 cos2  sen cos2 +  2 sen3 sen2 ]
= −  2 [sen 3  ( sen 2  + cos 2  ) + cos 2  ( sen  cos 2 + sen  sen 2  )]

= −  2 [sen 3  + cos 2  sen  (cos 2 + sen 2  )]

= −  2 [sen (cos2 + sen2 )] = −  2 sen , obteniendo finalmente

J (  ,  , ) =  2 sen

Si J(,,)  0 , entonces 2 sen   0 , por lo que sen   0 y   0.

Pero si  = 0 implica que son los puntos sobre el eje z (+) que es un conjunto de medida nula en 3
por lo que la formula de transformación es:

 f (x, y, z) dx dy dz =  f [ sen cos ,  sen sen,  cos ] 2 sen d d d
R R

Ejemplo 11:
a) Calcular el volumen de una esfera x 2 + y 2 + z 2 = a 2
utilizando coordenadas esféricas.

Solución:

28
 INTEGRALES MÚLTIPLES

z
2 a 
V =  dx dy dz =     2 sen  d d d
R 0 0 0 
2 a  2 a 
V =   −  2 cos  d d =   2   2 d d
0 0 0 0 0
 a y
a
2 2 2 2 4
V =  3 d = a 3  d =  a 3
30 0
3 0 3 x

b) Hallar el volumen de la región encima del plano xy limitado por el paraboloide

z = x² + y² y el cilindro x² + y² = a²

Solución: se trabaja con coordenadas cilíndricas, así la ecuación del paraboloide se transforma
en z = 2 y la ecuación del cilindro en 2 = a². z

2 a 2

V=
    dz
0 0 0
d d
z = x² + y² x² + y² = a²

y
2 a 2
2 a a
a4
V=

0 0
 z d d =

0 0
ρ dρ d = 2π
3
4
x

 4
V= a
2

La gráfica utilizando Maple es:

29
 INTEGRALES MÚLTIPLES

3.3. Aplicaciones de las integrales triples

3.3.1. Masa total de un sólido.

Las integrales triples pueden emplearse para calcular volúmenes, masas, centros de gravedad,
momentos de inercia y otros conceptos físicos asociados a sólidos. Si R es un sólido su volumen V
está definido por la integral triple:

V=
 dx dy dz
R

Si al sólido se le asigna una densidad  (x,y,z) en cada uno de sus puntos (x,y,z), su masa M es:

M=
  (x, y, z) dx dy dz
R

3.3.2. Centro de gravedad de un sólido.

Las coordenadas ( x , y , z ) del centro de gravedad de un sólido R se determinan por definición de la

siguiente forma:

1 1
M  M 
x= x  (x, y, z) dx dy dz y= y  (x, y, z) dx dy dz
R R
1
M 
z= z  (x, y, z) dx dy dz
R

3.3.3. Momento de inercia de un sólido.

El momento de inercia respecto al plano xy, Ixy , se determina por:

 z
2
I xy =  ( x, y, z) dx dy dz
R

existiendo fórmulas similares para I yz , I zx.

El momento de inercia respecto de una recta L, IL , se define como:

I L =  d 2 ( x, y, z)   (x, y, z) dx dy dz
R
donde d(x,y,z) representa la distancia de un punto genérico de R a la recta L.
30
 INTEGRALES MÚLTIPLES

Ejemplo 12: Hallar el momento de inercia con respecto al eje z, del Ejemplo 16 b),
suponiendo que la densidad es constante ,  (x,y,z) = cte.

Solución:
Iz =
 (x
R
2
+ y 2 ) ( x, y, z) dx dy dz

Como se trabaja en coordenadas cilíndricas

2 a 2
D a 6
Iz = D

0 0 0
 3 dz d d =
3

Resumen: En este curso hemos estudiado diferenciales de área y de volumen para diferentes tipos
de coordenadas.

Para trabajar con integrales dobles:

Xu Yu
Coordenadas curvilíneas dA = du dv
Xv Yv

Coordenadas Cartesianas dA = dx dy

Coordenadas Polares dA =  d d

Coordenadas Polares Generalizadas dA =  a bd d

Para trabajar con integrales triples:

Xu Yu Zu
Coordenadas curvilíneas dV = X v Yv Z v du dv dw
Xw Yw Zw

Coordenadas cilíndricas dV =  d d dz

Coordenadas esféricas dV =  2 sen  d d d

31