Integrales Multiples-IV

Vista previa del material en texto

INTEGRALES MULTIPLES
EMPEZAREMOS CON LAS INTEGRALES DOBLES
Las integrales dobles se aplican a las funciones de dos variables, z = f(x, y) con la finalidad de calcular el volumen que está debajo de ella. Pero iniciaremos nuestro estudio considerando algunas situaciones particulares que facilitarán su aprendizaje, como z>0 y de dominio rectangular [a, b]x[c, d], como se muestra en la figura siguiente:
Dividimos el intervalo [a, b] en “n” partes iguales y el intervalo [c, d] en m partes iguales, de modo que: y luego tenemos una malla de mxn sub-intervalos de área Aij= 
Tomando el sub-intervalo genérico ij, en donde el punto medio es 
 determinamos el valor de z en dicho punto será f
El área de dicho sub-intervalo es , de modo que dicho rectángulo de dimensiones y f definen un prisma recto de volumen f como se ve en la figura
 Y
d
.
.
yj
.
Y1 
C
 a x1 x2……… …….xi ………………… b x
Δy
 Δx
Como tenemos mxn sub-intervalos, tendremos mxn prismas rectos. El volumen del prisma genérico Vij será:
Ahora vamos a sumar ordenadamente los volúmenes de estos mxn prismas.
Primero sumando los volúmenes
en cada fila j:
J =1
J =2 
 .
.
J = m
Aij= 
Sumando los volúmenes de todas las m filas:
Ahora vamos a sumar los volúmenes en cada columna i:
Si i = 1
Si i = 2 
 . 
 .
Si i = n 
Sumando los volúmenes de todas las n columnas:
El volumen total es el mismo, solo se ha variado el orden de la suma, por tanto:
 =
Llevando al límite cuando n→ y m →
Como y entonces →y→
 
Teorema de Fubini
 
Integrales Iteradas
 
Hay que tener en cuenta que el orden de integración “es de dentro hacia fuera”.
Esto es se integra la integral que está 
dentro del corchete respecto a x manteniendo constante y, y luego el resultado, se integra respecto a y.
EJEMPLO Evaluar 
PROPIEDADES DE LAS INTEGRALES DOBLES
1)
2)
3)
4)
5)
Ejemplo
Calcular el volumen del sólido que se encuentra arriba del cuadrado [0, 2]x[0, 2] y debajo del paraboloide z=16- .
Luego el volumen encerrado V = 48 
 y
 Paraboloide
 z = 16- 
 0 2 y
 2 (2, 2)
 x
 
GENERALIZACIÓN DEL DOMINIO
Dada la función z = f(x, y) > 0, cuyo dominio es una región R simple y cerrada.
Definimos una función F(x, y) cuyo
dominio D =[a,b]x[c, d] que contiene
a R. 
Integrando esta función F(x, y):
R
 y
d
 D
C
 0 a b
En la siguiente figura se muestra el gráfico de F(x, y)
Como F(x, y) = 0 para cuando (x, y) se
encuentra en la región D-R,
el volumen en esa región es cero.
El gráfico de:
Por tanto:
 z
 z=f(x, y)
 c d y
 a
 R D
 b
TIPOS DE REGIONES
Región Tipo I: D1={ (x, y)/ a ≤ x ≤ b y u1(x) ≤ y ≤ u2(x) }
Cuando la región es Tipo I, x varía entre a y b, mientras que y va, lo vemos en el gráfico, que para un valor x є [a, b] y va de u1 a u2
La integral doble para esa región será:
 y y y
 u2(x)
 u2(x) u2(x)
 
 D1
 u1(x) u1(x) u1(x)
 0 a x b 0 a b 0 a b
Región Tipo II: D2 = {(x, y) / c ≤ y ≤ d y h1(y) ≤ x ≤ h2(y) }
Cuando la región es Tipo I, vemos en el gráfico, que para un valor y є [c,d] y x va de h1 a h2.
La integral doble para esa región será:
 y y y
 d d d 
 
 y
 h1(y) h2(y) h1(y) h2(y) h1(y) h2(y)
 D2 D2 D2
 c c c
 0 x 0 x 0 x
 
Teorema de Fubini para regiones generales
Cálculo de A(x):
Tomando el plano x = x corta a la
superficie S dada por z =f(x,y) en
una curva en donde x permanece 
Constante x = x, luego el área A(x):
A(x)= 
Para calcular el volumen consideramos un dx que define con A(x) un dV = A(x)dx, integrando entre a y b tenemos:
VOLUMEN = 
Ahora vamos a verlo desde el otro punto de vista
Cálculo de A(y):
Tomando el plano y = y corta a la
superficie S dada por z =f(x,y) en
una curva en donde y permanece 
constante y = y, luego el área A(y):
A(y) = 
Para calcular el volumen consideramos un dy que define con A(y) un dV = A(y)dy, integrando entre c y d tenemos:
VOLUMEN = 
Ejercicio
Evalue donde D es la región limitada por las 
curvas 
En donde
Hacemos el gráfico, vemos que la 
región es tipo I:
D={(x, y)/ -1 ≤ x ≤ 1 y 2
Luego la integral será: 
D
 y
 2
 1
 -1 0 1 x
Ejemplo 
Evalúe donde D es la región limitada por la recta 
y=x-1 y la parábola = 2x+6. 
18
INTRGRALES DOBLES EN CC. POLARES
Relación entre coordenadas polares (r, θ) y las rectangulares (x, y) de un punto:
Se usan cuando las regiones están limitadas por funciones que limitan las regiones son circulares. Como circunferencias, cardiodes, rosa de cuatro pétalos, trebol, etc.
Dada la función z = f(x, y) > 0, cuyo dominio R es un sector de corona circular, 
como se muestra en la figura:
Dividimos el intervalo de r en 
“n” partes y el intervalo de θ
En “m” partes, de modo que:
De modo que tendríamos una malla de mxn sub-intervalos, si nos detenemos en el sub-intervalo genérico
El punto medio de dicho sub-intervalo será: 
Θ=α
Θ=β
a
b
α
En dondey 
El área de un sector circular viene dada por:
El área del sub-intervalo genérico será:
Determinando el valor de f(x, y) en el punto genérico:
Estamos frente a un prisma cuya sección recta es el área del sub-intervalo genérico y la altura el valor de f(x, y) en el punto genérico, cuyo volumen estará dado por:
Como tenemos mxn sub-intervalos, tendremos mxn prismas, sumando los volúmenes, que se obtienen sumándolas ordenadamente (suma de Riemann):
Llevando al límite cuando m y n tienden a infinito:
Luego: 
EJEMPLO
Halle el volumen del sólido que se encuentra encima de z=0, dentro del cilindro 
z=
Como la región es un círculo
conviene usar C.C. Polares
x= rcosθ e y=rsenθ
La circunferencia:
La región de integración es el círculo, que en C.C. Polares queda definido de la siguiente manera;
 z
 4
 D
 x 2 1 0
 y 
z=
D = {(r, θ)/ - y 0 ≤ r ≤ 2cos θ} 
z = (x, y) = 
z = f(r, θ) = reemplazando:
 y
 r=2cosθ
 r
 θ
 D
 0 1 2 x
EJEMPLO
Determinar el volumen del sólido que está encima del plano z=0, limitado por los planos x=4 ; y = 6-x, y= 0; x=0 y debajo de 
 como se muestra en la figura siguiente.
 La región es tipo I:
			 D={(x, y)/ 0≤x≤4 y 0≤y≤6-x}
			 El volumen estará dado por:
			 
APLICACIONES
Valor medio, es el valor “z = f(x, y)” que multiplicada por el área de la base, sea igual al volumen que se encuentra debajo de la superficie S definida por z = f(x, y).
Área de la región D, se determina de modo indirecto, ya que en realidad lo que se calcula es el volumen del sólido cuyo volumen coincide con el valor del área de la región D.
 Sabemos que: 
	
1
EJEMPLO
Con el ejemplo anterior vamos a:
a) Calcular el área de D:
b) Calcular el valor medio:
Como el volumen esta dado por V = 48 
El valor medio z = = 
Cálculo de la masa de la región D
 SI f(x, y) es la densidad superficial = ρ(x, y) = 
 Entonces dm = ρ(x, y) dA integrando:
 En el ejemplo anterior si ρ(x, y) = 2x
Cálculo de los momentos de primer orden y centro de masa
Sabemos que la masa viene dada por:
Y sus momentos alrededor de los 
dos ejes coordenados son: 
El centro de masa se localiza en el punto
Donde: 
x
y
0 dm
D
MOMENTOS DE INERCIA
 El Momento de Inercia es una medida de
 la materia a resistirse a cambios en el 
 movimiento de rotación.
 Si consideramos un diferencial de masa
 en un punto cualquiera de la lámina, el
 
 dIy = y
 dIx = Integrando:
 
 El momento de Inercia respecto al origen: Io = Ix+Iy
x x
Y
y
EJEMPLO
Una lámina de densidad ρ(x,y)=xy está limitada por el eje de las x, la recta x = 8 y la curva y=x 2/3 . Calcular: la masa, el centro de masa y momentos de inercia.
Solución: En la figura, se aprecia la región correspondiente a la lámina
y=x 2/3 
D
Calculando los momentos de primer orden, se tiene:
Por lo tanto, el Centro de Masa de la lamina es:
 CM=
