Calculo Vectorial Ernesto Acorta
Calculo Vectorial Ernesto Acorta
Calculo Vectorial Ernesto Acorta
VECTORIAL
ii
Indice general
1. Vectores en R2 y R3
13
3. Funciones vectoriales
29
45
5. Diferenciaci
on Parcial
57
6. Regla de la cadena
71
7. Extremos de funciones
83
8. Integrales m
ultiples
93
9. Integrales m
ultiples
111
10.Campos vectoriales
125
11.Divergencia y rotacional
137
12.Teoremas Fundamentales
149
iii
iv
INDICE GENERAL
LECCION
Vectores en R2 y R3
1. VECTORES EN R2 Y R3
LECCION
Y
y
(x, y)
Gr
afica 1.1
Dada una pareja de n
umeros reales (x, y), se ubica en el eje X el punto Px de coordenada x y en el eje Y el punto Py de coordenada y. La recta perpendicular al eje
X que contiene al puntoPx y la recta perpendicular al eje Y que contiene al punto
Py se intersectan en un punto P del plano. Este punto tiene coordenadas (x, y), ya
que las proyecciones de P sobre los ejes numericos X y Y son los puntos Px y Py
con coordenadas x y y respectivamente.
Y
Py
Px
P (x, y)
Gr
afica 1.2
El mismo procedimiento se puede emplear para introducir un sistema de coordenadas
cartesianas rectangulares en el espacio. En este caso tendremos tres ejes numericos
X, Y , Z, perpendiculares entre si. Las coordenadas (x, y, z) de un punto P en el
espacio, se obtienen con el mismo procedimiento.
3
Z
z
(x, y, z)
y Y
x
X
Gr
afica 1.3
Dos puntos P y Q definen los segmentos de recta dirigidos PQ y QP. Estos dos segmentos de recta se distinguen u
nicamente por su orientaci
on. Es decir, el segmento
dirigido PQ lo pensamos recorrido desde P hasta Q, y el segmento dirigido QP lo
pensamos recorrido desde Q hasta P .
Si introducimos un sistema de coordenadas en el espacio y (x1 , y1 , z1 ), (x2 , y2 , z2 )
son las coordenadas de P y Q, respectivamente, los tres n
umeros x2 x1 , y2 y1 y
z2 z1 , son las componentes escalares del segmento dirigido PQ. Observese que cada
una de las componentes del segmento dirigido QP es el opuesto de la componente
correspondiente del segmento dirigido PQ, de aqu la importancia de la orientaci
on
de los segmentos dirigidos. Dos segmentos dirigidos son equivalentes si sus componentes escalares correspondientes son iguales. No es difcil mostrar que teniendo que
los puntos P y R son diferentes, los segmentos dirigidos PQ y RS son equivalentes
si, y solamente si, el cuadrilatero P QSR es un paralelogramo.
La igualdad de las componentes escalares correspondientes entre segmentos dirigidos
es una relaci
on de equivalencia y cada clase de equivalencia es un vector geometrico.
Es decir, un vector (en este caso tridimensional) esta representado por un segmento dirigido PQ, o por cualquier otro segmento dirigido PQ equivalente. As, un
vector esta completamente determinado por las componentes escalares de cualquier
segmento dirigido que lo representa. Se usa frecuentemente la notaci
on ha, b, ci para
denotar el vector con componentes escalares a, b y c y para notar los vectores las
letras con tipo negrita, por ejemplo, u , v , etc. En adelante para decir que el vector
1. VECTORES EN R2 Y R3
LECCION
5
Estas operaciones tienen propiedades bien conocidas, que se deducen de las propiedades de la adicion y multiplicaci
on de los n
umeros reales, que son las que las
definen. Por ejemplo, es claro que la adicion de vectores es asociativa, conmutativa, que hay un elemento neutro (vector nulo) y que cada vector tiene su opuesto.
Adem
as, la multiplicaci
on por escalar distribuye sobre la adicion de vectores, entre
otras. El producto escalar de vectores es conmutativo y distribuye sobre la adicion
de vectores. En cambio, aunque el producto vectorial distribuye sobre la adicion de
vectores no es conmutativo.
Dividir un vector por un escalar no nulo consiste en multiplicar el vector por el
inverso multiplicativo de dicho escalar. Por ejemplo, podemos dividir a un vector no
nulo por su magnitud para obtener un vector de magnitud 1, con la misma direccion
y sentido. Los vectores de longitud 1 se denominan vectores unitarios.
Otra forma de definir el producto escalar de dos vectores, que normalmente se usa
cuando no estan referidos a un sistema de coordenadas es:
v
uv
v,
=
kvk
kvk2
1. VECTORES EN R2 Y R3
LECCION
u
u
v
v
proyv u
proyv u
Gr
afica 1.4
al escalar kuk cos se le denomina la componente escalar de u en la direcci
on de v
compv u = kuk cos =
uv
.
kvk
7
EJERCICIOS
1. Para cada par de puntos determine las componentes escalares de los segmentos
dirigidos PQ y QP.
a) P (1, 2, 1) y Q(3, 1, 2).
b) P (1, 3, 4) y Q(6, 4, 7)
c) P (2, 5) y Q(1, 2)
2. Si las coordenadas de los puntos P y Q son (3, 5) y (2, 1) respectivamente,dibuje en un plano cartesiano:
a) El segmento dirigido PQ.
b) Un segmento dirigido equivalente a PQ con punto inicial en el origen.
c) Un segmento dirigido equivalente a QP con punto inicial en el origen.
d ) Un segmento dirigido equivalente a PQ con punto inicial en R(3, 5).
3. Si P = (1, 1, 1), Q = (2, 2, 2), R = (2, 0, 1), S = (3, 1, 2), se tiene que
PQ = RS? justifique su respuesta.
4. Si P = (0, 0, 0), Q = (1, 3, 2), R = (1, 0, 1), S = (2, 3, 4), se tiene que PQ =
RS? justifique su respuesta.
5. Dados los puntos P (1, 3), Q(2, 1) y R(6, 5), determine todos los puntos
S del plano tales que P, Q, R y S sean los vertices de un paralelogramo.
6. Dados los puntos P (p1 , p2 , p3 ), Q(q1 , q2 , q3 ) y R(r1 , r2 , r3 ) respectivamente, determine todos los puntos S del plano tales que P, Q, R y S sean los vertices de
un paralelogramo.
7. Muestre que la suma de vectores es conmutativa.
8. Muestre que el producto escalar de dos vectores es conmutativo.
9. Muestre que el producto vectorial no es conmutativo. Que relaci
on existe
entre los vectores u v y v u?
10. Calcule la magnitud de cada uno de los siguientes vectores:
(a) v = h2, 1i (b) v = h2, 1, 0i (c) v = h3, 2, 2i
(d) v = h0, 0, 1i (e) v = h6, 4, 4i
1. VECTORES EN R2 Y R3
LECCION
(d) 21 (v w).
v
kvk
(b) v + w.
(c)
(e) 12 (v + w).
(f) 2 v + 4 w.
(g) v 2 w.
12. El producto vectorial de dos vectores tambien se define como: Un vector que
es perpendicular a los dos vectores dados, que tiene la orientaci
on dada por
la regla de la mano derecha y cuya norma es el producto de las normas de
los vectores y el seno del menor angulo formado por ellos. Muestre que esta
definicion es equivalente a la dada en este texto.
13. Dados los puntos P (x1 , y1 , z1 ), Q(x2 , y2 , z2 ), R(a, b, c), determine las condiciones que deben satisfacer las coordenadas de todos los puntos S(x, y, z) del espacio para que los vectores PQ y RS sean paralelos. Describa el lugar geometrico
de dichos puntos.
14. Dados los puntos P (x1 , y1 , z1 ), Q(x2 , y2 , z2 ), R(a, b, c), determine las condiciones que deben satisfacer las coordenadas de todos los puntos S(x, y, z) del
espacio para que los vectores PQ y RS sean perpendiculares. Describa el lugar
geometrico de dichos puntos.
15. Determine el
area del tri
angulo que tiene vertices en los puntos (0, 7, 4), (2, 3, 6)
y (1, 4, 5) .
16. Determine el volumen del tetraedro que tiene sus vertices en los puntos (1, 2, 1),
(5, 2, 2), (3, 1, 4) y (1, 1, 3).
17. Cuales rectas en R2 son subespacios de R2 ? Y de R3 ? (Las soluciones de
las ecuaciones son parejas de n
umeros reales que representan las componentes
escalares de vectores).
18. Describa todos los subespacios vectoriales de R3 .
19. Encuentre la matriz de la transformacion lineal de R3 en R3 , que corresponde
a la reflexi
on con respecto al plano de ecuaci
on x = y. Encuentre la que
corresponde a la reflexi
on con respecto al plano de ecuaci
on ax + by + cz = 0.
Es la reflexi
on con respecto al plano de ecuaci
on ax + by + cz = d, d 6= 0, una
transformacion lineal de R3 en R3 ?
9
20. Si u, v son vectores y A es una matriz, cu
ales de los siguientes son n
umeros
y cu
ales son vectores: u v, kuk, u v, det A, Au.
21. Sea P el paraleleppedo en R3 definido por 0 x a, 0 y b y 0 z c.
a) Encuentre el
angulo entre la diagonal de P y sus lados.
b) Encuentre el
angulo entre la diagonal de P y sus caras.
c) Que tiene que ver el
area de las caras de P y sus angulos con la diagonal
de P ?
22. Ejemplos tpicos de vectores unitarios son los vectores i =h1, 0, 0i, j =h0, 1, 0i
y k =h0, 0, 1i, que forman una base ortonormal del espacio de vectores tridimensionales, conocida como la base canonica. Claramente, los vectores i , j y
k son vectores ortogonales entre s. Sabemos que todo vector de R3 se puede
escribir como una combinacion lineal de los vectores i, j y k. Se puede escribir
todo vector de R3 como una combinacion de i y j u
nicamente? Justifique su
respuesta.
23. Calcule v w:
a) v = h5, 1, 2i y w = h4, 4, 3i
b) v = 3 i 2 j k y w = 6 i + 4 j + 2 k
24. Encuentre el
angulo entre los vectores v y w:
a) v = h5, 1, 2i, w = h4, 4, 3i
d ) v = i, w = 3 i + 2 j + 4k
10
1. VECTORES EN R2 Y R3
LECCION
d ) v = i, w = 3 i + 2 j + 4k
29. Calcule el
area del tri
angulo P QR:
a) P = (5, 1, 2), Q = (4, 4, 3), R = (2, 4, 0)
b) P = (4, 0, 2), Q = (2, 1, 5), R = (1, 0, 1)
30. Calcule el
area del paralelogramo P QRS:
a) P = (2, 1, 3), Q = (1, 4, 5), R = (2, 5, 3), S = (3, 2, 1)
b) P = (2, 2), Q = (1, 4), R = (6, 6), S = (3, 0)
31. Encuentre el volumen del paraleleppedo con lados adyacentes u, v, w:
a) u = h1, 1, 3i, v = h2, 1, 4i, w = h5, 1, 2i
b) u = h1, 3, 2i, v = h7, 2, 10i, w = h1, 0, 1i
32. Calcule u (v w) y u (v w):
a) u = h1, 1, 1i, v = h3, 0, 2i, w = h2, 2, 2i
b) u = h1, 0, 2i, v = h1, 0, 3i, w = h2, 0, 2i
33. Calcule (u v) (w z) si u = h1, 1, 1i, v = h3, 0, 2i, w = h2, 2, 2i, z = h2, 1, 4i.
34. Si v y w son vectores unitarios en R3 , bajo que condiciones se tendra que
v w fuera tambien un vector unitario en R3 ? Justifique su respuesta.
11
35. Muestre que si v w = 0 para todo w en R3 , entonces v = 0.
36. Demuestre que, para todo v y w en R3 se tiene que:
a) kv wk2 + |v w|2 = kvk2 kwk2
b) Si v w = 0 y v w = 0, entonces v = 0 o w = 0.
37. Encuentre todas las soluciones de la ecuaci
on vectorial a x = b en R3 ,
donde a 6= 0.
38. Demuestre la identidad de Jacobi :
u (v w) + v (w u) + w (u v) = 0
39. Muestre que u, v, w estan en un mismo plano en R3 si y solo si u (v w) = 0.
40. Para todo los vectores u y v en R3 , muestre que
u (v u) = (u u)v (u v)u.
41. Use el ejercicio anterior para mostrar que para todo los vectores u, v, w, z en
R3 ,
(u v) (w z) = (z (u v))w (w (u v))z
y que
(u v) (w z) = (u (w z))v (v (w z))u
Tienen significado geometrico estas ecuaciones?
42. Calcule las magnitudes de los siguientes vectores:
(a) v = h2, 1i (b) v = h2, 1, 0i (c) v = h3, 2, 2i
43. Para los siguientes puntos, se tiene que P Q = RS?
a) P = (1, 1, 1), Q = (2, 2, 2), R = (2, 0, 1), S = (3, 1, 2)
b) P = (0, 0, 0), Q = (1, 3, 2), R = (1, 0, 1), S = (2, 3, 4)
12
1. VECTORES EN R2 Y R3
LECCION
(b) Calcule v + w.
(e) Calcule
12 (v + w)
.
(c) Calcule
v
.
kvk
LECCION
La introducci
on de un sistema de coordenadas en el espacio permite representar
algebraicamente subconjuntos del mismo . Este procedimiento consiste en encontrar
las condiciones algebraicas que satisfacen las coordenadas de todos los puntos que
pertenecen a dicho subconjunto. En general, la representaci
on algebraica de un subconjunto de puntos del espacio se puede hacer mediante ecuaciones, inecuaciones o
un sistema de ecuaciones e inecuaciones con tres incognitas (las coordenadas de los
puntos). M
as precisamente, una ecuaci
on con tres incognitas representa un conjunto
de puntos del espacio si y s
olo si las coordenadas de cada punto del conjunto satisfacen la ecuaci
on y cada tripla de n
umeros reales que satisface la ecuaci
on son las
coordenadas de un punto del conjunto.
Por ejemplo, dados un punto Q y un vector v , el lugar geometrico de todos los
puntos P del espacio para los que los vectores QP y v son ortogonales es un plano.
Si las coordenadas de Q son (x0 , y0 , z0 ), las de P son (x, y, z) y v = ha, b, ci, la
condici
on de ortogonalidad, es decir, el que QP v = 0, se escribe en coordenadas
as:
a(x x0 ) + b(y y0 ) + c(z z0 ) = 0
que es la ecuaci
on del plano que contiene a Q y tiene vector normal v =ha, b, ci.
13
14
El lugar geometrico de todos los puntos P (x, y, z) del espacio cuya distancia a un
punto fijo Q(x0 , y0 , z0 ) es un n
umero real positivo r, es una esfera de radio r y centro
Q. Como lo vimos en la leccion anterior, la condici
on geometrica la distancia de P
a Q es r, se escribe algebraicamente
p
(x x0 )2 + (y y0 )2 + (z z0 )2 = r
o,
(x x0 )2 + (y y0 )2 + (z z0 )2 = r2
que es precisamente la ecuaci
on de dicha esfera.
Sea C una curva en un plano , el conjunto de todas las rectas perpendiculares
al plano que intersectan a la curva C es una superficie cilndrica (o cilndro).
Por ejemplo, dada la circunferencia C en el plano XY con centro en el origen y
radio r, el conjunto de todas las rectas perpendiculares al plano XY que intersectan
la circunferencia C es una superficie cilndrica circular recta. Un punto P (x, y, z)
pertenece a esta superficie si, y solamente si, su proyecci
on ortogonal sobre el plano
XY es un punto de la circunferencia, es decir, si las coordenadas x y y del punto
P satisfacen la ecuaci
on de la circunferencia C. As, el punto P (x, y, z) esta en
la superficie si y solamente si, x2 + y 2 = r2 , que es precisamente la ecuaci
on de
la superficie . Observese que la ecuaci
on no impone ninguna restriccion sobre la
coordenada z.
Z
Y
C
15
Gr
afica 2.1
o
Ax2 + By 2 + Iz = 0.
Por ejemplo, xy = 1 es la ecuaci
on de una hiperbola, pero rotando 45 grados los ejes
X y Y en el sentido contrario al movimiento de las manecillas del reloj, la ecuaci
on
2
2
de la misma hiperbola es x y = 2, en la nuevas coordenadas.
Para representar gr
aficamente una superficie a partir de su ecuaci
on, se busca informacion haciendo cortes de la superficie con planos (frecuentemente paralelos a los
planos coordenados). Por ejemplo, si se intersecta la superficie cuadratica de ecuaci
on z = x2 + y 2 con un plano de ecuaci
on z = k, k > 0, se obtiene la ecuaci
on
2
2
x + y = k que para cada valor de k es la ecuaci
on de una circunferencia en el plano
z = k con centro en el punto (0, 0, k). Por otro lado, si se intersecta la superficie con
un plano de ecuaci
on x = k, se obtiene z = y 2 + k 2 que para cada valor de k es la
ecuaci
on de una par
abola en el plano x = k. Es decir, la interseccion de la superficie
con el plano x = k es una par
abola que abre hacia arriba, al intersectar la superficie con un plano de ecuaci
on y = k tambien se obtienen una par
abola. Dibujando
simult
aneamente estas circunferencias y par
abolas para diferentes valores de k, se
obtiene un bosquejo de la superficie en el espacio, que es un paraboloide circular con
eje de simetra el eje Z.
16
Las curvas de interseccion obtenidas con este procedimiento se conocen como trazas.
En este caso las trazas obtenidas son las siguientes:
x2 + y 2 = k
z
= k
Y
k
k
k
k
=
=
=
=
4
3
2
1
k= 1
4
Gr
afica 2.2
Que son trazas paralelas al plano XY por estar sobre planos paralelos al mismo, en
algunos contextos se denominan trazas horizontales.
z = k2 + y2
x = k
z = x2 + k 2
y = k
Que son trazas perpendiculares al plano XY por estar sobre planos perpendiculares
al mismo, en algunos contextos se denominan trazas verticales.
Z
k=
k=
k=
k=
k=1
k=0
Gr
afica 2.3
k=1
k=0
Grfica 2.4
17
Gr
afica 2.5
Y
C
18
EJERCICIOS
1. Determine la ecuaci
on de la esfera con centro en el origen y radio 2.
2. Determine la ecuaci
on de la esfera con centro en el punto C(3, 2, 4) y radio 5.
3. Determine la ecuaci
on de la esfera con centro en el punto C(1, 2, 5) y que
contiene al origen.
4. Determine el centro y el radio de la esfera de ecuaci
on x2 +y 2 +z 2 4x+6z = 3.
5. Determine el centro y el radio de la esfera de ecuaci
on 2x2 + 2y 2 + 2z 2 + 4x
y + 3z = 0.
6. Determine la ecuaci
on de una esfera con centro en el punto (1, 3, 2) y tangente
a la esfera de ecuaci
on x2 2x + y 2 + 4y + z 2 + 2z + 2 = 0.
7. Determine la ecuaci
on del plano que contiene el punto (3, 5, 1) y es perpendicular al vector h2, 7, 4i.
8. Determine la ecuaci
on del plano que contiene los puntos (2, 3, 5), (1, 1, 3),
(2, 4, 7).
9. Determine la ecuaci
on del plano paralelo al plano de ecuaci
on x 2y + 7z = 4
y que contiene al punto (7, 3, 4).
10. z Es posible determinar la ecuaci
on de un plano que contiene al punto (2, 1, 3),
es perpendicular al plano de ecuaci
on x y + 2z = 5 y al vector h1, 8, 4i?.
11. Determine la ecuaci
on del plano que contiene al punto (2, 1, 3), es perpendicular al plano de ecuaci
on x y + 2z = 5 y al vector h0, 8, 4i.
12. Determine la ecuaci
on del plano que contiene al punto (2, 1, 3), es perpendicular al plano de ecuaci
on x y + 2z = 5 y paralelo al vector h1, 8, 4i.
13. Las siguientes ecuaciones representan superficies cuadraticas, para cada una
dibuje las trazas y a partir de ellas haga un bosquejo de la gr
afica. De los
siguientes nombres deduzca cu
al le corresponde cada ecuaci
on: Elipsoide, Paraboloide, Paraboloide Hiperb
olico, Hiperboloide de una hoja,Hiperboloide de
dos hojas, Cono de doble hoja.
19
a) z = x2 + y 2 .
b) x2 + y 2 = z 2
c) x2 + y 2 = z 2 1
d ) x2 + y 2 = z 2 + 1
e) x2 y 2 = z
f)
x2 y 2 z 2
+ 2 + 2 =1
a2
b
c
Adem
as del sistema de coordenadas cartesianas para el espacio, com
unmente se usan
el sistema de coordenadas cilndricas y el sistema de coordenadas esfericas.
Un sistema de coordenadas cilndricas se introduce fijando un plano en el espacio
en el que se considera un sistema de coordenadas polares y un eje numerico W
perpendicular a ese plano de tal forma que el origen del eje coincida con el polo de
las coordenadas polares. Las coordenadas cilndricas de un punto P seran (r, , w)
donde (r, ) son las coordenadas polares de la proyecci
on ortogonal de P sobre el
plano, y w es la coordenada de la proyecci
on ortogonal de P sobre el eje numerico.
W
w
(r, , w)
r
Gr
afica 3.1
20
Obviamente, la representaci
on algebraica de un conjunto de puntos en el espacio,
depende del sistema de coordenadas que se escoge, por ejemplo si se ha introducido
un sistema de coordenadas cartesianas en el espacio para representar un conjunto de
puntos, y se quiere determinar la representaci
on usando un sistema de coordenadas
cilndricas se puede hacer escogiendo cualquiera de los planos XY , XZ o Y Z para
poner coordenadas polares (de manera conveniente) y el eje numerico restante, Z, Y
o X, respectivamente, para la tercera coordenada. En el primer caso, por ejemplo, se
puede establecer la relaci
on existente entre las coordenadas (x, y, z) y las coordenadas
(r, , w) de un punto P . En efecto, basta recordar la relaci
on que existe entre las
coordenadas cartesianas y las coordenadas polares de un punto en el plano. As, se
tiene que
p
r =
x2 + y 2
x = r cos ,
r [0, )
y
= arctan
y
y = r sen ,
[0, 2)
x
z = w,
w (, +)
w = z.
Teniendo en cuenta esta equivalencia entre las coordenadas de los puntos en los dos
sistemas , la ecuaci
on cartesiana de la esfera con centro en el punto (0, 0, 0) (origen
del sistema) y radio a, que es x2 + y 2 + z 2 = a2 corresponde a la ecuaci
on cilndrica
2
2
2
r +z =a .
Por otro lado, si se escoge poner el sistema polar en el plano XZ, se tiene que
x = r sen
y = w
z = r cos .
En este caso la ecuaci
on cilndrica de la esfera es y 2 + r2 = a2 .
Un sistema de coordenadas esfericas en el espacio se introduce fijando un semiplano,
el semiplano esferico, delimitado por una recta numerica, eje esferico, cuyo origen
es el punto O. Las coordenadas esfericas (, , ) de un punto P , se obtienen determinando la distancia del origen O al punto P , el angulo diedro (en sentido
dextrogiro) formado por el semiplano esferico y el semiplano delimitado por el eje
esferico y que contiene al punto P , y el angulo que forma el segmento dirigido OP
con el semieje esferico positivo.
21
P (, , )
semiplano esferico
Al igual que en el caso anterior, se puede establecer una equivalencia entre las coordenadas cilndricas y las coordenadas esferica de un punto en el espacio. Escogiendo,
por ejemplo, el eje W como eje esferico, el semiplano delimitado por el eje W que
contiene el eje polar como semiplano esferico y el origen del sistema como el polo O,
la relaci
on entre las coordenadas (r, , ) y las coordenadas (, , ) de un punto P
se obtiene de la siguiente manera. De acuerdo con la siguiente gr
afica
W
Q
r
P (, , )
w
O
el tri
angulo OQP es rectangulo y el angulo QOP es por lo tanto w = |OQ| =
cos , y r = |QP | = sen . Observese que en estas circunstancias el angulo es
22
[0, )
[0, 2)
[0, )
=
r2 + w2
=
r
= arctan .
w
Teniendo en cuenta lo establecido antes entre los sistemas de coordenadas cilndricas y cartesianas obtenemos la siguiente relaci
on entre las coordenadas esfericas y
cartesianas:
p
=
x2 + y 2 + z 2
x = sen cos ,
y
=
arctan
y
y = sen sen ,
x
p
x2 + y 2
z = cos ,
= arctan
.
z
23
EJERCICIOS
1. En un sistema de coordenadas rectangulares describa el lugar geometrico representado por cada una de las siguientes ecuaciones y haga un bosquejo de la
gr
afica:
a) x = 2
b) y = /4
c) z = 1
d) x = y
2. En un sistema de coordenadas cilndricas describa el lugar geometrico representado por cada una de las siguientes ecuaciones y haga un bosquejo de la
gr
afica:
a) r = 2
b) = /4
c) z = 1
d) z = r
e) r = cos
f ) r = sen
3. En un sistema de coordenadas esfericas describa el lugar geometrico representado por cada una de las siguientes ecuaciones y haga un bosquejo de la
gr
afica:
a) = 2
b) = /4
c) = 1
d ) = cos
e) = sen
f ) = cos
g) = sen
24
25
d ) Como se determina la distancia de un punto a un plano?
e) C
omo se determina la distancia entre dos planos paralelos?
f ) Como se determina la interseccion entre dos rectas?
g) C
omo se determina la interseccion entre dos planos?
h) Como se determina la distancia entre dos puntos, dadas sus coordenadas
cilndricas (sin pasarlos a coordenadas rectangulares)?
i ) Como se determina la distancia entre dos puntos dadas sus coordenadas
esfericas (sin pasarlos a coordenadas rectangulares)?
6. Describa la recta L que pasa por el punto P y es paralela al vector v en la
forma: (a) vectorial, (b) parametrica, y (c) simetrica si:
a) P = (2, 3, 2), v = (5, 4, 3)
26
x6
= y+3 = z
4
c) 2x2 + 2y 2 + 2z 2 + 4x + 4y + 4z 44 = 0
d ) x2 + y 2 z 2 + 12x + 2y 4z + 32 = 0
x2
a2
x2
a2
y2
b2
y2
b2
z
c
z2
c2
= 1 en el plano
en el plano xy.
27
21. Describa un procedimiento para encontrar la ecuaci
on de una esfera que pasa
por cuatro puntos no coplanares dados. Encuentre la ecuaci
on de la esfera que
pasa por los puntos (0, 0, 0), (0, 0, 2), (1, 4, 3) y (0, 1, 3).
22. Muestre que el hiperboloide de una sola hoja es una superficie doblemente
reglada, es decir, cada punto de la superficie esta en dos rectas contenidads
completamente en la superficie.
23. Muestre que el paraboloide hiperb
olico es una superficie doblemente reglada.
24. Sea S la esfera con radio 1 y centro (0, 0, 1), y sea S = S N , donde N es el
polo norte de S: (0, 0, 2). Sea P (a, b, c) un punto arbitrario en S . La recta que
pasa por N y por P intersecta el plano xy en alg
un punto Q(x, y, 0). Encuentre
este punto Q en terminos de a, b y c. El punto Q es la proyecci
on estereogr
afica
de P . Esta correspondencia identifica S con R2 .
25. Encuentre las coordenadas (a) cilndricas y (b) esfericas de los puntos cuyas
coordenadas cartesianas se dan:
a) (2, 2 3, 1)
b) (5, 5, 6)
c) ( 21, 7, 0)
d ) (0,
2, 2)
26. Escriba las ecuaciones dadas en coordenadas (a) cilndricas y (b) esfericas:
a) x2 + y 2 + z 2 = 25
b) x2 + y 2 = 2y
c) x2 + y 2 + 9z 2 = 36
28
LECCION
Funciones vectoriales
Las funciones de las que trata esta leccion son funciones definidas en alg
un intervalo
I de n
umeros reales y que toman valores en el espacio de vectores bidimensionales
o tridimensionales. M
as precisamente, una funcion vectorial (tridimesional) es una
funcion
r : I R3 , dada por r(t) = hx(t), y(t), z(t)i
(3.1)
Por ser una funcion que toma valores en R3 , tiene tres funciones componentes, a
saber, x, y, z. Cada una de estas es una funcion escalar definida en I y que toma
valores en R. El intervalo I, que es el dominio de la funcion vectorial puede ser
de cualquiera de las siguientes formas: (a, b), (a, b], [a, b), [a, b], [a, +), (a, +),
(, b], (, b) o (, +).
Siendo x, y, z funciones definidas en un intervalo real y que toman valores en R, se
extienden naturalmente los conceptos estudiados en los cursos de calculo diferencial
e integral, para funciones reales de valor real, a las funciones vectoriales. Por ejemplo,
si los lmites lmta x(t), lmta y(t) y lmta z(t) existen, el lmite lmta r (t) existe
y
D
E
lm r (t) = lm x(t), lm y(t), lm z(t) .
(3.2)
ta
ta
ta
ta
30
y
3. FUNCIONES VECTORIALES
LECCION
dz
(a) existen, la derivada de la funcion vectorial r =r (t) en t = a es
dt
dx
dy
dz
dr
(a) =
(a), (a), (a) .
dt
dt
dt
dt
(3.3)
(r s)(t) = hy(t)e
z (t) z(t)e
y (t), z(t)e
x(t) x(t)e
z (t), x(t)e
y (t) y(t)e
x(t)i ,
la funcion escalar
(r s)(t) = x(t)e
x(t) + y(t)e
y (t) + z(t)e
z (t),
si f es una funcion escalar definida en el mismo intervalo, la funcion vectorial
f r(t) = hf (t)x(t), f (t)y(t), f (t)z(t)i ,
si r es derivable en I, la funcion vectorial
r (t) = x (t), y (t), z (t) ,
31
posicion de un punto P de coordenadas (x(t), y(t), z(t)). As, cuando t recorre el intervalo I, el punto P describe una curva en el espacio, cuyas ecuaciones parametricas
estan dadas por
x = x(t),
y = y(t),
z = z(t).
(3.5)
En otras palabras, el conjunto
(x, y, z) R3 : x = x(t),
y = y(t),
z = z(t),
tI ,
(3.6)
y b = tv2 ,
z c = tv3 ,
o
x = a + tv1 ,
y = b + tv2 ,
z = c + tv3 ,
t R.
(3.7)
32
3. FUNCIONES VECTORIALES
LECCION
Si adem
as z = (2 cos t) + 3, el punto (x, y, z) pertenece al plano z = x + 3. As, una
parametrizacion de la curva de interseccion de la superficie cilndrica y el plano es
x = 2 cos t,
y = 2 sen t,
z = 2 cos t + 3,
t [0, 2)
33
EJERCICIOS
dr
r(a + h) r(a)
(a) = lm
, utilice las propiedades de la suma
h0
dt
h
de vectores, el producto por escalar y (3.2) para demostrar (3.3).
1. Sabiendo que
a)
a) r(t) = t2 , ln t, t .
D
E
b) r(t) =
t2 4, t, 1/t .
4. Determine la integral de la funcion r(t) = t2 , 2 sen t, ln t en el intervalo [1, ]
5. Determine una funcion vectorial r tal que r (t) = hsen t, cos t, ti y r(/3) =
h2, 5, 1i.
D
E
34
3. FUNCIONES VECTORIALES
LECCION
8a3
.
x2 + 4a2
35
Geometra de las curvas y movimiento en el espacio
Considere la curva en el espacio, cuyas ecuaciones parametricas son
x = x(t),
y = y(t),
z = z(t),
t I.
dr
Supongamos, que r es derivable en I y que
(t) = r (t) 6= 0, para todo t I .
dt
Entonces, el vector
dy
dz
dx
es un vector tangente a la curva en el punto P (x(t), y(t), z(t)). As, el vector unitario
T(t) =
r (t)
,
kr (t)k
es un vector unitario tangente a la curva en el punto P (x(t), y(t), z(t)). Una propiedad interesante de los vectores T(t) y T (t) es que son perpendiculares para todo
t. Este hecho se deduce facilmente de calcular la derivada de la funcion constante
(T T)(t) = 1. En efecto,
0=
d
[T(t) T(t)] = 2T(t) T (t)
dt
36
3. FUNCIONES VECTORIALES
LECCION
Recordemos que la longitud de arco entre P (x(a), y(a)) y Q(x(t0 ), y(t0 )) de una
curva C en el plano, parametrizada por
x = x(t),
es
t0
a
y = y(t),
t I,
p
x ( )2 + y ( )2 d,
t0
a
p
x ( )2 + y ( )2 + z ( )2 d.
Observese que, siendo la funcion vectorial r (t) = hx(t), y(t), z(t)i una parametrizaci
on de la curva C, la integral anterior se puede escribir en la forma
Z
t0
a
r ( )
d.
r ( )
d.
En el caso en que esta funcion fuese invertible (por lo general lo es), podemos escribir
las coordenadas del punto Q en terminos de la longitud de arco as: Q(x(t(s)), y(t(s)), z(t(s))),
donde t = t(s) es la inversa de s. Se obtienen as la parametrizacion por longitud de
arco de la curva C:
x=x
e(s) = x(t(s)),
y = ye(s) = y(t(s))
z = ze(s) = z(t(s))
t I.
37
Tenemos entonces que r (t) = e
r(s(t)) y por lo tanto,
r (t) =
d
de
r
ds
de
r
[e
r(s(t))] = (s(t)) (t) = (s(t))
r (t)
dt
ds
dt
ds
La direccion de la curva C en cada punto Q(x(t), y(t), z(t)) esta representada por
el vector T (t). La curvatura de C es una medida de la rapidez con que cambia la
direccion de la curva con respecto al desplazamiento sobre ella. Esta medida viene
dada por
dT(t(s))
=
ds
donde s es el par
ametro de longitud de arco. Si desarrollamos la derivada
usando la regla de la cadena obtenemos.
dT(t(s))
ds
dT
(t(s))
dT(t(s))
dT
dt
dT
1
=
(t(s)) (s) =
(t(s))
.
= dt
ds
ds
dt
ds
dt
kr (t(s))k
(t(s))
dt
De aqu se obtiene la relaci
on
dT
dt (t(s))
= r (t(s))
que necesitaremos m
as adelante.
38
3. FUNCIONES VECTORIALES
LECCION
y por lo tanto,
2
ds 2
) T T =
r
T T .
dt
Calculando la norma de este producto, teniendo en cuenta que T y T son perpendiculares y que kT k = kr k se obtiene
2
2
3
r r
=
r
T T
=
r
T
=
r
.
r r = (
T (t)
, que es un vector unitario perpendicukT (t)k
lar a T (t). Es decir, N es un vector unitario perpendicular a la curva C en el punto
P (x(t), y(t), z(t)) que se conoce como vector normal unitario. Este vector indica la
direccion en la que se dobla curva. El plano que contiene al punto P (x(t), y(t), z(t))
y es paralelo a T (t) y a N(t), se denomina plano osculador de la curva en el punto
P (x(t), y(t), z(t)). Un vector normal a este plano es
Consideremos ahora el vector N(t) =
B = T N,
que tambien es unitario y que se conoce como vector binormal de C en el punto P .
En el plano osculador, la circunferencia que pasa por el punto P , tiene los mismo
vector tangente unitario y normal unitario que C en P y tiene radio r = 1/ se
denomina circunferencia osculadora y es la que mejor aproxima a la curva C en el
punto P , esto ultimo quiere decir que la curva C y la circunferencia tienen la misma
curvatura en el punto P .
En el contexto de la cinematica la posicion de un objeto que se mueve en el espacio se
puede modelar mediante una funcion vectorial r = r(t), t I, tomando el par
ametro
t como el transcurrir del tiempo.
As, la derivada r (t) = v(t) representa la velocidad del objeto y contiene la informacion de como cambia su posicion con respecto al transcurrir del tiempo en cada
punto P (x(t), y(t), z(t)) en que se encuentra el objeto, es decir, la direccion de su
39
movimiento y la rapidez kv(t)k con que lo hace.
La derivada v (t) = a(t), representa la aceleracion del objeto en el instante t y contiene la informaci
on de como cambia su velocidad con respecto al tiempo en cada
punto P (x(t), y(t), z(t)) en la que se encuentra.
El cambio en la velocidad con respecto al tiempo, involucra la razon de cambio en
la direccion del movimiento y la razon de cambio en la rapidez del movimiento (las
dos caractersticas de la velocidad). La primera se denomina aceleraci
on normal aN
y la segunda se denomina aceleraci
on tangencial aT .
Para encontrar la descomposicion de la aceleracion en terminos de la aceleracion
ds
normal y la tancencial calculemos la derivada de v(t) = (t)T(t). Obtenemos:
dt
a(t) = v (t) =
ds
d2 s
(t)T(t) + (t)T (t)
2
dt
dt
ds
d2 s
(t)k T (t)
(t)T(t)
+
(t)
kT
dt2
dt
kT (t)k
2
ds
d2 s
(t)T(t) +
(t) N(t)
dt2
dt
=
=
40
3. FUNCIONES VECTORIALES
LECCION
EJERCICIOS
1. Dada la elipse de ecuaci
on cartesiana
x2 y 2
+
=1
4
9
2
a) r(t) = t + 1, t2 + 1, t3 + 1
b) r(t) = et + 1, e2t + 1, et + 1
c) r(t) = hcos 2t, sin 2t, ti
d ) r(t) = sin 2t, 2 sin2 t, 2 cos t
41
8. Considere una partcula de masa constante m con vector posicion r(t), velocidadv(t),
aceleracion a(t) y momento p(t) en el instante t. El momento angular L(t) de
la partcula con respecto al origen en el instante t se define como L(t) =
r(t) p(t). Si F(t) es una fuerza actuando sobre la partcula en el instante
t, se define el torque N(t) con respecto al origen, actuando sobre la partcula
por N(t) = r(t) F(t). Muestre que L (t) = N(t).
dg
dr
dh
d
(g h) + r
h +r g
9. Muestre que (r (g h)) =
.
dt
dt
dt
dt
10. El Teorema del Valor Medio no se satisface en general para funciones vectoriales. Muestre que para r(t) = hcos t, sin t, ti, no existe un t en el intervalo
(0, 2) tal que
r(2) r(0)
.
r (t) =
2 0
11. La curva de Bezier b30 (t) para cuatro puntos no colineales B0 , B1 , B2 y B3
en R3 , con vectores posicion b0 , b1 , b2 , b3 , esta definida por el siguiente
procedimiento:
Se definen las curvas
1) b10 (t) = (1 t)b0 + tb1
Por u
ltimo, a partir de las anteriores se define la curva
b30 (t) = (1 t)b20 (t) + tb21 (t)
a) Muestre que b30 (t) = (1 t)3 b0 + 3t(1 t)2 b1 + 3t2 (1 t)b2 + t3 b3 .
b) Escriba la expresi
on explcita de la curva de Bezier para los puntos B0 =
(0, 0, 0), B1 = (0, 1, 1), B2 = (2, 3, 0), B3 = (4, 5, 2), y representela usando
un programa graficador.
42
3. FUNCIONES VECTORIALES
LECCION
12. Sea r(t) el vector posicion de una partcula que se mueve en R3 . Muestre que
d
(r (v r)) = krk2 a + (r v)v (kvk2 + r a)r.
dt
13. Sea r(t) el vector posicion en R3 de una partcula que se mueve con rapidez
constante c > 0 en una circunferencia de radio a > 0 en el plano XY . Muestre
que a(t) tiene sentido contrario al de r(t) para todo t.
14. Calcule la longitud de arco de la curva parametrizada por r(t) sobre el intervalo
dado:
a) r(t) = h3 cos 2t, 3 sin 2t, 3ti en [0, /2]
c) r(t) = 2 cos 3t, 2 sin 3t, 2t3/2 en [0, 1]
43
20. Si T, N y B son respectivamente los vectores tangente unitario, normal y
binormal, muestre que
dB
B=0
dt
dB
b)
T=0
dt
dB
c)
es m
ultiplo escalar de N
dt
a)
22. Dada la trayectoria r(t) = t4 i + cos t2 j + (t2 1)k determine las ecuaciones
de la recta tangente, el plano normal y el plano osculador, as como el radio y
centro de la circunferencia osculadora, en el punto ( 2 , 1, 1).
23. Sea r(t) una curva diferenciable tal que r(t) 6= 0 para todo t. Muestre que
d
r(t)
r(t) (r (t) r(t))
.
=
dt kr(t)k
kr(t)k3
d
Ayuda: Para calcular dt
(kr(t)k), calcule de dos formas la derivada de kr(t)k2
y luego use el ejercicio 41 de la Lecci
on 1.
24. Sea r(t) una curva suave tal que r (t) 6= 0 para todo t y T es el vector tangente
unitario. Muestre que
T (t) =
44
3. FUNCIONES VECTORIALES
LECCION
25. Suponga que r (t) y r (t) no son paralelos, es decir T (t) 6= 0. Entonces,
existe el vector normal unitario N. Muestre que
N(t) =
r (t) r (t)
.
kr (t) r (t)k
Z bq
(t)2 + ((t)2 sin2 (t)) (t)2 + (t)2 (t)2 dt.
a
LECCION
Las funciones de las que trata esta leccion son funciones definidas en alg
un subn
conjunto D de R y que toman valores en R. M
as precisamente, una funcion de n
variables es una funcion
f : D R,
El dominio de una funcion de varias variables puede ser de diferente naturaleza. Sin
embargo, hay ciertos dominios que juegan un papel muy importante en el estudio
de propiedades de la funcion f . El prop
osito de los proximos parrafos es el de definir
algunas carctersticas de dichos conjuntos en R2 (las definiciones se extienden de
manera natural a R3 ).
Dado un subconjunto D de R2 podemos establecer tres conjuntos a partir de el: el
conjunto de puntos del interior de D, el conjunto de puntos del exterior de D y el
conjunto de puntos de la frontera de D. Decimos que un punto P es exterior a D, si
existe una vecindad con centro en P y radio , V (P, ) = {Q : d(P, Q) < }, tal que
D V (P, ) = . Un punto P es interior de D, si existe una vecindad V (P, ) tal
que V (P, ) D. Un punto P es punto frontera de D si no es ni exterior ni interior
a D, en otras palabras, si para toda vecindad V (P, ), se tiene que D V (P, ) 6=
45
46
z = f (x, y)} .
47
1. todos los puntos (x, y) tales que x2 + y 2 < 1 son puntos interiores del dominio,
2. todos los puntos (x, y) tales que x2 + y 2 > 1 son puntos exteriores del dominio
y
3. todos los puntos (x, y) tales que x2 + y 2 = 1 son puntos frontera del dominio.
Como todos los puntos frontera del dominio de f pertenecen al dominio, se tiene
que este es un conjunto cerrado.
Por otro lado, el dominio de f esta contenido en la vecindad con centro en (0, 0) y
radio 2, por lo tanto es un conjunto acotado. Las curvas de nivel de f ,
p
1 x2 y 2 = k} = {(x, y) : x2 + y 2 = 1 k 2 },
k=0
k = 2/3
k = 5/6
k = 13/14
48
Lmites
Para extender el concepto de lmite, usaremos la idea de los lmites laterales de
funciones de una sola variable. En otras palabras, recordemos que el lmite de una
funcion en un punto existe si los lmites laterales (por la izquierda y por la derecha)
existen y son iguales. Es decir, para estudiar el lmite de la funcion en un punto, nos
acercamos al punto por dos caminos posibles: por la izquierda y por la derecha. Esta
idea es un poco m
as complicada para funciones de dos (y m
as) variables, ya que hay
muchsimas formas de acercarse a un punto en el plano. Sea f : D R y P (a, b)
un punto en el plano, que no necesariamente pertenece a D. Un acercamiento a P
es una funcion vectorial de la forma aP : [1, 0] R2 tal que
i) aP es continua en 0,
ii) aP (t) D, para todo t [1, 0),
iii) P = aP (t) si y s
olo si t = 0.
Observese que un acercamiento no es otra cosa que un camino que conduce a P
usando puntos de D. Ahora, si existe por lo menos un acercamiento a P podemos
definir el lmite de f cuando nos acercamos a P (a, b). Diremos que L es el lmite de
f cuando (x, y) se acerca a (a, b) y escribiremos lm(x,y)(a,b) f (x, y) = L, si para
todo acercamiento aP se tiene que
lm f (aP (t)) = L
t0
t0
t0
3t3
3t
= lm
=0
2
2t
t0 2
Este calculo no demuestra la existencia del lmite. Sin embargo, nos dice que si
existiese el lmite debera ser 0. Veamos que efectivamente existe y es 0. Sea aP (t) =
hx(t), y(t)i cualquier acercamiento a P . Entonces,
3x(t)2 y(t)
x(t)2
=3
0 |f (aP (t))| =
|y(t)| 3 |y(t)|
x(t)2 + y(t)2
x(t)2 + y(t)2
49
Usando el Teorema del Emparedado y el hecho de que lmt0 y(t) = 0, porque aP
es un acercamiento a P (0, 0), se tiene que lmt0 f (aP (t)) = 0 y por consiguiente
lm(x,y)(0,0) f (x, y) = 0.
De acuerdo con la definicion de lmite, este no existir
a si encontramos alg
un acercamiento aP para el que lmt0 f (aP (t)) no exista, o si encontramos por lo menos
dos acercamientos a1P y a2P para los que
lm f (a1P (t)) 6= lm f (a2P (t))
t0
t0
xy
y tomamos
x2 + y 2
los acercamientos a P (0, 0), a1P (t) = ht, 0i y a2P (t) = ht, ti, se tiene que
t0
t0
t0
t0
1
t2
= ,
2
2t
2
(x,y)(a,b)
f (x, y) = a.
En efecto, sea aP (t) = hx(t), y(t)i un acercamiento a P (a, b). Es decir, lmt0 x(t) =
a y lmt0 y(t) = b. Entonces,
lm f (aP (t)) = lm f (x(t), y(t)) = lm x(t) = a
t0
t0
t0
50
(x,y)(a,b)
51
EJERCICIOS
1. Sea f (x, y) = ln(x + y 1)
a) Eval
ue f (1, 1) y f (e, 1).
b) Encuentre el dominio de f
c) Encuentre la imagen de f
d ) Represente gr
aficamente el dominio D de f : es abierto D?, es cerrado
D?, no es ni abierto ni cerrado D?, es acotado D? Por que?
2. Sea f (x, y) = ex
2 y
a) Eval
ue f (2, 4) y f (0, 1).
b) Encuentre el dominio de f
c) Encuentre la imagen de f
d ) Represente gr
aficamente el dominio D de f : es abierto D?, es cerrado
D?, no es ni abierto ni cerrado D?, es acotado D? Por que?
3. Sea g(x, y) =
p
36 9x2 4y 2
a) Eval
ue g(1, 2) y g(0, 1).
b) Encuentre el dominio de g
c) Encuentre la imagen de g
d ) Represente gr
aficamente el dominio D de g: es abierto D?, es cerrado
D?, no es ni abierto ni cerrado D?, es acotado D? Por que?
4. Sea f (x, y, z) = x2 ln(x y + z)
a) Eval
ue f (3, 6, 4)
b) Encuentre el dominio D de f .
c) Encuentre la imagen de f .
d ) Es abierto D?, es cerrado D?, no es ni abierto ni cerrado D?, es
acotado D? Por que?
52
1
5. Sea f (x, y, z) = p
2
2
x + y + z2 1
a) Eval
ue f (1, 3, 4)
b) Encuentre el dominio D de f .
c) Encuentre la imagen de f .
d ) Es abierto D?, es cerrado D?, no es ni abierto ni cerrado D?, es
acotado D? Por que?
6. Encuentre el dominio de f y diga si es cerrado, abierto, ni cerrado ni abierto,
acotado o no acotado, justificando su respuesta:
a) f (x, y) =
b) f (x, y) =
x+y
x+ y
c) f (x, y) = ln(9 x2 9y 2 )
x 3y
d ) f (x, y) =
x + 3y
3x + 5y
e) f (x, y) = 2
x + y2 4
f ) f (x, y) = y x ln(y + x)
p
g) f (x, y) = xy x2 + y
p
h) f (x, y) = x2 + y 2 1 + ln(4 x2 y 2 )
p
i ) f (x, y, z) = 1 x2 y 2 z 2
53
7. De por lo menos dos acercamientos convenientes, calcule los lmites al restringir
la funcion a los acercamientos que di
o, encuentre el lmite de las funciones si
existen, o demuestre que el lmite no existe
a) lm(x,y)(5,2) (x5 + 4x3 y 5xy 2 )
b) lm(x,y)(6,3) xy cos(x 2y)
c) lm(x,y)(0,0)
x2
x2 + y 2
d ) lm(x,y)(0,0)
(x + y)2
x2 + y 2
e) lm(x,y)(0,0)
8x2 y 2
x4 + y 4
f ) lm(x,y)(0,0)
g) lm(x,y)(0,0)
h) lm(x,y)(0,0)
x3 + xy 2
x2 + y 2
xy
p
x2 + y 2
xy + 1
2
x + y2 + 1
i ) lm(x,y)(0,0)
2x2 y
x4 + y 2
j ) lm(x,y)(0,0)
x3 y 2
x2 + y 2
x2 + y 2
k ) lm(x,y)(0,0) p
x2 + y 2 + 1 1
xy 2y
l ) lm(x,y)(0,0) 2
x + y 2 4x + 4
m) lm(x,y,z)(3,0,1) exy sin(z/2)
n) lm(x,y,z)(0,0,0)
n
) lm(x,y,z)(0,0,0)
x2 + 2y 2 + 3z 2
x2 + y 2 + z 2
xy + yz + xz
x2 + y 2 + z 2
1
x2 y
54
xy
+ y2 + 1
c) F (x, y) = arctan(x + y)
b) F (x, y) =
x2
h) f (x, y, z) = x + y + z
2 3
x y
si
i ) f (x, y) = 2x2 + y 2
1
si
xy
2
2
j ) f (x, y) = x + xy + y
0
9. Suponga que
lm
(x,y)(a,b)
f (x, y) y
(x, y) 6= (0, 0)
(x, y) = (0, 0)
si
(x, y) 6= (0, 0)
si
(x, y) = (0, 0)
lm
(x,y)(a,b)
Demuestre que:
[f (x, y) + g(x, y)] =
lm f (x, y) +
lm g(x, y)
(x,y)(a,b)
(x,y)(a,b)
(x,y)(a,b)
(b)
lm k f (x, y) = k
lm f (x, y)
(x,y)(a,b)
(x,y)(a,b)
(c)
lm [f (x, y)g(x, y)] =
lm f (x, y)
lm g(x, y)
(a)
lm
(x,y)(a,b)
(x,y)(a,b)
f (x, y)
(d)
lm
=
(x,y)(a,b) g(x, y)
lm
f (x, y)
lm
g(x, y)
(x,y)(a,b)
(x,y)(a,b)
(x,y)(a,b)
si
lm
(x,y)(a,b)
lm
(x,y)(a,b)
f (x, y) = L.
g(x, y) 6= 0
lm
(x,y)(a,b)
g(x, y) = 0, enton-
55
a)
c)
e)
g)
i)
k)
lm
cos(xy)
b)
lm
x2 y 2
x2 + y 2
d)
(x,y)(0,0)
(x,y)(0,0)
x2 2xy + y 2
xy
lm
(x,y)(1,1)
f)
lm
x2 y 2
xy
h)
lm
y 4 sin(xy)
x2 + y 2
j)
lm
x
y
l)
(x,y)(1,1)
(x,y)(0,0)
(x,y)(0,0)
lm
exy
lm
xy 2
x2 + y 4
lm
xy 2
x2 + y 2
(x,y)(0,0)
(x,y)(0,0)
(x,y)(0,0)
lm
(x,y)(0,0)
lm
(x,y)(0,0)
lm
(x,y)(0,0)
x2 2xy + y 2
xy
1
2
2
(x + y ) cos
xy
1
cos
xy
1
(x +y )/2 , para > 0, es constante en la circun11. Muestre que f (x, y) = 2
2e
ferencia de radio r > 0 centrado en el origen. Esta funcion se llama Gaussian
blur, y se usa como filtro en los programas de procesamiento de imagenes para
producir un efecto de desenfoque.
12. Suponga que f (x, y) f (y, x) para todo (x, y) en R2 . Muestre que f (x, y) =
f (y, x) para todo (x, y) en R2 .
p
13. Use la sustituci
on r = x2 + y 2 para mostrar que
p
sin x2 + y 2
p
=1.
lm
(x,y)(0,0)
x2 + y 2
56
LECCION
Diferenciacion Parcial
Para definir lmite y continuidad de una funcion en un punto P introdujimos el concepto de acercamiento al punto P . Usaremos un procedimiento similar para estudiar
la variaci
on de una funcion en un punto P . Introduciremos el concepto de paso por
el punto P .
Sea f : D R una funcion definida en un subconjunto abierto D del plano y P un
punto de D. Un paso por el punto P (x, y) es una funcion vectorial pvP : (, ) D
de la forma pvP (t) = hx + v1 t, y + v2 ti, donde v = hv1 , v2 i es un vector unitario y
> 0 es tal que pvP (t) D para todo t (, ).
La derivada de f en el punto P (x, y) en la direccion de v que se escribe como
o como Dv f (x, y) es:
d
d
df
v
(x, y) = f (pP (t)) = f (x + v1 t, y + v2 t)
dv
dt t=0
dt t=0
df
(x, y)
dv
58
5. DIFERENCIACION
PARCIAL
LECCION
df
dg
d(f + g)
(x, y) =
(x, y) +
(x, y). Adem
as,
cursos anteriores, como por ejemplo
dv
dv
dv
entre otras propiedades, se satisface la regla del producto para derivadas.
Las derivadas de f en P (x, y) en la direccion de los vectores i y j se conocen como
derivadas parciales de f en P (x, y) con respecto a x y a y respectivamente y se suele
escribir:
df
f
df
f
(x, y) =
(x, y),
(x, y) =
(x, y)
di
x
dj
y
Como veremos m
as adelante todas las derivadas direccionales se van a poder expresar
en terminos de las derivadas parciales y de aqu su importancia.
Veamos un ejemplo. Sea f : R2 R definida por f (x, y) = x3 2y 2 y calculemos
df
2/2, 2/2 , v =i y v =j . Tenemos
la derivada
(x, y) en los casos en que v =
dv
que:
2/2, 2/2
En el caso en que v =
d
df
2 3
(x
+
(x, y) =
2 t) 2(y +
dv
dt t=0
=
= 3
2
2 (x
2 2
2 t)
2
2 (y
2 2y,
2 2
2 x
en el caso en que v
=i
df
d
(x + t)3 2y 2
(x, y) =
di
dt t=0
=
3(x + t)2 t=0
= 3x2 =
f
(x, y),
x
y
en el caso en que v =j
2 2
2 t)
2
2 t)t=0
59
df
(x, y) =
dj
=
d
x3 2(y + t)2
dt t=0
4(y + t)|t=0
= 4y =
f
(x, y).
y
y = b + v2 t,
z = f (a + v1 t, b + v2 t),
t (, ).
df
(a, b) ,
v1 , v2 ,
dv
que tambien es tangente a la gr
afica de la funcion f en el punto Q(a, b, f (a, b)). La
df
pendiente de este vector con respecto al plano XY es precisamente
(a, b).
dv
e=
El vector tangente a esta curva en el punto Q(a, b, f (a, b)) es v
Dentro de estos vectores tangentes hay dos que nos interesan que son ei = h1, 0, fx (a, b)i
y ej = h0, 1, fy (a, b)i que nos dan las pendientes en las direcciones i y j de la gr
afica de
f en el punto Q(a, b, f (a, b)), respectivamente. A partir de estos vectores podemos
encontrar la ecuaci
on del plano tangente a la g
arfica de f en el punto Q(a, b, f (a, b)).
En efecto, el vector ei ej = hfx (a, b), fy (a, b), 1i es un vector normal al plano tangente, por lo que dicha ecuaci
on sera
fx (a, b)(x a) fy (a, b)(y b) + (z f (a, b)) = 0
o, lo que es equivalente
z = f (a, b) + fx (a, b)(x a) + fy (a, b)(y b).
60
5. DIFERENCIACION
PARCIAL
LECCION
Esta u
ltima ecuaci
on corresponde a la ecuaci
on de la funcion afn
(a,b)
Lf
que aproxima a f en una vecindad del punto Q. Esta funcion se usa frecuentemente
para aproximar valores de f al hacer incrementos peque
nos x y y en a y en b
(a,b)
respectivamente. Es decir, Lf (a + x, b + y) se toma como valor aproximado de
f (a+x, b+y). En general, al hacer esto se comete un error E = E(f, a, b, x, y)
que no es otra cosa que la diferencia de los dos anteriores. Es decir,
(a,b)
E = f (a + x, b + y) Lf (a + x, b + y)
= f (a + x, b + y) f (a, b) fx (a, b)(x) fy (a, b)(y).
Este error es peque
no cuando los incrementos x y y son peque
nos, pero se requiere, para que la aproximacion sea de utilidad, que este sea mucho m
as peque
no
en el siguiente sentido. Diremos que E es bien peque
no si existe h1 , 2 i tal que
1. E = hx, yi h1 , 2 i.
2. Si hx, yi h0, 0i se tiene que h1 , 2 i h0, 0i.
No todas las funciones se pueden aproximar mediante su funcion afn en (a, b) con
un error bien peque
no. Las que satisfacen esta condici
on se llaman funciones diferenciables en (a, b). En otras palabras, la funcion f es diferenciable en (a, b) si E
es bien peque
no. En general, no es facil encontrar h1 , 2 i, por lo que verificar si
una funcion es diferenciable en (a, b) se vuelve una tarea dispendiosa. Sin embargo,
contamos con el siguiente resultado que usaremos sin demostraci
on(ver el final de la
leccion):
61
buena aproximacion de f = f (a + x, b + y) f (a, b) si f es diferenciable en
(a, b). En efecto,
(a,b)
E = f (a + x, b + y) Lf (a + x, b + y)
= f (a + x, b + y) f (a, b) fx (a, b)x fy (a, b)y
= f df.
Es decir, f = df + E.
Como ilustraci
on del uso de las diferenciales, estimemos el error m
aximo al calcular
el volumen de un cono cuyas dimensiones, radio r de 10cm y base h de 25 cm, se han
obtenido con un posible error de medici
on de 0.1cm. El error cometido en la medici
on
2
del volumen V = r h/3 es V , que por lo que hemos visto es aproximadamente
igual a
dV = Vr r + Vh h = (2rh/3)r + (r2 /3)h.
Como |r| 0,1 y |h| 0,1. Entonces,
500
100
(0,1) +
(0,1) = 20
|dV | = (2rh/3)r + (r2 /3)h
3
3
Es decir, el m
aximo error cometido al calcular el volumen a partir de las mediciones
es de aproximadamente 63 cm3 .
Veamos ahora que efectivamente:
Para esto, demostraremos que existe h1 , 2 i tal que E = h1 , 2 i hx, yi y tal que
62
5. DIFERENCIACION
PARCIAL
LECCION
63
de modo que el rectangulo de vertices (a, b), (a, b + k), (a + h, b), (a + h, b + k) este
contenido en D. Consideramos las funciones F (x) = f (x, b + k) f (x, b) definida en
el intervalo [a, a + h] y L(y) = f (a + h, y) f (a, y) definida en el intervalo [b, b + k].
Claramente, F (a + h) F (a) = L(b + k) L(b), valor com
un que denotaremos .
Adem
as, F (x) = fx (x, b + k) fx (x, b) y L (y) = fy (a + h, y) fy (a, y). Por el
Teorema del Valor Medio, existen un puntos 1 en el intervalo abierto (a, a + h) y
2 en el intervalo abierto (b, b + k) tales que
F (a + h) F (a) = F (1 )h = (fx (1 , b + k) fx (1 , b))h
y
L(b + k) L(b) = L (2 )k = (fy (a + h, 2 ) fy (a, 2 ))k.
Consideremos ahora las funcion G(y) = fx (1 , y) definida en el intervalo [b, b + k] y
M (x) = fy (x, 2 ) definida el intervalo [a, a + h]. Por el Teorema del Valor Medio,
existen puntos 1 en el intervalo abierto (b, b+k) y 2 en el intervalo abierto (a, a+h)
tales que
G(b + k) G(b) = G (1 )k = fxy (1 , 1 )k.
y
M (a + h) M (a) = M (2 )h = fyx (2 , 2 )h
Por lo anterior y teniendo en cuenta que = F (a + h) F (a) = L(b + k) L(b) =
(G(b + k) G(b))h = (M (a + h) M (a))k, se tiene que
= fxy (1 , 1 )kh = fyx (2 , 2 )hk,
y, por consiguiente,
fxy (1 , 1 ) = fyx (2 , 2 ).
Teniendo en cuenta que fxy y fyx son continuas y que (1 , 1 ) y (2 , 2 ) tienden a
(a, b) cuando (h, k) tiende a (0, 0), se concluye que fxy (a, b) = fyx (a, b), que es lo
que queramos mostrar.
64
5. DIFERENCIACION
PARCIAL
LECCION
EJERCICIOS
1. Calcule las derivadas direccionales de las siguientes funciones en el punto dado
y la direccion dada. No olvide que la direccion la da un vector unitario.
a) f (x, y) = x2 + y 2 en el punto (1, 1) en la direccion del vector v= h1, 1i
b) f (x, y) = cos(x + y) en el punto (0, ) en la direccion del vector v= h2, 1i
p
c) f (x, y) = x2 + y + 4 en el punto (1, 0) en la direccion del vector v=
h1, 1i
x+1
d ) f (x, y) =
en el punto (1, 5) en la direccion del vector v= h1, 2i
y+1
2. Calcule las derivadas parciales
f
f
y
de las siguientes funciones:
x y
a) f (x, y) = exy + xy
b) f (x, y) = x2 y 2 + 6xy + 4x 8y + 2
c) f (x, y) = x4
d ) f (x, y) = x + 2y
e) f (x, y) =
g) f (x, y) =
p
x2 + y 2
f ) f (x, y) = sin(x + y)
p
3
x2 + y + 4
i ) f (x, y) = e(x
h) f (x, y) =
2 +y 2 )
j ) f (x, y) = ln(xy)
k ) f (x, y) = sin(xy)
3. Calcule las derivadas parciales
a) f (x, y, z) = x
xy + 1
x+y
l ) f (x, y) = tan(x + y)
f f
f
,
y
de las siguientes funciones:
x y
z
p
y2 + z2
2 +y 3 )/z
4. Encuentre la ecuaci
on del plano tangente a la superficie de ecuaci
on z = f (x, y)
en el punto P :
65
a) f (x, y) = x2 + y 3 , P = (1, 1, 2)
c) f (x, y) = x2 y, P = (1, 1, 1)
f ) f (x, y) =
x2 + y 2 , P = (3, 4, 5)
66
5. DIFERENCIACION
PARCIAL
LECCION
x+1
.
y+1
2
3
3x y y
x2 + y
f (x, y) =
si
(x, y) 6= (0, 0)
si
(x, y) = (0, 0)
tiene derivadas direccionales en todas las direcciones en (0, 0), pero no es diferenciable en (0, 0).
12. Sea f definida por
f (x, y) =
x2
xy
+ y2
si
(x, y) 6= (0, 0)
si
(x, y) = (0, 0)
2 2
sin(x y )
x2 + y 2
f (x, y) =
si
(x, y) 6= (0, 0)
si
67
a) Explique lo que significara que f sea diferenciable en (0, 0).
b) Muestre que las derivadas parciales fx (0,0) y fy (0,0) existen y calculelas.
c) z Es f diferenciable en (0, 0)?
14. Considere la funcion f definida por
sin(xyz)
f (x, y, z) = x4 + y 4 + z 4
si
(x, y, z) 6= (0, 0, 0)
si
68
5. DIFERENCIACION
PARCIAL
LECCION
d ) El
area de un lote triangular es (1/2)ab sin C, donde a y b son las longitudes de dos de sus lados y C es la medida del angulo formado por
esos dos lados. Al efectuar medidas topogr
aficas sobre el terreno del lote
triangular, se han obtenido medidas de a, b y C de 150 m, 200 m, y 60
respectivamente. Estime el error en el calculo del area del lote si la medida a, b de sus lados tienen errores de 15 cm cada una y C un error de
2 .
16. Calcuar todas las derivadas parciales de segundo orden de las siguientes funciones.
a) f (x, y) = sin(x2 3xy)
b) f (x, y) = x2 y 2 e2xy
c) f (x, y) = arctan xy
17. Una funcion u = f (x, y) es una funcion arm
onica si sus derivadas parciales
2
2
segundas existen, son continuas y satisfacen la ecuaci
on de Laplace xu2 + yu2 =
0. Muestre que la funcion f (x, y) = x3 3xy 2 es armonica.
xy(x2 y 2 )
x2 +y 2
f
x
f
y
si (x, y) 6= (0, 0)
b) Verificar que
f
x (0, 0)
c) Verifique que
2f
xy (0, 0)
f
y (0, 0)
= 0 (usando la definicion).
= 1 y que
2f
yx (0, 0)
= 1.
2f 2f
,
x2 y 2
2f
y x :
a) f (x, y) = x2 + y 2
c) f (x, y) =
p
x2 + y + 4
e) f (x, y) = exy + xy
b) f (x, y) = cos(x + y)
d ) f (x, y) =
x+1
y+1
f ) f (x, y) = x2 y 2 + 6xy + 4x 8y + 2
69
g) f (x, y) = x4
h) f (x, y) = x + 2y
i ) f (x, y) = ln(xy)
j ) f (x, y) = sin(xy)
=0.
x2
y 2
21. Sean u y v funciones dos veces diferenciables de una sola variable, y sea c 6= 0
una constanteee. Muestre que f (x, y) = u(x + cy) + v(x cy) es una solucion
de la ecuaci
on general de onda unidimensional.
2f
1 2f
=0.
x2 c2 y 2
70
5. DIFERENCIACION
PARCIAL
LECCION
LECCION
Regla de la cadena
72
6. REGLA DE LA CADENA
LECCION
z(t + t) z(t)
t
= fx (x(t), y(t))
y
E
x
+ fy (x(t), y(t))
+
,
t
t
t
dx
dy
dz
= fx
+ fy ,
dt
dt
dt
que es la Regla de la Cadena en esta situacion particular.
Si definimos el vector gradiente de f para cada (x, y) D como f (x, y) =
hfx (x, y), fy (x, y)i, la regla de la cadena se puede expresar de la siguiente forma
dz
dx dy
= hfx , fy i
,
dt
dt dt
= f r
o, m
as explcitamente,
dz
(t) = f (r(t)) r (t).
dt
Teniendo en cuenta la regla de la cadena, podemos reformular la derivada direccional
de f en la direccion de un vector v as:
df
d
f (x + v1 t, y + v2 t)
(x, y) =
dv
dt t=0
=
f (x + v1 t, y + v2 t) hv1 , v2 i|t=0
= f (x, y) v
73
A partir de esta u
ltima, se deducen varias propiedades del vector gradiente. En
efecto,
df
(x, y) = f (x, y) v
dv
= kf (x, y)k kvk cos
= kf (x, y)k cos ,
siendo el
angulo entre f (x, y) y v. De aqu que, la direccion v en la que la derivada direccional es m
as grande, es aquella en la que = 0, ya que en este valor
df
del angulo, cos toma su m
aximo valor, que es 1. Es decir,
(x, y) es m
axima
dv
en la direccion y sentido de f (x, y), que es en la que f (x, y) y v son paraledf
(x, y) = kf (x, y)k. De igual manera, en el sentido opuesto
los. En este caso,
dv
al del gradiente se obtiene el m
aximo decrecimiento de la funcion en el punto y
df
(x, y) = kf (x, y)k.
dv
Por otro lado, si r = r(t) es una parametrizacion de la curva de nivel Nkf , se tiene
que k = f (r(t)) y por consiguiente, 0 = f (r(t)) r (t), lo que quiere decir que
f (r(t)) es perpendicular a la curva de nivel Nkf en el punto (x(t), y(t)).
Ahora bien, todo lo que hemos hecho con funciones de dos variable hasta el momento,
se puede generalizar a funciones de tres y m
as variables. Por ejemplo, mediante un
razonamiento muy similar al u
ltimo aqu presentado, se puede concluir que para una
funcion f : R3 R, el gradiente f (x0 , y0 , z0 ) es perpendicular a la superficie de
nivel Nkf en el punto (x0 , y0 , z0 ) Nkf . Este u
ltimo resultado permite encontrar la
f
ecuaci
on del plano tangente a Nk en el punto (x0 , y0 , z0 ) Nkf , ya que un vector
normal al plano tangente en ese punto es f (x0 , y0 , z0 ). Explcitamente, la ecuaci
on
de dicho plano es
fx (x0 , y0 , z0 )(x x0 ) + fy (x0 , y0 , z0 )(y y0 ) + fz (x0 , y0 , z0 )(z z0 ) = 0
Supongase ahora que f es una funcion diferenciable de dos variables x, y, que a su vez
dependen de dos variables independientes s y t, es decir, z(s, t) = f (x(s, t), y(s, t)).
Suponiendo que x y y tienen derivadas parciales con respecto a s y a t, queremos
calcular las derivadas parciales con respecto a estas u
ltimas variables. Al hacer un
74
6. REGLA DE LA CADENA
LECCION
f
(x(s + s, t) x(s, t)) f
(y(s + s, t) y(s, t))
E
(x, y)
+
(x, y)
+
x
s
x
s
s
+ h
que tiende a
f
x f
y
(x, y)
+
(x, y) , cuando s tiende a cero. As,
x
s
y
s
z
f x f y
=
+
.
s
x s
y s
f dx f dy
+
x dt
y dt
f
f
+
g (t)
x y
f
(t, g(t))
g (t) = x
f
(t, g(t))
y
75
siempre y cuando
f
(t, g(t)) 6= 0.
y
76
6. REGLA DE LA CADENA
LECCION
EJERCICIOS
1. Suponga que el movimiento de un objeto en el espacio esta dado por r(t) =
cos ti + sin tj + tk y que la temperatura en un punto cualquiera del espacio
esta dada por T (x, y, z) = x2 + y 2 + z 2 2x y.
a) Determine si existe un punto del recorrido del objeto en el cual la temperatura es 0.
b) Determine la razon de cambio de la temperatura del objeto en t =
7
3 .
77
a) En que direccion debe moverse para que la temperatura baje lo m
as
r
apidamente posible?
b) Si la nave se mueve con una rapidez de e6 metros por segundo, a que velocidad bajar
a la temperatura cuando se desplace en esa direccion?
c) Desafortunadamente, el metal del casco se fracturar si se entra a una
+3z
y la
4. Comprobar la regla de la cadena para la funcion f (x, y, z) = ln 1+2x
1+y 2
2
trayectoria r(t) = 1 ti + t j + cos tk.
5. Encuentre la ecuaci
on del plano tangente a la superficie z = f (x, y) en el punto
P:
a) f (x, y) = x2 + y 3 , P (1, 1, 2)
c) f (x, y) = x2 y, P (1, 1, 1)
f ) f (x, y) =
x2 + y 2 , P (3, 4, 5)
6. Encuentre la ecuaci
on del plano tangente a la superficie dada en el punto P :
a)
x2
4
y2
9
z2
16
b) x2 + y 2 + z 2 = 9, P (0, 0, 3)
= 1, P 1, 2, 2 311
c) x2 + y 2 z 2 = 0, P (3, 4, 5)
d ) x2 + y 2 = 4, P ( 3, 1, 0)
7. Calcule el gradiente de f , f :
a) f (x, y) = x2 + y 2 1
c) f (x, y) =
p
x2 + y 2 + 4
e) f (x, y) = ln(xy)
b) f (x, y) =
1
x2 + y 2
d ) f (x, y) = x2 ey
f ) f (x, y) = 2x + 5y
78
6. REGLA DE LA CADENA
LECCION
g) f (x, y, z) = sin(xyz)
h) f (x, y, z) = x2 eyz
i ) f (x, y, z) = x2 + y 2 + z 2
j ) f (x, y, z) =
p
x2 + y 2 + z 2
8. Calcule
la
derivada direccional de f en el punto P en la direccion de v =
1 , 1
.
2
2
a) f (x, y) = x2 + y 2 1, P (1, 1)
c) f (x, y) =
p
x2 + y 2 + 4, P (1, 1)
b) f (x, y) =
x2
1
, P (1, 1)
+ y2
d ) f (x, y) = x2 ey , P (1, 1)
9.
Calcule la derivada
direccional de f en el punto P en la direcc
on de v =
1
1
1
, ,
.
3
3
3
a) f (x, y, z) = sin(xyz), P (1, 1, 1)
b) (f + g) = f + g
c) (f g) = f g + g f
d ) (f /g) =
e)
df
df
=
d(v)
dv
f)
g f f g
, g(x, y) 6= 0
g2
df
d(cf )
=c
dv
dv
79
df
dg
d(f g)
dg
df
d(f + g)
=
+
h)
=f
+g
dv
dv dv
dv
dv
dv
p
13. La funcion r(x, y) = x2 + y 2 es la longitud del vector posicion r = x i + y j
1
para cada punto (x, y) en R2 . muestre que r = r cuando (x, y) 6= (0, 0), y
r
que (r2 ) = 2 r.
g)
xy
si
(x, y) 6= (0, 0)
f (x, y) = x2 + y 2
0
si
(x, y) = (0, 0).
a) Muestre que fx y fy existen en (0, 0).
+v
15. Para la funcion z(x, y) = f (u(x, y), v(x, y)), con f (u, v) = uu2 v
2 , u(x, y) =
xy
xy
e
y v(x, y) = e . Calcule cada una de las siguientes derivadas usando la
regla de la cadena.
a)
b)
z
x
z
y
c)
d)
2z
xy
2z
y 2
16. Para la funcion z(x, y) = f (u(x, y), v(x, y)). Determine cada una de las siguientes derivadas usando la regla de la cadena.
a)
b)
z
x
z
y
c)
d)
2z
xy
2z
y 2
17. Dada la funcion z(u, v) = f (x(u, v), y(u, v)), con f (x, y) = x2 + 2xy + 3y 3 ,
z z
u
y
cuando
x(u, v) = y y(u, v) = ln(u v). Determine los valores de
v
u v
a) u = 4 y v = 2.
b) x = 2 y y = 0.
80
6. REGLA DE LA CADENA
LECCION
18. Dada la funcion z(u, v) = f (x(u, v), y(u, v)), con f (x, y) = arctan(
x(u, v) = uv y y(u, v) = eu/v . Determine los valores de
xy
), con
x+y
z z
y
cuando
u v
a) u = 4 y v = 2.
b) x = 4 y y = e
19. Las siguientes ecuaciones definen implcitamente a y como funcion de x. Calcule dy/dx en el punto dado.
a) 2x3 3y 2 + xy = 0 en el punto (1, 1).
81
a) Calcule w/r y w/ en terminos de f /x y f /y.
b) Escriba f /x y f /y en terminos de w/r y w/.
c) Use los items anteriores para demostrar que
f
x
2
f
y
2
w
r
2
1
r2
2
INTODUCIR EL LAPLACIANO
23. Encuentre el Laplaciano de la funcion f en coordenadas cartesianas:
a) f (x, y, z) = x + y + z b) f (x, y, z) = x5
d ) f (x, y, z) = ex+y+z
c) f (x, y, z) = (x2 + y 2 +
z 2 )3/2
2 y 2 z 2
82
6. REGLA DE LA CADENA
LECCION
LECCION
Extremos de funciones
84
7. EXTREMOS DE FUNCIONES
LECCION
En algunos contextos se usan las palabras relativo en lugar de local, global en lugar de absoluto y a un valor m
aximo o un valor mnimo se le denomina valor extremo.
No es difcil demostrar que si f tiene un valor m
aximo local en un punto P (a, b)
interior de D y f tiene derivada en el punto P (a, b) en la direccion de v, entonces
df
(a, b) = 0. En efecto, sea p : (, ) D un paso por P (a, b) en la direccion de
dv
v, como f tiene un m
aximo local en P (a, b) entonces la funcion f p, que es de una
sola variable, tiene un m
aximo local en 0. As,
df
d
(a, b) = f (p(t)) = 0.
dv
dt t=0
Recordemos que si f es diferenciable en P (a, b) entonces f tiene derivada en el punto
df
P (a, b) en la direccion de v, para todo v y
(a, b) = f (a, b) v, por lo tanto podv
demos concluir que si f tiene un m
aximo local en P (a, b) entonces f (a, b) v = 0,
para todo v, por lo tanto f (a, b) = hfx (a, b), fy (a, b)i = h0, 0i.
Un an
alisis completamente an
alogo al anterior se puede hacer para el caso en que f
tenga un mnimo local en (a, b).
N
otese que si f es una funcion de una sola variable que tiene un m
aximo local o
un mnimo local en un punto a y f es diferenciable en a, entonces f (a) = 0. Por
lo tanto, la recta tangente a la gr
afica en el punto (a, f (a)) es paralela al eje X y
tiene ecuaci
on y = f (a). De igual manera, si f es una funcion de dos variables que
tiene un valor m
aximo local o un mnimo local en (a, b) y f es diferenciable en (a, b)
entonces f (a, b) = h0, 0i. Por lo tanto, la ecuaci
on del plano tangente a la gr
afica
de f en el punto (a, b, f (a, b)) es z = f (a, b), que es la ecuaci
on de un plano
paralelo al plano XY .
Un punto P (a, b) D es un punto crtico de f si f es diferenciable en P (a, b) y
fx (a, b) = 0 = fy (a, b), o si f no es diferenciable en P (a, b).
Claramente si se quieren encontrar los puntos en los que una funcion alcanza sus
valores m
aximos y mnimos relativos, se debe buscar entre los puntos crticos. Una
vez identificados estos puntos, viene la tarea de su clasificacion. Es decir, determinar
85
en cu
ales de estos puntos f tiene un valor mnimo local, en cu
ales f tiene un valor
m
aximo local, y en cu
ales no tiene m
aximo ni mnimo local.
A continuaci
on se har
a la deduccion del criterio de clasificacion de puntos crticos
para funciones de dos variables, herramienta muy u
til en muchos casos en los que la
funcion f tiene segundas derivadas parciales continuas en el punto crtico (a, b).
Supongase que (a, b) es un punto crtico de la funcion f en el que f es diferenciable
y tiene segundas derivadas parciales continuas. Sea p un paso por el punto (a, b) en
la direccion de v = hv1 , v2 i y ponga g(t) = f (p(t)). Aplicando la regla de la cadena,
se tiene
g (t) = f (p(t)) v
= v1 fx (p(t)) + v2 fy (p(t)),
y adem
as se tiene que g (0) = f (p(0)) v = hfx (a, b), fy (a, b)i v = 0, por lo tanto
0 es un punto crtico de g. Aplicando nuevamente la regla de la cadena se obtiene:
g (t) =
d
{v1 fx (p(t)) + v2 fy (p(t))}
dt
= v1 {v1 fxx (p(t)) + v2 fxy (p(t))} + v2 {v1 fxy (p(t)) + v2 fyy (p(t))}
= v12 fxx (p(t)) + 2v1 v2 fxy (p(t)) + v22 fyy (p(t)).
Si fxx (a, b) 6= 0, podemos multiplicar y dividir por fxx (a, b) la u
ltima expresi
on y
completar el cuadrado, obteniendo
2 (p(t)) v 2
{v1 fxx (p(t)) + v2 fxy (p(t))}2 + fyy (p(t))fxx (p(t)) fxy
2
,
g (t) =
fxx (p(t))
y evaluando en t = 0, podemos escribir:
2 (a, b) v 2
{v1 fxx (a, b) + v2 fxy (a, b)}2 + fyy (a, b)fxx (a, b) fxy
2
g (0) =
.
fxx (a, b)
Haciendo
(a, b) = fxx (a, b)fyy (a, b) (fxy (a, b))2 ,
de la u
ltima expresi
on se deduce que si > 0 y fxx (a, b) > 0 entonces g (0) > 0,
independientemente del vector v, y que si > 0 y fxx (a, b) < 0 entonces g (0) < 0,
86
7. EXTREMOS DE FUNCIONES
LECCION
87
1. Encuentre los puntos crticos de f en D,
2. Encuentre los puntos de la frontera de D en los que f tiene valores extremos,
3. Compare los valores de f en los puntos anteriores y escoja aquel en el que f
alcanza el mayor valor y aquel en el que alcanza el menor valor.
En general, la identificaci
on de los puntos de la frontera en los que la funcion tiene
valores extremos no es tarea facil. Sin embargo, contamos con la tecnica de los
multiplicadores de Lagrange que en muchos casos es de gran ayuda, sobre todo
cuando la frontera tiene una ecuaci
on de la forma g(x, y) = k, es decir cuando esta
frontera es una curva de nivel Nkg de una cierta funcion g.
Supongamos que la funcion f : D R, definida en el conjunto cerrado D, tiene un
valor extremo en el punto P (x0 , y0 ) de la frontera de D. Supongamos tambien que
la frontera tiene una ecuaci
on de la forma g(x, y) = k. Si parametrizamos la frontera
mediante una funcion vectorial r : (, ) D de tal forma que r(0) = hx0 , y0 i, se
tiene, por un lado, que 0 es un valor extremo de la funcion z = f r y, por otro, que
g(r(t)) = k para todo t (, ). Entonces, usando la regla de la cadena, se tiene
que
dz
1. 0 =
= f (x0 , y0 ) r (0) y
dt t=0
2. 0 = g(x0 , y0 ) r (0).
Esto quiere decir que en los puntos (x0 , y0 ) de la frontera en los que f tiene un valor
extremo local, los gradientes f (x0 , y0 ) y g(x0 , y0 ) son paralelos, debido a que
ambos son perpendiculares a r (0). En resumen, si f tiene un valor extremo local
en el punto (x0 , y0 ) de la frontera g(x, y) = k de D, existe R (multiplicador de
Lagrange) tal que:
f (x0 , y0 ) = g(x0 , y0 ).
En otras palabras, si estamos interesados en encontrar los puntos de la frontera en
los que f tiene un valor extremo local, debemos encontrar soluciones de la ecuaci
on
anterior.
Lo que se ha hecho en el contexto de funciones de dos variables se puede extender
de manera natural a funciones de m
as variables.
88
7. EXTREMOS DE FUNCIONES
LECCION
La identificaci
on de puntos crticos en la frontera del dominio de una funcion es un
problema de optimizaci
on con restricciones. En otras palabras, lo que hemos hecho
es restringir una funcion f al conjunto R = {(x, y); g(x, y) = k} y encontrar puntos
crticos en R . Esta situacion se puede generalizar, por ejemplo a funciones de tres
variables y con dos restricciones. Explcitamente, supongamos que hemos restringido
la funcion al conjunto de puntos
R = {(x, y, z) : g1 (x, y, z) = k1
g2 (x, y, z) = k2 } ,
89
EJERCICIOS
1. Para cada una de las siguientes funciones determine los puntos crticos y clasifquelos
a) f (x, y) = x3 3x + y 2
b) f (x, y) = x3 12x + y 2 + 8y
c) f (x, y) = x3 3x + y 3 3y
d ) f (x, y) = x3 + 3x2 + y 3 3y 2
g) f (x, y) =
p
x2 + y 2
i ) f (x, y) = x + 2y
b) Determine el m
aximo y el mnimo de f en B = (x, y, z) : x2 + y 2 + z 2 = 1
5. Para un s
olido rectangular de 1000 metros c
ubicos de volumen, encuentre las
dimensiones que maximizan el area de su superficie.
90
7. EXTREMOS DE FUNCIONES
LECCION
and
f2 (x, y) = 0 ,
yn+1 = yn +
f1 (xn , yn )
f1
(xn , yn )
x
f2 (xn , yn )
f2
(xn , yn )
x
D(xn , yn )
f2
f1
f2
1
donde D(xn , yn ) = f
x (xn , yn ) y (xn , yn ) y (xn , yn ) x (xn , yn ) . Entonces
la sucesion de puntos (xn , yn )
on. Escriba un programa
n=1 converge a una soluci
de computador que use este algoritmo para encontrar soluciones aproximadas
del sistema de ecuaciones:
sin(xy) x y = 0
e2x 2x + 3y = 0 .
91
14. Un alambre de longitud A se corta en dos pedazos. Uno se dobla en forma
de cuadrado y el otro en forma de crculo, encuentre las dimensiones de los
dos pedazos de manera que la suma de las areas del cuadrado y el crculo sea
m
axima y las dimensiones de los dos pedazos de manera que la suma de las
areas sea mnima.
15. Resuelva el ejercicio anterior en el caso en que los alambres se doblan en forma
de tri
angulo equil
atero y cuadrado.
16. Resuelva el ejercicio anterior en el caso en que los alambres se doblan en forma
de tri
angulo rectangulo isosceles y cuadrado.
17. Una caja rectangular sin tapa tiene base cuadrada, encontrar las dimensiones
de la caja de volumen constante y de mnima area superficial.
18. La longitud de generatriz de un cono circular recto es a, determine la medida
del
angulo formado por la generatriz y el eje de simetra del cono de volumen
m
aximo.
19. Muestre que el producto de los senos de los angulos de un tri
angulo es m
aximo
cuando el tri
angulo es equil
atero.
20. Para las siguientes funciones determine todos los puntos crticos y clasifquelos.
3x4 4x3 12x2 + 18
12(1 + 4y 2 )
b) f (x, y) = (2 + cos x) sen y
a) f (x, y) =
c) f (x, y) = y sen x
21. Determine los m
aximos y mnimos absolutos de la funci
on f (x, y) = 2x2 4x+
y 2 4y + 1 en la regi
on triangular de vertices (0, 0), (0, 2) y (1, 2)
22. Determine la distancia m
as corta del punto (a, 0), a 6= 0, a la curva de ecuaci
on
2
x=y .
23. Determine el punto de la curva de interseccion del plano 2y + 4z = 5 y el cono
z 2 = 4x2 + 4y 2 m
as cercano al origen.
24. Descomponer el n
umero 97 en 4 sumandos de manera que la suma de sus
cuadrados sea mnima.
92
7. EXTREMOS DE FUNCIONES
LECCION
25. Descomponer el n
umero positivo p en n factores positivos, de manera que la
suma de sus inversos sea mnima.
LECCION
Integrales multiples
94
8. INTEGRALES MULTIPLES
LECCION
f (xk )x.
k=1
R
Rb
La integral de f sobre el intervalo I = [a, b], que se escribe I f (x)dx o a f (x)dx,
es el lmite de estas sumas de Riemann cuando n tiende a infinito,
Z
b
a
f (x)dx = lm R(f, n, ) = lm
n
n
X
f (xk )x.
k=1
c y d} .
95
una de las mn 2-celdas [xk , xk+1 ] [yj , yj+1 ], y definimos la suma de Riemann de
una funcion continua f : R R2 R, correspondiente a esta particion y esta
escogencia de punto de muestra as:
R(f, n, m, )
)A + f (x , y )A + + f (x , y )A
= f (x11 , y11
21 21
n1 n1
)A + f (x , y )A + + f (x , y )A
+ f (x12 , y12
n2 n2
22 22
...
)A + f (x , y )A + + f (x , y )A
+ f (x1m , y1m
nm nm
2m 2m
m
n X
X
f (xkj , ykj
)A
k=1 j=1
R
La integral de f sobre R que se escribe R f (x, y)dA es el lmite de estas sumas de
Riemann cuando n y m tienden a infinito:
Z
f (x, y)dA = lmn,m R(f, n, m, )
R
= lmn,m
m
n X
X
f (xkj , ykj
)A.
k=1 j=1
= {(x, y, z) : a x b,
= [a, b] [c, d] [e, h].
c y d,
e z h}
96
8. INTEGRALES MULTIPLES
LECCION
l
m X
n X
X
f (xkji , ykji
zkji
)V.
R
La integral de f sobre P que se escribe P f (x, y, z)dV es el lmite de estas sumas
de Riemann cuando n, m y l tienden a infinito:
Z
f (x, y, z)dV =
P
= lmn,m,l R(f, n, m, l, )
= lmn,m,l
l
m X
n X
X
, zkji
)V
f (xkji , ykji
Esta
se conoce como integral iterada. El Teorema de Fubini establece que la integral
de f sobre R es igual a esta integral iterada. Explcitamente, si f es continua en R
97
se tiene que
R b R d
f
(x,
y)dy
dx
a
c
R d R b
= c
a f (x, y)dx dy
R f (x, y)dA =
los puntos suspensivos son para las otras tres integrales iteradas que faltan (c
uales?).
Para dar cuenta de integrales de funciones definidas en regiones m
as generales que
las celdas, recurriremos a una estrategia de parametrizacion, concepto que no es muy
utilizado en los libros de calculo pero que mata varios p
ajaros de un solo tiro.
Sea f : D R una funcion continua, definida en una regi
on D del plano. Una
parametrizacion de D es una funcion vectorial
r : [a, b] [c, d] R2 ,
98
8. INTEGRALES MULTIPLES
LECCION
= (1 v) u, u2 + v u, u2 + 2u + 4
= (1 v)u + vu, u2 + v(u2 + 2u + 4)
P m Pn
i=1
rv k uv,
99
donde ru y rv son los vectores de derivadas parciales de x = x(u, v) y y = y(u, v)
con respecto a u y a v en Pij respectivamente. La aparicion de ku ru v rv k en
la suma de Riemann se debe a que esta es el area del paralelogramo generado por los
vectores u ru (Pij ) y v rv (Pij ) , que es una aproximacion del area de la subregi
on
Dij .
La integral de f sobre D es el lmite de R(f, m, n, r) cuando m y n tienden a infinito:
ZZ
n
m X
X
f dA = lmm,n
f (Pi,j ) kru rv k uv
D
i=1 j=1
ZZ
[ab][cd]
[a,b][c,d]
Es decir, las integrales sobre rectangulos son un caso particular de integrales sobre
regiones m
as generales.
Para ilustrar lo anterior, calculemos la integral de f (x, y) = x2 + y 2 sobre el dis
co D = (x, y) : x2 + y 2 4 con respecto a la parametrizacion de D que dimos
anteriormente. Tenemos que
kru rv k = khcos v, sen vi hu sen v, u cos vik = u,
y por lo tanto,
ZZ
f dA =
ZZ
u ududv =
[0,2][0,2]
2 Z 2
0
u dv du = 8
y tal que
100
8. INTEGRALES MULTIPLES
LECCION
D = {(x(u, v, w), y(u, v, w), z(u, v, w)) : (u, v, w) [a, b] [c, d] [e, h]} .
Supondremos adem
as que las funciones x = x(u, v, w), y = y(u, v, w),
z(u, v, w) tiene derivadas parciales continuas en (a, b) (c, d) (e, h).
z =
es una parametrizacion de la esfera D = (x, y, z) : x2 + y 2 + z 2 4 . Es claro que
D = {(u cos v sen w, u sen v sen w, u cos w) : (u, v, w) [0, 2] [0, 2] [0, ]} .
Una particion de la regi
on D en subregiones Dijk se obtiene a partir de una particion
de [a, b] [c, d] [e, h] en 3-celdas Rijk = [ui1 , ui ] [vj1 , vj ] [wk1 , wk ]:
Dijk = {(x(u, v, w), y(u, v, w), z(u, v, w)) : (u, v, w) [ui1 , ui ] [vj1 , vj ] [wk1 , wk ]} .
Escogemos en cada subregi
on Dijk el punto
Pijk (x(ui1 , vj1 , wk1 ), y(ui1 , vj1 , wk1 ), z(ui1 , vj1 , wk1 )).
Una suma de Riemann de f sobre D con respecto a la parametrizacion r y a una
particion de D en mnl subregiones, esta definida por
R(f, m, n, l, r) =
=
P m Pn
i=1
j=1
P m Pn
i=1
j=1
Pl
Pl
v rv ) w rw |
rv ) rw | uvw,
donde ru , rv y rw son los vectores cuyas componentes son las derivadas parciales
con respecto a u, a v y a w de las componentes de r en Pijk , respectivamente. La
aparicion de |(u ru v rv ) w rw | en la suma de Riemann se debe a que este
es el volumen del paraleleppedo generado por u ru , v rv y w rw , que es una
aproximacion del volumen de la subregi
on Dijk .
3
101
La integral de f sobre D con respecto a r es el lmite de R(f, m, n, l, r) cuando m,
n y l tienden a infinito:
ZZZ
f dV =
D
= lmm,n,l
=
ZZZ
Pm P n
i=1
[ab][cd][e,h]
j=1
Pl
rv ) rw | uvw
f dr =
[a,b][c,d][e,h]
ZZZ
f (u, v, w)dudvdw.
[a,b][c,d][e,h]
Es decir, las integrales sobre 3-celdas son un caso particular de integrales sobre regiones m
as generales.
Como ejemplo calculemos la integral de f (x, y, z) = x2 + y 2 + z 2 sobre la esfera
D = (x, y, z) : x2 + y 2 + z 2 4 con respecto a la parametrizacion de D que dimos
anteriormente. Tenemos que
|(ru rv ) rw | =
= |(hcos v sen w, sen v sen w, cos wi hu sen v sen w, u cos v sen w, 0i)
hu cos v cos w, u sen v cos w, u sen wi|
= u2 sen w,
102
8. INTEGRALES MULTIPLES
LECCION
y por lo tanto,
ZZZ
f dV
ZZZ
Z
u2 u2 sen wdudvdw
[0,2][0,2][o,]
2 Z 2
0
= 128/5
Z
u sen wdw dv du
0
103
EJERCICIOS
1. Calcule la suma de Riemman R (f, n, m, ) para la funcion f (x, y) = x2 y + y,
la regi
on D = [2, 4] [1, 5], m = 4, n = 5 tomando como punto muestra el
centro de cada subcelda.
2. Calcule las siguientes integrales iteradas.
Z 1Z 3Z 2
x2 y 2 z 2 dzdydx
a)
b)
c)
d)
e)
2
2Z 3
(x2 y + xy 3 )dydx
1Z 4
x
3
+ x2 )dydx
y
2
1 Z /4
0
bZ dp
s2 + t2 dsdt
3. Parametrizar la regi
on del plano limitada por las curvas dadas, de manera que
el dominio de la parametrizacion sea una 2-celda.
a) y = x, y = 0, x = 2
b) y = x2 , x = 0, y = 4
c) y = 0, y = 1 x2
p
d ) x = 1 y2 , x = 0
e) y = 3x + 2, y = x + 4, x = 0
f ) y = x2 , y = x
g) x + y = 3, x + y = 7, x y = 4, x y = 1.
4. Sea S la esfera con centro en el origen y radio 2 y P la superficie conica
p
de ecuaci
on z = x2 + y 2 , parametrice cada una de las regiones del espacio
de manera que el dominio de los par
ametros sea una 1-celda (intervalo), una
2-celda o una 3-celda, dependiendo de la naturaleza de la regi
on.
a) La superficie S.
104
8. INTEGRALES MULTIPLES
LECCION
b) El s
olido limitado por S.
c) La parte de la superficie S que esta dentro de la superficie conica P .
d ) La parte de la superficie conica que se encuentra dentro de la esfera.
e) La curva de interseccion de las dos superficies.
f ) El s
olido interior a la esfera y exterior al cono
5. Dado que S es la parte del plano de ecuaci
on 2x + 4y + 3z = 6 que se encuentra
en el primer octante, parametrice cada una de las regiones del espacio de
manera que el dominio de la parametrizacion sea una 1-celda, una 2-celda o
una 3-celda, dependiendo de la naturaleza de la regi
on.
a) La superficie S.
b) El s
olido limitado por la superficie S y los planos coordenados.
c) El borde de S.
6. Escriba cada integral como una integral sobre una 1-celda, una integral iterada
sobre una 2-celda o una integral iterada sobre una 3-celda, dependiendo de la
naturaleza de la regi
on.
R
a) C (x + y + z) dl donde C es el borde de la regi
on del numeral 3c).
R
b) C (x + y + z) dl donde C es la curva descrita en el numeral 3e).
R
c) C (x + y + z) dl donde C es la curva descrita en el numeral 5c).
R
d ) (x + y + z)dA donde es la superficie descrita en el numeral 3c).
R
e) (x + y + z)dA donde es la superficie descrita en el numeral 3d).
R
f ) (x + y + z)dA donde es la superficie descrita en el numeral 3f).
R
g) (x + y + z)dA donde es la superficie descrita en el numeral 3g).
R
h) (x + y + z)dA donde es la superficie descrita en el numeral 4a).
R
i ) (x + y + z)dA donde es la superficie descrita en el numeral 4c).
R
j ) (x + y + z)dA donde es la superficie descrita en el numeral 4d).
R
k ) (x + y + z)dA donde es la superficie descrita en el numeral 5a).
R
l ) S (x + y + z)dV donde S es s
olido descrito en el numeral 4b).
R
m) S (x + y + z)dV donde S es s
olido descrito en el numeral 4f).
105
n)
S (x
+ y + z)dV donde S es s
olido descrito en el numeral 5b).
a) Sea f : R2 R continua en R = [a, b] [c, d]; para a < x < b y c < x <
RxRy
2F
2F
=
.
d, se define F (x, y) = a c f (u, v)dudv. Demostrar que
xy
yx
Use este ejemplo para establecer la relaci
on que existe entre el teorema
de Clairaut y el de Fubini.
R 1 R 1 x2 y 2
dxdy
b) Determine el valor de cada una de las siguientes integrales 0 0 2
(x + y 2 )
R 1 R 1 x2 y 2
dydx Contradice esto el teorema de Fubini?
y 0 0 2
(x + y 2 )
a) f (x, y) = 4xy,
b) f (x, y) = ex+y ,
c) f (x, y) = x3 + y 2 ,
R = [0, 1] [1, 1]
d ) f (x, y) = x4 + xy + y 3 ,
R = [0, 1] [0, 1]
R = [1, 2] [0, 2]
10. Eval
ue la integral doble dada:
a)
c)
1Z 2
1
2Z 1
0
(1 y)x dx dy
(x + 2) dx dy
b)
d)
1Z 2
x(x + y) dx dy
2
1
x(xy + sin x) dx dy
1
106
8. INTEGRALES MULTIPLES
LECCION
e)
g)
/2 Z 1
xy cos(x2 y) dx dy
f)
2Z 4
xy dx dy
h)
c)
e)
g)
1Z 1
RdRb
b)
24x y dy dx
2 Z ln x
d)
4x dy dx
/2 Z y
f)
cos x sin y dx dy
2Z y
1 dx dy
h)
0
1
1
/2
sin x cos(y ) dx dy
1 dx dy
1
M dx dy = M (d c)(b a).
Z
Z
sin x dx dy
0
2 Z 2y
ey dx dy
xye(x
2 +y 2 )
dx dy
1 Z x2
2 dy dx
3Z 2Z 1
0
xyz dx dy dz
b)
x Z xy
0
x2 sin z dz dy dx
d)
1Z xZ y
1Z z
0
xyz dz dy dx
y
0
zey dx dy dz
107
e)
g)
eZ y
0
1/y
f)
2Z 4Z 3
2
x2 z dx dz dy
1 dx dy dz
h)
2 Z y2
0
y1
z2
yz dx dz dy
0
1 Z 1x Z 1xy
0
R z 2 R y2 R x 2
z1
x1
1 dz dy dx
M dx dy dz = M (z2 z1 )(y2
21. Escriba un programa que use el metodo de Monte Carlo para aproximar la
RR xy
e dA,donde R = [0, 1] [0, 1]. Muestre el resultado del
integral doble
R
programa para N = 10, 100, 1000, 10000, 100000 y 1000000 puntos aleatorios.
22. Escriba un programa que use el metodo de Monte Carlo para aproximar la
RRR xyz
e dV , donde S = [0, 1] [0, 1] [0, 1]. Muestre el resulintegral triple
S
tado del programa para N = 10, 100, 1000, 10000, 100000 y 1000000 puntos
aleatorios.
x2
a2
y2
b2
=1
108
8. INTEGRALES MULTIPLES
LECCION
x2
a2
+ yb2 + zc2 = 1
1
0
tx1 (1 t)y1 dt ,
satisface la relaci
on B(y, x) = B(x, y) para x > 0, y > 0.
32. Use la sustituci
on t = u/(u + 1), para mostrar que la funcion Beta se puede
escribir como
Z
ux1
B(x, y) =
du , for x > 0, y > 0.
(u + 1)x+y
0
33. Encuentre el centro de masa de la regi
on R con funcion de densidad dada
(x, y):
a) R = {(x, y) : 0 x 2, 0 y 4 }, (x, y) = 2y
b) R = {(x, y) : 0 x 1, 0 y x2 }, (x, y) = x + y
c) R = {(x, y) : y 0, x2 + y 2 a2 }, (x, y) = 1
d ) R = {(x, y) : y 0, x 0, 1 x2 + y 2 4 }, (x, y) =
e) R = {(x, y) : y 0, x2 + y 2 1 }, (x, y) = y
x2 + y 2
x2
dx.
(x)2 /2 2
1
2 e
dx.
109
37. Para variables aleatorias continuas X, Y con funcion de distribucion de probabilidad conjunta F , defina los segundos momentos E(X 2 ) y E(Y 2 ) por
Z Z
Z Z
2
2
2
y 2 f (x, y) dx dy ,
x f (x, y) dx dy y E(Y ) =
E(X ) =
and
110
8. INTEGRALES MULTIPLES
LECCION
LECCION
Integrales multiples
Estudiaremos dos tipos de dominios del plano denominados regiones del tipo I y
regiones del tipo II. Las regiones del tipo I son aquellas limitadas por rectas x = a
y x = b y por las gr
aficas de dos funciones y = g1 (x) y y = g2 (x) diferenciables en
[a, b]. Las regiones del tipo II son aquellas limitadas por las rectas y = c y y = d y
por las gr
aficas de dos funciones x = h1 (y) y x = h2 (y) diferenciables en [c, d].
Las regiones del tipo I son de la forma:
= {(x, y) : a x b,
111
g1 (x) y g2 (x)} .
112
9. INTEGRALES MULTIPLES
LECCION
y = g2 (x)
y = g1 (x)
a
rv
y por consiguiente que
ZZ
f dA =
ZZ
[a,b][0,1]
Z b Z
a
1
0
du.
113
Haciendo el cambio de variable y = (1v)g1 (u)+vg2 (u) y tomando x = u se obtiene
!
Z
Z
ZZ
b
g2 (x)
f dA =
g1 (x)
c y d} ,
Y
d
x = h1 (y)
x = h2 (y)
c
X
h2 (y)
f dA =
f (x, y)dx dy
h1 (y)
114
9. INTEGRALES MULTIPLES
LECCION
f dV =
ZZ
g2 (x,y)
f (x, y, z)dz
g1 (x,y)
dA,
si es del tipo I,
ZZZ
f dV =
ZZ
l2 (x,z)
f dV =
ZZ
h2 (y,z)
115
Si el dominio es una superficie en el espacio podemos proceder exactamente de la
misma manera como lo hicimos en el caso de regiones del plano.
Sea f : R una funcion continua, definida en una superficie acotada del
espacio. Una parametrizacion de es una funcion vectorial
r : [a, b] [c, d] R3 ,
continua
y tal que
= {(x(u, v), y(u, v), z(u, v)) : (u, v) [a, b] [c, d]} .
Supondremos adem
as que r tiene derivadas parciales continuas en (a, b) (c, d).
Por ejemplo, una parametrizacion de la superficie de la esfera = (x, y, z) : x2 + y 2 + z 2 = 4
es r : [0, 2] [0, ] : R3 definida por r(u, v) = h2 cos u sen v, 2 sen u sen v, 2 cos vi.
Es claro que
= {(2 cos u sen v, 2 sen u sen v, 2 cos v) : (u, v) [0, 2] [0, ]} .
Una particion de la superficie en subregiones ij se obtiene a partir de una partici
on de [a, b] [c, d] en 2-celdas Rij = [ui1 , ui ] [vj1 , vj ]:
ij = {(x(u, v), y(u, v), z(u, v)) : (u, v) [ui1 , ui ] [vj1 , vj ]} .
Escogemos en cada subregi
on ij el punto
Pij (x(ui1 , vj1 ), y(ui1 , vj1 ), z(ui1 , vj1 )).
Una suma de Riemann de f sobre con respecto a la parametrizacion r y a una
particion de en m n subregiones, esta definida por
Pm Pn
R(f, m, n, r) =
i=1
j=1 f (Pij ) ku ru v rv k ,
=
f
(P
)
kr
r
k
uv,
ij
u
v
j=1
P P
m
n
i=1
donde ru y rv son las derivadas parciales de r con respecto a u y a v en Pij respectivamente. La aparicion de ku ru v rv k en la suma de Riemann se debe a que esta
1
116
9. INTEGRALES MULTIPLES
LECCION
ZZ
n
m X
X
f dA = lmm,n
f (Pi,j ) kru rv k uv
i=1 j=1
ZZ
[a,b][c,d]
4 sen vdudv =
[0,][0,2]
2
0
Z
Nos queda estudiar el caso de una funcion definida en una curva. Sea f : R
una funcion continua, definida en una curva acotada del espacio (o del plano).
Una parametrizacion de es una funcion vectorial
r : [a, b] R3 ,
117
Por ejemplo, una parametrizacion de la circunferencia
= (x, y, z) : x2 + y 2 = 4, z = 5
esta dada por la funcion vectorial r : [0, 2] : R3 definida por r(u) = h2 cos u, 2 sen u, 5i.
Es claro que
= {(2 cos u, 2 sen u, 5) : u [0, 2]} .
Una particion de en subarcos i se obtiene a partir de una particion de [a, b] en
subintervalos Ii = [ui1 , ui ]:
i = {(x(u), y(u), z(u)) : u [ui1 , ui ]} .
Escogemos en cada subarco i el punto Pi (x(ui1 ), y(ui1 ), z(ui1 )). Una suma de
Riemann de f sobre con respecto a la parametrizacion r y a una particion de
en m subarcos, esta definida por
R(f, m, r) =
m
X
i=1
f (Pi )
u r (ui1 )
,
m
X
i=1
r (u
!
f (Pi )
r (ui1 )
u,
donde
on de ku r (ui1 )k en la suma de
i1 ) es la derivada r en Pi . La aparici
Riemann se debe a que esta es la longitud del vector tangente u r (ui1 ) que es
una aproximacion de la longitud del subarco i .
La integral de f sobre con respecto a r es el lmite de R(f, m, r) cuando m tiende
a infinito:
!
Z
m
X
f dL = lmm
f (Pi )
r (ui1 )
u
i=1
[a,b]
f (x(u), y(u), z(u))
r (u)
du
118
9. INTEGRALES MULTIPLES
LECCION
y por lo tanto,
Z
f dL =
[0,2]
29 2 du =
58 du = 116.
0
119
EJERCICIOS
1. Suponiendo que la densidad en cada punto (x, y) de la una placa delgada
esta dada por (x, y) = x + y
7
6
5
4
3
2
1
1
1
II
III
1 2 3 4 5 6 7 8 9
|y|
y
2
4Z 2
ex dxdy
120
9. INTEGRALES MULTIPLES
LECCION
c)
d)
1
3
Z 9y2
1Z
xdxdy
9y 2
arctan y
x2 y 2
dxdy
(x2 + y 2 )
ZZ
(x, y) dA
ZZ
G(x, y)(x, y) dA M
ZZ
(x, y) dA,
R
RR
G(x, y)(x, y) dA
RRR
R (x, y) dA
M.
RR
G(x, y)(x, y) dA
RRR
R (x, y) dA.
Este resultado se conoce como el Teorema del Valor Medio para Integrales
dobles.
Use la desigualdad del valor medio para demostrar que
a) Si f (x, y) = esin(x+y) y R = [, ] [, ] entonces
1
1
2
e
4
b) Si f (x, y) =
ZZ
f (x, y)dA e.
dA
y R = [1, 1] [1, 2] entonces
x2 + y 2 + 1
ZZ
dA
1
6
2
2
R x +y +1
121
RRR
R1RxRy
olido y plantee
4. Si
f
(x,
y,
z)dV
=
p
2
2
2
olido limitado por las
a) f (x, y, z) = x2 + y 2 + z 2 e(x +y +z ) y es el s
esferas con centro en el origen y radios a y b respectivamente, tomando
a < b.
1
b) f (x, y, z) = p
y es la bola unidad (esfera con centro en el
x2 + y 2 + z 2
origen y radio 1).
p
c) f (x, y, z) = x2 + y 2 + z 2 y = (x, y, z) : 21 z 1, x2 + y 2 + z 2 1
6. Sea el elipsoide s
olido
x2 y 2 z 2
+ 2 + 2 1 donde a, b, c son positivos.
a2
b
c
a) Determine el volumen de .
2
RRR
y2 z2
x
dV
b) Calcule
a 2 + b2 + c 2
7. Sea la regi
on en el plano XY interior a la gr
afica de x2 + y 2 = 1 y exterior
a la gr
afica de x2 + y 2 = 2y con x, y 0.
a) Dibujar la regi
on .
b) Si u = x2 + y 2 y v = x2 + y 2 2y, dibujar la regi
on que corresponde
a bajo esta transformacion.
RR
c) Calcular xey dA usando la transformacion.
122
9. INTEGRALES MULTIPLES
LECCION
RRR m n p
12. Calcular
olida con centro en el origen y radio
x y z dV si es la esfera s
1 y m, n, p son n
umeros enteros no negativos.
RRR p q r
s
13. Calcular la integral de Dirichlet
x y z (1 x y z) dV con p, q, r, s
n
umeros positivos si es el s
olido limitado por las superficies x + y + z = 1,
x = 0, y = 0 y z = 0, haciendo x + y + z = , y + z = , z =
14. En que razon divide la superficie x2 + y 2 + az = 4a2 el volumen de la esfera
x2 + y 2 + z 2 4az?
15. Suponga que el plano z = intersecta la esfera s
olida con centro en el origen
y radio 1, determinando dos s
olidos. Encuentre el valor de para el cual la
razon entre los vol
umenes de los dos s
olidos es 3.
16. Dada la esfera de ecuaci
on x2 + y 2 + z 2 = a2 y el cono de ecuaci
on z =
p
k x2 + y 2 , determine el valor de la constante k para el cual la razon entre las
areas de la partes de la esfera interior al cono y la parte de la esfera exterior
al cono es 2/5.
17. Determine la longitud de cada una de las siguientes curvas:
a) La parte de la par
abola en el plano XY de ecuaci
on y = x2 desde el punto
(1, 1) hasta el punto (2, 4).
b) La parte de la elice de ecuaci
on r(t) = cos ti + sen tj + tk desde el punto
(1, 0, 0) hasta el punto (1, 0, 2).
c) La elpse en el plano XY de ecuaci
on
x2
a2
y2
b2
= 1.
18. Determine el
area de cada una de las siguientes superficies:
a) La parte de la esfera con centro en el origen y radio a que se encuentra
p
dentro del cono de ecuaci
on z = x2 + y 2 .
b) La parte de la superficie cilndrica de ecuaci
on z = 4 y 2 cuyos puntos
satisfacen la desigualdad x2 + y 2 4.
x2
a2
y2
b2
=1
123
19. Determine el centro de masa de:
a) Un alambre delgado modelado mediante la curva r(t) = ti+2tj+( 23 )t3/2 k,
0 t 2, si la densidad es = 3 5 + t. (Thomas)
p
b) Un s
olido acotado por el cono de ecuaci
on z = x2 + y 2 y por el plano
z = 1 si la densidad en cada punto esta dada por = z 2 .(Thomas)
c) Una l
amina delgada modelada por la superficie que esta acotada por las
curvas de ecuaciones x = y 2 y x = 2y y 2 si la densidad en cada punto
esta dada por = y + 1. (Thomas)
p
d ) Una l
amina delgada modelada por la superficie del cono z = x2 + y 2
cortada por el cilndro circular x2 + y 2 = 2x. (Thomas)
e) Un s
olido acotado por el el elipsoide
cada punto esta dada por = z
x2
a2
y2
b2
z2
c2
= 1 si la densidad en
124
9. INTEGRALES MULTIPLES
LECCION
LECCION
10
Campos vectoriales
(x, y) D.
(x, y, z) D.
126
El dominio D de un campo vectorial es la interseccion de los dominios de sus funciones componentes. El rango es el conjunto de las imagenes de los puntos del dominio.
La gr
afica es el conjunto de todas las parejas formadas por cada elemento del dominio y su respectiva imagen.
El lmite de un campo vectorial en un punto de su dominio existe si existen los lmites de sus funciones componentes y las derivadas parciales de un campo vectorial
existen si existen las derivadas de sus funciones componentes, los calculos se hacen
componente a componente. Por u
ltimo, un campo vectorial es continuo o diferenciable si sus funciones componentes lo son.
La mejor forma de representar gr
aficamente un campo vectorial bidimensional es
dibujando a partir de unos cuantos puntos (x, y) D los vectores F(x, y) correspondientes. Para ilustrar la representaci
on gr
afica de un campo vectorial, veamos el
siguiente ejemplo:
Consideremos el campo vectorial F(x, y) = h1, x + yi y dibujemos en el plano los
vectores F(0, 0), F(1, 1), F(1, 1), F(2, 0), F(0, 2), F(2, 0), y F(0, 2).
Y
La circulaci
on y el flujo.
Los dos conceptos m
as importantes, desde el punto de vista fsico, relacionados con
campos vectoriales son el flujo y la circulaci
on, y que se hacen evidentes al estudiar
127
la interacci
on de los campos con curvas y superficies. Estudiaremos primero estos
dos conceptos en campos vectoriales bidimensionales para los cuales nos interesa
estudiar su interacci
on con curvas contenidas en sus dominios.
Consideremos un campo vectorial bidimensional F : D R2 definido en un subconjunto D de R2 y una curva en D parametrizada por una funcion vectorial
r : [a, b] R2 . Definimos la densidad de circulaci
on, c (x(t), y(t)), del campo F
en el punto (x(t), y(t)) de la curva como la componente tangencial de F(x(t), y(t)),
esto es
c (x(t), y(t)) = F(x(t), y(t)) T(t),
donde T(t) es el vector unitario tangente a la curva en el punto (x(t), y(t)).
Y
F
N
Definimos la densidad de flujo del campo F en el punto (x(t), y(t)) de la curva como
la componente normal de F(x(t), y(t)), esto es
f (x(t), y(t)) = F(x(t), y(t)) N(t),
donde N(t) es el vector unitario normal a la curva en el punto (x(t), y(t)), que se
obtiene al girar T(t) 90 grados en sentido positivo. En otras palabras, la densidad de
flujo de F en el punto (x(t), y(t)), f (x(t), y(t)), es la componente normal del campo
128
N=
x (t)
y (t)
,p
p
x (t)2 + y (t)2
x (t)2 + y (t)2
c dL =
Z
Z
Z
Z
b
a
b
F(r(t)) T(t)
r (t)
dt
F(r(t))
a
b
a
b
r (t)
r (t)
dt
kr (t)k
F(r(t)) r (t)dt
(P (r(t))x (t) + Q(r(t))y (t))dt.
Debido a la forma de la u
ltima integral, la circulaci
on de F a lo largo de se suele
R
R
escribir como P dx + Qdy o como F dr.
El flujo de F a traves de es
f dL,
129
ponentes escalares de F as:
Z
f dL =
Z
Z
Z
Z
b
a
b
F(r(t))
a
b
a
b
a
F(r(t)) N(t)
r (t)
dt
*
y (t)
x (t)
p
,p
x (t)2 + y (t)2
x (t)2 + y (t)2
r (t)
dt
hP (r(t)), Q(r(t))i y (t), x (t) dt
(Q(r(t))x (t) P (r(t))y (t))dt.
Debido a la forma de la u
ltima integral, el flujo de F a traves de se suele escribir
R
R
e dr, donde F
e = hQ, P i se conoce como el campo
como Qdx P dy o como F
dual del campo F. Lo que se muestra en la cadena de igualdades arriba, es que el
flujo de un campo a traves de una curva es igual a la circulaci
on de su campo dual
a lo largo de la misma curva.
Densidad rotacional y densidad de expansi
on.
Los conceptos de densidad rotacional y de densidad de expansion de un campo vectorial 1 , son fundamentales para entender los conceptos de rotacional y divergencia
de un campo vectorial y comprender los teoremas fundamentales: el Teorema de
Stokes y el Teorema de la Divergencia. Estas densidades no son otra cosa que la
circulaci
on y el flujo por unidad de
area, respectivamente, en cada uno de los puntos
del dominio del campo vectorial.
Calculemos la densidad rotacional de un campo vectorial bidimensionale F(x, y) =
hP (x, y), Q(x, y)i en el punto (0, 0) y que estamos suponiendo esta en el dominio de
F. Supondremos adem
as que las componentes escalares del campo, P y Q, tienen
derivadas parciales continuas en una vecindad de (0, 0) y consideremos la circulaci
on
de F a lo largo del borde del rectangulo de vertices V1 (x, y), V2 (x, y),
1
130
Circulacion de F a lo largo de
area de
4xy{Qx (, 0) Py (0, )}
4xy
= Qx (0, 0) Py (0, 0)
Calculemos ahora la densidad de expansion de F en el punto (0, 0), con las mismas
suposiciones hechas anteriormente. Para esto, consideremos el flujo de F a traves
del borde del rectangulo de vertices V1 (x, y), V2 (x, y), V3 (x, y)
y V4 (x, y), orientado negativamente. El flujo de F a traves del lado V1 V2 , que
2
131
escribiremos FV1 V2 F es, aproximadamente, el producto de la componente normal de
F en el punto medio de V1 V2 y de la longitud de V1 V2 , es decir FV1 V2 F = {F(0, y)
j}{2x}. El flujo del campo a traves de los demas lados se puede expresar de la
misma forma, y as, se tiene que el flujo de F a traves de que se nota como f F es:
f F = FV1 V2 F + FV2 V3 F + FV3 V4 F + FV4 V1 F
{F(0, y) j}{2x} + {F(x, 0) (i)}{2y}
+{F(0, y) (j)}{2x} + {F(x, 0) i}{2y}
= 2y{P (x, 0) P (x, 0)} + 2x{Q(0, y) Q(0, y)}
= 4xyPx (, 0) + 4xyQy (0, )
= 4xy{Px (, 0) + Qy (0, )},
donde (x, x) y donde (y, y). En el pen
ultimo paso se hizo uso del
Teorema del Valor medio para derivadas. Si usamos la continuidad de las derivadas
parciales de P y Q obtenemos la densidad de expansi
on de F en (0, 0), que es el
flujo por unidad de
area de F en (0, 0) y denotaremos por E (0, 0), es:
E (0, 0) = lm(x,y)(0,0)
= lm(x,y)(0,0)
Flujo de F a traves de
area de
4xy{Px (, 0)) + Qy (0, )}
4xy
= Px (0, 0) + Qy (0, 0)
Las densidades rotacional y de expansion se pueden calcular en todos los puntos
del dominio del campo F de la misma forma en que se calcul
o en el punto (0, 0),
para encontrar que la densidad rotacional de F en (x, y), R (x, y) y la densidad de
espansi
on de F en (x, y), E (x, y), estan dadas por:
R (x, y) = Qx (x, y) Py (x, y),
respectivamente.
132
Debemos se
nalar que el concepto de circulaci
on a lo largo de una curva se puede
extender de forma natural a campos vectoriales tridimensionales, pero el concepto
de flujo a traves de una curva no (por que?). En la siguiente leccion discutiremos los
conceptos, flujo, densidad de expansion y densidad de rotaci
on de campos vectoriales
tridimensionales.
Los campos gradiente
Un ejemplo de campo vectorial es el que proviene del gradiente de una funcion. El
argumento que presentaremos a continuaci
on es valido para funciones de dos o de
tres variables, y demuestra que la circulaci
on de un campo gradiente de una funcion
f depende exclusivamente de los valores de f en los puntos inicial y final de la
curva. Considere una funcion f de varias variables con dominio D y una curva en
D parametrizada por una funcion vectorial r : [a, b] D continua. Calculemos la
circulaci
on del campo vectorial f a lo largo de :
Z
f dL =
b
a
b
a
f (r(t)) r (t)dt
d(f r)
(t)dt
dt
= f (r(b)) f (r(a))
Este resultado significa que si F es un campo gradiente, la circulaci
on de F a lo
largo de cualquier curva dentro del dominio s
olamente depende de los puntos inicial
y final de la curva, recprocamente, si la circulaci
on de un campo F a lo largo de
cualquier curva dentro del dominio s
olamente depende de los puntos inicial y final
de la curva, entonces el campo vectorial es un campo gradiente. En efecto, debido
R
a la suposicion hecha, esta bien definida la funcion f (x, y, z) = [(a,b,c),(x,y,z)] c dl,
donde la integral se calcula sobre cualquier curva que parta del punto (a, b, c) y
termine en el punto (x, y, z), y se puede demostrar que F(x, y, z) = f (x, y, z). Los
campos vectoriales F para los que existe una funcion f tal que F = f , se llaman
campos conservativos y, en este caso, a la funcion f se le llama funcion potencial de f .
133
Para terminar esta secci
on anotamos que si el campo F es un campo de fuerzas, la
R
Rb
integral c dL = a F(r(t)) T(t) kr (t)k dt representa el trabajo realizado por el
campo al desplazar un objeto sobre la curva desde el punto inicial hasta el punto
final.
134
EJERCICIOS
1. Suponga que (0, 0, 0) es un punto interior del dominio del campo vectorial
F(x, y, z) = P (x, y, z)i + Q(x, y, z)j + R(x, y, z)k y que las funciones P, Q y
R tienen derivadas parciales de primer orden, continuas en (0, 0, 0). Calcule la
densidad de expansion de F en (0, 0, 0). Ayuda: imite el calculo de la densidad
de expansion de un campo bidimensional considerando una 3-celda conveniente
con centro en (0, 0, 0).
R
2. Muestre que si F dr depende u
nicamente de los puntos inicial y final de ,
entonces F = f para alguna funcion escalar f .
3. Sea una curva suave parametrizada por r = r(t) y F un campo vectorial
continuo definido sobre .
R
a) Si F(r(t)) es perpendicular a r (t), muestre que F dr = 0.
R
R
b) Si F(r(t)) es paralelo a r (t), muestre que F dr = kFk dL.
5. Calcule la circulaci
on de F(x, y, z) = (z 3 + 2xy)i + x2 j + 3xz 2 k a lo largo del
borde de la cara superior del cubo de vertices (1, 1, 1).
2
135
8. Para cada uno de los siguientes campos vectoriales haga una representaci
on
R
gr
afica del campo y la curva, y calcule F dr.
a) F(x, y) = i j;
: x = 3t, y = 2t, 0 t 1
b) F(x, y) = y i x j;
: x = cos t, y = sin t, 0 t 2
c) F(x, y) = x i + y j;
: x = cos t, y = sin t, 0 t 2
d ) F(x, y) = (x2 y) i + (x y 2 ) j;
e) F(x, y) = xy 2 i + xy 3 j;
hasta (0, 1) hasta (0, 0)
f ) F(x, y) = (x2 + y 2 ) i;
9. Demuestre que
f (x, y) dL
: x = cos t, y = sin t, 0 t 2
f (x, y) dL.
10. Sea una curva suave con longitud de arco L, y suponga que F(x, y) =
P (x, y) i + Q(x, y) j es un campo vectorial tal que kF(x, y)k M para todo
R
(x, y) en . Muestre que | F dr| M L.
11. Calcule
12. Calcule
13. En cada caso, si hay una funcion potencial f para F encuentre una:
a) F(x, y) = y i x j.
b) F(x, y) = x i y j.
c) F(x, y) = xy 2 i + x3 y j.
14. Sean F y G campos vectoriales, sean a y b constantes, y sea una curva en
R2 . Muestre que
Z
Z
Z
(a F b G) dr = a F dr b G dr .
1 dL
= .
136
x
i + x2 +y
17. Suponga que F(x, y) = x2y
2 j para todo (x, y) 6= (0, 0) es un campo
+y 2
de fuerzas, y sea : x = cos t, y = sin t, 0 t 2.
LECCION
11
Divergencia y rotacional
En esta secci
on presentaremos la generalizaci
on de los conceptos de densidad de expansion y densidad rotacional a campos vectoriales tridimensionales y enunciaremos
los Teoremas de la Divergencia y del Rotacional.
Densidad de expansi
on para campos tridimensionales
As como la densidad de expansion para un campo vectorial bidimensional es el flujo
del campo por unidad de
area, la densidad de expansion de un campo tridimensional
es el flujo del campo por unidad de volumen. Para hacer esta generalizaci
on debemos
precisar el concepto de flujo de un campo a traves de una superficie:
Sea F(x, y, z) = P (x, y, z)i + Q(x, y, z)j + R(x, y, z)k un campo vectorial tridimensional y una superficie contenida en el dominio del campo. Supongamos que es
una superficie orientable, esto es, se puede parametrizar por una funcion vectorial
r(u, v) = x(u, v)i + y(u, v)j + z(u, v)k, con (u, v) [a, b] [c, d], ru , rv continuas y
ru rv 6= 0 para todo (u, v). Bajo estas condiciones, en cada punto de podemos
ru rv
considerar el vector normal unitario N =
y por consiguiente definir la
kru rv k
densidad de flujo f de F en cada punto de , como la componente normal de F,
esto es f = F N . As, el flujo de F a traves de la superficie es la integral de f
137
138
sobre la superficie:
R
f dr
=
=
=
RR
F NdA
RbRd
a
RbRd
a
F(r(u, v))
ru rv
kru rv k dvdu
kru rv k
Con este concepto de flujo de un campo a traves de una superficie podemos generalizar el concepto de densidad de expansion de un campo tridimensional en los
puntos de su dominio. Supongamos que (0, 0, 0) pertenece al dominio del campo
F(x, y, z) = P (x, y, z)i + Q(x, y, z)j + R(x, y, z)k y que las componentes del campo
tienen derivadas parciales continuas en una vecindad de (0, 0, 0). Consideremos la
3-celda [x, x] [y, y] [z, z], que esta centrada en (0, 0, 0), con x,
y y z lo suficientemente peque
nos para que la 3-celda este contenida en el dominio del campo. Procediendo de manera similar a como se hizo en la leccion anterior,
podemos hacer una aproximacion del flujo del campo F a traves de la superficie de
la 3-celda, dividir esta cantidad por el volumen de la 3-celda y hacer tender las longitudes de las aristas a cero. Se Obtiene as la densidad de expansion E del campo
F en el punto (0, 0, 0):
E (0, 0, 0) = Px (0, 0, 0) + Qy (0, 0, 0) + Rz (0, 0, 0).
Al hacer estos mismos calculos para cualquier punto (x, y, z) del dominio de F se
obtiene que E (x, y, z) = Px (x, y, z) + Qy (x, y, z) + Rz (x, y, z). Es frecuente en la
literatura sobre el tema llamar a la densidad de expansion divergencia del campo y
escribirlo como:
divF(x, y, z) = Px (x, y, z) + Qy (x, y, z) + Rz (x, y, z).
As, podemos escribir
E (x, y, z) = divF(x, y, z)
139
La circulaci
on de un campo tiene sentido s
olamente a lo largo de curvas, debido a que
para cada punto en una curva se puede determinar un u
nico vector tangente unitario, mientras que si el punto se considera en una superficie, existen infinitos vectores
unitarios tangentes a la superficie en el punto. Por lo tanto, al igual que para campos
bidimensionales, definimos la densidad rotacional de un campo tridimensional como
la circulaci
on por unidad de
area. Sin embargo, note que si se tienen una curva y
un punto en el espacio, existen muchas superficies que contienen el punto y tienen
como borde la curva, por lo tanto, para precisar la curva y el area determinaremos
la densidad de circulaci
on referida a un plano que contenga al punto y haremos los
calculos de la misma manera que se hizo en el caso bidimensional.
Supongamos que (0, 0, 0) pertenece al dominio del campo F(x, y, z) = P (x, y, z)i +
Q(x, y, z)j + R(x, y, z)k y que el campo tiene derivadas parciales continuas en una
vecindad de (0, 0, 0). Tomemos el plano de ecuaci
on ax + by z = 0, que se puede
parametrizar con la funcion r(x, y) = hx, y, ax + byi, tomemos el cuadrilatero de
vertices V1 (x, y, ax+by), V2 (x, y, ax+by), V3 (x, y, ax
by) y V4 (x, y, ax by). Sea la curva cerrada formada por los segmentos
V1 V2 , V2 V3 , V3 V4 y V4 V1 , recorrida en el sentido descrito y la regi
on del plano
limitada por . Calculemos ahora una aproximacion de la circulaci
on del campo F
a lo largo , de manera identica a como se hizo en el caso bidimensional. Si C F es
la circulaci
on de F a lo largo de :
C F
1
h0, 1, bi
2y 1 + b2 F(x, 0, ax)
1 + b2
1
+2x 1 + a2 F(0, y, by)
h1, 0, ai
1 + a2
1
h0, 1, bi
+2y 1 + b2 F(x, 0, a x)
1 + b2
1
h1, 0, ai
+2x 1 + a2 F(0, y, b)
1 + a2
140
R (0, 0, 0)
lm(x,y)(0,0)
lm(x,y)(0,0)
C F
A()
1 + a 2 + b2
1 + a 2 + b2
Procediendo de manera similar para cualquier punto (x, y, z) del dominio del campo F
obtenemos que la densidad rotacional del campo, en el punto, referida al plano de ecuacion
ax + by z=0, es:
R (x, y, z) =
1
ha, b, 1i hRy Qz , (Rx Pz ), Qx Py i |(x,y,z)
1 + a 2 + b2
1
ha, b, 1i
1 + a 2 + b2
141
Para recordar las componentes del vector rotF se usa la siguiente formula mnemotecnica
i
j
k
rotF = x y z
P Q R
Podemos ahora calcular la densidad rotacional de un campo F referida a una superficie orientable , calculando su densidad rotacional en cada punto de la superficie referida al plano
tangente a la misma en el punto, es decir calculando la componente normal del rotacional
de F. Explcitamente si (x, y, z) es un punto de la superficie orientable , parametrizada por
r(u, v) se tiene que:
ru rv
.
R (x, y, z) = rotF(r(u, v))
kru rv k
Antes de enunciar los teoremas fundamentales, haremos algunas aclaraciones sobre la notacion. Primero que todo, recordemos que para definir el concepto de integral m
ultiple sobre
curvas, superficies y s
olidos, introdujimos el concepto de parametrizaci
on mediante funciones vectoriales r de una, dos y tres variables respectivamente. Para poder diferenciar la
dimension de la integral introducimos aqu los siguientes elementos de integraci
on:
1. Elemento de longitud de arco:
dL = r (t)dt
dL = kr (t)k dt
2. Elemento de
area de superficie:
dA = (ru rv )dudv
dA = kru rv k dudv
3. Elemento de volumen de s
olido:
dV = |(ru (rv rw )| dudvdw
Con esta terminologa podemos escribir:
1. Si es una curva parametrizada por la funcion vectorial derivable r = r(t), a t b
Z
y
F dL =
f dL =
b
a
F(r(t)) r (t)dt
b
a
142
f dA =
b
a
d
c
3. Si es una s
olido parametrizado por la funcion vectorial diferenciable r = r(u, v, w),
a u b, c v d, e w h
ZZZ
Z bZ dZ h
f (r(u, v, w)) |ru (ru rw )| dwdvdu.
f dV =
e
R
=
=
=
RbRd
a
RbRd
a
RR
rotF(r(u, v))
ru rv
kru rv k dvdu
kru rv k
rotF dA
La expansi
on total de un campo F sobre el s
olido :
RRR
RRR
divFdV
dV =
E
=
RbRdRh
a
Rb
a
Rb
F(r(t))
r (t)
kr (t)k dt
kr (t)k
F(r(t)) r (t)dt
F dL
143
Teorema de la divergencia: Suponga que F= hP, Q, Ri es un campo vectorial y una
superficie cerrada que encierra un s
olido contenidos en el dominio del campo y que
se puede parametrizar mediante una funcion vectorial diferenciable r : P R3 . Si las
componentes escalares de F tienen derivadas parciales continuas, entonces, la expansi
on
total del campo en es igual al flujo de F a traves de la frontera de . Mas precisamente
se tiene que:
ZZZ
ZZ
divF dV =
F dA
es una superficie con borde contenidos en el dominio del campo y que se puede
parametrizar mediante una funcion vectorial diferenciable r : R R3 . Si las componentes
escalares de F tienen derivadas parciales continuas, entonces, la rotacion total del campo
sobre es igual a la circulacion de F a lo largo del borde de . Mas precisamente se tiene
que:
ZZ
Z
rotF dA =
F dL
144
EJERCICIOS
: x = 3t, y = 2t, 0 t 1
b) F(x, y) = y i x j;
: x = cos t, y = sin t, 0 t 2
c) F(x, y) = x i + y j;
: x = cos t, y = sin t, 0 t 2
d ) F(x, y) = (x y) i + (x y 2 ) j;
e) F(x, y) = xy 2 i + xy 3 j;
(0, 0).
f ) F(x, y) = (x2 + y 2 ) i;
: x = cos t, y = sin t, 0 t 2
2. Sea una curva con longitud de arco l y suponga que F(x, y) = P (x, y) i + Q(x, y) j
es un campo vectorial tal que kF(x, y)k M para todo (x, y) en . Muestre que
Z
F dr M l.
3. Calcule
4. Calcule
5. En cada uno de los numerales determine si el campo F tiene una funcion potencial y
si la tiene determnela
a) F(x, y) = y i x j.
b) F(x, y) = x i y j.
c) F(x, y) = xy 2 i + x3 y j.
d ) F(x, y, z) = y i x j + z k
e) F(x, y, z) = a i + b j + c k (a, b, c constantes)
f ) F(x, y, z) = (x + y) i + x j + z 2 k
g) F(x, y, z) = xy i (x yz 2 ) j + y 2 z k
6. Use el Teorema del Rotacional1 para calcular las siguientes integrales. Se considera
orientada positivamente.
1
Para superficies planas el Teorema del Rotacional se conoce como Teorema de Green
145
a)
2x }
Z
x2 y dx + 2xy dy; es la frontera de = { (x, y) : 0 x 1, x2 y x }
b)
Z
c)
2y dx 3x dy; es la circunferencia x2 + y 2 = 1
Z
2
2
d)
(ex + y 2 ) dx + (ey + x2 ) dy; es la frontera del triangulo de vertices (0, 0),
a dx +
7. Use el Teorema de Green para demostrar que si es una curva plana, cerrada, simple
(que no se cruza ella misma), orientada positivamente y es la regi
on del plano
acotada por entonces el
area de , A() est
a dada por
Z
Z
Z
1
xdy ydx.
A() =
y dx =
x dy =
2
C
C
8. Calcule
vertices (1, 1), (1, 1), (1, 1) y (1, 1), orientada positivamente.
9. Calcule el flujo del campo F a traves de la superficie
a) F(x, y, z) = xi + yj + zk, : borde del cubo s
olido = { (x, y, z) : 0 x, y, z
1 }.
b) F(x, y, z) = x2 i + xyj + zk, es la parte del plano 6x + 3y + 2z = 6 con x 0,
y 0, and z 0, y vector normal unitario con tercera componente positiva.
c) F(x, y, z) = x2 i + xyj + xzk, es la semiesfera x2 + y 2 + z 2 = 4, z 0.
146
11. Muestre que el flujo de un campo vectorial constante a traves de cualquier superficie
cerrada es cero.
12. Encuentre un campo vectorial tal que su flujo a traves de la superficie de una esfera
de radio r sea el
area de esa esfera.
13. El elipsoide
x2
a2
: x = t sin t, y = t cos t, z =
2 2 3/2
,
3 t
0t1
15. En cada numeral calcule la circulacion total del campo vectorial F a lo largo de la
curva dada.
a) F(x, y, z) = i j + k;
: x = 3t, y = 2t, z = t, 0 t 1
b) F(x, y, z) = y i x j + z k;
c) F(x, y, z) = x i + y j + z k;
: x = cos t, y = sin t, z = t, 0 t 2
: x = cos t, y = sin t, z = 2, 0 t 2
: x = t, y = 2t, z = t2 1,
f ) F(x, y, z) = xy i + (z x) j + 2yz k;
(1, 0, 0) a (1, 2, 0) y hasta (1, 2, 2).
16. Verifique el Teorema del Rotacional para el campo vectorial F y la superficie dados.
a) F(x, y, z) = 2y i x j + z k;
b) F(x, y, z) = xy i + xz j + yz k;
: x2 + y 2 + z 2 = 1, z 0
: z = x2 + y 2 , z 1
147
18. Para F(r, , z) = r er + z sin e + rz ez en coordenadas cilndicas, encuentre div F y
rot F.
19. Para F(, , ) = e + cos e + e en coordenadas esfericas, encuentre div F y
rot F.
20. Demuestre las siguientes formulas siendo r = krk la norma del campo vectorial de
posicion r(x, y, z) = x i + y j + z k.
a) (1/r) = r/r3
b) (1/r) = 0
c) (r/r3 ) = 0
d ) (ln r) = r/r2
21. Demuestre las siguientes propiedades de la divergencia y el rotacional
a) div (F + G) = div F + div G
b) rot (F + G) = rot F + rot G
c) div (f F) = f div F + F f
d ) div (F G) = G rot F F rot G
e) div (f g) = 0
f ) rot (f F) = f rot F + (f ) F
g) rot (rot F) = (div F) F
h) (f g) = f g + g f + 2(f g)
22. Muestre que si un campo vectorial es conservativo, su rotacional es el vector cero.
23. Muestre que si F es un campo vectorial tridimensional, con segundas derivadas parciales continuas entonces la divergencia del rotacional de F es cero.
24. Para cada uno de los siguientes campos F muestre que no existe ning
un campo G
cuyo rotacional sea F.
a) F(x, y, z) = (2x + y)i + (y 2z)3/2 j + y 3 k
b) F(x, y, z) = cos(x 2z)i + sen( z)
x
3/2
j + y3 k
25. Si F es un campo vectorial tridimensional, con derivadas parciales continuas, determine la densidad rotacional del campo para:
a) Un punto arbitrario del plano XY
b) Un punto arbitrario del plano XZ
148
f
r
er + 1r
e +
RRR
(u v v u) dV =
RR
(u v v u) dA
LECCION
12
Teoremas Fundamentales
Dedicaremos esta u
ltima seccion a la demostraci
on de los teoremas fundamentales.
=
=
Pm Pn
i=1
j=1
Pm Pn
i=1
149
rv k uv,
150
donde ru y rv est
an evaluadas en Pij , es una aproximacion de la rotacion total de F sobre
, que no es otra cosa que la integral de R sobre :
Z bZ d
ZZ
R (x(u, v), y(u, v), z(u, v)) kru rv k dvdu
R dA =
a
b
a
Cij
que es la circulacion de F a lo largo de ij por
kru rv k uv
unidad de area. Entonces,
Tenemos que R (Pij )
R(R , m, n, r)
Pm P n
i=1
Pm P n
i=1
j=1
Pm P n
i=1
j=1
rv k uv
Cij
kru rv k uv
kru rv k uv
Cij .
1
2
3
La circulacion de F a lo largo de ij es igual a la suma de las circulaciones Cij
, Cij
, Cij
4
y Cij
de F a lo largo de cada una de las curvas correspondientes a los lados de Rij =
1
3
2
4
[ui1 , ui ] [vj1 , vj ]. Pero como Cij
= C(i+1)j
y Cij
= Ci(j+1)
, se tiene que
R(R , m, n, r)
=
=
=
Pm P n
j=1
i=1
Pm P n
1
j=1 (Cij
i=1
Pm
2
i=1 (Cim
Pm
i=1
Cij .
2
3
4
+ Cij
+ Cij
+ Cij
)
4
)+
+ Ci0
2
+
Cim
Pm
i=1
Pn
1
j=1 (Cnj
4
Ci0
+
Pn
j=1
3
+ C0j
)
1
Cnj
+
Pn
j=1
3
C0j
Pm
Pm
Pn
Pn
1
2
3
4
donde i=1 Cim
, i=1 Ci0
, j=1 Cnj
y j=1 C0j
son las sumas de Riemann de las circulaciones de F a lo largo de cada uno de las curvas correspondientes a los lados de la 2-celda
R = [a, b] [c, d] y que conforman el borde de la superficie . Por consiguiente:
RR
rotF dA
lmm,n0 R(R , m, n, r)
lmm,n0
F dL
P
m
i=1
2
Cim
+
Pm
4
i=1 Ci0 +
Pn
1
j=1 Cnj +
Pn
3
j=1 C0j
151
El teorema de la divergencia
Sea F un campo vectorial tridimensional definido en un conjunto abierto D
del espacio. Sea un solido encerrado por una superficie cerrada contenidos
en D. Suponga que esta parametrizado por una funcion vectorial
r : P R3 ,
R(divF, m, n, l, r) =
Pm Pn
j=1
i=1
Pl
k=1
divFdV
= lmm,n,l
=
Z bZ
a
d
c
h
e
Pm Pn
i=1
j=1
Pl
k=1
152
R(E , m, n, l, r) =
=
Pm Pn
j=1
i=1
Pm Pn
j=1
i=1
Pm Pn
j=1
i=1
Pl
Pl
k=1
Pl
k=1
rv ) rw | uvw
ijk
|(ru rv ) rw | uvw
|(ru rv ) rw | uvw
ijk .
R(E , m, n, l, r)
=
=
=
=
Pm Pn
i=1
j=1
Pm Pn
i=1
j=1
Pm Pn
i=1
Pl
k=1
Pl
3
j=1 (ijl
Pm Pn
i=1
k=1
P6
s=1
sijk
+ 4ij0 ) +
3
j=1 ijl +
Pm Pl
i=1
k=1
ijk
Pm Pl
i=1
Pm Pn
5i0k +
2
k=1 (ink
4
j=1 ij0 +
i=1
Pn
j=1
Pl
k=1
+ 5i0k ) +
Pm Pl
k=1
i=1
1mjk +
Pn
j=1
2ink
Pl
k=1
Pn
j=1
Pl
1
k=1 (mjk
+ 60jk )
60jk
153
ZZZ
divFdV
= lmm,n,l0 R(E , m, n, l, r)
= lmm,n,l0
P
+
=
m
i=1
Pn
j=1
3ijl +
Pm Pl
ZZ
i=1
5
k=1 i0k
Pm Pn
i=1
Pn
j=1
j=1
4ij0 +
Pl
1
k=1 mjk
P m Pl
k=1
i=1
Pn
j=1
2ink
Pl
6
k=1 0jk
F dA
Para terminar esta leccion, consideremos el Teorem de Stokes en el caso particular en el que la superficie este contenida en el plano XY . Obsrvese que en
este caso, el vector normal unitario en cada uno de los puntos de la superficie,
no es otro que el vector k. Por consiguiente, la densidad de rotacion del campo
F = hP, Q, Ri es rotF k = Qx Py . Por consiguiente, en este caso:
ZZ
ZZ
rotF kdA =
(Qx Py ) dA
y
F dL =
P dx + Qdy.
154
EJERCICIOS
1. Aplicando el teorema de la divergencia, muestre que si un campo es conservativo en todo el espacio, el flujo del campo a traves de cualquier
superficie cerrada orientable es cero.
2. Muestre que si un campo es conservativo en el espacio, excepto en un
punto, El flujo del campo atraves de cualquier par de superficies cerradas
que contengan al punto en su interior es igual.
3. Para cada uno de los siguientes campos vectoriales determine si existe
una funcion potencial, si existe encuentrela
a) F(x, y) = (y 2 + 3x2 ) i + 2xy j.
b) F(x, y) = (x3 cos(xy) + 2x sin(xy)) i + x2 y cos(xy) j.
c) F(x, y) = (8xy + 3) i + 4(x2 + y) j.
4. Use una integral de superficie para mostrar que el area de la superficie
de una esfera de radio r es 4r2 .
5. Use una integral para mostrar que el area lateral de un cono circular recto
de radio R y altura h es R h2 + R2 .
6. Si la superficie esta dada por z = f (x, y), con (x, y) en la region R2 ,
muestre que el area de esta dada por
ZZ r
2
f 2
dA .
S =
1 + x + f
y
7. Construya una banda de Mobius con una tira de papel, dibuje una lnea
a lo largo del centro y cortela siguiendo esta lnea. Cuantas superficies
resultan? Como las describira? Son orientables?
8. Muestre que el Teorema de Green es un caso especial del Teorema de
Stokes.
9. Suponga que u = 0 (es decir, u es arm
onica) sobre R3 . Defina la
u
derivada normal n de u sobre una superficie cerrada con vector normal
RR u
u
= Dn u = n u. Muestre que
d = 0.
unitario exterior n por n
n