Resolución Detallada:
Área de la esfera dentro del cilindro
Enunciado
Halla el área de la parte de la esfera
x2 + y 2 + z 2 = 4
que queda dentro del cilindro
x2
+ y 2 = 1.
4
Observa que el cilindro atraviesa toda la esfera, por lo que la región buscada consta de dos tapas
gemelas (la parte superior con z≥0 y la parte inferior con z ≤ 0).
1. Denición de supercies
Esfera. La ecuación de la esfera de radio R=2 es
x2 + y 2 + z 2 = 4.
Si la describimos como supercie z = f (x, y), tenemos
p
z = ± 4 − x2 − y 2 .
p
Para la parte superior (z = + 4 − x2 − y 2 ), el elemento de área sobre la esfera (proyectado en el
plano xy ) se obtiene de
q
dS = 1 + (fx )2 + (fy )2 dx dy.
Aquí
∂ p x y
fx = 4 − x2 − y 2 = − p , fy = − p .
∂x 4 − x2 − y 2 4 − x2 − y 2
Entonces
x2 y2 4
1 + fx2 + fy2 = 1 + + = ,
4 − x2 − y 2 4 − x2 − y 2 4 − x2 − y 2
y se concluye
r
4 2
dS = dx dy = p dx dy.
4 − x2 − y 2 4 − x2 − y 2
p q
Para la parte inferior (z = − 4 − x2 − y 2 ) es exactamente la misma expresión de dS , porque fx2 + fy2
no cambia signo.
1
Cilindro. El cilindro está dado por
x2
+ y 2 = 1.
4
Proyectado en el plano xy , dene la región
x2
D = (x, y) ∈ R2 : + y2 ≤ 1 ,
4
que es una elipse de semiejes 2 en x-dirección y 1 en y -dirección.
2. Integral de área
Queremos calcular ZZ ZZ
2
A= dS = p dx dy.
D D 4 − x2 − y 2
Como D es la elipse { x2 /4 + y 2 ≤ 1}, haremos un cambio que la convierta en el círculo unidad.
3. Cambio de variable: elipse → círculo
Denimos
x = 2 u, y = v,
entonces
x2
+ y 2 = u2 + v 2 ≤ 1,
4
es decir, el dominio en (u, v) es el círculo unidad
(u, v) : u2 + v 2 ≤ 1 .
Duv =
El jacobiano es
2 0
dx dy = det D(x, y)/D(u, v) du dv = det du dv = 2 du dv.
0 1
| {z }
=2
Además,
x2 + y 2 = (2u)2 + v 2 = 4u2 + v 2 , 4 − x2 − y 2 = 4 − (4u2 + v 2 ) = 4 − 4u2 − v 2 .
2
Por tanto, el integrando p dx dy se transforma en
4 − x2 − y 2
2 4
√ × (2 du dv) = √ du dv.
4 − 4u2 − v 2 4 − 4u2 − v 2
Así, para la tapa superior (z ≥ 0) la integral queda:
ZZ
4
Ssup = √ du dv.
u2 +v 2 ≤1 4 − 4u2 − v 2
(Al nal multiplicaremos por 2 para sumar tapa inferior).
2
4. Paso a coordenadas polares en (u, v)
En (u, v) usamos polares:
u = ρ cos θ, v = ρ sin θ, 0 ≤ ρ ≤ 1, 0 ≤ θ < 2π,
con
du dv = ρ dρ dθ.
Además,
4 − 4u2 − v 2 = 4 − 4(ρ2 cos2 θ) − (ρ2 sin2 θ) = 4 − ρ2 4 cos2 θ + sin2 θ .
Denamos
a(θ) = 4 cos2 θ + sin2 θ = 1 + 3 cos2 θ.
Entonces
4 − 4u2 − v 2 = 4 − ρ2 a(θ),
y el elemento de área se convierte en
4
dS = p (ρ dρ dθ).
4 − ρ2 a(θ)
Por tanto,
Z 2π Z 1
4ρ
Ssup = p dρ dθ.
θ=0 ρ=0 4 − a(θ) ρ2
5. Integral interior en ρ
Sea Z 1
4ρ
I a(θ) = p dρ,
0 4 − a(θ) ρ2
Z 1
4ρ
donde a(θ) = 1 + 3 cos2 θ. Para resolver I(a) = p dρ, hacemos el cambio
0 4 − a ρ2
dt
t = 4 − a ρ2 =⇒ dt = −2a ρ dρ =⇒ ρ dρ = − .
2a
Cuando ρ = 0, t = 4. Cuando ρ = 1, t = 4 − a. Entonces
dt
1 4−a
4−a
4 − 2a
Z Z Z
4ρ 2 1
I(a) = p dρ = √ = − t− 2 dt.
ρ=0 4 − a ρ2 t=4 t a 4
Z √
Como t−1/2 dt = 2 t,
2 h √ i 4−a 4 √ √ 4 √
I(a) = − 2 t =− 4−a − 4 = 2− 4−a .
a 4 a a
Por tanto,
Z 2π Z 2π
4 p
Ssup = I a(θ) dθ = 2− 4 − a(θ) dθ,
0 0 a(θ)
2
p p √
con a(θ) = 1 + 3 cos θ, y 4 − a(θ) = 4 − (1 + 3 cos2 θ) = 3 |sin θ|.
Entonces el integrando se separa como
√
4 √ 8 4 3 |sin θ|
2− 4−a = 2
T 1(θ) − 2
T 2(θ).
a |1 + 3{zcos θ} |1 + 3{zcos θ}
3
Denimos √
2π 2π
4 3 |sin θ|
Z Z
8
P = dθ, Q = dθ.
0 1 + 3 cos2 θ 0 1 + 3 cos2 θ
Así
Ssup = P − Q.
6. Cálculo de P
Usamos la fórmula estándar (ver, por ejemplo, libros de cálculo multivariable):
Z 2π
dθ 2π
=p , a > 0, a + b > 0.
0 a + b cos2 θ a (a + b)
p √
Aquí a = 1, b = 3. Entonces a(a + b) = 1 · 4 = 2,
Z 2π
dθ 2π
2
= = π.
0 1 + 3 cos θ 2
Por ende,
Z 2π
8
P = dθ = 8 × π = 8π.
0 1 + 3 cos2 θ
7. Cálculo de Q
√
Z 2π
4 3 |sin θ| √ Z 2π |sin θ|
Q= 2
dθ = 4 3 dθ.
0 1 + 3 cos θ 0 1 + 3 cos2 θ
|sin θ|
Notemos que f (θ) = tiene periodo π y es simétrica en cada cuadrante. Entonces
1 + 3 cos2 θ
Z 2π π/2
|sin θ|
Z
sin θ
dθ = 4 dθ.
0 1 + 3 cos2 θ 0 1 + 3 cos2 θ
Z π/2
sin θ
Calculemos dθ. Hacemos el cambio
0 1 + 3 cos2 θ
u = cos θ =⇒ du = − sin θ dθ.
Cuando θ = 0, u = 1. Cuando θ = π/2, u = 0. Por tanto,
π/2 0 1
− du
Z Z Z
sin θ du
dθ = = .
0 1 + 3 cos2 θ u=1 1 + 3u2 0 1 + 3u2
Y Z
du
Z
du 1 √
= √ = √ arctan 3 u .
1 + 3u2 1 + ( 3 u)2 3
Evaluando de 0 a 1:
Z 1 √ 1 √
du 1 1 π
= √ arctan 3u = √ arctan( 3) = √ .
0 1 + 3u2 3 0 3 3 3
Luego
2π
|sin θ|
Z
π 4π
2
dθ = 4 × √ = √ .
0 1 + 3 cos θ 3 3 3 3
Finalmente,
√ 4π 16π
Q=4 3 × √ = .
3 3 3
4
8. Conclusión para la tapa superior
Tenemos
16π
P = 8π, Q= ,
3
así que
16π 24π − 16π 8π
Ssup = P − Q = 8π − = = .
3 3 3
Eso corresponde al área de la parte de la esfera con z ≥ 0 dentro del cilindro.
9. Tapa inferior por simetría
Por simetría respecto al plano z = 0, la tapa inferior (z ≤ 0) tiene la misma área:
8π
Sinf = Ssup = .
3
10. Área total dentro del cilindro
Sumamos ambas tapas:
8π 8π 16π
Stotal = Ssup + Sinf = + = .
3 3 3
11. Vericación numérica
Una integración numérica (por ejemplo, con SciPy o Monte Carlo) arroja aproximadamente:
(num) (num)
Ssup ≈ 8,3775804, Stotal ≈ 16,7551608,
mientras que
16π
≈ 16,7551608,
3
lo cual coincide dentro del margen de error de redondeo.
Resultado Final
16π 2
Stotal = unidades .
3