Para hallar los momentos de inercia aplicamos:
x
y
y=x 2/3 
D
AREA DE UNA SUPERFICIE
Sea S una superficie con ecuación z = f(x, y), donde F tiene derivadas parciales continuas. Para facilitar el análisis consideramos un dominio rectangular D=[a, b]x[c, d]
Dividimos en el intervalo [a, b] en
n partes iguales y el [c, d] en m, de
modo que:
 y 
Tendremos una malla de mxn sub-
intervalos iguales. Tomamos el sub-
Intervalo genérico, ij y definimos el
punto genérico (xi, yj) y por ese 
punto trazamos el plano tangente.
c
d
a
 b
x
Dicho plano queda dividido por 
las caras del prisma correspon-
diente al sub-intervalo genérico
en un paralelogramo, como se
muestra en la figura en azul.
Si hacemos un zoom en esta
zona tendremos:
Vamos aproximar el área ΔSij
con el área del plano tagente ΔTij
a
b
ΔSij
ΔTij
Como ΔTij es igual //a x b// vamos a calcular a y b:
Como a es paralelo al plano y=0 no tiene componente “j” luego a=(Δx, 0, Δx fx(xi, yj)) y b por ser paralelo al plano x=0 no tiene componente “i” luego b = (o, Δy, Δy fy(xi, yj)) 
 
Como tenemos mxn sub-intervalos, tenemos que sumarlos, aplicando la suma de Riemann:
 
Llevando al límite Tij cuando m y n tienden a infinito
Expresión que nos permite calcular el área de la superficie S
EJEMPLOS ILUSTRATIVOS
Hallar el área de la porción del plano 
z = 2 – x – y que está interceptada por
El cilindro 1 en el primer 
cuadrante 
La región de integración D, al ser parte
De un círculo, conviene usar C.CPolares
Luego D={(r, θ)/ 0≤ θ ≤ π/2 y 0 ≤ r ≤ 1}
Como el área viene dada por:
Tenemos que calcular fx y fy:
fx(x, y) = -1 y fy(x, y)= -1
Remplazando en la integral los datos:
EJEMPLO 2
Hallar el área de la superficie de sobre el círculo de radio 1.
Calculamos las derivadas parciales
 y
 
Como R es un círculo, usamos CC. Polares
D ={(r, θ)/ 0 ≤ r ≤ 1 y 0 ≤ θ ≤ 2π}
INTEGRALES TRIPLES
Las integrales triples se aplican para funciones de tres variables como w = f(x, y, z), esto nos presenta una gran dificultad debido a que no sabemos graficar dicha función debido a que se encuentra en un espacio tetra-dimensional, pero por inducción, aunque no sepamos como es, intentaremos imaginarlo. Las integrales triples se logran entender mejor en las aplicaciones.
Tomaremos como dominio un prisma recto, pero para facilitar el análisis nuestro dominio D =[0, a]x[0, b]x[0, c], y dividiremos 
El [0, a] en n partes iguales; el [0, b] en m partes iguales y el [0, c] en ñ partes iguales, de modo que tendremos mxnxñ sub-intervalos iguales, cuyas dimensiones son:
Tomaremos el elemento genérico ijk
De coordenadas 
El prisma genérico tiene como 
punto medio a 
y su volumen 
Definimos el valor de f en ese
punto w = f
c
a
b
Asociando este valor con su correspondiente volumen, tendríamos:
Como tenemos mxnxñ sub-intervalos y que cada uno lleva asociado el producto 
Sumándolos ordenadamente tendríamos, según la suma de Riemann:
Llevando al límite cuando n, m y ñ → :
 
TEOREMA DE FUBINI
Si “f” es continua sobre su dominio D =[0, a]x[0, b]x[0, c], entonces:
El orden de integración es de adentro hacia afuera, en la integral de la derecha, empezamos integrando respecto a z, manteniendo constantes a y y a z. luego integramos respecto a y manteniendo constante a x y por último integramos respecto a x. Hay otras posibles órdenes de integración:
Integrales Triples para Regiones más Generales
Tomemos una función f(x, y, z)>0 continua en una región general E, que es un sólido simple y cerrada situada en el primer octante.
Definimos una función F(x, y, z) de
dominio D=[0, a]x[0, b]x[0, c], que
contiene a E, del modo siguiente:
Como F satisface la definición: la 
Integral será:
 
0 E
D
c
b
a
Expresando la integral de modo explícito:
Por tanto:
TIPOS DE REGIONES DEINTEGRACION
Región Tipo I: E1={(x, y, z)/ (x, y) є R y u1(x, y) ≤ z ≤ u2(x, y)}
La integral será:
Para un punto (x, y)єR vemos que z va de u1 a u2
 z
 z = u2(x, y)
 z
 z=u1(x, y)
 0 y
 R
 x
(x, y)
Esta integral presenta dos alternativas dependiendo de como es la región R.
a) Si R es tipo I: R1 ={ (x, y)/ a ≤ x ≤ b y h1(x) ≤ y ≤ h2(x) }
b) Si R es tipo II: R2 ={ (x, y)/ c ≤ y ≤ d y g1(y) ≤ x ≤ g2(y) }
2) Región Tipo II: E2={(x, y, z)/ (y, z) є R y u1(y, z) ≤ z ≤ u2(y, z)}
La integral será:
Para un punto (y, z)єR vemos que x va de u1 a u2
Si R es tipo I: 
 R1={ (y, z)/ c ≤ y ≤ d y h1(y) ≤ z ≤ h2(y) }
 
b) Si R es tipo II: R2={ (y, z)/ e ≤ z ≤ f y r1(z) ≤ y ≤ r2(z) }
 z
 R
 E u1(y, z)
 u2(y, z) y
 x
3) Región Tipo III: E3={(x, y, z)/ (x, z) є R y u1(x, z) ≤ y ≤ u2(x, z)}
Para un punto (y, z)єR vemos que x va de u1 a u2
Si R es tipo I: 
R1={ (y, z)/ a ≤ x ≤ b y h1(x) ≤ z ≤ h2(x) }
b) Si R es tipo II: R2={ (y, z)/ e ≤ z ≤ f y r1(z) ≤ x ≤ r2(z) }
 z
 u1(x, z) u2(x, z) 
 R 
 0 y
x 
EJEMPLO Evaluar la integral en donde Q es:
 La región E será: E={(x, y, z) / (x, y)ϵ R y o ≤ z ≤ }
 
 
 z
 
 
 x=0
 0 E
 y=0
 z=0
 R
Como la región R es un sector de corona circular, usaremos CC. Polares: x =rcos. y=rsen:
Luego R={(r, )/ 1
EJEMPLO
Evaluar la integral en el recinto señalado
Es una región Tipo I: graficando R Y E:
La región es Tipo I: R = { } 
Luego E = 
Con esa información podemos construir los gráficos
R
0 1 2 x
E
Luego:
PROPIEDADES DE LAS INTEGRALES TRIPLES
APLICACIONES DE LAS INTEGRALES TRIPLES
Cálculo de Volúmenes
De modo análogo, el Volumen se calcula haciendo f(x, y, z) = 1 en la integral:
 
EJEMPLO
Calcular el volumen del sólido limitado por las superficies: 
Plano ; Cilindro parabólico y los planos coordenados
Según el gráfico la región E es Tipo I:
E={(x, y, z)/ (x,y)ϵR y o≤ z ≤ }
Donde R=[(x, y) / 0 ≤ x ≤ 4; 0 ≤ y ≤6-x}
Luego la integral :
planteadas. Luego; 
se puede visualizar el sólido que se genera de la intersección de 
Cálculo de la Masa del sólido definido por el Dominio
dm = ρ(x, y, z) dV Integrando:
EJEMPLO: Calcularemos la masa del sólido visto en el ejemplo anterior si ρ(x, y, z) = 2x
Remplazando en la ecuación:
Cálculo de Momentos de Primer Orden y el centro de masa
Se calculan respecto a los planos coor-
denados:
dMxz= dm = y 
dMxy= dm = z 
dMyz= dm = x 
Integrando:
x
z
y
dm
El Centro de Masa se calcula:
La masa viene dada por:
Los Momentos de Primer Orden:
Luego las coordenadas del centro de masa es una media ponderada:
62
EJEMPLO:
Calcular el centro de masa del cubo unidad cuyos vértices son (0,0,0); (0,1,0); (1,1,0), (0,1,0); (0,0,1); (1,0,1); (1,1,1); (0,1,1). Si la densidad en el punto (x,y,z) es proporcional al cuadrado de su distancia al origen
Cálculo de Myz:
Por simetría y como la densidad es constante 
por tanto
 
 CM = 
EJEMPLO
Encuentre el centro de masa de un sólido de densidad constante que está limitado por el cilindro parabólico x = y los planos x = z, z = 0 y x = 1.
Densidad 
E={(x, y,z)/ (x, y)єR; 0 ≤ z ≤ x}
R={(x, y)/ -1≤ y ≤ 1; ≤ x ≤ 1}
 z
 -1
 x = 
 0
 E
 1 
 1 R y 
 x 
Haciendo lo mismo y por simetría y = 0 
5. MOMENTOS DE INERCIA RESPECTO A LOS EJES
Momento de inercia respecto a x:
Momento de inercia respecto a y:
Momento de inercia respecto a y:
 y
 d3 
 d2 d2 
 d1 z
 x
MOMENTOS DE INERCIA RESPECTO A LOS PLANOS COORDENADOS
Momento respecto al plano z=0
Momento respecto al plano y=0
Momento respecto al plano x=0
Estos momentos se relacionan de esta manera:
CAMBIO DE VARIABLES EN INTEGRALES MÚLTIPLES
Un cambio de variables está dado por una Transformación Lineal T del plano “uv” al plano “xy”
Se denota así: T(u, v) = (x, y)
Donde x e y están relacionados con u y v por las ecuaciones:
x = g(u, v) e y = h(u, v)
Y g y h tienen derivadas parciales continuas de primer orden.
Una transformación T es una función cuyo dominio (o imagen), son subconjuntos de Si T(u, v) = (x, y) entonces el punto (u, v) se llama imagen de (x, y). 
Si no hay dos puntos que tengan la misma imagen , T se llama biunívoca. En el siguiente gráfico se muestra el efecto de una transformación T de una región S del plano uv. T transforma a S en una región R del plano xy llamado imagen de S, formado por las imágenes de todos los puntos de S.
 v y
 T
 
 
 0 u 0 x
EJEMPLO
Una transformación está definida por las ecuaciones:
	x= y = 2uv
S 
(u, v)
 S
 (x, y)
Encuentre la imagen de S={(u, v)/ 0
x= y = 2uv
Teniendo la frontera de S vamos a definir 
La frontera de R;
1) Camino 1: v=0 0 u 1
 Luego: x= y = 0
 Como 0 u 1 → 0 1 → 0 x 1 
 Luego y = 0 si 0 x 1 
2) Camino 2: u = 1 si 0 v 1
 Luego : x= ; y = 2v → v = y/2
 Reemplazando: x = 1 - → = -2(x-1)
 Como 0 v 1 → 0 2v 2 entonces 0 y 2
 v
 1 C3 (1, 1)
 S
 C4 C2
 0C1 1 u
