Integrales Multiples
Integrales Multiples
Integrales Multiples
Geraldine Cisneros
1. INTEGRALES DOBLES
En este trabajo se extiende el concepto de la integral de una funcin real de variable
real a funciones de varias variables, comenzando en este captulo con integrales de
funciones de dos variables; es decir, funciones del tipo f : D
. La integral
, xi 1 , xi ,
, x n 1 , x n }.
Sus contribuciones
destacaron en las reas de
anlisis y geometra
diferencial, la
fisicomatemtica y la
teora de funciones de
variables complejas.
Su nombre tambin est
relacionado con la
funcin zeta.
La
longitud
del
subintervalo genrico se
calcula de la siguiente
manera:
xi = xi xi 1
RP = f ( xi* ) xi
(I.1)
i =1
Geraldine Cisneros
Significado geomtrico
de la suma de Riemann
Si la funcin
es
[ ]
positiva
x a, b ,
entonces la suma de
Riemann corresponde a
un valor aproximado del
rea de la regin
comprendida bajo la
grfica de la funcin f ,
sobre el eje x, y entre las
rectas x = a y x = b .
Figura 1.1
Significado geomtrico de la Suma de Riemann para una funcin
positiva en el intervalo cerrado [a, b ] .
de la Integral Definida.
El smbolo lo introdujo
el matemtico alemn
Gottfried Wilhelm
f en [a ,b ]
von Leibniz
(1646, 1716).
b
a
f ( x ) dx = Lm f ( xi* )x
si el lmite existe.
p 0
i =1
(I.2)
Geraldine Cisneros
La
Integral
f (x )dx
b
Definida
es un nmero
real
que
puede
interpretarse como el rea
bajo la grfica de la
funcin f , sobre el eje
x y entre las rectas x = a
y x = b , si la funcin es
positiva.
f ( x ) es el
Una particin
Px
del
intervalo [a, b] es un
conjunto
finito
de
elementos, donde se
cumple:
{( x, y )
a x b c y d (I.3)
Px = {x 0 , x1 , x 2 , , xi 1 , xi , , x n 1 , x n }
(I.4)
Py = {y 0 , y1 , y 2 , , y j 1 , y j , , y m 1 , y m }
(I.5)
P = Px Py
(I.6)
entonces
m subintervalos y j 1, y j
de longitud y j = y j y j 1 , entonces la
Geraldine Cisneros
En la figura 1.2, se
aprecia que:
*
xi 1 xi xi
*
y j 1 y j y j
Figura 1.2
Particin
Figura 1.3
Subrectngulo
El
punto
(x
Dij
, y j ) Dij
Esquina inferior
izquierda
*
xi , y j * = xi 1 , y j 1
) (
Esquina inferior
derecha
xi * , y j * = xi , y j 1
Esquina superior
izquierda
xi * , y j * = xi 1 , y j
) (
) (
) (
(I.7)
x + x y j 1 + y j
, y j * ) = i 1 i ,
2
2
S D = f xi , y j Aij
i =1 j =1
Punto medio
(x
Aij = xi y j
Esquina superior
derecha
*
*
xi , y j = x i , y j
(I.8)
Geraldine Cisneros
)
)A
(
+ f (x
S D = f x1 , y1 A11 + f x1 , y 2 A12 +
*
1
f x 2 , y A21 + f x 2 , y 2
22
)
)A
+ f x1 , y m A1m +
f x n , y1 An1 + f x n , y 2 An 2 +
, ym
2m
(I.9)
+ f x n , y m Anm
P 0
i =1 j =1
(I.10)
f (x , y )dxdy
D
i =1 j =1
(I.11)
si el lmite existe.
f (x , y )dydx
D
f (x , y )dA = L
D
donde L
(I.12)
Geraldine Cisneros
Definicin del lmite de
una funcin:
es
El lmite
Lim f ( x ) = L
x x0
f (x
n
i =1 j =1
*
i
(I.13)
P <
(I.14)
, y j Aij L <
*
siempre que
0 < x x0 <
Siempre que:
0 < P no se coloca ya
que la norma
particin
de la
es
una
(x
*
i
, yj
) en el subrectngulo D
ij
es positiva.
Geraldine Cisneros
z = f (x , y )
(I.15)
z = f ( x, y )
Figura 1.4
Grfica de una funcin f :
10
Geraldine Cisneros
Figura 1.5
Slido S definido sobre la regin D y bajo la grfica de f
V Vij
(I.16)
i =1 j =1
Vij = f xi , y j Aij
(I.17)
f (x
n
i =1 j =1
*
i
11
Geraldine Cisneros
V f xi , y j Aij
i =1 j =1
(I.18)
(x
*
i
, y j* , f ( xi* , y j* )
z = f ( x, y )
D
y = y j 1
y = yj
x = xi 1
x = xi
Figura 1.6
Paraleleppedo de base
12
Geraldine Cisneros
Figura 1.7
Paraleleppedos empleados para aproximar el volumen del slido S
definido sobre la regin D
i =1 j =1
S , es decir:
n
i =1 j =1
(I.19)
13
Geraldine Cisneros
Figura 1.8
Paraleleppedos empleados en la aproximacin del volumen del slido S
con una particin refinada sobre
EJEMPLO 1.1
D.
z = x2 + 4y
arriba
del
rectngulo
Solucin:
a) Sea V el volumen del slido debajo de la superficie z = x 2 + 4 y
y arriba del rectngulo D . Entonces se desea estimar a V de la
siguiente manera:
2
i =1 j =1
14
Geraldine Cisneros
(x
*
i
m=3
n = 2,
entonces
x =
ba 20
=
=1
n
2
y =
d c 30
=
=1
m
3
) (
Como xi* , y j* = xi , y j
),
xi = x0 + ix = 0 + i(1) = i
y j = y0 + jy = 0 + j(1) = j
Figura 1.10
Particin empleada para el ejemplo1.1
Y adems:
Aij = xy = (1)(1) = 1
Recuerde:
n
k = kn
i =1 j =1
si k
i =1
i=
i =1
n
i
i =1
n ( n + 1)
2
n ( n + 1 )( 2 n + 1 )
6
V f ( i, j )(1) = ( i 2 + 4 j ) = ( 3i 2 + 24 ) = 15 + 48 = 63
i =1 j =1
i =1 j =1
i =1
V = ( x 2 + 4 y )dA 63
D
15
Geraldine Cisneros
Cuando
(x
se
selecciona
, y j = (xi , y j ) ,
*
en el
parte
aproximacin
a,
la
del
es inferior
Figura 1.11
Paraleleppedos empleados para aproximar el volumen del slido S
descrito en el ejemplo 1.1 parte a
(x
x + x y + y j i 1+ i j 1+ j 1
1
, y j * ) = i 1 i , j 1
,
=
= i , j
2
2 2
2
2
2
Luego:
2
3
2
3
1
1
V f ( xi* , y j* )Aij = f i , j (1)
2
2
i =1 j =1
i =1 j =1
(x
, y j * ) en la funcin
16
Geraldine Cisneros
(x
1
1
, y j* ) = i , j
2
2
f ( xi* , y j * ) = ( xi* ) + 4 y j *
2
i
j
1 1
,
2 2
3 1
,
2 2
9
4
17
4
1 3
,
2 2
3 3
,
2 2
25
4
33
4
1 5
,
2 2
3 5
,
2 2
41
4
49
4
Cuadro 1.1
Valores de f xi , y j
En el ejemplo 1.1 parte
b, cuando se selecciona
(x
, y j * ) como el punto
medio
de
cada
subrectngulo se puede
Por lo tanto
9 25 41 17 33 49
+ + + + +
= 43,5
4 4
4 4 4
4
V = ( x 2 + 4 y ) dA 43,5
D
atraviesa a
si
la
aproximacin
volumen del slido
del
es
Figura 1.12
Paraleleppedos empleados para aproximar el volumen del slido S
descrito en el ejemplo 1.1 parte b
17
Geraldine Cisneros
EJEMPLO 1.2
se
divide
subcuadrados,
en
n=m=2
x =
ba 40
=
=2
n
2
y =
d c 40
=
=2
m
3
V = ( 36 x 2 y 2 ) dA f ( xi* , y j * )Aij
entonces
i =1 j =1
donde xi , y j = (xi , y j )
*
Figura 1.14
Particin empleada para el ejemplo 1.2
18
Geraldine Cisneros
Aij = xy = ( 2 )( 2 ) = 4
) (
Como xi* , y j* = xi , y j
),
entonces
2
2
V f ( xi* , y j * )Aij = f ( 2i, 2 j )( 4 ) = 4 16 ( 2i ) ( 2 j )
i =1 j =1
i =1 j =1
i =1 j =1
xi = x0 + ix = 0 + i ( 2) = 2i
yj = y0 + jy = 0 + j ( 2) = 2 j
S.
Figura 1.15
Volumen aproximado en el ejemplo 1.2 parte a
19
Geraldine Cisneros
V = ( 36 x 2 y 2 ) dA f ( xi* , y j * )Aij
D
i =1 j =1
Diez mil
100
100
402, 7648
Un milln
1.000
1.000
405, 077248
f (x
i =1 j =1
, y j * )Aij
Cuadro 1.2
Aproximaciones del volumen del slido planteado en el ejemplo 1.2
Con
b, se aprecia que la
aproximaciones:
aproximacin
del
la
ayuda
del
software
se
obtuvo
V = ( 36 x 2 y 2 ) dA 402 , 7648
D
V = ( 36 x 2 y 2 ) dA 405, 077248
D
las
siguientes
20
Geraldine Cisneros
EJEMPLO 1.3
cada subrectngulo.
Solucin:
a) Sea V el volumen del slido debajo de la superficie z = 1 y y
arriba del rectngulo D .
La regin D y su particin se muestran en la siguiente figura
Figura 1.17
Particin empleada para el ejemplo 1.3 parte a
21
Geraldine Cisneros
Si
m=2
n = 3,
entonces
b a 30
x =
=
=1
n
3
y =
d c 1 ( 1)
=
=1
m
2
Aij = (x )(y ) = 1
i =1 j =1
y j = y0 + jy = 1+ j
, y j * ) es el punto
f xi , y j
y j 1 = y0 + ( j 1)y = j 2
(x
)= 1 y
*
j
2 j 3 5
= 1
= j
2 2
V j (1) = 2 = 6
i =1 j =1 2
i =1
V = (1 y ) dA 6
D
atraviesa
los
paraleleppedos, por lo
cual no se puede asegurar
si la aproximacin es por
exceso o por defecto.
Figura 1.18
Aproximacin del volumen para el ejemplo 1.3 parte a
22
Geraldine Cisneros
esta particin.
Figura 1.19
Particin empleada para el ejemplo 1.3 parte b
Si
n=6
m=8
entonces
x =
ba 1
=
n
2
y =
d c 1
=
m
4
A ij =
i =1 j =1
2 j 9 17 j
*
*
*
f xi , y j = 1 y j = 1
=
8 8 4
1
8
y j 1 = y0 + ( j 1)y =
y j = y0 + jy =
j 4
4
j 5
4
V = (1 y ) dA 6
D
23
Geraldine Cisneros
Figura 1.20
Aproximacin del volumen para el ejemplo 1.3 parte b
24
Geraldine Cisneros
Segundo Teorema
Fundamental del
Clculo
Si f es una funcin
continua en el intervalo
cerrado [a,b] y si F es
una antiderivada de
entonces:
de
la
funcin
f ( x, y )dxdy
sobre
D,
denotada
por
d
b
f ( x, y ) dxdy = f ( x, y ) dx dy
a
c
(I.20)
b
d
f ( x, y ) dydx = f ( x, y ) dy dx
c
a
(I.21)
O tambin
25
Geraldine Cisneros
Recuerde
integral
que
en
f ( x, y) dx , la
a
d
c
f ( x,y) dx = A( y)
f ( x, y ) dy se integra
parcialmente respecto a
la variable y; es decir x es
considerada constante.
(I.22)
Finalmente:
f ( x, y ) dxdy = f ( x, y ) dx dy =
En
forma
anloga,
en
la
expresin
(I.21),
A ( y ) dy (I.23)
la
integral
f ( x,y) dy = A( x)
(I.24)
EJEMPLO 1.4
b
d
b
f ( x, y ) dydx = f ( x, y ) dy dx = A ( x ) dx
a
c
a
(I.25)
a)
c)
e)
(x
2
+ 4 y ) dxdy
( 36 x
1
1
y 2 )dxdy
(1 y ) dxdy
0
b)
d)
f)
3
0
(x
3
+ 4 y ) dydx
( 36 x
4
y 2 )dydx
(1 y ) dydx
1
26
Geraldine Cisneros
Recuerde
que
en
el
V = ( x + 4 y ) dA 63
D
Al
comparar
estas
puede
8
2
0 3 + 8 y dy = 3 y + 4 y
8
= + 8y
3
= 8 + 36 = 44
0
Por lo tanto:
(x
valor de la integral, en
se
aproximaciones con el
efecto
x3
x
+
y
dx
=
4
) 3 + 4 xy
0 (
V = ( x + 4 y ) dA 43,5
(x
parcialmente respecto a x,
maneras diferentes:
3
0
mejor
+ 4 y ) dxdy = 44
aproximacin.
b) Se desea resolver
(x
3
( 3x
2
2
0
(x
3
+ 4 y ) dydx :
+ 4 y ) dy = x 2 y + 2 y 2
+ 18 ) dx = x 3 + 18 x
(x
2
2
0
3
0
= 3 x 2 + 18
= 8 + 36 = 44
+ 4 y ) dydx = 44
4
0
( 36 x
4
y 2 )dxdy , primero se
x3
y
dx
=
x
y2 x
36
36
(
)
0
3
= 144
64
368
4 y2 =
4 y2
3
3
27
Geraldine Cisneros
Recuerde
que
en
4 3
368
368
2
0 3 4 y dy = 3 y 3 y
4
=
0
=
405,3
3
3
3
Por lo tanto:
( 36 x
b, se obtuvo:
y 2 )dxdy = 405,3
V 402,7648 y
V 405,077248
la
integral
( 36 x
4
doble,
d) Resolviendo
y 2 )dxdy = 405,3
( 36 x
4
y 2 )dydx
1 3
2
0 ( 36 x y ) dy = 36 y x y 3 y
4
= 144 4 x 2
64 368
=
4 x2
3
3
la estimacin de la parte
a.
4 3
368
368
2
0 3 4 x dx = 3 x 3 x
4
( 36 x
En el ejemplo 1.3, se
aproxim
doble
la
integral
mediante
dos
diferentes,
=
0
1472 256
= 405,3
3
3
y 2 )dydx = 405,3
(1 y ) dxdy ,
1
primero se integra
una
0 (1 y ) dx = (1 y ) x
particiones
donde
3
0
= 3 (1 y )
en
V = (1 y ) dA 6
D
3
2
1 3 (1 y )dy = 2 [1 y ] 1 = 6
1
Es decir:
28
Geraldine Cisneros
(1 y ) dxdy = 6
f) Ahora se resuelve
(1 y ) dydx
en el orden de integracin
1 3
y2
= = 2
1
y
dx
y
(
)
1
2 1 2 2
3
0
2dx = 6
(1 y ) dydx = 6
29
Geraldine Cisneros
El nombre de Teorema
de Fubini se debe al
matemtico italiano:
Guido Fubini
(1879, 1943).
f ( x, y ) dA =
D
f ( x, y ) dxdy =
b
a
d
c
f ( x, y ) dydx (I.25)
Demostracin intuitiva:
Considere que la funcin f es positiva, es decir, f (x, y ) 0 , por lo
El principio de Cavalieri
se debe al matemtico
italiano
Bonaventura
FrancescoCavalieri
(1598, 1647).
f ( x, y ) dA
D
V = A ( x ) dx
a
(I.26)
30
Geraldine Cisneros
z = f(x,y)
A(x)
a
x = x0
b
Figura 1.21
Interpretacin geomtrica del Teorema de Fubini
A ( x ) = f ( x, y ) dy
(I.27)
b
a
d
c
f ( x, y )dydx
(I.28)
V = A ( y ) dy
c
(I.29)
31
Geraldine Cisneros
Figura 1.22
A ( y ) = f ( x, y ) dx
A( y)
Interpretacin
geomtrica del
teorema de Fubini
A( y)
seccin transversal del
slido
S que es
perpendicular al eje y y al
plano xy .
(I.30)
es el rea de la
d
c
b
a
f ( x, y )dxdy
(I.31)
f ( x, y ) dA =
D
f ( x, y ) dydx =
d
c
b
a
f ( x, y ) dxdy (I.32)
32
Geraldine Cisneros
Figura 1.23
Regin D con una forma ms general
{ ( x, y )
a xb
g ( x) y h ( x)
(I.33)
33
Geraldine Cisneros
Figura 1.24
Regiones de tipo 1
{ ( x, y )
g ( y) x h( y) c y d
(I.34)
respectivamente.
En la figura 1.25 se aprecian algunas regiones de tipo 2.
34
Geraldine Cisneros
Figura 1.25
Regiones de tipo 2
y = 1 x2
la siguiente definicin:
Sea f :
Figura 1.27
El crculo unitario
como una regin
tipo 2
si ( x, y ) D
f ( x, y )
F ( x, y ) =
0 si ( x, y ) D ( x, y ) R
La integral doble de f sobre D , denotada
(I.35)
f ( x, y ) dA , est
D
dada por:
f ( x, y )dA = F ( x, y )dA
D
f ( x, y )dA ,
D
(I.36)
se debe
35
Geraldine Cisneros
Figura 1.28
Rectngulo R que contiene a la regin D de tipo 1
resulta:
b
F ( x, y )dA =
R
F ( x, y ) dydx
(I.37)
d
c
F ( x, y ) dy =
h( x )
g ( x)
F ( x, y ) dy =
h( x )
g ( x)
f ( x, y ) dy
(I.38)
36
Geraldine Cisneros
D=
{ ( x, y )
g ( x) y h ( x)
a xb
(I.39)
h( x )
g ( x)
f ( x, y )dA =
D
f ( x, y )dydx
(I.40)
x = h(y)
D
R
c
x = g(y)
a
Figura 1.29
Rectngulo
tiene:
d
F ( x, y )dA =
R
F ( x, y ) dxdy
(I.41)
37
Geraldine Cisneros
b
a
F ( x, y ) dx =
h( y )
g( y)
F ( x, y ) dx =
h( y )
g( y)
f ( x, y ) dx
(I.42)
D=
{ ( x, y )
g ( y) x h( y) c y d
(I.43)
h( y )
g( y)
f ( x, y )dA =
D
COMENTARIO
f ( x, y )dxdy
(I.44)
h( x )
g ( x)
38
Geraldine Cisneros
: Indica cual es el
valor de la variable
y a la salida de la
regin D (lmite
superior).
: Indica cual es el
valor de la variable
y a la entrada de la
regin D (lmite
inferior).
Figura 1.30
Orden de integracin para la integral doble de
Por otra parte, la ecuacin (I.44) seala que en una regin de tipo
2, la integral doble de la funcin
d
h( y )
g( y)
se obtiene como
y
: Indica cual es el
valor de la variable
x a la salida de la
regin D (lmite
superior).
: Indica cual es el
valor de la variable
x a la entrada de la
regin D (lmite
inferior).
x = h(y)
D
R
c
x = g(y)
a
Figura 1.31
Orden de integracin para la integral doble de
39
Geraldine Cisneros
EJEMPLO 1.5
1
0
2
0
x2
0
b)
dydx
4 y2
4 y
d)
dxdy
3 x +1
2x
ey
dydx
xdxdy
Solucin:
a) Para resolver la integral
1
0
x2
0
Figura 1.32
Slido del ejemplo 1.5
parte a
x2
1
x
dydx = dy dx = y
0
0
x2
x2
0
dydx =
3
dx = 1 x 2 dx = x
0
=
0
1
3
1
3
Regin tipo 1: D =
{( x, y )
Regin tipo 2: D =
{( x, y )
0 x 1 0 y x2
y x 1 0 y 1
40
Geraldine Cisneros
y = x2
En el ejemplo 1.5 a la
integral 1 x dydx = 1 , lo
unidades de longitud.
Valor de y a
la entrada de D
y=0
Figura 1.34
3 x +1
2x
2
1
3 x +1
2x
dydx
2
3 x +1
2
2
3 x +1
dydx =
dy dx = y 2 x dx = ( x + 1)dx
2x
1
1
1
( x + 1)
1 ( x + 1)dx =
Figura 1.35
Slido del ejemplo 1.5
parte b
3 x +1
2x
dydx =
5
2
5
2
D=
{( x, y )
1 x 2
2 x y 3x + 1
41
Geraldine Cisneros
y = 3x + 1
x=2
x =1
D
Valor de y a
la entrada de D
y = 2x
Figura 1.37
Regin
D del ejemplo1.5 b
Figura 1.38
Slido del ejemplo 1.5
parte c
Esta
4 y2
4 y2
2
4 y 2
2
dxdy
=
0 4 y2 dx dy = 0 x
4 y2
4 y2
integral
se
resuelve
dxdy , se tiene:
4 y2
4 y2
dy = 2 2 4 y 2 dy
0
empleando
una
sustitucin
trigonomtrica:
2
4 y 2
4 y
Figura 1.39
Funcin f definida en
dxdy = 2
2
4 y 2 dy = 4
la regin D del
ejemplo 1.5 parte c
2
0
y
1 dy
2
42
Geraldine Cisneros
4 y2
4 y2
dxdy = 4
2
0
dy = 2 cos d
y = 2 = arcsen1 =
1 + cos ( 2 )
sen ( 2 )
dxdy = 8 2
d = 4 +
2
4 y
0
2
2
4 y2
4 y2
4 y 2
= 2
dxdy = 2
4 y2
4 y2
Regin tipo 1: D =
{( x, y )
2 x 2 0 y 4 x2
Regin tipo 2: D =
{( x, y )
4 y2 x 4 y2
}
}
0 y2
dxdy = 2
Valor de x a
la entrada de D
Valor de x a
la salida de D
x = 4 y2
x = 4 y2
Figura 1.40
Regin D del ejemplo 1.5 c
y =0
43
Geraldine Cisneros
ey
d) La integral
ey
ey
Figura 1.41
Slido del ejemplo 1.5
parte d
e
xdxdy =
0
y
1
12 3
xdx dy = x 2
0 3
2 2 3y
2 5
xdxdy = e 2 y 2
3 3
5
ey
Figura 1.42
Funcin f definida en
=
0
ey
y
1 2 3y
3
dy = 0 e 2 y 2 dy
3
2 2 3 2 2 2 4 3 2 32
e = e
3 3
5 3 9
45
4 3 32
xdxdy = e 2
9
45
la regin D del
ejemplo 1.5 parte d
D=
{( x, y )
y x ey 0 y 1
Valor de x a
la entrada de D
Valor de x a
la salida de D
x= y
x = ey
y=0
Figura 1.43
Regin D del ejemplo 1.5 d
44
Geraldine Cisneros
general D.
1.7.1 Propiedad de linealidad
2
Sean f :
y g:
f ( x, y ) + g ( x, y ) dA = f ( x, y ) dA + g ( x, y ) dA
D
D
(I.45)
Sean f :
y g:
f ( x, y )dA g ( x, y )dA
D
(I.46)
f ( x, y )dA =
D
D1
f ( x, y )dA + f ( x, y )dA
D2
(I.47)
45
Geraldine Cisneros
EJEMPLO 1.6
( x + y + 1) dA ,
D
D=
{ ( x, y ) y x
+ 2x
y3
y 3x + 6
Solucin:
El primer paso para resolver este ejercicio es identificar si la regin
D es tipo 1 o tipo 2. En la siguiente figura se muestra la regin D .
Figura 1.44
Funcin f definida en
la regin D del
ejemplo 1.6
y=3
y = 3x+6
y = x2+2x
Figura 1.45
Regin D del ejemplo 1.6
{( x, y )
D2 =
2 x 1 x 2 + 2 x y 3 x + 6 y
{( x, y )
1 x 1 x2 + 2 x y 3
46
Geraldine Cisneros
Valor de y a
la salida de D2
y =3
Valor de y a
la salida de D1
y = 3x + 6
D2
D1
Valor de y a
la entrada de D1
y = x2 + 2 x
Valor de y a
la entrada de D2
y = x2 + 2 x
Figura 1.46
Por lo tanto:
I = ( x + y + 1) dA =
D
3 x+6
( x + y + 1)dydx + 1 x + 2 x ( x + y + 1)dydx
2 x + 2 x
2
1
1 15
x4
5x2
x4
I = 24 3 x3 +
+ 25 x dx + 3 x3 + 5 x 2 + x dx
2
1 2
2
2
2
I=
29 172
+
60 15
I = ( x + y + 1) dA =
D
239
20
{( x, y )
1 1 + y x 1 + 1 + y
DB = ( x, y )
y6
x 1 + 1 + y 0 y 3
3
1 y 0 y
47
Geraldine Cisneros
Valor de x a
la entrada de DB
Valor de x a
la salida de DB
y 6
x=
3
x = 1 + 1 + y
DB
y=0
Valor de x a
la entrada de DA
x = 1 1 + y
Valor de x a
la salida de DA
DA
x = 1 + 1 + y
Figura 1.47
Regin D del ejemplo 1.6 dividida en dos regiones tipo 2
I =
0
1
1+ 1+ y
1+ 1+ y
( x + y + 1)dxdy + 0 y 6
1 1+ y
( x + y + 1)dxdy
Resolviendo se obtiene I =
8 749
+
, luego:
15 60
I = ( x + y + 1) dA =
D
239
20
48
Geraldine Cisneros
EJEMPLO 1.7
(10 + 4 x
y ) dA , D =
{ ( x, y )
yx 2
x2 + y 2 4
Solucin:
Tal como se explic en los ejemplos anteriores, el primer paso
para resolver la integral doble planteada consiste en clasificar a la
Figura 1.48
Funcin f definida en
la regin D del
ejemplo 1.7
inecuaciones:
yx =
x y si
yx
y<x
( x, y )
que
se
encuentran
dentro
sobre
la
y = x+2
x2 + y 2 = 4
y = x2
Figura 1.49
Regin
49
Geraldine Cisneros
x=0
y = 4 x2
Valor de y a
la salida de D1
y = x+2
D2
D1
Valor de y a
la entrada de D1
Valor de y a
la entrada de D2
y = x2
y = 4 x2
Figura 1.50
Regin D del ejemplo 1.7 dividida en dos regiones tipo 1
{( x, y ) 2 x 0 4 x y x + 2} y
D = {( x, y ) 0 x 2 x 2 y 4 x }
Donde: D1 =
Por lo tanto:
I = (10 + 4 x 2 y ) dA = (10 + 4 x 2 y ) dA +
D
D1
I =
0
2
x+2
4 x
D2
10 + 4 x 2 y )dydx +
2 (
2
0
4 x2
x2
(10 + 4 x
(10 + 4 x
y ) dA
y )dydx
(8x + 20 + 10 4 x + 4 x + 7 x + 4 x 4 x ) dx +
+ ( 12 x + 20 + 10 4 x 4 x + 9 x + 4 x 4 x ) dx
I =
50
Geraldine Cisneros
I = 14 + + (14 + 24 )
3
(10 + 4 x
EJEMPLO 1.8
y ) dA = 28 +
152
3
x + y dA , D =
{ ( x, y )
y0
x2 + y 2 9
Solucin:
La regin D es una regin tipo 1 tal como se muestra en la
siguiente figura.
y = 9 x2
Figura 1.51
Funcin f definida en
la regin D del
ejemplo 1.8
y=0
Figura 1.52
Regin
f ( x, y ) = x + y =
( x + y ) si x + y < 0
51
Geraldine Cisneros
y = 9 x2
D1
y = -x
x+ y 0
D2
x+ y <0
y=0
Figura 1.53
Interseccin de la regin
Entonces se tiene:
x + y si ( x, y ) D1
f ( x, y ) =
( x + y ) si ( x, y ) D
2
Donde: D1 =
{ ( x, y )
D2 =
{ ( x, y )
y0
y0
x2 + y2 9
x+ y 0 y
x2 + y 2 9
x+ y <0
D1
( x + y ) dA + D
( x + y ) dA
y = 9 x
D1 = D1.A D1.B
x=0
Valor de y a
la salida de D1 .B
y = 9 x2
3 3
,
2 2
D1.A
D1.B
Valor de y a
la entrada de D1.A
y = x
Valor de y a
la entrada de D1.B
Figura 1.54
y=0
52
Geraldine Cisneros
D1
( x + y ) dA =
3
2
( x + y )dydx +
3
0
9 x2
0
( x + y )dydx
3
9
9 x2
9 x 2 + dx + x 9 x 2 + dx
0
2
2 2
D ( x + y ) dA = 3 x
1
9 x2
( x + y ) dA = ( 9
D1
2 9 + 18
2 2
D2
Valor de x a
la entrada de D2
Valor de x a
la salida de D2
x = y
x = 9 y2
y=0
Figura 1.54
Regin D2 del ejemplo 1.8
As:
D2 ( x + y ) dA = 0
3
2
D2
y
9 y
( x + y )dxdy = 0 2
( x + y ) dA = 9 2 9
Por lo tanto
I = x + y dA = 18 2
D
y 9 y 2 dy
2
53
Geraldine Cisneros
2. INTEGRALES TRIPLES
En esta seccin se presenta la integral triple para funciones de tres variables, funciones
del tipo f : B
dobles. As como se define la integral triple a partir de una triple suma de Riemann y se
ilustra el proceso de resolucin de la misma, de manera similar se puede esbozar la
definicin y el clculo de integrales mltiples de funciones del tipo f : Q
2.1 INTEGRAL
TRIPLE
SOBRE
UNA
CAJA
RECTANGULAR
Sea f una funcin definida sobre la caja rectangular B , esto es
f :B
(II.1)
o tambin:
B=
La caja rectangular B ,
tambin recibe el nombre
3
, O
de rectngulo en
intervalo tridimensional,
aunque el nombre ms
apropiado para B es
paraleleppedo.
{( x, y,z )
a xb c yd r zs
(II.2)
(II.3)
Py = {y 0 , y1 , y 2 , , y j 1 , y j , , y m 1 , y m }
(II.4)
(II.5)
entonces
54
Geraldine Cisneros
(II.6)
Figura 2.1
Particin
[ z k 1 ,zk ]
de longitud
volumen
de
cada
una
de
estas
pequeas
cajas
55
Geraldine Cisneros
(II.7)
(x
*
i
, y j* ,zk* ) del
(II.8)
i =1 j =1 k =1
En la figura 2.2, se
aprecia que:
xi 1 xi* xi
y j 1 y j * y j
zk 1 zk * zk
Figura 2.2
Paraleleppedo genrico Bijk de la particin P .
56
Geraldine Cisneros
P 0
(II.9)
i =1 j =1 k =1
en un
paraleleppedo B .
si el lmite existe.
i =1 j =1 k =1
(II.10)
57
Geraldine Cisneros
una
funcin
continua
en
el
paraleleppedo
f ( x, y,z ) dV =
B
Al igual que en el
captulo anterior; para la
resolucin de integrales
triples, se emplearn los
siguientes smbolos para
identificar los lmites de
integracin:
: Valor de la variable
a la salida de la
regin B (lmite
superior).
: Valor de la variable
a la entrada de la
regin B (lmite
superior).
f ( x, y,z )dxdydz
(II.11)
f ( x, y,z ) dV =
f ( x, y,z )dxdzdy
(II.12)
f ( x, y,z )dydxdz
(II.13)
f ( x, y,z )dydzdx
(II.14)
(II.15)
f ( x, y,z ) dV =
B
f ( x, y,z ) dV =
B
(II.16)
58
Geraldine Cisneros
EJEMPLO 2.1
f ( x, y,z )dV
y dibuje el paraleleppedo
z= 2
Es comn llamar I a la
integral triple que desea
resolverse.
Valor de z a
la entrada de B
z=0
Figura 2.3
Paraleleppedo B del ejemplo 2.1.
I = f ( x, y, z )dV =
B
plano xy
La proyeccin de B
mostrada en la figura 2.4,
indica que la segunda
integracin se realiza
respecto a la variable x.
I =
I =
1
-2
3
2
1
xz 4 (1 y )
4
1
-2
2
0
2
0
xz 3 (1 y )dzdxdy
dxdy =
1
-2
x (1 y )dxdy
2
3
1 2
5 1
5
2
x (1 y ) dy = (1 y )dy = (1 y )
-2 2
2
2 -2
4
1
1
-2
45
4
59
Geraldine Cisneros
-2
EJEMPLO 2.2
2
0
xz 3 (1 y )dzdxdy =
f ( x, y,z )dV
45
4
y dibuje el paraleleppedo
B , donde f ( x, y,z ) = x + y cos z y B = [ 1, 2] [ 0,1] 0 , .
2
Solucin:
En la figura 2.5 se muestra el paraleleppedo B y se indica,
adems, que la primera integracin parcial se realiza respecto a la
variable x.
Valor de x a
la entrada de B
x = 1
Figura 2.6
Proyeccin del
paraleleppedo B del
ejemplo 2.2 en el
plano yz
Valor de x a
la salida de B
La proyeccin de B en
el plano yz, muestra que
la segunda integracin
parcial se realiza respecto
a la variable z.
x=2
Figura 2.5
Paraleleppedo
I = f ( x, y, z )dV =
B
I =
1
0
2
0
x2
2 + xy cos z
1
0
( x + y cos z )dxdzdy
2
0
dzdy =
1
0
2
0
2 + 3 y cos z dzdy
60
Geraldine Cisneros
1 3
3
3
2
I = z + 3 ysenz dy = + 3 y dy = + y
0 2
0
24
0
4
1 3
3
3
2
I = z + 3 ysenz dy = + 3 y dy = + y
0 2
0
4
2
4
3 3
( x + y cos z )dxdzdy = 4 + 2
1
2
0
2 1
2 1
61
Geraldine Cisneros
Figura 2.6
Regin general B del espacio tridimensional
sobre la
62
Geraldine Cisneros
z = 2 ( x, y )
D
Valor de z a
la entrada de B
z = 1 ( x, y )
Figura 2.7
Primer orden de integracin para una regin general B
B=
{( x, y,z ) ( x, y ) D
1 ( x, y ) z 2 ( x, y )
(II.17)
sobre la regin
f ( x, y,z ) dV =
2 ( x ,y )
(
1
x ,y )
f ( x, y,z )dzdA
(II.18)
63
Geraldine Cisneros
y = g2 ( x )
x=b
x=a
Valor de y a
la entrada de D
y = g1 ( x )
Figura 2.8
Proyeccin de la regin general B sobre el plano xy
B=
{( x, y,z ) a x b,
g1 ( x ) y g 2 ( x ) , 1 ( x, y ) z 2 ( x, y )
(II.19)
Se tiene que:
f ( x, y,z ) dV =
b
a
g2 ( x )
g1 ( x )
2 ( x ,y )
(
1
x ,y )
f ( x, y,z )dzdydx
(II.20)
B=
{( x, y,z ) h ( z ) x h ( z ) ,
1
1 ( x,z ) y 2 ( x,z ) , r z s}
(II.21)
Entonces:
s
h2 ( z )
h1 z
f ( x, y,z ) dV =
B
2 ( x ,z )
( ) (
1
x ,z )
f ( x, y,z )dydxdz
(II.22)
B=
{( x, y,z ) ( y,z ) x
1
( y,z ) ,
c y d,
j1 ( y ) z j2 ( y )
64
Geraldine Cisneros
j2 ( y )
j1 ( y )
f ( x, y,z ) dV =
B
EJEMPLO 2.3
2 ( y ,z )
(
1
y ,z )
f ( x, y,z )dxdzdy
(II.24)
B=
{( x, y,z ) 0 x 2,
x y x2 , x + y z x2 + y 2
Solucin:
Para evaluar
Valor de z a
la salida de B
z = x2 + y 2
En la figura 2.9, se
aprecia que el recinto B
est
limitado
superiormente por la
superficie z = x 2 + y 2 e
inferiormente por la
superficie z = x + y .
Valor de z a
la entrada de B
z = x+ y
Figura 2.9
Regin
65
Geraldine Cisneros
z, considerando a x y a y constantes.
En la figura 2.10 se muestra la proyeccin de la regin B sobre el
plano xy. Adicionalmente se ilustra el segundo orden de
integracin seleccionado.
Valor de y a
la salida de D
y = 4x
x=2
Cuando se proyecta la
regin B sobre el plano
xy, tal como se muestra
en la figura 2.10, se
obtiene una regin D
bidimensional de tipo 1.
Valor de y a
la entrada de D
y = x2
Figura 2.10
Proyeccin de la regin B del ejemplo 2.3 en el plano xy
I = dV =
B
2
0
x2
x2 + y 2
4x
x+ y
dzdydx =
2
0
(y
x2
4x
+ x 2 x y )dydx
x 6 x 4 79 x3
I = +
12 x 2 dx
0
3
3 2
x2
x2 + y 2
4x
x+ y
UC. Facultad de Ingeniera. Departamento de Matemtica.
dzdydx =
6724
105
66
Geraldine Cisneros
EJEMPLO 2.4
dV , donde
es la regin
x + y + z = 10 .
Solucin:
dV , es necesario ilustrar el
En la figura 2.11, se
aprecia que el recinto B
est limitado
por la
izquierda por el plano
cartesiano xz y por la
derecha por el plano de
ecuacin y = 10 x z
Valor de y a
la salida de B
Valor de y a
la entrada de B
y = 10 x z
y=0
Figura 2.11
Regin
dV .
67
Geraldine Cisneros
Valor de z a
la salida de D
z = 10 x
x=0
En la figura 2.12 se
observa
que
la
proyeccin de la regin B
sobre el plano xz, es una
regin D de tipo 1 o tipo
2; sin embargo, se trabaja
como una regin tipo 1.
D
Valor de z a
la entrada de D
z=0
Figura 2.12
Proyeccin de la regin
10
0
10 x
0
10 x z
0
dydzdx =
10
0
10 x
0
10
x2
I = 50 + 10 x dx
0
2
10
10 x
10 x z
0
dydzdx =
500
3
(10 x z )dzdx
68
Geraldine Cisneros
Sean f :
y g:
f ( x, y, z ) + g ( x, y, z ) dV
=
+
f ( x, y, z )dV
g ( x, y, z )dV
B
(II.25)
Sean f :
definidas
en
y g:
una
regin
B,
tales
que
f ( x, y, z )dV g ( x, y, z )dV
B
(II.26)
f ( x, y, z )dv = f ( x, y, z )dV +
B
B1
B2
f ( x, y, z )dV
(II.27)
69
Geraldine Cisneros
EJEMPLO 2.5
y B
{( x, y,z ) x
+ y 2 + z 2 4 , x 2 + y 2 1, x 0 , y 0 , z 0
Solucin:
El recinto B es la regin del primer octante que se encuentra
dentro de la esfera, de radio 4 y centro en el origen del sistema de
coordenadas; y fuera del cilindro circular recto de radio 1 y que
tiene como eje directriz al eje z. En la figura 2.13 se muestra el
recinto B .
Valor de z a
la salida de B
z = 4 x2 y 2
En la figura 2.13, se
muestra el recinto B, del
ejemplo 2.5. Esta regin
est
acotada
superiormente por la
superficie de ecuacin
e
z = 4 x2 y 2
inferiormente por
plano
cartesiano
( z = 0 ).
el
xy
B
Valor de z a
la entrada de B
z=0
Figura 2.13
Regin
4 x2 y 2
0
xyzdzdA
70
Geraldine Cisneros
y = 4 x2
Cuando se proyecta la
regin B sobre el plano
xy, se obtiene una regin
D que no es de tipo 1 ni
de tipo 2, por lo que se
divide en dos subregiones
tipo 1.
x =1
Valor de y a
la salida de D2
D1
y = 4 x2
D2
Valor de y a
la entrada de D1
y = 1 x2
Valor de y a
la entrada de D2
y=0
Figura 2.14
Luego, D = D1 D2 , donde:
{( x, y ) 0 x 1
D = {( x, y ) 1 x 2
D1 =
1 x2 y 4 x2
0 y 4 x2
I =
1
0
4 -x 2
1-x
1
0
4 -x 2
1-x 2
4 x2 y 2
0
xyzdzdydx +
2
xy
2
2
4
y
dydx
+
(
)
1 0
2
4 -x 2
2
1
4 -x 2
0
4 x2 y 2
0
xy
2
2
2 ( 4 x y ) dydx
5
1 9x
2 x
9 9
I = dx + x 3 + 2 x dx = +
0
1
16 16
8
8
xyzdV =
9
8
xyzdzdydx
71
Geraldine Cisneros
EJEMPLO 2.6
( xyz ) dV ,
B
Valor de z a
la salida de B
z = 2 ( x2 + y 2 )
z=4
En la figura 2.15, se
muestra el recinto B.
Observa que la flecha
que se encuentra a la
izquierda sale de la
regin en el plano de
ecuacin z = 4 , mientras
que la flecha de la
derecha sale de la regin
por la
superficie del
cono z = 2 ( x 2 + y 2 ) .
Valor de z a
la entrada de B
Valor de z a
la entrada de B
z = x2 + y 2
z = x2 + y2
Figura 2.15
Regin
( xyz ) dV =
D1
4
x2 + y
( xyz ) dzdA +
2
D2
2 x2 + y 2
x2 + y 2
( xyz ) dzdA
72
Geraldine Cisneros
y = 16 x 2
La regin D1, de la figura
2.16, no es de tipo 1 ni de
tipo 2, por lo que se
divide en dos subregiones
tipo 1.
x= 8
D1.A
Valor de y a
la salida de D1.B
y = 16 x 2
Valor de y a
la entrada de D1.A
D1.B
y = 8 x2
Valor de y a
la entrada de D1.B
y=0
Figura 2.16
Primera proyeccin de la regin
{( x, y ) 0 x 8
D = {( x, y )
8x4
D1.A =
8 x 2 y 16 x 2
0 y 16 x 2
1.B
4
x2 + y 2
D1
( xyz ) dzdA
16 x 2
8 x 2
16 x 2
4
x2 + y 2
( xyz ) dzdydx +
4
x2 + y 2
( xyz ) dzdydx
D2
2 x2 + y2
x2 + y2
integral
en la siguiente figura.
73
Geraldine Cisneros
y = 8 x2
x=0
D2
8 ,0
Valor de y a
la entrada de D2
y=0
Figura 2.17
Donde:
D2 =
{( x, y )
0 x 8 0 y 8 x2
Resultando:
D2
2 x2 + y 2
x2 + y 2
8 x2
( xyz ) dzdA =
2 x2 + y 2
x2 + y 2
( xyz ) dzdydx
Por lo tanto:
I = ( xyz ) dV
B
+
Resolviendo, se tiene:
16 x 2
8 x
16 x 2
0
8 x 2
x2 + y 2
x2 + y 2
( xyz ) dzdydx
( xyz ) dzdydx
2 x2 + y 2
x2 + y 2
( xyz ) dzdydx
74
Geraldine Cisneros
I=
16 x 2
8 x 2
8 x 2
8
0
xy
(16 y 2 x 2 )dydx +
2
xy 2
( y + x 2 ) dydx
2
8 xdx +
16 x 2
x5
4 x 3 + 32 x dx +
8
8
I = 32 +
32 64
+
= 64
3
3
I = ( xyz ) dV = 64
B
8
0
xy
(16 y 2 x2 ) dydx +
2
x5
+ 8 x dx
8
75
Geraldine Cisneros
f ( x, y ) dA ,
D
f ( x, y ) dA =
D
dA
(III.1)
que:
i =1 j =1
dA = Lim Aij
P 0
i =1 j =1
(III.2)
76
Geraldine Cisneros
xi
(xi*,yj*)
d = ym
yj
yj-1
Dij
yj
c = y0
a = x0
xi-1
xi
xn= b
Figura 3.1
Regin D dividida en subrectngulos Dij
dA representa el volumen de un
. Sea A el
(III.3)
77
Geraldine Cisneros
Recuerde que una regin
D es de tipo 1 si se
cumple:
( x, y ) a x b
D=
f ( x ) y g ( x )
A=
b
a
g( x)
f ( x)
dydx =
b
a
g ( x)
[ y ] f ( x ) dx
(III.3)
A = g ( x ) f ( x ) dx
a
(III.4)
EJEMPLO 3.1
dxdy y
dydx , D =
{ ( x, y ) x y
2y
x 4 y2
Solucin:
La regin D se encuentra acotada por las grficas de las
parbolas horizontales x = y 2 2 y y x = 4 y 2 , tal como se puede
observar en la siguiente figura.
x = y2 2 y
x = ay 2 + by + c
Es una parbola
horizontal
x = 4 y2
Figura 3.2
Regin
78
Geraldine Cisneros
D
Observe que la regin D
es una regin tipo 2, por
lo cual el rea se obtiene
empleando una sola
integral doble
de la
.
forma
dxdy
Valor de x a
la entrada de D
Valor de x a
la salida de D
x = 4 y2
x = y2 2 y
Figura 3.3
Regin
2
1
4 y 2
y
dxdy = 4 2 y 2 + 2 y dy = 9
2 y
1
A = dxdy = 9
D
y = 1 1+ x
Para la segunda curva:
x = 4 y2
entonces:
y = 4 x
79
Geraldine Cisneros
Valor de y a
la salida de D2
x=0
y = 1+ 1+ x
y = 4 x
x=3
Valor de y a
la salida de D 3
y = 4 x
D1
D2
Valor de y a
la entrada de D1
D3
y = 1 1+ x
Valor de y a
la entrada de D2
Valor de y a
la entrada de D3
y = 4 x
y = 1 1+ x
Figura 3.4
Entonces D = D1 D2 D3 , donde:
{( x, y ) 1 x 0 1 1 + x y 1 + 1 + x }
D = {( x, y ) 0 x 3 1 1 + x y 4 x }
D = {( x, y ) 3 x 4 4 x y 4 x }
D1 =
2
D2
A=
Al comparar los dos
clculos de rea de la
regin D del ejemplo 3.1,
resulta
mucho
ms
sencillo
emplear
la
que
integral
dxdy
1+ 1+ x
1 1 1+ x
D3
dydx +
A = 2 1 + xdx +
A=
3
0
4 x
1 1+ x
dydx +
4 x
4 x
dydx
4 x 1 + 1 + x dx + 2 4 xdx
4 19 4
+ + =9
3 3 3
A = dydx = 9
D
80
Geraldine Cisneros
EJEMPLO 3.2
dxdy
dydx .
C2
C1
C3
D
Figura 3.5
Regin
Solucin:
Las ecuaciones de las curvas que limitan a la regin D son:
Las ecuaciones de las
curvas en funcin de la
variable y son:
y 20
C1 : x =
16
20 y
C2 : x =
2
y
C1 : x =
2
C1 : y = 16 x + 20
C2 : y = 2 x + 20
C3 : y = 4 x 2
a) Para el clculo del rea de la regin D por medio de la integral
doble
81
Geraldine Cisneros
La regin D no es una
regin tipo 2, sin
embargo se puede dividir
en tres regiones: D1, D2
y D3., que s lo son. Por
esta razn, para resolver
la
integral
doble
dxdy
se
debe
emplear la propiedad
aditiva respecto a la
regin de integracin.
x=
Valor de x a
la salida de D3
y 20
16
D3
x=
20 y
2
y = 16
Valor de x a
la entrada de D2
y 20
x=
16
y=4
Valor de x a
la entrada de D1
x=
Valor de x a
la salida de D2
D2
x=
y
2
Valor de x a
la salida de D1
y
x=
2
D1
y
2
Figura 3.6
Regin
D2
donde:
D1 = ( x, y )
D2 = ( x, y )
D3 = ( x, y )
A=
A=
4
0
4
0
y
2
y
y
y
x
2
2
y
y 20
x
16
2
y 20
20 y
x
16
2
dxdy +
16
4
y
2
y 20
16
0 y 4
4 y 16
16 y 20
dxdy +
20
16
20 y
2
y 20
16
dxdy
16
20 45
y y 20
9y
ydy +
dy + dy
4 2
16
16
4 16
A=
16 157 9
+
+ = 36
3
6
2
D3
82
Geraldine Cisneros
Valor de y a
la salida de D2
y = 2 x + 20
Valor de y a
la salida de D1
y = 16 x + 20
La regin D puede
dividirse en dos regiones
tipo 1, identificadas
como: D1 y D2 ; es decir:
D = D1 D2
D2
Valor de y a
la entrada de D2
D1
y = 4x2
Valor de y a
la entrada de D1
y = 4x2
x=0
Figura 3.7
Regin D del ejemplo 3.2 como dos regiones tipo 1
D2
{( x, y )
= {( x, y )
}
y 2 x + 20}
D1 =
1 x 0 4 x 2 y 16 x + 20
D2
0 x 2 4 x2
A=
0
1
16 x + 20
4x
dydx +
A = (16 x + 20 4 x 2 ) dx +
0
2
0
2 x + 20
4 x2
( 2 x + 20 4 x ) dx = 323 + 763 = 36
A = dydx = 36
D
dydx
2
83
Geraldine Cisneros
EJEMPLO 3.3
La regin D planteada en
el ejemplo 3.3 recibe el
nombre
de
corona
circular, y su rea es:
A = R2 r 2
x2 + y2 = 4
donde
R: Radio externo
r: radio interno
D
x 2 + y 2 = 16
Figura 3.8
Regin
A = 4 A1
La regin denotada como D1, se muestra en la figura 3.9.
84
Geraldine Cisneros
y = 16 x 2
x=2
Valor de y a
la salida de D1.B
D1.A
y = 16 x 2
D1 = D1.A D1.B
D1.B
Valor de y a
la entrada de D1.A
y = 4 x2
Valor de y a
la entrada de D1.B
y=0
Figura 3.9
Regin D1 del ejemplo 3.3
Luego: A1 =
D1. A
dydx +
D1. B
dydx , donde:
{( x, y ) 0 x 2
D = {( x, y ) 2 x 4
D1.A =
1.B
A1 =
A1 =
2
0
2
0
16 x 2
4 x2
4 x 2 y 16 x 2
0 y 16 x 2
dydx +
4
2
16 x 2
0
16 x 2 4 x 2 dx +
4
2
A1 = 2 3 + + 2 3 +
3
3
A = dydx = 12
D
dydx
16 x 2 dx
= 3
85
Geraldine Cisneros
f ( x, y ) dA
D
y g:
( x, y ) D .
Sea
el
volumen
del
slido
acotado
V = g ( x, y ) f ( x, y ) dA
D
EJEMPLO 3.4
(III.5)
z = 20 x 2 y 2
dobles.
Solucin:
En la figura 3.10 se muestra el slido S de este ejemplo, donde la
superficie superior es z = 20 x 2 y 2 y la superficie inferior viene
dada por la ecuacin z = 2 x 2 + y 2 .
86
Geraldine Cisneros
z = 20 x2 y 2
Es una semiesfera (parte
superior).
Valor de z a
la entrada de S
z = 2 x2 + y 2
z = 2 x2 + y2
Es un cono .
Figura 3.10
Slido S del ejemplo 3.4
V = 20 x 2 y 2 2 x 2 + y 2 dA
D
z = 2 x 2 + y 2
z = 20 x 2 y 2
2 x 2 + y 2 = 20 x 2 y 2
4 ( x 2 + y 2 ) = 20 x 2 y 2
x2 + y 2 = 4
Entonces:
D=
{( x, y )
x2 + y 2 4
87
Geraldine Cisneros
y = 4 x2
Valor de y a
la entrada de D
y = 4 x2
Figura 3.11
Regin
Es decir, D =
En el siguiente captulo,
se
mostrar
como
resolver una integral de
este tipo, empleando un
cambio
de
variable
apropiado.
{( x, y )
2 x 2 4 x2 y 4 x2
2
2
4 x2
4 x2
20 x 2 y 2 2 x 2 + y 2 dydx
88
Geraldine Cisneros
EJEMPLO 3.5
x2 + y 2 = 1
Valor de z a
la entrada de S
z =1
Figura 3.12
Slido S del ejemplo 3.5
89
Geraldine Cisneros
D
En este ejemplo, la
regin D es de tipo 1 y
tambin tipo 2, pero se
trabaja como una regin
tipo 2.
Valor de x a
la entrada de D
Valor de x a
la salida de D
x = 1 y2
x = 1 y2
Figura 3.13
Regin D del ejemplo 3.5
{( x, y )
1 y2 x 1 y2
1 y 1
1 y 2
1 1 y
[3 + xy ] dxdy = 1 6
1
1 y 2 dy = 3
V = [3 + xy ] dA = 3
D
EJEMPLO 3.6
Solucin:
En la figura 3.14 se observa el slido S, acotado superiormente por
z = 1 + x 3 y + xy 3
e inferiormente por
90
Geraldine Cisneros
Valor de z a
la entrada de S
z=0
Figura 3.14
Slido S del ejemplo 3.6
y = x3 x
En la figura 3.15, se
observa que la regin D
del ejemplo 3.6 es una
regin de tipo 1.
Valor de y a
la entrada de D
y = x2 + x
Figura 3.15
Regin
91
Geraldine Cisneros
{( x, y )
1 x 0 x 2 + x y x3 x
0
1
x3 x
x2 + x
1 + x3 y + xy 3 dydx
13
0 x
7 x9
517
x11 +
x8 4 x 7 2 x 6 + x3 x 2 2 x dx =
V =
1
4
1260
4
517
V = 1 + x3 y + xy 3 dA =
D
1260
cada punto ( x, y ) D .
Adicionalmente:
( x, y ) = 0 ( x, y ) D
La
densidad
tiene
unidades de masa por
rea unitaria.
Para esta aplicacin,
considere que la funcin
densidad es continua
en la regin D .
Figura 3.16
Regin D no homognea
92
Geraldine Cisneros
(III.6)
m ( xi* , y j * ) Aij
(III.7)
i =1 j =1
(III.8)
i =1 j =1
i =1 j =1
(III.9)
El clculo de masa de
una regin D , tambin
puede emplearse para
calcular
la
carga
elctrica, Q, distribuida
sobre una regin D .
Q = ( x, y ) dA
D
Donde es la funcin
densidad de carga.
m = ( x, y ) dA
D
(III.10)
93
Geraldine Cisneros
EJEMPLO 3.7
Solucin:
Recuerde que la densidad se calcula como m = ( x, y ) dA , por
D
Valor de x a
la entrada de D
Valor de x a
la salida de D
x = 2 y2 2
x = y2 1
Figura 3.17
Regin
{( x, y )
2 y2 2 x y2 1 1 y 1
Por lo tanto:
m=
y 2 1
2
1 2 y 2
dxdy = (1 y 2 ) dy =
m=
y 2 1
1 2 y 2 2
dxdy =
4
3
4
3
94
Geraldine Cisneros
EJEMPLO 3.8
y=
3 2
x 6 x + 4 y y = 2 x 2 , cuya densidad vara de acuerdo a la
2
funcin ( x, y ) = 1 + 2 x .
Solucin:
Segn la definicin del
valor absoluto
x 2 si x 2 0
x2 =
2 x si x 2 < 0
m = (1 + 2 x ) dA
D
entonces
2 x 4 si
y=
4 2 x si
x2
x<2
y = 2x 4
y = 2 x + 4
La regin D debe
dividirse en dos regiones
tipo 1, tal que:
D = D1 D2
D
y=
3 2
x 6x + 4
2
Figura 3.18
Regin
Entonces:
m = (1 + 2 x ) dA = (1 + 2 x ) dA +
D
D1
D2
(1 + 2 x ) dA
Donde
3 2
D1 = ( x, y ) 0 x 2
x 6 x + 4 y 2 x + 4
2
3 2
D2 = ( x, y ) 2 x 4
x 6 x + 4 y 2 x 4
2
95
Geraldine Cisneros
Valor de y a
la salida de D1
y = 4 2x
x=2
D2
D1
Valor de y a
la entrada de D2
3
y = x2 6 x + 4
2
Valor de y a
la entrada de D1
y=
3 2
x 6x + 4
2
Figura 3.19
Entonces:
m=
2
0
4 2 x
3 2
x 6 x + 4
2
2 x4
(1 + 2 x ) dydx + 2 3 x 6 x+ 4 (1 + 2 x ) dydx
2
2
4
13
29 2
m = 3x 3 + x 2 + 4 x dx + 8 3x3 +
x 8 x dx
0
2
2
2
m=
40 80
+
= 40
3
3
m = (1 + 2 x ) dA = 40
D
96
Geraldine Cisneros
una
lmina
placa
plana
D,
dividida
en
M x es una medida de la
tendencia a girar en torno
al eje x, anlogamente,
M y es una medida de la
tendencia
a
girar
alrededor del eje y.
Figura 3.20
Regin general
D no homognea
M x ij = y j * ( xi* , y j * ) Aij
(III.11)
M x y j * ( xi* , y j * ) Aij
i =1 j =1
(III.13)
97
Geraldine Cisneros
(III.14)
i =1 j =1
i =1 j =1
(III.15)
i =1 j =1
(III.16)
, la cual es continua
M x = y ( x, y ) dA
D
(III.17)
M y = x ( x, y ) dA
D
(III.18)
98
Geraldine Cisneros
EJEMPLO 3.9
Solucin:
Los momentos estticos se calculan de la siguiente manera:
M x = y ( x, y ) dA y M y = x ( x, y ) dA .
D
Entonces:
Mx =
Y se encuentra acotada
por las curvas x = y 2 1
y x = 2 y2 2 .
La densidad es :
( x, y ) = 1
2
2
( x, y ) 2 y 2 x y 1
D=
1 y 1
My =
y 2 1
2
1 2 y 2
ydxdy = y (1 y 2 ) dy = 0
1
1
3 3
8
xdxdy = y 4 + 3 y 2 dy =
2
1
5
2 2
y 2 1
1 2 y
Por lo tanto, los momentos estticos para una lmina con la forma
de la regin D del ejemplo 3.7 son:
M x = ydA = 0
D
M y = xdA =
D
EJEMPLO 3.10
8
5
Solucin:
Los momentos estticos se calculan como: M x = y ( x, y ) dA y
D
M y = x ( x, y ) dA .
D
La densidad:
( x, y ) = 1 + 2 x
Donde D = D1 D2
( x, y ) 0 x 2
D1 =
3 2
x 6 x + 4 y 2 x + 4
( x, y ) 2 x 4
D2 =
3 2
x 6 x + 4 y 2 x 4
Mx =
2
0
42 x
3 2
x 6 x + 4
2
y (1 + 2 x ) dydx +
4
2
2 x4
3 2
x 6 x+ 4
2
y (1 + 2 x ) dydx
2
9
135 4
M x = x5 +
x 35 x 3 + 10 x 2 + 16 x dx +
0
8
4
4
9
135 4
+ x5 +
x 35 x3 + 10 x 2 + 16 x dx
2
8
4
99
Geraldine Cisneros
2
0
4 2 x
3 2
x 6 x + 4
2
x (1 + 2 x ) dydx +
My =
4
2
2 x4
3 2
x 6 x + 4
2
x (1 + 2 x ) dydx
64
3
1424
M y = x (1 + 2 x ) dA =
D
15
M x = y (1 + 2 x ) dA =
D
coordenadas
(x,y) D ,
x=
y=
El centro de gravedad
recibe el nombre de
centroide cuando la
densidad es constante.
My
m
Mx
m
(III.19)
(III.20)
100
Geraldine Cisneros
, la cual es continua
x=
1
x ( x, y ) dA
m D
(III.21)
y=
1
y ( x, y ) dA
m D
(III.22)
( x, y ) dA .
D
EJEMPLO 3.11
x = 2 y2 2 .
Su densidad es :
( x, y ) = 1
Y
adicionalmente
obtuvo:
se
4
m = 2 dxdy =
1 2 y 2
3
1
Solucin:
y 2 1
M x = ydA = 0
D
M y = xdA =
D
8
5
8
6
x=
=5 =
4
5
m
3
My
y=
Mx 0
= =0
4
m
3
101
Geraldine Cisneros
Entonces:
6
P ( x, y ) = ,0
5
En la siguiente figura se observa el centro de masa o de gravedad
de la placa D descrita en el ejemplo 3.7
6
,0
5
Figura 3.21
Centro de masa de la regin
EJEMPLO 3.12
m = (1 + 2 x ) dA = 40
D
64
D
3
1424
M y = x (1 + 2 x ) dA =
D
15
M x = y (1 + 2 x ) dA =
Solucin:
Sustituyendo el valor de la masa y los momentos estticos en las
ecuaciones que permiten calcular las coordenadas del centro de
masa, se tiene:
1424
178
x=
= 15 =
40
75
m
My
102
Geraldine Cisneros
178 8
P ( x, y ) =
,
75 15
En la figura 3.22 se aprecia la regin D y su centro de masa:
178 8
,
75 15
Figura 3.22
Centro de masa de la regin
103
Geraldine Cisneros
i =1 j =1
El momento de inercia
alrededor del origen
tambin es conocido
como momento polar de
inercia.
I0 = I x + I y
i =1 j =1
(III.24)
2
2
I 0 = Lim ( xi* ) + ( y j * ) ( xi* , y j * ) Aij = ( x 2 + y 2 ) ( x, y ) dA
D
P 0
i =1 j =1
(III.25)
MOMENTOS DE INERCIA DE FIGURAS PLANAS
Sea D una regin del plano xy, tal que su densidad viene
dada por la funcin
, la cual es continua
(III.26)
I y = x 2 ( x, y ) dA
(III.27)
I 0 = ( x 2 + y 2 ) ( x, y ) dA
D
(III.28)
104
Geraldine Cisneros
EJEMPLO 3.13
La grfica de la regin D
del ejemplo 3.7 se
muestra a continuacin:
Solucin:
Los momentos de inercia alrededor de los ejes coordenados se
calculan
de
la
siguiente
manera:
I x = y 2 ( x, y ) dA ,
D
I y = x 2 ( x, y ) dA y I 0 = ( x 2 + y 2 ) ( x, y ) dA .
D
Ix =
Cuya densidad es :
( x, y ) = 1
2
2
( x, y ) 2 y 2 x y 1
D=
1 y 1
Iy =
I0 =
y 2 1
1 2 y 2 2
y 2 dxdy = y 2 (1 y 2 ) dy =
1
4
15
1 7
7 6
32
2
4
2
x
dxdy
=
1 3 3 y + 7 y 2 y dy = 15
1 2 y 2 2
y 2 1
y 2 1
2
1 2 y 2
(x
1 7
7
12
+ y 2 ) dxdy = y 6 + 6 y 4 6 y 2 dy =
1 3
3
5
I0 = I x + I y =
4 32 36 12
+
=
=
15 15 15 5
4
15
32
I y = x 2 dA =
D
15
I x = y 2 dA =
D
I 0 = ( x 2 + y 2 ) dA =
D
12
5
105
Geraldine Cisneros
EJEMPLO 3.14
La grfica de la regin D
del ejemplo 3.8 se
observa a continuacin:
Ix =
2
0
4 2 x
3 2
x 6 x + 4
2
Ix =
(1 + 2 x )
4
2
2 x4
3 2
x 6 x + 4
2
y 2 (1 + 2 x ) dydx
3
3 2
( 4 2 x ) x 6 x + 4 dx +
2
(1 + 2 x )
( 2 x 4 )
x2 6x + 4
2
dx
Ix =
Donde D = D1 D2
( x, y ) 0 x 2
D1 =
3 2
x 6 x + 4 y 2 x + 4
( x, y ) 2 x 4
D2 =
3 2
x 6 x + 4 y 2 x 4
y 2 (1 + 2 x ) dydx +
2
0
42 x
3 2
x 6 x + 4
2
x 2 (1 + 2 x ) dydx +
4
2
2 x4
3 2
x 6 x+ 4
2
x 2 (1 + 2 x ) dydx
2
4
13
29 4
I y = 3x5 + x 4 + 4 x3 dx + 3 x5 +
x 8 x3 8 x 2 dx
0
2
2
2
Iy =
2
0
42 x
3 2
x 6 x + 4
2
(x
+ y 2 ) (1 + 2 x ) dydx +
O tambin como: I 0 = I x + I y =
4
2
2 x4
3 2
x 6 x+ 4
2
(x
+ y 2 ) (1 + 2 x ) dydx
576
D
7
3856
I y = x 2 (1 + 2 x ) dA =
D
15
35632
I 0 = ( x 2 + y 2 ) (1 + 2 x ) dA =
D
105
I x = y 2 (1 + 2 x ) dA =
106
Geraldine Cisneros
f ( x, y,z ) dV ,
B
como el lmite de
i =1 j =1 k =1
(III.29)
(III.30)
107
Geraldine Cisneros
EJEMPLO 3.15
Solucin:
Para calcular el volumen del slido B, se emplea la integral triple
Valor de z a
la salida de B
z = 5 x2 y 2
y=x
B
y = 2 x
Valor de z a
la entrada de B
z =1
Figura 3.23
Slido
5 x 2 y 2
D 1
dzdA
108
Geraldine Cisneros
Valor de y a
la entrada de D
y=x
Figura 3.24
Proyeccin del slido B del ejemplo 3.15 en el plano xy
{( x, y )
0 x 1
x y 2 x
Luego:
V =
1
0
2 x
5 x 2 y 2
dzdydx =
1
0
(4 x
2 x
y 2 )dydx
1 16
8
8
V = + x 3 4 x 2 4 x dx =
0
3
3 3
V = dV =
B
EJEMPLO 3.16
8
3
y = 4 , y = x2 , z = 0 y z = 4 y .
Solucin:
El clculo de volumen del slido B, se realiza por medio de la
integral triple
109
Geraldine Cisneros
Valor de z a
la salida de B
z = 4 y
Valor de z a
la entrada de B
z=0
Figura 3.25
Slido B del ejemplo 3.16
4 y
0
dzdA
Valor de y a
la entrada de D
y = x2
Figura 3.26
Proyeccin del slido B del ejemplo 3.16 en el plano xy
110
Geraldine Cisneros
{( x, y )
2 x 2
x2 y 4
Luego:
V =
2
2
4 y
x2
dzdydx =
V = dV =
B
EJEMPLO 3.17
( 4 y ) dydx = 2 8 4 x 2 +
2 x
x4
256
dx =
2
15
256
15
La regin tridimensional
comprendida entre las
dos esferas concntricas
es simtrica respecto al
origen, razn por la cual,
dicha regin se divide en
8 partes correspondientes
a cada cuadrante.
x2 + y 2 + z 2 = 1
Figura 3.27
Slido
111
Geraldine Cisneros
Valor de z a
la salida de B1
z = 16 x 2 y 2
Valor de z a
la salida de B1
z = 16 x2 y 2
B1
Valor de z a
la entrada de B1
Valor de z a
la entrada de B1
z=0
z = 1 x2 y 2
Figura 3.28
Slido B1 del ejemplo 3.17
Entonces:
V = dV = 8 dV
B
B1
D2
Figura 3.29
Proyeccin del slido B1
sobre el plano xy
V = 8 dV = 8
dzdA +
B1
D2
D1 0
16 x 2 y 2
2
1 x y
dzdA
112
Geraldine Cisneros
La regin bidimensional
D1 se divide en dos
regiones tipo 1; es decir:
Valor de y a
la salida de D1.2
D1 = D1.1 D1.2
y = 16 x2
D1.1
D1.2
Valor de y a
la entrada de D1.1
y = 1 x
Valor de y a
la entrada de D1.2
y=0
Figura 3.30
Regin D1 del ejemplo 3.17
{( x, y ) 0 x 1
D = {( x, y ) 1 x 4
D1.1 =
1.2
1 x 2 y 16 x 2
0 y 16 x 2
Valor de y a
la salida de D2
y = 1 x2
La regin bidimensional
D2 es una regin tipo 1.
D2
Valor de y a
la entrada de D2
Figura 3.31
y=0
113
Geraldine Cisneros
D2 =
{( x, y )
0 x 1 0 y 1 x2
V = 8 ( 32.98672287 )
V =8
1
0
+8
16 x 2
1
0
1 x 2
1 x 2
0
16 x 2 y 2
0
dzdydx + 8
16 x 2 y 2
1 x 2 y 2
16 x 2
0
16 x 2 y 2
0
dzdydx +
dzdydx
(III.31)
114
Geraldine Cisneros
EJEMPLO 3.18
Valor de z a
la entrada de B
z=0
Figura 3.32
m =
1 y
0
2 zdzdA
115
Geraldine Cisneros
La grfica de la ecuacin:
4 x2
2
y=
x2 + 4 y2 = 4
Es una elipse horizontal.
D
La regin bidimensional
D del ejemplo 3.18 es
una regin tipo 1 y
tambin una regin tipo
2.
Valor de y a
la entrada de D
y=
4 x2
2
Figura 3.33
Regin D del ejemplo 3.18
D = ( x, y )
2 x 2
4 x2
y
2
4 x2
2
m=
2
2
4 x 2
2
4 x 2
1 y
0
2 zdzdydx =
1
4 x2
m=
1+
2 3
2
2
2
4 x2
1
2
m = 2 zdV =
B
EJEMPLO 3.19
4 x2
2
4 x2
(1 y )
dydx
dx = 5
2
5
2
( x, y , z ) = x + y + z + 1 .
116
Geraldine Cisneros
integracin parcial.
Valor de z a
la salida de B
z = 8 4 x2 4 y2
Valor de z a
la entrada de B
z = 4 x2 + 4 y 2
Figura 3.34
Slido
8 4 x 2 4 y 2
4 x2 +4 y 2
( x + y + z + 1) dzdA
2
2
z = 8 4x 4 y
117
Geraldine Cisneros
obtiene la ecuacin x 2 + y 2 = 1 .
Valor de y a
la salida de D
x + y =1
2
y = 1 x2
Es una circunferencia.
Valor de y a
la entrada de D
y = 1 x2
Figura 3.35
Regin D del ejemplo 3.19
D=
{( x, y )
1 x 1 1 x2 y 1 x2
1 x 2
1 1 x 2
1 x 2
1 1 x 2
( 40 8 x
8 4 x 2 4 y 2
4 x2 + 4 y 2
( x + y + z + 1) dzdydx
40 x 2 + 8 x 8 xy 2 8 x 2 y + 8 y 40 y 2 8 y 3 ) dydx
1 160
32
160 2
32
m=
1 x2 + x 1 x2
x 1 x 2 x 3 1 x 2 dx = 20
1
3
3
3
3
m = ( x + y + z + 1) dV = 20
B
118
Geraldine Cisneros
(III.32)
M yz = x ( x, y, z ) dV
(III.33)
M xz = y ( x, y, z ) dV
(III.34)
EJEMPLO 3.20
4 x2
4 x2
B = ( x, y,z ) 2 x 2
y
2
2
0 z 1 y
119
Geraldine Cisneros
2
2
4 x2
2
4 x2
1 y
z ( 2 z ) dzdydx =
2
2
4 x2
2
4 x2
2 (1 y )
dydx
3
3
3
2 2
1
7
M xy = 4 x 2 + ( 4 x 2 ) 2 dx =
2 3
6
3
2
2
4 x2
2
4 x2
1 y
x ( 2 z ) dzdydx =
2
2
4 x2
2
4 x2
x (1 y ) dydx
2
3
2
1
M yz = x 4 x 2 + x ( 4 x 2 ) 2 dx = 0
2
12
2
2
4 x2
2
4 x2
1 y
0
y ( 2 z ) dzdydx =
2
2
4 x2
2
4 x2
y (1 y ) dydx
2
3
2 1
M xz = ( 4 x 2 ) 2 dx =
2
6
M xy = z ( 2 z ) dV =
B
7
3
M yz = x ( 2 z ) dV = 0
B
M xz = y ( 2 z ) dV =
B
EJEMPLO 3.21
B=
{( x, y,z ) 1 x 1
1 x2 y 1 x2 4 x2 + 4 y 2 z 8 4x2 4 y 2
120
Geraldine Cisneros
M xy =
M xy =
M xy =
1
1
1 x 2
1 x 2
8 4 x 2 4 y 2
4 x2 + 4 y2
z ( x + y + z + 1) dzdydx
32 1 1 x2
( 4 x 4 8x 2 + 8 x2 y 2 + 3x 8 y 2 + 3 y + 4 y 4 + 19 )( x 2 + y 2 1) dx
3 1 1 x2
1
1
1
1 x 2
1 x 2
M yz = 8
M yz =
1
1
8 4 x 2 4 y 2
4 x2 + 4 y2
1 x 2
1 1 x 2
x ( x + y + z + 1) dzdydx
x ( x + 5 + y ) ( x 2 + y 2 1) dx
32
160 3 32 4 2
160
x + x2
x x =
1 x2
3
3
3 3
3
1 x 2
1 1 x 2
M xz = 8
1
1
8 4 x 2 4 y 2
4 x2 + 4 y 2
1 x 2
1 x 2
y ( x + y + z + 1) dzdydx
y ( x + 5 + y ) ( x 2 + y 2 1) dx
32 1 x 2
2
M xz =
1 2 x2 + x4 ) =
(
1
15
3
1
272
B
3
2
M yz = x ( x + y + z + 1) dV =
B
3
2
M xz = y ( x + y + z + 1) dV =
B
3
M xy = z ( x + y + z + 1) dV =
121
Geraldine Cisneros
M yz
x=
y=
z=
M xz
m
M xy
m
(III.35)
(III.36)
(III.37)
Entonces:
CENTRO DE MASA DE UN SLIDO DEL ESPACIO
Sea B un recinto del espacio, tal que su densidad viene dada
por la funcin :
x=
1
x ( x, y, z ) dV
m B
(III.38)
y=
1
y ( x, y, z ) dV
m B
(III.39)
z=
1
z ( x, y, z ) dV
m B
(III.40)
m = ( x, y, z ) dV .
B
122
Geraldine Cisneros
EJEMPLO 3.22
5
2
M xy = z ( 2 z ) dV =
B
sin embargo,
M yz = x ( 2 z ) dV = 0
M xz = y ( 2 z ) dV =
x=
y=
M yz
m
0
=0
5
2
M xz
2
=
=
5
5
m
2
7
14
= 3 =
z=
5
15
m
2
M xy
Entonces:
2 14
P ( x , y ,z ) = 0 , ,
5 15
2 14
0, ,
5 15
Figura 3.36
Centro de gravedad del slido
123
Geraldine Cisneros
EJEMPLO 3.23
m = ( x + y + z + 1) dV = 20
272
3
2
M yz =
3
2
M xz =
3
M xy =
III.36 y III.37:
M yz
x=
3= 1
20 30
2
M xz
3= 1
=
20 30
m
y=
z=
M xy
m
272
3 = 68
20
15
Entonces:
1 1 68
P ( x , y ,z ) = , ,
30 30 15
En la siguiente figura, se aprecia el centro de masa del slido B.
1 1 68
, ,
30 30 15
Figura 3.37
Centro de gravedad del slido
124
Geraldine Cisneros
(III.41)
I y = ( x 2 + z 2 ) ( x, y, z ) dV
(III.42)
I z = ( x 2 + y 2 ) ( x, y, z ) dV
(III.43)
EJEMPLO 3.24
(III.44)
4 x2
4 x2
B = ( x, y,z ) 2 x 2
y
2
2
0 z 1 y
125
Geraldine Cisneros
Ix =
2
2
4 x2
2
4 x2
4 x2
2
4 x2
(y
1 y
+ z 2 ) ( 2 z ) dzdydx
4
2
1
2
2 (1 y ) + y (1 y ) dydx
5
3
2 1
3
1
27
I x = 4 x2 +
4 x 2 ) 2 + ( 4 x 2 ) 2 dx =
(
2 2
160
3
8
Respecto al eje y:
Iy =
Iy =
2
2
2
2
4 x2
2
4 x2
4 x2
2
4 x2
(x
1 y
+ z 2 ) ( 2 z ) dzdydx
4
2
1
2
2 (1 y ) + x (1 y ) dydx
4 x2 5 2
(
) + ( 3 + x 2 ) 4 x2 3 2 + 4 x 2 1 + x 2 dx = 119
Iy =
(
)
2
160
12
24
2
Respecto al eje z:
Iz =
2
2
4 x2
2
4 x2
(x
1 y
+y
) ( 2 z ) dzdydx =
2
4 x2
2
4 x2
(y
+ x 2 ) (1 y ) dydx
2
5
3
2 1
1
37
I z = ( 4 x 2 ) 2 + (1 + x 2 )( 4 x 2 ) 2 + x 2 4 x 2 dx =
2 80
12
12
I0 =
2
2
2
2
4 x2
2
4 x2
4 x2
2
4 x2
(x
1 y
+ y 2 + z 2 ) ( 2 z ) dzdydx
4
2
1
2
2
2 (1 y ) + ( y + x ) (1 y ) dydx
126
Geraldine Cisneros
3
4 + x2 )
(
)
(
4 x2
137
2
2
2
2
4 x ) +
I0 =
+
+ x 4 x dx =
(
2
160
12
2
24
27
8
119
I y = ( x 2 + z 2 ) ( 2 z ) dV =
B
24
37
I z = ( x 2 + y 2 + z 2 ) ( 2 z ) dV =
B
12
137
I 0 = ( x 2 + y 2 + z 2 ) ( 2 z ) dV =
B
24
I x = ( y 2 + z 2 ) ( 2 z ) dV =
B
EJEMPLO 3.25
{( x, y,z ) 1 x 1
1 x2 y 1 x2 4 x2 + 4 y 2 z 8 4x2 4 y 2
1 x 2
1 1 x 2
8 4 x 2 4 y 2
4 x2 + 4 y 2
(y
+ z 2 ) ( x + y + z + 1) dzdydx
8 1 1 x2 2
( x + y 2 1) 16 ( x5 + y5 ) + 16 x 4 (13 + y ) +
3 1 1 x2
+ 32 x3 ( y 2 1) + 32 x 2 (13 y 2 + y 3 y 13) + 16 x ( 3 + y 4 ) +
Ix =
127
Geraldine Cisneros
1 x2
(143968 + 22240 x 251584 x 2 30656 x3 +
1 105
+ 160864 x 4 + 12512 x5 53248 x 6 4096 x 7 ) dx
1
I x = 462
Los momentos de inercia respecto a los ejes coordenados y al
origen para el slido del ejemplo 3.19, se muestran a continuacin:
I x = ( y 2 + z 2 ) ( x + y + z + 1) dV = 462
B
I y = ( x 2 + z 2 ) ( x + y + z + 1) dV = 462
B
20
3
1396
I 0 = ( x 2 + y 2 + z 2 ) ( x + y + z + 1) dV =
B
3
I z = ( x 2 + y 2 + z 2 ) ( x + y + z + 1) dV =
B
128
Geraldine Cisneros
4.1 INTRODUCCIN
Recuerde que emplear un
cambio de variable de un
integral definida implica
que el cambio afecta: el
intervalo de integracin,
el integrando y la
diferencial.
La expresin:
g ([ c,d ]) [ a,b ]
Significa
que
las
imgenes de la funcin g
son un subconjunto de
[ a,b] .
f ( x )dx =
x = g (t )
t=d
(IV.1)
CV
dx = g ( t ) dt
t=c
f g ( t ) g ( t ) dt
CLI
x = g (c) = a
x = g (d ) = b
129
Geraldine Cisneros
(IV.2)
Donde:
En otras palabras, la
funcin T transforma
todo punto ( u,v ) D
T1 ( u,v ) = x
(IV.3)
T2 ( u,v ) = y
(IV.4)
en un punto ( x, y ) D .
(IV.5)
T ( u,v ) =
=
T2 ( u,v ) y
(IV.6)
130
Geraldine Cisneros
f ( x, y ) dA
D
f (T ( u,v ) ) = f ( u,v )
(IV.7)
Figura 4.1
Transformacin geomtrica de la regin D en la regin D
D = T ( D ) D = T 1 ( D )
por lo cual
inyectiva.
debe ser
La expresin:
f (T ( u,v ) )
suele escribirse:
tambin
( u ,v ) D \ 2 en ( x, y ) D \ 2 ,
( )
f ( x, y ) dA = f (T ( u,v ) ) ( u,v ) dudv
x, y
(IV.8)
131
Geraldine Cisneros
Al determinante
jacobiano:
del
se obtiene como:
( x, y )
( u,v )
tambin se
jacobiano.
le
( x, y )
se conoce como determinante del jacobiano y
( u,v )
x
u
( x, y )
= det
( u,v )
y
llama
x
v
(IV.9)
Sin
embargo,
en
algunas
xv
yv
ocasiones,
(IV.10)
se
desconoce
la
( x, y )
se obtiene mediante la propiedad:
( u,v )
( x, y ) ( u,v )
=1
( u,v ) ( x, y )
(IV.11)
u x
= det
( x, y )
vx
(IV.12)
En donde:
( u,v )
uy
v y
132
Geraldine Cisneros
(IV.13)
dudv
( x, y )
dA
representa el rea de la
regin D.
dA =
( x, y )
dudv
( u,v )
(IV.14)
D = ( u,v ) u0 u u0 + u v0 v v0 + v
(IV.15)
Figura 4.2
Una regin
D en el plano uv
A ( u0 ,v0 ) ,
D ( u0 + u,v0 + v ) .
B ( u0 + u,v0 ) ,
C ( u0 ,v0 + v )
133
Geraldine Cisneros
(IV.16)
Figura 4.3
Regin
u
u i =
0
0
v j =
v
Por otra parte,
x
u
x
y
u = u , u
u
u
u
x
v x
y
v = v , v
v
v
x
u
T ( ui ) =
y
u
x
x
u
v u
= u
y 0
y
v
u
(IV.17)
x
u
T ( v j ) =
y
u
x
x
v 0
v
= v
y v
y
v
v
(IV.18)
134
Geraldine Cisneros
El
rea
de
un
paralelogramo
cuyos
lados estn definidos por
los vectores:
( a,b ) y ( c,d )
Se obtiene como el valor
absoluto
del
determinante:
a b a c
=
c d b d
x
u u
det
y
u u
x
x
v
u
v
= det
y
y
v
u
v
x
x
u
v
u v = det
y
y
u
v
x
v
u v
y
v
(IV.19)
Recuerde que u y v
son longitudes, por lo
tanto:
u = u
v = v
x
u u
det
y
u u
x
v
( x, y )
v ( x, y )
=
uv =
uv (IV.20)
( u,v )
y ( u,v )
v
v
( x, y )
uv , entonces el rea de
( u,v )
AT ( D)
( x, y )
uv
( u,v )
(IV.21)
AT ( D) =
( x, y )
dudv
( u,v )
(IV.22)
135
Geraldine Cisneros
dA =
( x, y )
dudv
( u,v )
Figura 4.4
Transformacin
EJEMPLO 4.1
T en una regin D
1
dA , empleando un cambio de
1 + xy
En este ejemplo,
transformacin
T ( u , v ) = ( x, y )
la
(1, 2 )
2 , 2
(1,1)
no
1
y=
x
T 1 ( x , y ) = ( u , v )
Figura 4.5
Regin D del ejemplo 4.1
y=
2
x
2, 2
y=x
136
Geraldine Cisneros
, xy = ( u, v )
x
y= x
y
=1 u =1
x
Es decir:
y
y = 2x = 2 u = 2
x
xy = 1 v = 1
T 1 ( x, y ) = ( u , v )
xy = 2 v = 2
D = ( u,v ) 1 u 2 1 v 2
Por
medio
de
T 1
v=2
la
tranformacin T , la
nueva
regin
de
integracin D es una
regin rectangular.
Valor de u a
la entrada de D
u =1
Valor de u a
la salida de D
u=2
v =1
Figura 4.6
Transformacin de la regin
( x, y )
, para lo
( u,v )
y
x = u
v
xy
y
2
( u,v )
= det x
( x, y )
1
x = y y = 2 y = 2u
x x
x
137
Geraldine Cisneros
I =
I =
EJEMPLO 4.2
ln 2
1
2
dv = ln ( 3) ln ( 2 ) ln ( 2 )
2 (1 + v )
2
1
1
2
dA = ln ( 3) ln ( 2 ) ln ( 2 )
1 + xy
2
yx
cos
dA , empleando un cambio
y+x
C2
C1
C4
Figura 4.7
Regin D del ejemplo 4.2
C3
138
Geraldine Cisneros
y+x=2
C3 : y = 0
C4 : y = 1 x
Con el cambio propuesto
se obtiene la regin D
y + x =1 v =1
y+x=2 v=2
y + x =1
yx
A partir de la funcin integrando f ( x, y ) = cos
, se propone
y+x
u = x
y =0
v=x
y x u
y + x = v
y = 0 u = v
u = y
x =0
v = y
Entonces:
x =0u =v
D = ( u,v )
v u v 1 v 2
D por medio de T 1 .
Por
medio
de
T 1
la
tranformacin T , la
nueva
regin
de
integracin D es una
regin tipo 2.
v=2
Valor de u a
la salida de D
u=v
Valor de u a
la entrada de D
u = v
v =1
Figura 4.8
Transformacin de la regin
139
Geraldine Cisneros
Recuerde que:
y x u
T 1 ( x, y ) =
=
y + x v
( x, y )
, se tiene que:
( u,v )
1 1
( u,v )
= 1 1 = 2
= det
( x, y )
1 1
Empleando la ecuacin IV.12 se tiene que:
2 v
2
yx
3
u 1
I = cos
dA = cos
dudv = sen (1) vdv = sen (1)
D
1
v
1
2
v 2
y+x
yx
3
cos
dA = sen (1)
2
y+x
f ( x, y ) dA ,
D
x 2 + y 2 = r2 2
x 2 + y 2 = r12
Figura 4.9
Una regin general D
140
Geraldine Cisneros
D=
{( x, y )
r12 x 2 + y 2 r2 2
tg (1 ) x y tg ( 2 ) x
(IV.23)
(IV.24)
D =
{( r, )
r1 r r2
1 2
(IV.25)
0 < 2
xy .
Figura 4.10
Transformacin de la regin
D en la regin D a travs de T ( r , ) = ( x, y )
( )
f ( x, y ) dA = f ( r cos ,rsen ) ( r, ) drd
x, y
(IV.26)
141
Geraldine Cisneros
Recuerde que:
r cos x
=
T ( r, ) =
rsen y
Y que la identidad
fundamental es:
x
r
( x, y )
= det
( r, )
y
x
cos
= det
y
sen
rsen
= r cos 2 + rsen 2
r cos
( x, y )
= r ( cos 2 + sen 2 ) = r
( r, )
cos 2 + sen 2 = 1
(IV.27)
por
D =
{( r, )
r1 r r2
1 2 ,
donde
0 2 1 < 2 , entonces:
(IV.28)
Figura 4.11
Una regin ms general D
142
Geraldine Cisneros
D =
{( r, )
r1 ( ) r r2 ( ) 1 2
(IV.29)
r2 ( )
(IV.30)
r1
D =
{( r, )
(IV.31)
( ) f ( r cos ,rsen ) rd dr
(IV.32)
1 ( r ) 2 ( r )
r1 r r2
En estos casos:
EJEMPLO 4.3
f ( x, y ) dA =
r2
r1
2 ( r )
1
4 y2
4 y 2
4 y2
4 y 2
dxdy = 2
D=
{( x, y )
4 y2 x 4 y2
0 y2
f T ( r, ) = 1
143
Geraldine Cisneros
r=2
Valor de r a
la entrada de D
Valor de a
la salida de D
r =0
Valor de a
la entrada de D
=0
Figura 4.12
Valores de
{( r, )
0 r 2 0
4 y2
4 y
dxdy =
2
0
rdrd = 2d = 2
rdrd = 2
144
Geraldine Cisneros
EJEMPLO 4.4
A = dydx = 12
D
{( x, y )
4 x 2 + y 2 16
r=4
D
Valor de r a
la entrada de D
= 2
r=2
=0
Figura 4.14
Valores de
{( r, )
2 r 4 0 2
A=
2
0
4
2
rdrd =
UC. Facultad de Ingeniera. Departamento de Matemtica.
2
0
6d = 12
rdrd = 12
145
Geraldine Cisneros
EJEMPLO 4.5
coordenadas polares.
Solucin:
En el ejemplo 3.4 del
captulo 3, y se obtuvo
que:
V = 19,77678464
V = 20 x 2 y 2 2 x 2 + y 2 dA
D
donde D =
{( x, y )
x2 + y 2 4
Figura 4.16
Slido S del ejemplo 4.5
Valor de r a
la salida de D
r=2
= 2
Valor de r a
la entrada de D
=0
r =0
Figura 4.18
Regin D ejemplo 4.5
Figura 4.17
Valores de r y
Como:
r cos x
=
T ( r, ) =
rsen y
Entonces:
2
2
x2 + y 2 = ( r cos ) + ( rsen )
x +y =r
2
Donde D =
{( r, )
0 r 2 0 2
2
2
V = 20 x 2 y 2 2 x 2 + y 2 dA = 20 r 2 2r rdrd
0
0
D
146
Geraldine Cisneros
80
40
160
40
3 5 3 d = 3 5 3 19, 77678464
2
0
Finalmente:
2
20 r 2 2r rdrd = 40 5 160
3
3
3
3
4.3 TRANSFORMACIN GEOMTRICA DE \ \
uvw.
Sea T una funcin definida como T : B \ 3 B \ 3 , tal que:
T ( u,v,w ) = (T1 ( u,v,w ) ,T2 ( u,v,w ) ,T3 ( u,v,w ) )
(IV.33)
Donde:
T1 ( u,v,w ) = x
(IV.34)
T2 ( u,v,w ) = y
(IV.35)
T3 ( u,v,w ) = z
(IV.36)
La funcin T tambin
suele escribirse como:
T1 ( u,v,w) x
T ( u,v,w) = T2 ( u,v,w) = y
T ( u,v,w) z
3
T ( u,v,w ) = ( x, y,z )
(IV.37)
147
Geraldine Cisneros
( u,v,w ) B \ 3 en ( x, y, z ) B \ 3 , mediante la
( u,v,w ) B , entonces:
(
)
f ( x, y,z ) dV = f (T ( u,v,w) ) ( u,v,w) dudvdw
x, y,z
El jacobiano
( x, y,z )
se obtiene como:
( u,v,w )
x
u
( x, y,z )
y
= det
( u,v,w )
u
z
El jacobiano tambin se
denota como:
xu xv xw
( x,y,z )
= det yu yv yw
( u,v,w)
zu zv zw
(IV.38)
x
v
y
v
z
v
x
w
y
w
z
(IV.39)
de
rectangular
coordenadas esfricas.
a:
coordenadas
cilndricas
148
Geraldine Cisneros
En el APNDICE A, se
presenta un repaso del
sistema de coordenadas
cilndricas.
que
se
desea
f ( x, y,z ) dV , donde
B
calcular
una
integral
triple
siguiente figura.
z = z2 ( x, y )
z = z1 ( x, y )
Figura 4.19
Una regin general
B=
{( x, y,z )
( x, y ) D
z1 ( x, y ) z z2 ( x, y )
(IV.40)
(IV.41)
149
Geraldine Cisneros
{( r, ,z ) r r r
1
1 2 z1 ( r, ) z z2 ( r, ) (IV.42)
T ( r, ,z ) = rsen = y
z z
x
r
( x, y,z )
y
= det
( r, ,z )
r
z
x
cos
z
y
= det sen
z
z
0
rsen
r cos
0
0 = r cos 2 + rsen 2
( x, y,z )
= r ( cos 2 + sen 2 ) = r
( r, ,z )
(IV.43)
f : \3 \
Sea
en una
B =
{( r, ,z ) r r r
1
1 2 z1 ( r, ) z z2 ( r, ) ,
(IV.44)
150
Geraldine Cisneros
EJEMPLO 4.6
B=
{( x, y,z ) x
+ y 2 + z 2 4 , x 2 + y 2 1, x 0 , y 0 , z 0
Solucin:
En el ejemplo 2.5 del
captulo 2, y se obtuvo
que:
9
B xyzdV = 8
Valor de z a
la salida de B
z = 4 r2
Cambiando la ecuacin
de
la
esfera
a
x2 + y2 + z2 = 4
coordenadas cilndricas
se tiene:
z = 4 x2 y2 = 4 r2
Valor de z a
la entrada de B
z=0
Figura 4.21
Regin B del ejemplo 4.6
4 x2 y 2
0
xyzdzdA
151
Geraldine Cisneros
r=2
Valor de r a
la entrada de D
r =1
=0
Figura 4.22
Regin
D = ( r, )
1 r 2 0
B ( xyz ) dV = 02 1 0
( xyz ) dV =
2
4 r 2
4r 2
2
0
r3 (4 r2 )
( xyz ) dV =
(r
2
0
9
9
cos sen d =
4
8
9
8
152
Geraldine Cisneros
EJEMPLO 4.7
( xyz ) dV ,
B
( xyz ) dV = 64
B
Solucin:
El slido B , junto con su proyeccin en el plano , xy se muestran
a continuacin, en la figura 4.23
Valor de z a
la salida de B
Valor de z a
la salida de B
z=4
Recuerde
que
las
funciones
del
tipo
deben
z = f ( x, y )
z = 2r
expresarse en funcin de
r y , por lo tanto:
z = 2 ( x 2 + y 2 ) = 2r
Valor de z a
la entrada de B
z=r
z = x2 + y2 = r
Valor de z a
la entrada de B
z=r
Figura 4.23
Slido
B ( xyz ) dV = D
1
4
2
x +y
( xyz ) dzdA
D2
2 x2 + y 2
2
x +y
( xyz ) dzdA
153
Geraldine Cisneros
Valor de r a
la salida de D
D1
r=4
Valor de r a
la entrada de D
r= 8
=0
Figura 4.24
Recuerde que al definir
una regin
D en
coordenadas
polares,
dicha regin se denota
Valor de r a
la salida de D
D2
r= 8
Valor de r a
la entrada de D
=0
r =0
Figura 4.25
Regin D2 del ejemplo 4.7
T ( r, ,z ) = rsen = y
z z
D1 = ( r, )
D2 = ( r, )
8 r 4 0
0r 8
154
Geraldine Cisneros
I = ( xyz ) dV = 2
B
+2
0
I =
2
0
8
0
(r
2r
r 3 (16 r 2 )
I =2
0
(r
4
8
0
r5
cos sen drd
2
256
128
cos sen d + 2
cos sen d
0
3
3
I=
128 64
+
= 64
3
3
Entonces:
2
0
4
r
2r
r
f ( x, y,z ) dV ,
B
donde B es una
155
Geraldine Cisneros
z = z1 ( x, y )
Figura 4.26
Una regin general B
(IV.45)
B =
{( , , )
1 2
1 2
(IV.46)
T ( , , ) = sen sen = y
cos z
y
( x, y,z )
= det
( , , )
z
x
cos sen
y
= det sen sen
z
cos
cos sen
( x, y,z )
= ( 2 sen 2 sen3 ) + ( 2 cos 2 cos 2 sen ) +
( , , )
( 2 sen 2 cos 2 sen ) ( 2 cos 2 sen3 )
156
Geraldine Cisneros
(IV.47)
f : \3 \
en una
B =
{( , , )
1 2
1 2 , donde
sen d d d
(IV.48)
B =
{( , , ) ( , )
1
( , )
1 2
1 2
(IV.49)
2 ( , )
sen d d d
(IV.50)
157
Geraldine Cisneros
EJEMPLO 4.8
El volumen pedido en
este ejercicio se plante
en el ejemplo 3.17 del
captulo 3; sin embargo,
ntese lo fcil que resulta
calcular dicho volumen
en coordenadas esfricas.
B=
{( x, y,z )
1 x 2 + y 2 + z 2 16
Y su volumen es V = dV
B
=4
Figura 4.27
Regin tridimensional
B del ejemplo 4.8
=0
=
Valor de a
la entrada de B
=1
Figura 4.28
Porcin de la regin tridimensional B del ejemplo 4.8
158
Geraldine Cisneros
{( , , ) 1 4
0 2
I = ( xyz ) dV =
4
V = ( R3 r 3 )
3
I =
4
V = ( 64 1) = 84
3
2
0
2 sen d d d
21 sen d d = 42 sen d d = 84
0
EJEMPLO 4.8
4
1
2 sen d d d = 84
La
funcin
T
de
transformacin
a
coordenadas esfricas es:
cos sen x
T ( , , ) = sen sen = y
cos z
por
B=
{( x, y,z ) x
+ y 2 + z 2 4 , x 2 + y 2 1, x 0 , y 0 , z 0
x2 + y 2 + z 2 = 4 = 2
x2 + y 2 = 1
( sen cos )
+ ( sen sen ) = 1
2
1
sen 2
= csc
=2
1
2 = csc 2 =
sen
= csc
sen =
1
2
159
Geraldine Cisneros
Recuerde que:
0 <
Valor de a
la salida de B
=2
Valor de a
la entrada de B
= csc
Figura 4.29
=0
Figura 4.30
Proyeccin del slido B
sobre el plano xy
Regin
B = ( , , )
csc 4 0
Se tiene:
160
Geraldine Cisneros
I = ( xyz ) dV = 2 2
B
I = 2 2
0
I = 2 2
0
2
csc
csc
1
sen3 cos sen cos ( 64 csc 6 ) d d
6
9
9
I = 2 cos sen d =
0 4
8
Finalmente:
2
0
EJEMPLO 4.9
csc
9
8
a) Coordenadas cartesianas.
b) Un cambio de variable adecuado.
Solucin:
El volumen de un slido B se obtiene mediante la integral
dV .
( x 1)
z = x 2 + y 2 es un paraboloide. A continuacin se
161
Geraldine Cisneros
Figura 4.30
Slido B del ejemplo 4.9
z = x2 + y2
Valor de z a
la entrada de B
z=0
Figura 4.31
Primer orden de integracin en coordenadas cartesianas
para el slido B del ejemplo 4.9
162
Geraldine Cisneros
x2 + y 2
dzdA
y = 2 x x2
Valor de y a
la entrada de D
y = 2x x2
Figura 4.32
Regin
D=
{( x, y )
0 x 2 2 x x2 y 2 x x2
Por lo cual:
V =
2
0
2 x x2
2 x x
x2 + y 2
0
dzdydx =
2
0
2 x x2
2 x x2
(x
+ y 2 )dydx
3
22
3
V = ( 2 x x 2 ) 2 + 2 x 2 2 x x 2 dx =
0 3
2
2 x x2
2 x x2
x2 + y 2
0
3
dzdydx =
2
163
Geraldine Cisneros
La transformacin es:
r cos x
T ( r, ,z ) = rsen = y
z z
donde:
r 2 = x2 + y 2
r ( r 2 cos ) = 0 , entonces:
x 2 + y 2 = 2 x r 2 = 2r cos
x 2 + y 2 = 2 x r = 2 cos
Para el paraboloide se tiene:
z = x2 + y 2
z = r2
r = 2 cos
La
grfica
de
r = 2 cos se obtiene
para [ 0, ] .
Cuando
0,
2
se
obtiene
la
semicircunferencia
superior, mientras que
para , , el radio
2
Valor de r a
la entrada de D
r =0
=0
=
Figura 4.33
Regin
D =
{( r, )
0 r 2 cos
V =
2 cos
0
r2
0
rdzdrd =
2 cos
r2
0
2 cos
0
3
r 3 drd = 4 cos 4 d =
0
2
3
rdzdrd =
2