3) Camino 3: v = 1 si 0 u 1
 x= y = 2u → u = y/2
 Reemplazando: x = → = 2(x+1)
 Como 0 u 1 → 0 2u 2 entonces 0 y 2 
Camino 4: u = 0 si 0 v 1
 x= y = 0
 Como 0 v 1 → 0 1 → -1 0 
 Luego: y = 0 si -1 0 
 
 Graficando lo calculado:
 y
 2
 = 2(x+1) = -2(x-1)
 R
 -1 0 1 x
 
CAMBIO DE VARIABLES EN UNA INTEGRAL DOBLE
Se tiene una función z = f(x, y) , en donde x e y son funciones de u y de v dadas por: x = g(u, v) e y = h(u, v) y se quiere expresar 
la integral en función de u y de v. Sabiendo que el dominio S de u y v, está dado por:
S = { (u, v) / a
De los datos se tiene la función vectorial:
 r(u, v) = g(u, v) i +h(u, v) j = (x, y)
 
Dividimos el intervalo [a, b] en “n” partes iguales y el intervalo [c, d] en m partes iguales, de modo que: y luego tenemos una malla de mxn sub-intervalos de área Aij= 
 v 
 d v=d
 
 vj
 u =a u=b
 c
 v=c
 0 a ui b u
S
CAMBIO DE VARIABLES EN UNA INTEGRAL DOBLE
Dada una función z = f(x, y) en donde x e y son funciones de dos parámetros u y v, definidas con las Transformación x=g(u, v) e y = h(u, v). Estando definidas los puntos (u, v) en una región S limitada por a
Dividimos el intervalo [a, b] en “n” partes iguales y el intervalo [c, d] en m partes iguales, de modo que: y y definimos las curvas r(u,v)
 y 
 r(u, d)
 v(a,v) 
 v(b,v) 
 
 r(u, c)
0 x
R
Definimos la frontera de R de modo que la imagen de u=a tiene como imagen la curva r(a, v), u = b a r(b, v), v=c a r(u, c) y v = d a r(u, d). Luego para cada u=ui y v = vj definimos las curvas reticulares r(ui, v) e r(u, vj)
Tomamos el sub.intervalo genérico ij que define en S un parche, trazamos las secantes a y b
En donde a = r(ui+Δu, vj) – r(ui, vj) y b=r(ui, vj+ Δv)- r(ui, vj) 
Una primera aproximación sería que el área del parche sea igual //axb//
(ui, vj)
u=ui
v=vj
b
Sij
Rij
b
r(ui+Δu. vj)
T
)
r(ui, vj)
 a
r(ui,vj+Δv) 
(ui, vj)
u=ui
v=vj
Por otro lado tenemos que las derivadas parciales:
Haciendo una aproximación tendríamos que a y b serían:
Luego: 
 Tij=
Por otro lado sabemos que r(u, v) = (x, y)
Luego como:
Luego:
Por tanto:
EJEMPLO
Evaluar la siguiente integral donde R es la
región trapezoidal de vértices (1,0), (2,0), (0,-2) y (0,-1)
Graficando la región:
Haciendo el cambio de variables:
Determinamos la frontera de S:
C1
C2
C3
C4
Graficando tenemos S:
S={(u, v)/ 1 ≤ v ≤ 2; -v ≤ u ≤ v }
Jacobiano: 
Cálculo de las derivadas parciales:
S={(u, v)/ 1 ≤ v ≤ 2; -v ≤ u ≤ v }
Haciendo el cambio de variables:
Como S = {(u, v)/ 1 ≤ v ≤ 2; -v ≤ u ≤ v }
Jacobiano: J = 2
CAMBIO DE VARIABLES PARA INTEGRALES TRIPLES
Si tenemos la función de tres variables w = f(x, y, z) definida en un dominio E, de modo que x = g(u, v, w); y = h(u, v, w) y z = r(u, v, w), entonces el Jacobiano o factor de escala viene dada por:
Luego el cambio de variable de una integral triple está dado:
INTEGRALES TRIPLES EN COORDENADAS CILINDRICAS
Algunos sólidos son difíciles de representar en coordenadas rectangulares y es preciso usar otro tipo de coordenadas-
Ecuaciones de conversión de rectangulares a cilíndricas
 
 En donde r = R es un cilindro
 θ = θ1 es un plano
 z = z1 es u plano
 
 R 1 
En este tipo de coordenadas el sólido más sencillo es un bloque cilíndrico, vienen a ser las coordenadas polares en el espacio.
Cálculo de la integral triple en coordenadas cilíndricas:
Tenemos la transformación dada por las siguientes ecuaciones:
El cambio de variables en una integral triple viene dada por:
En nuestro caso u = r; v = θ; z = z, cálculo del Jacobiano:
Como:
Remplazando:
Luego la integral triple en coordenadas cilíndricas viene dada por:
Para ilustrar su aplicación vamos a
calcular el volumen del sólido que
se muestra en la figura en donde E
será:
E = { (r, θ, z)/ r1≤ r ≤ r2; θ1≤ θ ≤ θ2 ;
z1 ≤ z ≤ z2 }
Luego como f(x, y, z) = 1 cuando se quiere calcular el volumen, entonces f(r, θ,, z) = 1, por tanto:
EJEMPLO
Calcular el volumen del sólido que
se encuentra dentro de la esfera y
 del cilindro dados por:
Graficando tendríamos :
 
 2
 0 1 2
Vamos a definir la región E, como el sólido se encuentra dentro de un cilindro conviene usar las coordenadas cilíndricas, ya que en en coordenadas cartesianas , la región asume una forma muy compleja, como mostraremos a continuación:
Como en polares es r =2senθ
Si vemos el sólido desde eje z tenemos:
R={(r, θ) / 0 ≤ θ ≤ π ; 0 ≤ r ≤ 2sen θ}
Mientras que z varía:
 
 r=2senθ
 R
 0 1 2 y
 
 x
En coordenadas polares será:
E={(r, θ) / 0 ≤ θ ≤ π ; 0 ≤ r ≤ 2sen θ ; }
Luego el volumen estará dado por:
INTEGRALES TRIPLES EN COORDENADAS ESFERICAS
Ecuaciones de conversión de rectangulares a cilíndricas
 φ = φ1 Semicono 
 θ = θ1 Semiplano 
 r = r1 Cilindro 
Cálculo de la integral triple en coordenadas esféricas:
Tenemos la transformación dada por las siguientes ecuaciones:
El cambio de variables en una integral triple viene dada por:
En nuestro caso u = r; v = θ; z = φ, cálculo del Jacobiano:
Como:
Determinamos las derivadas parciales:
Remplazando
Remplazando valores tenemos:
Sentido de variación de las variables
El radio varía: El ángulo φ varía: El ángulo θ varía:
 0 ≤ r ≤ R 0 ≤ φ ≤ π 0 ≤ θ ≤ 2π
 
 z
 0 R y
 x
 z
 
 0 φ1 y
 
 x
 
 z0 y
 θ
 x 
EJEMPLO
Hallar el volumen de la región sólida Q acotada por abajo por la hoja superior del cono y por la esfera de radio 3 y centrada en el origen
Graficando:
La esfera será: 
Interceptando con 
Tendríamos: 
En donde z = 
En el triángulo OPR rectángulo
Isósceles: φ = 
Luego la región E será:
E={(r, θ, φ)/ (r, θ)єR; 0 ≤ φ ≤ } y R={(r, θ)/ 0 ≤ r ≤ 3; 0 ≤ θ ≤ 2π}
 3 Esfera
 
 P R 
 E
Cono
 φ
 
 0 
 
 
Luego el volumen pedido estará dado por la integral que en coordenadas esféricas viene dada por:
EJEMPLO
Plantee la integral triple, en coordenadas cilíndricas, del volumen del sólido mostrado en la figura:
Del gráfico se deduce que el sólido
está limitado por las superficies:
Cilindro: 
Planos: z = 2; z = 0; x = 0; y = 0
De acuerdo a los datos la región
será:
E={(r, θ, z)/ (r, θ)єR; 0 ≤ z ≤ } y R={(r, θ)/ 0 ≤ r ≤ 3; 0 ≤ θ ≤ }
 
 z
 
 x=0
 z= 2
 Plano
 E
 y=0 0 3
 R y
 3 z=0
 Cilindro
 x
y
x
y
x
f
V
j
i
ij
D
D
=
)
,
(
*
*
å
=
D
D
=
n
i
i
i
y
x
y
x
f
V
1
*
1
*
1
)
,
(
å
=
D
D
=
n
i
i
i
y
x
y
x
f
V
1
*
2
*
2
)
,
(
å
=
D
D
=
n
i
m
i
im
y
x
y
x
f
V
1
*
*
)
,
(
)
,
(
*
*
j
i
y
x
)
,
(
*
*
j
i
y
x
f
å
å
=
=
D
D
=
m
j
n
i
j
i
s
y
x
y
x
f
V
1
1
*
*
)
,
(
å
=
D
D
=
m
j
j
j
y
x
y
x
f
V
1
*
*
1
1
)
,
(
å
=
D
D
=
m
j
j
j
y
x
y
x
f
V
1
*
*
2
2
)
,
(
å
=
D
D
=
m
j
j
n
nj
y
x
y
x
f
V
1
*
*
)
,
(
å
å
=
=
D
D
=
n
i
m
j
j
i
s
y
x
y
x
f
V
1
1
*
*
)
,
(
å
å
=
=
D
D
m
j
n
i
j
i
y
x
y
x
f
1
1
*
*
)
,
(
ò
ò
å
å
ò
ò
å
å
=
D
D
=
D
D
=
=
¥
®
¥
®
=
=
¥
®
¥
®
b
a
d
c
n
i
m
j
j
i
m
n
d
c
b
a
m
j
n
i
j
i
m
n
dy
dx
y
x
f
y
x
y
x
f
dy
dx
y
x
f
y
x
y
x
f
)
,
(
)
,
(
lim
)
,
(
)
,
(
lim
1
1
*
*
1
1
*
*
ò
ò
ò
ò
=
d
c
b
a
b
a
d
c
dy
dx
y
x
f
dx
dy
y
x
f
)
,
(
)
,
(
dy
dx
y
x
f
dx
dy
y
x
f
d
c
b
a
b
a
d
c
ò
ò
ò
ò
ú
û
ù
ê
ë
é
=
ú
û
ù
ê
ë
é
)
,
(
)
,
(
dy
dx
y
x
f
d
c
b
a
ò
ò
ú
û
ù
ê
ë
é
)
,
(
[
]
2
27
2
9
9
3
:
2
27
3
2
3
2
3
2
:
2
1
2
2
1
2
1
3
0
3
2
1
3
0
2
3
0
3
3
0
2
2
1
3
0
2
2
3
0
2
1
2
=
ú
û
ù
ê
ë
é
=
=
ú
û
ù
ê
ë
é
=
=
=
ú
û
ù
ê
ë
é
=
ú
û
ù
ê
ë
é
=
ú
û
ù
ê
ë
é
=
=
ò
ò
ò
ò
ò
ò
ò
ò
y
dy
y
dy
y
x
dy
ydx
x
I
x
a
respecto
primero
Integramos
x
dx
x
dx
y
x
dx
ydy
x
I
y
a
respecto
primero
Integramos
ò
ò
3
0
2
1
2
dx
ydy
x
3
2
0
3
2
0
2
0
2
2
0
2
0
2
3
2
0
2
0
2
2
2
0
2
0
48
3
144
3
32
176
3
4
3
88
4
3
88
2
3
16
)
2
16
(
)
,
(
u
y
y
V
dy
y
dy
xy
x
x
V
dxdy
y
x
dxdy
y
x
f
V
=
=
-
=
ú
û
ù
ê
ë
é
-
=
=
ú
û
ù
ê
ë
é
-
=
ú
û
ù
ê
ë
é
-
-
=
=
-
-
=
=
ò
ò
ò
ò
ò
ò
î
í
ì
-
Î
Î
=
R
D
y
x
si
R
y
x
si
y
x
f
y
x
F
)
,
(
0
)
,
(
)
,
(
)
,
(
òò
òò
òò
òò
=
+
=
-
R
R
D
R
D
dA
y
x
f
dA
dA
y
x
f
dA
y
x
F
)
,
(
0
)
,
(
)
,
(
òò
òò
òò
òò
òò
òò
=
=
+
=
-
R
D
R
R
D
R
D
dA
y
x
f
dA
y
x
F
tiene
se
y
x
F
de
definición
la
Por
dA
y
x
f
dA
y
x
F
dA
y
x
F
dA
y
x
F
)
,
(
)
,
(
:
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
ò
ò
d
c
y
h
y
h
dy
dx
y
x
f
)
(
)
(
2
1
)
,
(
òò
+
D
dA
y
x
)
2
(
2
2
2
1
1
2
x
y
x
y
e
x
y
=
-
+
=
=
2
1
x
y
=
-
2
2
x
y
=
[
]
(
)
ò
ò
ò
ò
òò
-
-
+
-
+
-
-
+
+
+
=
+
=
=
+
=
+
1
1
4
3
2
2
3
1
1
1
2
2
1
1
1
2
4
2
)
1
(
)
2
(
)
2
(
2
2
2
2
dx
x
x
x
x
x
dx
y
xy
I
I
dx
dy
y
x
dA
y
x
x
x
x
x
D
[
]
(
)
(
)
15
26
15
24
20
30
5
)
2
(
3
4
0
3
)
2
(
2
2
0
2
5
3
4
3
2
2
3
2
1
4
2
)
1
(
1
1
5
4
3
2
1
1
4
3
2
1
1
4
3
2
2
3
1
1
1
2
2
2
2
=
-
+
=
ú
û
ù
ê
ë
é
-
-
+
+
=
=
ú
û
ù
ê
ë
é
-
-
+
+
=
-
-
+
+
=
-
-
+
+
+
=
+
=
-
-
-
-
+
ò
ò
ò
I
x
x
x
x
x
dx
x
x
x
x
I
dx
x
x
x
x
x
dx
y
xy
I
x
x
òò
D
xydA
q
q
rsen
y
r
x
=
=
cos
x
y
tg
y
x
r
=
+
=
q
2
2
2
(
)
{
}
b
q
a
q
£
£
Ù
£
£
=
b
r
a
r
R
:
,
]
,
[
]
,
[
1
1
|
j
j
i
i
x
r
r
q
q
-
-
)
,
(
*
*
j
i
r
q
)
,
(
*
*
j
i
r
q
(
)
þ
ý
ü
î
í
ì
£
£
£
£
=
-
-
j
j
i
i
j
i
ij
r
r
r
r
R
q
q
q
q
1
1
:
,
(
)
*
*
,
j
i
r
q
)
(
2
1
1
*
i
i
i
r
r
r
+
=
-
)
(
2
1
1
*
j
j
j
q
q
q
+
=
-
2
2
1
r
A
sc
q
=
)
,
cos
(
*
*
*
*
j
i
j
i
sen
r
r
f
q
q
[
]
[
]
[
]
*
1
1
2
1
2
2
1
2
2
1
2
1
2
1
2
1
i
i
j
i
i
i
i
j
ij
i
i
j
i
j
i
j
ij
r
r
r
r
r
r
A
r
r
r
r
A
D
D
=
+
-
D
=
-
D
=
D
-
D
=
-
-
-
-
q
q
q
q
q
[
]
*
*
*
*
*
)
,
cos
(
i
i
j
j
i
j
i
ij
r
r
sen
r
r
f
V
D
D
=
q
q
q
å
å
=
=
D
D
=
n
i
m
j
i
i
j
j
i
j
i
s
r
r
sen
r
r
f
V
1
1
*
*
*
*
*
)
,
cos
(
q
q
q
ò
ò
å
å
==
=
D
D
=
=
=
¥
®
¥
®
b
a
n
i
m
j
i
i
j
j
i
j
i
m
n
dr
d
rsen
r
rf
I
r
r
sen
r
r
f
I
b
a
q
q
q
q
q
q
)
,
cos
(
)
,
cos
(
lim
1
1
*
*
*
*
*
ò
ò
ò
ò
=
=
b
a
b
a
dr
d
r
rf
dr
d
rsen
r
rf
V
b
a
b
a
q
q
q
q
q
)
,
(
)
,
cos
(
x
y
x
2
2
2
=
+
q
q
q
cos
2
0
)
cos
2
(
cos
2
2
=
®
=
-
®
=
r
r
r
r
r
[
]
2
3
2
2
cos
1
4
cos
4
4
)
(
)
,
(
2
2
2
2
2
4
2
2
cos
2
0
4
2
2
cos
2
0
2
p
q
q
q
q
q
q
q
q
p
p
p
p
p
p
q
p
p
q
b
a
ò
ò
ò
ò
ò
ò
ò
-
-
-
-
=
ú
û
ù
ê
ë
é
+
=
=
=
ú
û
ù
ê
ë
é
=
=
=
d
V
d
d
r
V
d
dr
r
r
dr
d
r
rf
V
b
a
4
4
2
x
z
-
=
x
y
-
=
6
4
=
x
4
4
2
x
z
-
=
(
)
3
4
0
4
3
2
3
2
2
48
16
1
2
1
2
24
4
1
2
3
4
24
4
4
6
4
0
4
0
u
x
x
x
x
V
dx
x
x
x
dx
x
x
V
==
ú
û
ù
ê
ë
é
==
=
÷
ø
ö
ç
è
æ
=
÷
÷
ø
ö
ç
ç
è
æ
=
+
-
+
-
-
-
-
-
ò
ò
dx
x
y
dx
dy
x
x
x
-
-
-
ò
ò
ò
ú
û
ù
ê
ë
é
÷
÷
ø
ö
ç
ç
è
æ
=
÷
÷
ø
ö
ç
ç
è
æ
-
6
0
2
2
4
0
4
0
6
0
4
4
4
4
òò
=
D
D
dA
y
x
f
A
z
)
,
(
1
D
D
D
A
dA
V
y
x
f
Haciendo
dA
y
x
f
V
=
=
®
=
=
òò
òò
1
)
,
(
)
,
(
[
]
[
]
2
4
0
2
2
6
0
16
2
6
2
6
6
4
0
4
0
4
0
4
0
6
0
u
x
x
dx
x
x
A
dx
x
dx
y
dx
dy
A
D
x
D
x
=
ú
û
ù
ê
ë
é
-
=
ú
û
ù
ê
ë
é
-
=
-
=
=
=
ò
ò
ò
ò
ò
-
-
òò
òò
=
=
D
D
D
dA
y
x
m
dm
)
,
(
r
[
]
[
]
m
x
u
x
x
dx
x
x
m
dx
xy
dx
xdy
m
x
3
160
3
2
6
)
2
12
2
2
4
0
3
2
2
6
0
4
0
4
0
4
0
6
0
=
ú
û
ù
ê
ë
é
-
=
-
=
=
=
=
ò
ò
ò
ò
-
-
(
)
dA
y
x
M
D
òò
=
,
r
(
)
y
x
,
(
)
(
)
dA
D
y
x
x
y
M
dA
D
y
x
y
x
M
Integrando
dA
y
x
x
xdm
dM
y
dA
y
x
y
ydm
dM
y
x
òò
=
òò
=
=
=
=
=
,
;
,
:
)
,
(
)
,
(
r
r
r
r
(
)
y
x
,
M
x
M
Y
M
y
M
X
=
=
;
dm
(
)
(
)
(
)
(
)
dA
y
x
x
I
y
dA
y
x
y
I
D
D
y
x
,
,
2
2
r
r
òò
òò
=
=
(
)
5
768
20
3
3
7
8
0
3
10
8
0
3
/
2
0
8
0
2
8
0
3
/
2
0
2
1
2
,
=
ú
û
ù
ê
ë
é
=
=
ú
û
ù
ê
ë
é
=
=
=
=
ò
ò
ò
ò
òò
x
dx
x
dx
y
x
M
dydx
xy
dA
y
x
M
x
x
D
r
(
)
(
)
13
12288
8
0
26
3
3
10
3
13
8
0
3
/
2
0
8
0
2
2
8
0
3
/
2
0
2
1
2
,
=
ú
û
ù
ê
ë
é
=
=
ú
û
ù
ê
ë
é
=
=
=
ò
ò
ò
ò
òò
x
dx
x
dx
y
x
M
dydx
xy
x
dA
y
x
x
M
x
x
y
D
y
r
(
)
(
)
[
]
3
1024
12
1
3
8
0
4
8
0
3
/
2
0
8
0
3
8
0
3
/
2
0
3
1
3
,
=
=
=
ú
û
ù
ê
ë
é
=
=
=
ò
ò
ò
ò
òò
x
dx
x
dx
y
x
M
dydx
xy
y
dA
y
x
y
M
x
x
x
D
x
r
M
x
M
Y
M
y
M
X
=
=
;
9
20
;
13
80
5
768
3
1024
;
5
768
13
12288
=
=
=
=
Y
X
Y
X
÷
÷
ø
ö
ç
ç
è
æ
9
20
,
13
80
(
)
(
)
(
)
(
)
dA
y
x
x
I
dA
y
x
y
I
D
D
y
x
,
,
2
2
r
r
òò
òò
=
=
(
)
(
)
6144
7
6144
8
0
0
2
8
0
0
2
3
/
2
3
/
2
=
=
=
=
ò
ò
ò
ò
dydx
xy
x
I
dydx
xy
y
I
x
y
x
x
1
)
,
(
)
,
(
0
0
2
2
+
+
D
D
=
D
D
D
D
=
j
i
y
j
i
x
y
j
j
x
i
i
y
x
f
y
x
f
y
x
f
y
y
f
x
x
k
j
i
axb
å
å
=
=
+
+
D
D
=
n
i
m
j
j
i
y
j
i
x
y
x
f
y
x
f
y
x
T
1
1
2
2
1
)
,
(
)
,
(
dx
dy
y
x
f
y
x
f
S
y
x
f
y
x
f
y
x
S
b
a
d
c
y
x
n
i
m
j
j
i
y
j
i
x
m
n
ò
ò
å
å
+
+
=
+
+
D
D
=
=
=
¥
®
¥
®
1
)
,
(
)
,
(
1
)
,
(
)
,
(
lim
2
2
1
1
2
2
4
3
2
3
2
3
3
1
1
1
2
0
1
0
2
0
2
2
0
1
0
2
0
1
0
p
q
q
q
q
p
p
p
p
=
=
ú
û
ù
ê
ë
é
=
=
=
+
+
=
ò
ò
ò
ò
ò
ò
d
d
r
S
d
rdr
d
rdr
S
dx
dy
y
x
f
y
x
f
S
b
a
d
c
y
x
ò
ò
+
+
=
1
)
,
(
)
,
(
2
2
2
2
1
y
x
z
+
+
=
x
y
x
f
x
2
)
,(
=
y
y
x
f
y
2
)
,
(
=
òò
×
+
+
=
R
dA
y
x
S
1
4
4
2
2
ò
ò
+
=
p
q
2
0
1
0
2
4
1
rdrd
r
S
(
)
(
)
(
)
(
)
6
1
5
5
6
1
5
5
1
5
5
12
1
4
1
12
1
4
1
2
0
2
0
1
0
2
3
2
2
0
1
0
2
p
p
q
q
q
p
p
p
-
=
=
-
=
-
=
=
ú
û
ù
ê
ë
é
+
=
+
=
ò
ò
ò
ò
S
d
S
d
r
rdrd
r
S
3
4
2
3
1
2
min
?
¿
)
,
,
(
min
)
,
(
min
)
(
V
en
io
Do
y
V
en
dV
z
y
x
f
V
en
io
Do
y
V
en
Volumen
dA
y
x
f
V
en
io
Do
y
V
en
Area
dx
x
f
b
a
f
e
d
c
b
a
d
c
b
a
ò
ò
ò
ò
ò
ò
=
=
=
ñ
c
z
m
b
y
n
a
x
=
D
=
D
=
D
;
;
i
i
i
ijk
z
y
x
V
D
D
D
=
D
)
,
,
(
*
*
*
k
j
i
z
y
x
i
x
D
j
y
D
k
z
D
)
,
,
(
k
j
i
z
y
x
)
,
,
(
*
*
*
k
j
i
z
y
x
k
j
i
k
j
i
z
y
x
z
y
x
f
D
D
D
)
,
,
(
*
*
*
å
å
å
=
=
=
D
D
D
n
i
m
j
ñ
k
k
j
i
k
j
i
z
y
x
z
y
x
f
1
1
1
*
*
*
)
,
,
(
ò
ò
ò
òòò
å
å
å
=
=
D
D
D
=
=
=
¥
®
¥
®
¥
®
a
b
c
D
n
i
m
j
ñ
k
k
j
i
k
j
i
ñ
m
n
dz
dy
dx
z
y
x
f
dV
z
y
x
f
z
y
x
z
y
x
f
0
0
0
1
1
1
*
*
*
)
,
,
(
)
,
,
(
)
,
,
(
lim
ò
ò
ò
òòò
=
a
b
c
D
dx
dy
dz
z
y
x
f
dV
z
y
x
f
0
0
0
)
,
,
(
)
,
,
(
î
í
ì
-
Î
Î
=
E
D
z
y
x
si
E
z
y
x
si
z
y
x
f
z
y
x
F
)
,
,
(
0
)
,
,
(
)
,
,
(
)
,
,
(
òòò
D
dV
z
y
x
F
)
,
,
(
òòò
òòò
òòò
-
+
=
E
D
E
D
dV
dV
z
y
x
f
dV
z
y
x
F
0
)
,
,
(
)
,
,
(
òò
ò
òòò
=
R
y
x
u
y
x
u
E
dA
dz
z
y
x
f
dV
z
y
x
f
)
,
(
)
,
(
2
1
)
,
,
(
)
,
,
(
òòò
òòò
=
E
D
dV
z
y
x
f
dV
z
y
x
F
)
,
,
(
)
,
,
(
ò
ò
ò
òò
ò
=
b
a
x
h
x
h
y
x
u
y
x
u
R
y
x
u
y
x
u
dzdydx
z
y
x
f
dA
dz
z
y
x
f
)
(
)
(
)
,
(
)
,
(
)
,
(
)
,
(
2
1
2
1
2
1
)
,
,
(
)
,
,
(
ò
ò
ò
òò
ò
=
d
c
y
g
y
g
y
x
u
y
x
u
R
y
x
u
y
x
u
dzdxdy
z
y
x
f
dA
dz
z
y
x
f
)
(
)
(
)
,
(
)
,
(
)
,
(
)
,
(
2
1
2
1
2
1
)
,
,
(
)
,
,
(
òò
ò
òòò
=
R
z
y
u
z
y
u
E
dA
dx
z
y
x
f
dV
z
y
x
f
)
,
(
)
,
(
2
1
)
,
,
(
)
,
,
(
ò
ò
ò
òò
ò
=
d
c
y
h
y
h
z
y
u
z
y
u
R
z
y
u
z
y
u
dxdzdy
z
y
x
f
dA
dx
z
y
x
f
)
(
)
(
)
,
(
)
,
(
)
,
(
)
,
(
2
1
2
1
2
1
)
,
,
(
)
,
,
(
ò
ò
ò
òò
ò
=
f
e
z
r
z
h
z
y
u
z
y
u
R
z
y
u
z
y
u
dxdydz
z
y
x
f
dA
dx
z
y
x
f
)
(
)
(
)
,
(
)
,
(
)
,
(
)
,
(
2
1
2
1
2
1
)
,
,
(
)
,
,
(
ò
ò
ò
òò
ò
=
b
a
x
h
x
h
z
x
u
z
x
u
R
z
x
u
z
x
u
dzdx
dy
z
y
x
f
dA
dy
z
y
x
f
)
(
)
(
)
,
(
)
,
(
)
,
(
)
,
(
2
1
2
1
2
1
)
,
,
(
)
,
,
(
òò
ò
òòò
=
R
z
x
u
z
x
u
E
dA
dy
z
y
x
f
dV
z
y
x
f
)
,
(
)
,
(
2
1
)
,
,
(
)
,
,
(
ò
ò
ò
òò
ò
=
f
e
z
r
z
r
z
x
u
z
x
u
R
z
x
u
z
x
u
dxdz
dy
z
y
x
f
dA
dy
z
y
x
f
)
(
)
(
)
,
(
)
,
(
)
,
(
)
,
(
2
1
2
1
2
1
)
,
,
(
)
,
,
(
òòò
Q
xyzdV
þ
ý
ü
î
í
ì
³
³
³
³
+
£
+
+
=
0
,
0
,
0
;
1
;
4
/
)
,
,
(
2
2
2
2
2
z
y
x
y
x
z
y
x
z
y
x
Q
1
2
2
³
+
y
x
4
2
2
2
£
+
+
z
y
x
=
ú
û
ù
ê
ë
é
-
-
=
=
ú
û
ù
ê
ë
é
=
=
ú
û
ù
ê
ë
é
=
òò
òò
ò
òòò
òò
-
-
-
-
R
R
y
x
o
Q
y
x
R
dA
xy
y
x
xy
I
dA
zdz
xy
xyzdV
dA
z
xy
I
2
4
2
3
3
4
4
0
2
2
2
2
2
q
q
q
q
q
q
q
q
q
q
q
q
p
p
p
d
sen
r
sen
r
I
drd
sen
r
sen
r
I
rdrd
sen
r
sen
r
sen
r
I
dA
xy
y
x
xy
I
R
ò
ò
ò
ò
ò
òò
÷
÷
ø
ö
ç
ç
è
æ
-
=
÷
÷
ø
ö
ç
ç
è
æ
-
=
ú
û
ù
ê
ë
é
-
-
=
=
ú
û
ù
ê
ë
é
-
-
=
2
0
2
1
6
4
2
0
2
1
5
3
2
0
2
1
3
4
3
4
2
3
3
24
4
2
2
2
1
4
4
2
2
2
1
2
cos
cos
2
2
2
4
4
15
96
)
1
1
(
63
4
)
1
1
(
15
2
1
96
4
cos
63
4
2
cos
15
2
1
24
4
63
2
2
15
2
1
2
0
2
0
=
ú
û
ù
ê
ë
é
-
+
-
-
-
=
=
ú
û
ù
ê
ë
é
+
-
=
÷
ø
ö
ç
è
æ
-
=
ò
I
d
sen
sen
I
p
p
q
q
q
q
q
òòò
Q
dV
{
}
2
2
2
;
4
;
2
0
/
)
,
,
(
y
x
z
y
x
x
y
x
x
z
y
x
Q
+
£
£
+
£
£
£
£
=
x
y
x
x
y
x
4
;
2
0
/
)
,
(
2
£
£
£
£
2
x
y
=
x
y
4
=
}
)
,
(
/
)
,
,
{(
2
2
y
x
z
y
x
y
R
y
x
z
y
x
+
£
£
+
Î
2
2
y
x
z
+
£
z
y
x
£
+
[
]
[
]
[
]
=
-
-
+
-
=
+
-
+
-
-
+
-
=
ú
û
ù
ê
ë
é
-
-
-
-
-
+
-
=
=
ú
û
ù
ê
ë
é
-
-
+
=
=
-
-
+
=
=
ò
ò
ò
ò
ò
ò
ò
ò
ò
+
+
dx
x
x
x
x
V
dx
x
x
x
x
x
x
x
x
V
dx
x
x
x
x
x
x
x
x
x
x
V
y
xy
y
y
x
V
dydx
y
x
y
x
dzdydx
V
x
x
x
x
x
x
y
x
y
x
2
0
2
6
3
4
2
0
4
2
3
2
6
3
4
3
2
0
4
2
2
6
3
2
2
2
0
4
2
3
2
2
0
4
2
2
2
0
4
72
2
158
3
6
1
3
48
6
24
2
128
6
24
6
1
2
16
)
4
(
3
64
)
4
(
2
3
2
2
2
2
2
105
6724
105
20160
3840
2016
66360
6
1
3
576
7
256
632
5
96
6
1
3
72
7
2
4
158
5
3
6
1
2
0
3
7
4
5
=
÷
ø
ö
ç
è
æ
-
-
-
=
ú
û
ù
ê
ë
é
-
-
+
-
=
ú
û
ù
ê
ë
é
-
-
+
-
=
V
V
x
x
x
x
V
òòò
òòò
=
×
Q
Q
dV
z
y
x
f
c
dV
z
y
x
f
c
)
,
,
(
)
,
,
(
.
1
[
]
=
±
òòò
dV
z
y
x
g
z
y
x
f
Q
)
,
,
(
)
,
,
(
.
2
òòò
òòò
±
=
Q
Q
dV
z
y
x
g
dV
z
y
x
f
)
,
,
(
)
,
,
(
òòò
òòò
òòò
±
=
2
1
)
,
,
(
)
,
,
(
)
,
,
(
.
3
Q
Q
Q
dV
z
y
x
f
dV
z
y
x
f
dV
z
y
x
f
2
1
Q
Q
Q
È
=
òòò
òòò
=
®
=
E
E
dV
V
dV
z
y
x
f
I
)
,
,
(
6
=
+
x
y
4
4
2
x
z
-
=
0
=
x
4
4
2
x
z
-
=
0
=
z
0
=
y
0
=
x
4
4
2
x
-
3
4
0
4
2
3
4
0
3
2
4
0
2
4
0
6
0
2
4
0
2
6
0
4
0
6
0
4
4
0
48
16
32
2
64
96
16
2
2
24
4
4
2
3
24
)
6
(
4
4
4
4
4
4
2
u
I
x
x
x
x
dx
x
x
x
I
dx
x
x
dx
yx
y
I
dydx
x
dzdydx
I
x
x
x
x
=
ú
û
ù
ê
ë
é
+
-
-
=
=
ú
û
ù
ê
ë
é
+
-
-
=
÷
÷
ø
ö
ç
ç
è
æ
+
-
-
=
=
-
÷
÷
ø
ö
ç
ç
è
æ
-
=
ú
û
ù
ê
ë
é
-
=
=
ú
û
ù
ê
ë
é
-
=
=
ò
ò
ò
ò
ò
ò
ò
ò
-
-
-
-
=
=
òòò
E
dV
z
y
x
m
)
,
,
(
r
15
1856
10
1024
3
512
192
384
10
3
8
4
3
2
48
2
8
3
48
)
6
(
2
8
2
8
2
8
2
4
0
5
3
4
2
4
0
4
2
3
4
0
3
4
0
6
0
3
4
0
3
6
0
4
0
6
0
4
4
0
2
=
+
-
-
=
÷
÷
ø
ö
ç
ç
è
æ
+
-
-
=
=
÷
÷
ø
ö
ç
ç
è
æ
+
-
-
=
=
-
÷
÷
ø
ö
ç
ç
è
æ
-
=
ú
û
ù
ê
ë
é
÷
÷
ø
ö
ç
ç
è
æ
-
=
=
ú
û
ù
ê
ë
é
-
=
=
ò
ò
ò
ò
ò
ò
ò
ò
-
-
-
-
x
x
x
x
I
dx
x
x
x
x
I
dx
x
x
x
dx
y
x
x
I
dydx
x
x
xdzdydx
m
x
x
x
x
òòò
òòò
òòò
=
=
=
E
yz
E
xy
E
xz
dV
z
y
x
x
M
dV
z
y
x
z
M
dV
z
y
x
y
M
)
,
,
(
)
,
,
(
)
,
,
(
r
r
r
òòò
òòò
òòò
=
=
=
E
yz
E
xy
E
xz
dV
z
y
x
x
M
y
dV
z
y
x
z
M
dV
z
y
x
y
M
)
,
,
(
)
,
,
(
;
)
,
,
(
r
r
r
m
M
z
m
M
y
m
M
x
xy
xz
yz
=
=
=
;
;
m
M
z
m
M
y
m
M
x
xy
xz
yz
=
=
=
;
;
)
(
2
2
2
z
y
x
k
+
+
=
r
òòò
=
Q
dV
z
y
x
masa
)
,
,
(
r
ò
ò
ò
+
+
=
1
0
1
0
1
0
2
2
2
)
(
dzdydx
z
y
x
k
masa
ò
ò
ò
+
+
=
1
0
1
0
1
0
2
2
2
)
(
dzdydx
z
y
x
x
k
M
yz
k
x
x
k
dx
x
k
m
dx
y
y
y
x
k
dydx
y
x
k
m
dydx
z
z
y
z
x
k
dzdydx
z
y
x
km
m
=
ú
û
ù
ê
ë
é
+
=
ú
û
ù
ê
ë
é
+
=
ú
û
ù
ê
ë
é
+
+
=
ú
û
ù
ê
ë
é
+
+
=
=
ú
û
ù
ê
ë
é
+
+
=
+
+
=
ò
ò
ò
ò
ò
ò
ò
ò
ò
1
0
3
1
0
2
1
0
1
0
3
2
1
0
1
0
2
2
1
0
1
0
1
0
3
2
2
1
0
1
0
1
0
2
2
2
3
2
3
3
2
3
3
3
1
3
)
(
12
7
3
4
3
2
3
3
3
)
(
)
(
1
0
2
4
1
0
3
1
0
1
0
3
3
1
0
1
0
2
3
1
0
1
0
1
0
2
2
3
1
0
1
0
1
0
2
2
2
k
x
x
k
dx
x
x
k
M
dx
xy
xy
y
x
k
dydx
x
xy
x
k
M
dzdydx
xz
xy
x
k
dzdydx
z
y
x
kx
M
yz
yz
yz
=
ú
û
ù
ê
ë
é
+
=
ú
û
ù
ê
ë
é
+
=
ú
û
ù
ê
ë
é
+
+
=
ú
û
ù
ê
ë
é
+
+
=
=
+
+
=
+
+
=
ò
ò
ò
ò
ò
ò
ò
ò
ò
ò
12
7
=
®
x
k
k
m
M
x
yz
12
/
7
=
=
12
/
7
=
=
=
z
y
x
(
)
÷
ø
ö
ç
è
æ
=
12
7
,
12
7
,
12
7
,
,
z
y
x
òòò
òòò
=
=
E
E
kdV
m
dV
z
y
x
m
)
,
,
(
r
5
4
5
4
5
2
2
2
5
2
5
2
2
1
2
1
1
5
1
1
5
1
1
4
1
1
1
2
1
1
1
1
1
1
0
2
2
2
k
m
k
k
y
y
k
m
y
y
k
dy
y
k
m
dy
kx
kxdxdy
m
kdzdxdy
kdV
m
y
y
E
y
x
=
=
÷
ø
ö
ç
è
æ
-
=
ú
û
ù
ê
ë
é
-
=
=
ú
û
ù
ê
ë
é
-
=
ú
û
ù
ê
ë
é
-
=
=
ú
û
ù
ê
ë
é
=
=
=
=
=
-
-
-
-
-
-
ò
ò
ò
ò
òòò
ò
ò
ò
[
]
7
5
:
7
4
21
2
14
21
2
3
2
3
1
3
1
1
6
1
1
1
3
1
1
1
2
1
1
1
0
1
1
1
0
2
2
2
2
=
=
=
-
=
ú
û
ù
ê
ë
é
-
=
=
ú
û
ù
ê
ë
é
-
=
ú
û
ù
ê
ë
é
=
=
=
=
=
=
=
ò
ò
ò
ò
ò
ò
òòò
ò
ò
ò
-
-
-
-
-
m
M
x
Luego
k
k
k
M
dy
y
k
dy
x
k
M
dxdy
x
k
dxdy
xz
k
M
xkdzdxdy
xkdV
M
yz
yz
y
yz
y
y
x
yz
E
y
x
yz
14
5
=
=
m
M
z
xy
òòò
+
=
Q
x
dV
z
y
x
z
y
I
)
,
,
(
)
(
2
2
r
òòò
+
=
Q
y
dV
z
y
x
z
x
I
)
,
,
(
)
(
2
2
r
òòò
+
=
Q
z
dV
z
y
x
y
x
I
)
,
,
(
)
(
2
2
r
òòò
=
Q
xy
dV
z
y
x
z
I
)
,
,
(
2
r
òòò
=
Q
xz
dV
z
y
x
y
I
)
,
,
(
2
r
òòò
=
Q
yz
dV
z
y
x
x
I
)
,
,
(
2
r
xy
xz
x
I
I
I
+
=
xy
yz
y
I
I
I
+
=
xz
yz
z
I
I
I
+
=
òò
R
dxdy
y
x
f
)
,
(
v
v
u
r
v
v
u
r
v
u
r
v
v
u
r
v
v
u
r
v
u
r
u
v
u
r
v
u
u
r
v
u
r
u
v
u
r
v
u
u
r
v
u
r
j
i
j
i
j
i
v
j
i
j
i
v
j
i
v
j
i
j
i
j
i
u
j
i
j
i
u
j
i
u
D
-
D
+
»
®
D
-
D
+
=
D
-
D
+
»
®
D
-
D
+
=
D
D
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
lim
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
lim
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
j
i
v
j
i
j
i
j
i
u
j
i
j
i
v
u
r
v
v
u
r
v
v
u
r
b
v
u
r
u
v
u
r
v
u
u
r
a
D
»
-
D
+
=
D
»
-
D
+
=
j
v
y
i
v
x
r
y
j
u
y
i
u
x
r
Entonces
v
u
h
y
e
v
u
g
x
y
y
x
v
u
r
v
u
¶
¶
+
¶
¶
=
¶
¶
+
¶
¶
=
=
=
=
:
)
,
(
)
,
(
)
,
(
)
,
(
[
]
(
)
(
)
v
u
j
i
v
j
i
u
r
x
r
v
u
v
u
r
v
x
v
u
r
u
axb
D
D
=
D
D
=
)
,
(
)
,
(
òò
òò
òò
òò
òò
å
å
å
å
òò
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
=
=
¶
¶
D
D
=
D
=
D
=
¶
¶
D
D
=
=
=
=
®
D
®
D
=
=
®
D
®
D
R
S
S
S
S
n
i
m
j
j
i
j
i
v
u
n
i
m
j
j
i
R
y
x
ij
dudv
v
u
y
x
v
u
f
dudv
v
u
y
x
v
u
h
v
u
g
f
dA
y
x
f
to
Por
dudv
v
u
y
x
v
u
f
dudv
v
u
y
x
v
u
h
v
u
g
f
v
u
y
x
v
u
v
u
h
v
u
g
f
A
y
x
f
dA
y
x
f
A
v
u
y
x
v
u
axb
T
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
))
,
(
),
,
(
(
)
,
(
:
tan
)
,
(
),
(
)
,
(
)
,
(
)
,
(
))
,
(
),
,
(
(
)
,
(
)
,
(
))
,
(
),
,
(
(
lim
)
,
(
lim
)
,
(
)
,
(
)
,
(
1
1
0
0
1
1
0
0
k
v
y
v
x
u
y
u
x
v
y
v
x
u
y
u
x
k
j
i
xr
r
v
u
¶
¶
¶
¶
¶
¶
¶
¶
=
¶
¶
¶
¶
¶
¶
¶
¶
=
0
0
î
í
ì
î
í
ì
£
£
=
=
=
®
=
=
®
=
£
£
=
2
1
0
2
1
2
u
v
u
v
u
y
v
y
x
u
y
x
C
î
í
ì
-
=
+
=
y
x
v
y
x
u
dA
e
R
y
x
y
x
òò
-
+
î
í
ì
£
£
-
=
=
î
í
ì
£
£
-
®
£
£
-
=
+
®
=
®
=
-
=
î
í
ì
-
£
£
-
-
=
=
î
í
ì
-
£
£
-
®
-
£
£
-
-
=
®
=
=
1
1
1
1
1
0
1
1
2
1
1
2
1
2
1
2
1
0
3
4
u
v
u
y
u
y
v
y
x
C
u
u
v
u
y
u
v
x
C
v
y
u
y
v
x
u
x
v
u
y
x
¶
¶
¶
¶
¶
¶
¶
¶
=
¶
¶
)
,
(
)
,
(
î
í
ì
£
£
-
=
=
î
í
ì
-
=
=
®
=
-
£
£
=
2
2
2
2
2
2
2
2
0
1
u
v
x
u
y
v
y
x
x
C
2
1
1
1
1
)
,
(
)
,
(
1
;
1
;
1
;
1
-
=
-
=
¶
¶
-
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
î
í
ì
-
=
+
=
v
u
y
x
v
y
u
y
v
x
u
x
y
x
v
y
x
u
òò
òò
òò
¶
¶
=
¶
¶
=
R
S
S
dudv
v
u
y
x
v
u
f
dudv
v
u
y
x
v
u
h
v
u
g
f
dA
y
x
f
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
))
,
(
),
,
(
(
)
,
(
dudv
v
u
y
x
e
dA
e
S
v
u
R
y
x
y
x
)
,
(
)
,
(
¶
¶
=
òò
òò
-
+
(
)
(
)
[
]
(
)
4
ln
ln
2
2
2
2
)
,
(
)
,
(
1
2
1
1
2
1
1
2
1
2
1
-
-
-
-
-
-
=
-
=
-
=
ú
ú
ú
û
ù
ê
ê
ê
ë
é
=
=
=
¶
¶
=
ò
ò
ò
ò
òò
e
e
I
v
e
e
v
dv
e
e
dv
v
e
I
dudv
e
dudv
v
u
y
x
e
I
v
v
v
u
v
v
v
u
S
v
u
w
z
v
z
u
z
w
y
v
y
u
y
w
x
v
x
u
x
w
v
u
z
y
x
J
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
=
¶
¶
=
)
,
,
(
)
,
,
(
dudvdw
w
v
u
z
y
x
w
v
u
f
dV
z
y
x
f
E
S
òòò
òòò
¶
¶
=
)
,
,
(
)
,
,
(
)
,
,
(
)
,
,
(
z
z
rsen
y
r
x
=
=
=
q
q
cos
2
1
r
r
r
£
£
2
1
q
q
q
£
£
2
1
z
z
z
£
£
z
z
rsen
y
r
x
=
=
=
;
;
cos
q
q
z
z
z
r
z
z
y
y
r
y
z
x
x
r
x
z
r
z
y
x
J
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
=
¶
¶
=
q
q
q
q
)
,
,
(
)
,
,
(
1
;
0
;
0
;
0
;
cos
;
0
;
;
cos
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
-
=
¶
¶
=
¶
¶
z
z
z
r
z
z
y
r
y
sen
r
y
z
x
rsen
x
r
x
q
q
q
q
q
q
q
r
rsen
r
z
r
z
y
x
J
r
sen
rsen
r
sen
rsen
z
r
z
y
x
J
=
+
=
¶
¶
=
=
-
-
=
-
=
¶
¶
=
q
q
q
q
q
q
q
q
q
q
q
q
2
2
cos
)
,
,
(
)
,
,
(
cos
cos
1
0
0
0
cos
0
cos
)
,
,
(
)
,
,
(
òòò
òòò
òòò
=
¶
¶
=
=
S
E
S
dz
rdrd
z
r
f
I
valores
emplazando
dz
drd
z
r
z
y
x
z
r
f
dV
z
y
x
f
I
q
q
q
q
q
)
,
,
(
:
Re
)
,
,
(
)
,
,
(
)
,
,
(
)
,
,
(
ò
ò
ò
òòò
=
=
2
1
2
1
2
1
r
r
z
z
E
dr
dzd
r
dz
rdrd
V
q
q
q
q
1
1
1
)
1
(
2
4
2
2
2
2
2
2
2
2
+
-
±
=
®
=
-
+
®
=
+
=
+
+
x
y
y
x
y
y
x
z
y
x
}
1
1
1
1
;
2
0
/
)
,
(
{
}
4
4
;
)
,
(
/
)
,
,
{(
2
2
2
2
2
2
+
-
£
£
+
-
-
£
£
=
-
-
£
£
-
-
Î
=
x
y
x
x
y
x
R
y
x
z
y
x
R
y
x
z
y
x
E
y
y
x
2
2
2
=
+
2
2
2
2
4
4
y
x
z
y
x
-
-
£
£
-
-
2
2
4
4
r
z
r
-
£
£
-
-
ò
ò
ò
-
-
-
=
p
q
q
0
2
0
4
4
2
2
sen
r
r
rdzdrd
V
ò
ò
-
=
2
/
0
2
0
2
4
2
2
p
q
q
sen
drd
r
r
q
q
p
d
sen
r
0
2
)
4
(
3
2
2
2
/
3
2
2
/
0
-
-
=
ò
ò
-
=
2
/
0
3
)
cos
8
8
(
3
4
p
q
q
d
ò
-
-
=
2
/
0
2
))
1
)(
(cos
1
(
3
32
p
q
q
q
d
sen
0
2
/
3
3
32
3
p
q
q
q
ú
û
ù
ê
ë
é
+
-
=
sen
sen
)
4
3
(
9
16
-
=
p
f
q
f
q
f
cos
cos
r
z
sen
rsen
y
rsen
x
=
=
=
f
q
f
q
f
cos
;
;
cos
r
z
sen
rsen
y
rsen
x
=
=
=
dudvdw
w
v
u
z
y
x
w
v
u
f
dV
z
y
x
f
E
S
òòò
òòò
¶
¶
=
)
,
,
(
)
,
,
(
)
,
,
(
)
,
,
(
f
q
f
q
f
q
f
q
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
¶
=
¶
¶
=
z
z
r
z
y
y
r
y
x
x
r
x
r
z
y
x
J
)
,
,
(
)
,
,
(
f
f
q
f
q
f
f
q
f
q
q
f
q
f
f
q
f
q
q
f
sen
r
z
z
r
z
sen
r
y
rsen
y
sen
sen
r
y
r
x
sen
rsen
x
sen
r
x
-
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
=
¶
¶
-
=
¶
¶
=
¶
¶
;
0
;
cos
cos
;
cos
;
cos
cos
;
;
cos
f
f
q
f
f
f
f
f
q
q
q
f
f
q
q
f
q
f
q
f
i
q
f
f
q
f
f
q
f
f
q
f
q
f
q
f
q
f
q
f
q
f
f
q
sen
r
r
z
y
x
J
sen
r
sen
sen
r
r
z
y
x
sen
sen
r
sen
sen
r
sen
sen
r
sen
r
sen
sen
r
sen
r
r
z
y
x
rsen
sen
r
rsen
sen
sen
r
sen
rsen
sen
r
z
y
x
J
2
2
2
2
2
2
2
2
2
2
2
3
2
2
3
2
2
2
2
2
2
2
2
3
2
)
,
,
(
)
,
,
(
)
cos
(
)
,
,
(
)
,
,
(
)
cos
(
cos
)
(cos
cos
cos
cos
cos
)
,
,
(
)
,
,
(
0
cos
cos
cos
cos
cos
cos
)
,
,
(
)
,
,
(
=
¶
¶
=
-
=
+
-
=
¶
¶
+
-
+
-
=
-
-
-
-
=
¶
¶
=
-
-
=
¶
¶
=
òòò
òòò
=
¶
¶
=
S
S
d
d
dr
sen
r
r
f
d
d
dr
r
z
y
x
r
f
I
f
q
f
f
q
f
q
f
q
f
q
2
)
,
,
(
)
,
,
(
)
,
,
(
)
,
,
(
2
2
2
y
x
z
+
=
ò
÷
÷
ø
ö
ç
ç
è
æ
-
=
p
q
2
0
2
2
1
9
d
)
2
2
(
9
-
=
p
ò
ò
=
p
p
q
f
f
2
0
4
/
0
9
d
d
sen
ò
-
=
p
q
p
f
2
0
0
4
/
cos
9
d
ò
ò
ò
=
p
p
q
f
r
f
2
0
4
/
0
3
0
2
d
d
d
sen
r
òòò
=
Q
dV
V
[
]
ò
ò
ò
ò
ò
=
=
=
2
0
3
0
2
0
2
0
3
0
2
0
p
p
q
q
drd
z
r
rdzdrd
2
9
2
2
2
2
0
3
0
2
2
0
3
0
p
q
q
p
p
=
ú
û
ù
ê
ë
é
=
=
ò
ò
ò
d
r
rdrd
V