Doce Lecciones de Calculo Vectorial Libr PDF

DOCE LECCIONES DE CÁLCULO VECTORIAL
Ernesto Acosta Gempeler

Bernarda Aldana Gómez
Escuela Colombiana de Ingenierı́a

Julio Garavito
ii
Índice general
1. Vectores en R2 y R3 1
2. Ecuaciones de superficies en el espacio 17
3. Funciones vectoriales 35
4. Funciones de varias variables 59
5. Diferenciación Parcial 73
6. Regla de la cadena 89
7. Extremos de funciones 101
8. Integrales múltiples 113
9. Integrales múltiples 129
10.Campos vectoriales 139
iii
iv ÍNDICE GENERAL
LECCIÓN 1
Vectores en R2 y R3
Hacemos un recuento muy rápido de los conceptos fundamentales de coordenadas de

puntos en el plano y en el espacio y de vectores bidimensionales y tridimensionales,
ası́ como de las operaciones entre vectores.
Para localizar puntos en el plano, se introducen dos ejes numéricos perpendiculares

entre sı́, cuyo punto O de intersección es el origen de los ejes numéricos, es decir, el
punto con coordenada cero en cada eje numérico. A uno de estos ejes lo designamos
eje X y al otro eje Y . Estos dos ejes numéricos definen en el plano un sistema de
coordenadas cartesianas rectangulares. Las coordenadas (x, y) de un punto P en el
plano se obtienen proyectándolo ortogonalmente sobre los ejes X y Y , siendo x la
coordenada de la proyección sobre el eje X y y la coordenada de la proyección sobre
el eje Y . La notación P (x, y) significa el punto P de coordenadas (x, y).
1
2 LECCIÓN 1. VECTORES EN R2 Y R3
y (x, y)
X
x
Figura 1.1
Dada una pareja de números reales (x, y), se ubica en el eje X el punto Px de coor-
denada x y en el eje Y el punto Py de coordenada y. La recta perpendicular al eje
X que contiene al punto Px 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 numéricos X y Y son los puntos Px y Py
con coordenadas x y y respectivamente.
Py P (x, y)
·
X
Px
Figura 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 numéricos
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
(x, y, z)
y Y
x
Figura 1.3
Dos puntos diferentes P y Q definen los segmentos de recta dirigidos PQ y QP.

Estos dos segmentos de recta se diferencian únicamente por su orientación. Es decir,
el segmento dirigido PQ se supone recorrido desde P hasta Q, y el segmento dirigido
QP se supone 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 números x2 − x1 , y2 − y1 y
z2 −z1 , son las componentes escalares del segmento dirigido PQ. Obsérvese que cada
una de las componentes del segmento dirigido QP es el opuesto de la componente
correspondiente del segmento dirigido PQ. Dos segmentos dirigidos son equivalentes
si sus componentes escalares correspondientes son iguales. No es difı́cil mostrar que
si P , Q, R y S son cuatro puntos no colineales, los segmentos dirigidos PQ y RS
son equivalentes si, y solamente si, el cuadrilátero P QSR es un paralelogramo.
La igualdad de las componentes escalares correspondientes entre segmentos dirigidos

es una relación de equivalencia y cada clase de equivalencia es un vector geométrico.
Es decir, un vector (en este caso tridimensional) está representado por un segmen-
to dirigido PQ, o por cualquier otro segmento dirigido equivalente. Ası́, un vector
está completamente determinado por las componentes escalares de cualquier seg-
mento dirigido que lo representa. Se usa frecuentemente la notación 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
v tiene componentes escalares hx, y, zi, escribiremos v = hx, y, zi.
Sea v el vector representado por el segmento dirigido OP, donde O es el origen del
sistema de coordenadas cartesianas. Si (x, y, z) son las coordenadas del punto P ,
entonces hx, y, zi son las componentes escalares del vector v , ya que las coordena-
das del punto O son (0, 0, 0). El punto P define un paralelepı́pedo rectangular de
caras paralelas a los planos coordenados, donde una de sus diagonales es el segmento
OP. Las dimensiones de esta caja rectangular son |x|, |y| y |z| y por consiguiente
p
la longitud del segmento OP es |OP| = x2 + y 2 + z 2 . Este último número repre-
senta además la magnitud kvk del vector v y la distancia d(O, P ) del punto O al
punto P . Ahora bien, si el vector v está representado por el segmento dirigido RS,
siendo R(x1 , y1 , z1 ) y S(x2 , y2 , z2 ), entonces el vector v tendrá componentes escalares
hx2 − x1 , y2 − y1 , z2 − z1 i y en consecuencia su magnitud será
p
kvk = (x2 − x1 )2 + (y2 − y1 )2 + (z2 − z1 )2 , (1.1)
este número representa también la longitud |RS| del segmento RS y la distancia

d(R, S) entre los puntos R y S.
Usando las componentes escalares de los vectores se definen las operaciones de adi-
ción, producto escalar y producto vectorial de vectores, ası́ como multiplicación por
escalar de la siguiente manera: Sean u =hu1 , u2 , u3 i y v =hv1 , v2 , v3 i dos vectores
tridimensionales y a un número real. La suma de u y v es
u + v = hu1 + v1 , u2 + v2 , u3 + v3 i , (1.2)
el producto escalar de u y v es
u · v = u1 v 1 + u2 v 2 + u3 v 3 , (1.3)
el producto vectorial de u y v es
u × v = hu2 v3 − u3 v2 , u3 v1 − u1 v3 , u1 v2 − u2 v1 i (1.4)
y la multiplicación del escalar a por el vector u es
au = hau1 , au2 , au3 i . (1.5)
Veremos más adelante una fórmula mnemotécnica para recordar el producto vecto-
rial.
5
Estas operaciones tienen propiedades bien conocidas, que se deducen de las pro-
piedades de la adición y multiplicación de los números reales, que son las que las
definen. Por ejemplo, es claro que la adición de vectores es asociativa, conmutativa,
que hay un elemento neutro (vector nulo) y que cada vector tiene su opuesto. Entre
otras propiedades se tiene que la multiplicación por escalar distribuye sobre la adi-
ción de vectores, el producto escalar de vectores es conmutativo y distribuye sobre
la adición de vectores. Sin embargo, aunque el producto vectorial distribuye sobre
la adición 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 dirección
y sentido. Los vectores de longitud 1 se denominan vectores unitarios. Los ejemplos
más conocidos de vectores unitarios en R3 son los vectores i = h1, 0, 0i, j = h0, 1, 0i,
y k = h0, 0, 1i.
Otra forma de definir el producto escalar de dos vectores, que normalmente se usa
cuando no están referidos a un sistema de coordenadas es:
u · v = kuk kvk cos θ (1.6)
donde θ es el menor ángulo formado por los dos vectores. De aquı́ podemos deducir
que si θ = π/2 se tiene que cos θ = 0 y por consiguiente u · v = 0, y que si u y v son
vectores no nulos y u · v = 0 entonces θ = π/2 . Es decir, dos vectores no nulos son
ortogonales si, y solo si, su producto escalar es cero.
Al proyectar ortogonalmente el vector u sobre el vector v se obtiene un vector de-

nominado el vector proyección, este vector tiene la misma dirección del vector v ,
su magnitud es |kuk cos θ| y el sentido depende del ángulo θ, si 0 < θ < π/2 tie-
ne el mismo sentido de v y si π/2 < θ < π tiene el sentido opuesto. Ası́ el vector
proyección de u sobre v es
v u·v
proyv u = kuk cos θ = v, (1.7)
kvk kvk2
u
u
v
v
proyv u
proyv u
Figura 1.4
al escalar kuk cos θ se le denomina la componente escalar de u en la dirección de v

u·v
compv u = kuk cos θ = . (1.8)
kvk
Otro hecho geométrico muy útil: Claramente el área del paralelogramo generado
por los vectores u y v es kuk kvk sen θ, donde θ es el ángulo entre u y v , y es fácil
mostrar que
ku × vk = kuk kvk sen θ, (1.9)
de aquı́ se deduce que dos vectores no nulos son paralelos si y sólo si su producto
vectorial es el vector nulo, que es el único vector cuya magnitud es cero. Además
si recordamos que u × v es ortogonal tanto a u como a v y que el volumen de un
paralelepı́pedo es igual al área de la base multiplicada por la altura , podemos
calcular el volumen del paralelepı́pedo generado por tres vectores u , v y w como:
±vol(u, v, w) = ku × vk kwk cos γ

(1.10)
= (u × v) · w
siendo γ el ángulo entre u × v y w .
Uno de los papeles más importantes de los sistemas de coordenadas es que nos
permiten escribir de manera algebraica las condiciones geométricas que definen con-
juntos de puntos en el plano o en el espacio. Un lugar geométrico de puntos es un
conjunto de puntos que está definido por una, o varias condiciones, por lo general de
carácter geométrico. Las condiciones geométricas se traducen en relaciones entre las
coordenadas de los puntos que las satisfacen y estas relaciones se escriben mediante
7
ecuaciones o inecuaciones en dos o tres variables.
Por ejemplo, la recta que contiene un punto dado P y es paralela a un vector v es el

lugar geométrico de puntos Q tales que el segmento dirigido PQ es paralelo a v. En
otras palabras, un punto Q está en esta recta si, y sólo si, su vector posición es la
suma de OP y un múltiplo escalar de v. Ası́, la recta que pasa por P y tiene vector
director v es el lugar geométrico de puntos Q tales que para algún número real t se
tiene que
OQ = OP + tv (1.11)
Si P tiene coordenadas (a, b, c) y v = hv1 , v2 , v3 i, un punto Q de coordenadas (x, y, z)

estará en la recta si, y sólo si, sus coordenadas satisfacen la ecuación
hx, y, zi = ha, b, ci + t hv1 , v2 , v3 i (1.12)
para algún número real t. Esta última es la ecuación vectorial de la recta que contiene
a P y tiene vector director v. Las ecuaciones paramétricas de esta recta son las tres
ecuaciones que se obtienen de la igualdad de los dos vectores en la ecuación anterior
x = a + tv1 , y = b + tv2 , z = c + tv3 . (1.13)
A partir de las ecuaciones paramétricas se pueden escribir las ecuaciones simétricas

al establecer las relaciones que se obtienen entre las coordenadas de los puntos de la
recta eliminando el parámetro. En el caso en que v1 , v2 y v3 son diferentes de cero,
estas son
x−a y−b z−c
= = . (1.14)
v1 v2 v3
En realidad, estas son tres ecuaciones (¿cuáles?) que corresponden a planos que
contienen a la recta, cada uno de los cuales es paralelo a uno de los ejes coordenados.
EJERCICIOS
1. Para cada par de puntos P y Q 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. En un plano cartesiano represente los puntos P (−3, −5) y Q(2, −1) y dibuje:
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(7, −2).
3. En cada numeral determine si los segmentos PQ y RS son equivalentes, jus-

tifique su respuesta.
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).
4. 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 vértices de un paralelogramo.
5. Dados los puntos P (p1 , p2 , p3 ), Q(q1 , q2 , q3 ) y R(r1 , r2 , r3 ), determine todos los

puntos S tales que P, Q, R y S sean los vértices de un paralelogramo.
6. Muestre que la suma de vectores es conmutativa.
7. Muestre que el producto escalar de dos vectores es conmutativo.
8. Muestre que el producto vectorial no es conmutativo ni asociativo. ¿Qué rela-

ción existe entre los vectores u × v y v × u?
9. 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 .
9
10. Sean v = h−1, 5, −2i y w = h3, 1, 1i. Calcule

v
(a) v − w. (b) v + w. (c) kvk .
(d) 12 (v − w). (e) 12 (v + w). (f) −2 v + 4 w.
(g) v − 2 w.
11. Abrir la mano derecha, doblar los dedos meñique y anular, poner el dedo pul-
gar perpendicular al ı́ndice y al corazón. Ası́ se ponen los dedos ı́ndice, corazón
y pulgar (en ese orden), orientados de acuerdo con la regla de la mano dere-
cha. Una tripla de vectores (u, v, w) en R3 están orientados de acuerdo con
la regla de la mano derecha, si sus posiciones relativas están de acuerdo con
las posiciones relativas de los dedos ı́ndice, corazón y pulgar de la mano de-
recha descritas anteriormente. Este comportamiento relativo de los vectores
se define desde el punto de vista matemático de la siguiente manera: Una tri-
pla de vectores (u, v, w) en R3 están orientados de acuerdo con la regla de la
mano derecha si y sólo si det(u, v, w) > 0. Por ejemplo, la tripla de vectores
(i, j, k), representados en un sistema cartesiano dextrógiro, están orientados de
acuerdo con la regla de la mano derecha porque det(i, j, k) = 1 > 0. Pero si
intercambiamos el orden de i y j tenemos que det(j, i, k) = −1 < 0 por lo tanto
estos vectores no están orientados de acuerdo con la regla de la mano derecha,
tendremos que cambiar a k por −k, para que la nueva tripla esté orientada
por la regla de la mano derecha, esto es det(j, i, −k) = 1 > 0.
El producto vectorial de dos vectores también se define como: Un vector que

es perpendicular a los dos vectores dados, que tiene la orientación 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 ángulo formado por ellos. Muestre que esta
definición es equivalente a la dada en (1.4).
12. Dados el punto C(1, 2, 5), determine las coordenadas de 3 puntos del espacio
tales que la distancia de cada uno de ellos al punto C sea 4 unidades. ¿Exis-
ten más puntos del espacio que cumplan la condición? ¿Existe un punto con
primera coordenada 3 que cumpla la condición? ¿Existe un punto con primera
componente 8 que cumpla la condición? ¿Cómo sabrı́a usted si un punto de
coordenadas (x, y, z) satisface la condición?
13. Dado el punto C(a, b, c), determine las condiciones que deben satisfacer las
coordenadas de todos los puntos S(x, y, z) del espacio para que kCSk = 4.
Describa el lugar geométrico de dichos puntos.
14. Dados los puntos P (3, −1, 4), Q(9, −2, 3), R(1, 4, 7), determine 4 puntos tales
que el vector PQ sea paralelo al vector cuyo punto inicial es R y punto final
el punto determinado. ¿Existe algún punto con primera coordenada 3 que
cumpla la condición? ¿Existe algún punto con primera coordenada 20 que
cumpla la condición? ¿Existe algún punto con primera coordenada 3 y segunda
coordenada 5 que cumpla la condición? ¿Existen más puntos del espacio que
cumplan la condición? ¿Cómo sabrı́a usted si un punto de coordenadas (x, y, z)
satisface la condición?
15. Dados 3 punto fijos P , Q y R, describa el lugar geométrico de todos los puntos
S del espacio tales que los vectores PQ y RS sean paralelos. Si P (x1 , y1 , z1 ),
Q(x2 , y2 , z2 ), R(a, b, c), determine las condiciones que deben satisfacer las coor-
denadas de todos los puntos S(x, y, z) del espacio para que los vectores PQ y
RS sean paralelos.
16. Dados 3 punto fijos P , Q y R, describa el lugar geométrico de todos los pun-
tos S del espacio tales que los vectores PQ y RS sean perpendiculares. Si
P (x1 , y1 , z1 ), Q(x2 , y2 , z2 ), R(a, b, c), determine las condiciones que deben sa-
tisfacer las coordenadas de todos los puntos S(x, y, z) del espacio para que los
vectores PQ y RS sean perpendiculares.
17. Determine el área del paralelogramo generado por los vectores u = h2, 1, −1i
y v = h1, 1, 4i.
18. Determine el área del paralelogramo que tiene vértices en los puntos (0, 7, 4),
(−2, 3, −6) y (1, −4, 5).
19. Determine el área del triángulo que tiene vértices en los puntos (0, 7, 4), (−2, 3, −6)
y (1, −4, 5) .
20. Determine el volumen del paralelepı́pedo generado por los vectores u = h2, 1, −1i,
v = h1, 1, 4i y w = h1, 3, 2i.
21. ¿Cuántos paralelepı́pedos diferentes pueden tener cuatro de sus vértices en los
puntos (1, 2, −1), (−5, 2, 2), (3, 1, 4) y (−1, 1, −3)? Determine el volumen de
11
por lo menos dos de ellos y justifique por qué se puede afirmar que todos tie-
nen el mismo volumen. ¿Es alguno de estos paralelepı́pedos un paralelepı́pedo
recto?
22. Sea P el paralelepı́pedo en R3 definido por
P = {(x, y, z) : 0 ≤ x ≤ a, 0 ≤ y ≤ b, 0 ≤ z ≤ c}.
a) Dibuje el paralelepı́pedo en un sistema de coordenadas cartesianas.

b) Determine los cosenos de los ángulos formados por la diagonal principal
de P y cada una de sus aristas y encuentre una identidad que los relacione.
c) Determine los ángulos entre la diagonal principal de P y cada una de sus
caras (¿Cómo se define el ángulo entre un vector y un plano?).
d ) Para cada uno de las caras del paralelepı́pedo, determine el punto medio,
tome el segmento dirigido que va del origen a dicho punto y determine
el coseno del ángulo formado por la diagonal principal y el segmento
dirigido.
e) Si el paralelepı́pedo no fuera recto, ¿cómo determinarı́a los cosenos de los
ángulos formados por cada una de las caras y la diagonal principal?
f ) ¿Cuántos paralelepı́pedos diferentes pueden tener la misma diagonal?
23. ¿En qué posición se encontrará un objeto si a partir del punto P (4, −7) se ha
desplazado 7 unidades en la dirección del vector h−1, 3i?
24. ¿En qué posición se encontrará un objeto si a partir del punto P (2, −6, 4) se
ha desplazado 4 unidades en la dirección del vector h2, −3, 6i?
25. ¿Qué significa que los vectores i =h1, 0, 0i, j =h0, 1, 0i y k =h0, 0, 1i forman una
base ortonormal del espacio de vectores tridimensionales? ¿Se puede escribir
todo vector de R3 como una combinación lineal de i y j únicamente? Justifique
su respuesta.
26. Sean v = h−1, 5, −2i y w = h3, 1, 1i.
a) Determine el vector u tal que u + v + w = i.

b) Determine el vector u tal que u + v + w = 2 j + k.
c) ¿Existe un escalar m tal que m(v + 2 w) = k?
27. Calcule v · w, v × w, la componente escalar de v en la dirección de w y la

proyección ortogonal de v sobre w, para cada uno de los pares de vectores v
y w siguientes:
a) v = h5, 1, −2i, w = h4, −4, 3i . b) v = − i+2 j+k, w = −3 i+6 j+3 k.
c) v = h1, 3, 2i, w = h7, 2, −10i . d ) v = i, w = 3 i + 2 j + 4k.
28. Para cada par de vectores determine el ángulo θ formado por ellos:
a) v = h5, 1, −2i, w = h4, −4, 3i . b) v = h7, 2, −10i, w = h2, 6, 4i .
c) v = h2, 1, 4i, w = h1, −2, 0i . d ) v = i, w = 3 i + 2 j + 4k
29. En cada numeral determine si v y w son o no perpendiculares. Justifique su

respuesta:
a) v = h8, 4, 3i y w = h−2, 1, 4i.

b) v = h6, 0, 4i y w = h0, 2, −1i.
30. Para cada una de las siguientes afirmaciones, demuéstrela si es verdadera o

dé un contraejemplo si es falsa:
a) Si u · v = u · w, entonces v = w.
b) Si v · w = 0 para todo v, entonces w = 0.
c) Si u · v = u · w para todo u, entonces v = w.
d ) kvk − kwk ≤ kv − wk para todo v y w.
e) Si v × w = 0 para todo w en R3 , entonces v = 0.
f ) Si u y v son unitarios, entonces ku × vk = 1
31. Sean α, β y γ los ángulos entre un vector no nulo v en R3 y los vectores i, j,

y k, respectivamente. Muestre que cos2 α + cos2 β + cos2 γ = 1.
(Nota: α, β, γ son los ángulos directores de v, y cos α, cos β, cos γ son los
cosenos directores de v.)
13
32. En cada numeral 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 = h−1, 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 tendrı́a que

v × w fuera también un vector unitario en R3 ? Justifique su respuesta.
35. 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.
c) Los vectores u = kvk w + kwk v y z = kvk w − kwk v son ortogonales.
d ) El vector u = kvk w+kwk v biseca al ángulo formado por los dos vectores.
36. Encuentre todas las soluciones de la ecuación vectorial a × x = b en R3 ,

donde a 6= 0.
37. Muestre que para todos los vectores u y v en R3 , se tiene que
u × (v × u) = (u · u)v − (u · v)u.
38. Demuestre la identidad de Jacobi :
u × (v × w) + v × (w × u) + w × (u × v) = 0.
39. Se dice que tres vectores u, v, w están en un mismo plano si los puntos P , Q y
R que tienen respectivamente a u, v, w como vectores posición, son coplanares
junto con el origen O. Demuestre que u, v, w están en un mismo plano si, y
sólo si, u, v, w son linealmente dependientes.
40. Muestre que u, v, w están en un mismo plano en R3 si y solo si u · (v × w) = 0.
41. Muestre que para todos 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 geométrico estas propiedades?
42. Muestre que si u, v, w, z son vectores no nulos en R3 que están en un mismo

plano entonces (u × v) × (w × z) = 0
43. ¿Es verdad que si u, v, w, z son vectores no nulos en R3 y (u×v)×(w×z) = 0

entonces están en el mismo plano? Justifique su respuesta.
44. Use vectores para demostrar que:
a) Los puntos medios de los lados de cualquier cuadrilátero son vértices de

un paralelogramo.
b) Si A, B y C son los vértices de un triángulo, A′ , B ′ y C ′ son los puntos
medios de los lados del triángulo y P es un punto exterior del triángulo
entonces PA + P B + PC = PA’ + PB’ + PC’.
c) Las diagonales de un rombo son perpendiculares.
45. Use vectores para demostrar la ley de los cosenos y la ley de los senos.
46. Determine el trabajo realizado por la fuerza F = 3i − j + 2k al mover un objeto

a lo largo del segmento dirigido OP = −i + 2j + ek
47. ¿Cuáles son los subespacios vectoriales propios de R2 ? ¿Cuáles son los de R3 ?
Nota. No siempre somos concientes de la dualidad presente en R2 y R3 , ya que
algunas veces son los conjuntos de puntos, del plano o del espacio, respectiva-
mente, y algunas veces son los espacios vectoriales de vectores bidimensionales
y tridimensionales, respectivamente. Podrı́amos en esta circunstancia, hacer
un paralelo con lo que acontece en el libro El extraño caso del Dr. Jekyll y Mr.
Hyde del escritor inglés Robert Louis Stevenson.
48. Encuentre la matriz asociada a la transformación lineal de R3 en R3 , que

corresponde a la reflexión con respecto al plano de ecuación x = y. Encuentre la
que corresponde a la reflexión con respecto al plano de ecuación ax+by+cz = 0.
¿Es la reflexión con respecto al plano de ecuación ax + by + cz = d, d 6= 0, una
transformación lineal de R3 en R3 ? Recuerde a Dr. Jekyll y a Mr. Hyde.
15
49. En cada numeral escriba la ecuación vectorial, las ecuaciones paramétricas y

las ecuaciones simétricas de la recta L que contiene el punto P y es paralela
al vector v:
a) P (2, 3, −2), v = (5, 4, −3)

b) P (3, −1, 2), v = (2, 8, 1)
c) P (2, 1, 3), v = (1, 0, 1)
d ) P (0, 0, 0), v = (7, 2, −10)
50. Escriba las ecuaciones paramétricas de la recta L que contiene los puntos P1
y P2 :
a) P1 (1, −2, −3), P2 (3, 5, 5)

b) P1 (4, 1, 5), P2 (−2, 1, 3)
51. Explique claramente el procedimiento solicitado en cada uno de los siguientes

numerales, suponiendo que se tiene las ecuaciones correspondientes:
a) ¿Cómo de determina la distancia de un punto a una recta?

b) ¿Cómo se determina la distancia entre dos rectas paralelas?
c) ¿Cómo se determina la distancia entre dos rectas que no son paralelas y
no se cortan?
d ) ¿Cómo se determina la intersección entre dos rectas?
52. Calcule la distancia del punto P a la recta L.
a) P (1, −1, −1), L : x = −2 − 2t, y = 4t, z = 7 + t

b) P (0, 0, 0), L : x = 3 + 2t, y = 4 + 3t, z = 5 + 4t
53. Determine las coordenadas del punto de intersección (si lo hay) de las rectas
dadas.
a) x = 7+3s, y = −4−3s, z = −7−5s y x = 1+6t, y = 2+t, z = 3−2t

x−6 x − 11 y − 14 z+9
b) =y+3=z y = =
4 3 −6 2
54. Determine la distancia entre las rectas dadas.
a) x = 6 + s, y = −2 − 3s, z = 4 − 5s y x = 1 + 2t, y = 2 + t, z = 3 − 4t
x−6 x−3 y−6 y−1z+9
b) =y+3=z y = =
4 8 2 2 2
LECCIÓN 2
Ecuaciones de superficies en el espacio
La introducción de un sistema de coordenadas permite representar algebraicamente

lugares geométricos de puntos mediante ecuaciones e inecuaciones. El procedimiento
consiste en encontrar las condiciones algebraicas que deben satisfacer las coordenadas
de un punto para pertenecer al lugar. La representación algebraica de un lugar
geométrico de puntos del espacio se puede hacer mediante ecuaciones, inecuaciones
o un sistema de ecuaciones e inecuaciones con tres incógnitas (las coordenadas de
los puntos). Más precisamente, una ecuación con tres incógnitas representa un lugar
geométrico de puntos del espacio si, y sólo si, las coordenadas de cada punto del
lugar satisfacen la ecuación y cada tripla de números reales que satisface la ecuación
son las coordenadas de un punto del lugar.
Por ejemplo, el plano que pasa por un punto Q y tiene vector normal v , es el
lugar geométrico de todos los puntos P del espacio para los que los vectores QP y
v son ortogonales. Si las coordenadas de Q son (x0 , y0 , z0 ), las de P son (x, y, z) y
v = ha, b, ci, la condición de ortogonalidad, es decir, el que QP · v= 0, se escribe en
coordenadas ası́:
a(x − x0 ) + b(y − y0 ) + c(z − z0 ) = 0 (2.1)
que es la ecuación que representa al plano que contiene a Q(x0 , y0 , z0 ) y tiene vector
normal v =ha, b, ci.
17
18 LECCIÓN 2. ECUACIONES DE SUPERFICIES EN EL ESPACIO
El lugar geométrico de todos los puntos P (x, y, z) del espacio cuya distancia a un
punto fijo Q(x0 , y0 , z0 ) es un número real positivo r, es la superficie de una esfera
de radio r y centro Q. Como lo vimos en la lección anterior, la condición geométrica
“la distancia de P a Q es r”, se escribe algebraicamente ası́:
p
(x − x0 )2 + (y − y0 )2 + (z − z0 )2 = r, (x − x0 )2 + (y − y0 )2 + (z − z0 )2 = r2
o,
(2.2)
Esta es precisamente la ecuación que representa dicha superficie.
Y
•
Γ
Figura 2.1
Sea Γ una curva en un plano Π, la unión de todas las rectas perpendiculares al

plano Π que intersectan a la curva Γ es la superficie cilı́ndrica recta(o cilı́ndro) con
directriz C. En otras palabras, la superficie cilı́ndica recta definida por Γ es el lugar
geométrico de puntos P del espacio cuya proyección ortogonal sobre el plano Π
está en la curva Γ. Por ejemplo, dada la circunferencia Γ en el plano XY con centro
en el origen y radio r, la unión de todas las rectas perpendiculares al plano XY que
intersectan la circunferencia Γ es una superficie cilı́ndrica circular recta. Un punto
P (x, y, z) pertenece a esta superficie si, y solamente si, su proyección ortogonal sobre
el plano XY es un punto de la circunferencia. Esto sucede únicamente cuando las
coordenadas x y y del punto P satisfacen la ecuación de la circunferencia Γ. Ası́, el
19
punto P (x, y, z) está en la superficie cilı́ndrica recta definida por Γ si y solamente

si, x2 + y 2 = r2 , que es precisamente la ecuación que representa a la superficie
en cuestión. Obsérvese que esta ecuación no impone ninguna restricción sobre la
coordenada z del punto P .
La superficie de la esfera y la superficie cilı́ndrica circular recta son superficies
cuadráticas (o cuádricas) ya que las ecuaciones que las representan son cuadráti-
cas. En general, si la ecuación de una superficie tiene forma
Ax2 + By 2 + Cz 2 + Dxy + Exz + F yz + Gx + Hy + Iz + J = 0, (2.3)
donde las letras mayúsculas son números reales fijos con la condición de que A, B,
C, D, E, F no sean todos nulos, se dice que la superficie es una superficie cuadrática.
Escogiendo convenientemente los ejes de coordenadas, es decir, transladando y ro-
tando los ejes de coordenadas originales, se puede demostrar que cualquier superficie
cuadrática admite una de las dos formas de ecuación siguientes:
ax2 + by 2 + cz 2 = d, con d = 1, −1, o 0 (2.4)
o
ax2 + by 2 = dz, con d = 1, −1, o 0. (2.5)
Por ejemplo, xy−1 = 0 es la ecuación de una superficie cilı́ndrica recta cuya directriz
es la hipérbola en el plano XY de igual ecuación, pero rotando 45 grados los ejes X
y Y en el sentido contrario al movimiento de las manecillas del reloj y dejando fijo
el eje Z, la ecuación de la misma superficie es x̃2 − ỹ 2 = 2, en la nuevas coordenadas
x̃2 − ỹ 2 x̃2 + ỹ 2
X̃ Ỹ Z. Se ha usado la sustitución x = √ ,y= √ , que nos dá las nuevas
2 2
coordenadas a partir de las primeras.
Estas superficies se clasifican según los signos de los coeficientes. Por ejempo, si los
coeficientes en (2.4) son todos distintos de cero, la superficie será un elipsoide si a, b,
c y d son positivos; un hiperboloide elı́ptico de una sola hoja, si a, b y d son positivos
y c es negativo; un hiperboloide elı́ptico de dos hojas, si a y b son positivos y c y d son
negativos; una superficie cónica si a y b son positivos, c es negativo y d es cero. Si los
coeficientes en (2.5) son distintos de cero, la superficie será un paraboloide elı́ptico si
a y b tienen el mismo signo y un paraboloide hiperólico si a y b son de signos contra-
rios. ¿Qué superficies serán aquellas en las que coeficientes en (2.4) y (2.5) se anulan?
Una forma de representar gráficamente una superficie a partir de su ecuación, es bus-

cando información de los cortes de la superficie con planos (frecuentemente paralelos
a los planos coordenados). Por ejemplo, si se intersecta la superficie cuadrática de
ecuación z = x2 + y 2 con un plano de ecuación z = k, k > 0, se obtiene la ecuación
x2 + y 2 = k que para cada valor de k es la ecuación 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 ecuación x = k, se obtiene z = y 2 + k 2 que para cada valor de k
es la ecuación de una parábola en el plano x = k. Es decir, la intersección de la
superficie con el plano x = k es una parbola que abre en el sentido positivo del eje
Z, al intersectar la superficie con un plano de ecuación y = k también se obtiene
una parábola. Dibujando simultaneamente estas circunferencias y parábolas para
diferentes valores de k, se obtiene un bosquejo de la superficie en el espacio, que es
un paraboloide circular con eje de simetrı́a el eje Z.
Las curvas de intersección obtenidas con este procedimiento se conocen como trazas.
Las trazas que se obtienen con planos paralelos a los planos coordenados se suelen
llamar trazas paralelas a los planos coordenados, las paralelas al plano XY , trazas
horizontales y las paralelas a los plano XZ o a Y Z, se suelen llamar trazas verticales.
Las trazas horizontales de la superficie de ecuación z = x2 + y 2 están determinadas
por el siguiente sistema de ecuaciones:
x2 + y 2 = k,
(2.6)
z = k.
Y
k = 4
k = 3
k = 2
k = 1
k= 1
4
Figura 2.2
21
Las trazas verticales están determinadas por los siguiente sistemas de ecuaciones:
z = k2 + y2 , z = x2 + k 2 ,
y y (2.7)
x = k, y = k,
respectivamente.
Z Z
√ √
k= 3 k= 3
√ √
k= 2 k= 2
k=1 k=1
k=0
Y k=0
X
Figura 2.3 Figura 2.4
Y
•
C
Figura 2.5
Coordenadas cilı́ndricas y coordenadas esféricas

Además del sistema de coordenadas cartesianas para el espacio, comúnmente se usan

el sistema de coordenadas cilı́ndricas y el sistema de coordenadas esféricas.
Un sistema de coordenadas cilı́ndricas se introduce fijando un plano en el espacio

en el que se considera un sistema de coordenadas polares y un eje numérico W
perpendicular a ese plano de tal forma que el origen del eje coincida con el polo de
las coordenadas polares. Las coordenadas cilı́ndricas de un punto P serán (r, θ, w)
donde (r, θ) son las coordenadas polares de la proyección ortogonal de P sobre el
plano, y w es la coordenada de la proyección ortogonal de P sobre el eje numérico.
W
eje cilı́ndrico
(r, θ, w)
r
ar
θ
l
po
e
ej
Figura 2.6
Obviamente, la representación 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 representación usando un sistema de coordenadas
cilı́ndricas se puede hacer escogiendo cualquiera de los planos XY , XZ o Y Z para
poner coordenadas polares (de manera conveniente) y el eje numérico restante, Z, Y
o X, respectivamente, para la tercera coordenada. En el primer caso, por ejemplo,
haciendo coincidir el eje polar con el eje positivo X, se puede establecer la relación
existente entre las coordenadas (x, y, z) y las coordenadas (r, θ, w) de un punto P .
En efecto, basta recordar la relación que existe entre las coordenadas cartesianas y
las coordenadas polares de un punto en el plano, bajo las condiciones establecidas.
23
Ası́, se tiene que
p
x = r cos θ, r ∈ [0, ∞) r = x2 + y 2
y
y = r sen θ, θ ∈ [0, 2π) y θ = arctan (2.8)
x
z = w, w ∈ (−∞, +∞) w = z.
Teniendo en cuenta esta equivalencia entre las coordenadas de los puntos en los dos
sistemas , la ecuación cartesiana de la superficie 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
ecuación cilı́ndrica de la superficie de la esfera r2 + z 2 = a2 .
Por otro lado, si se escoge poner el sistema polar en el plano XZ y el eje polar con
el semieje positivo X, se tiene que
x = r sen θ
y = w (2.9)
z = r cos θ.
En este caso la ecuación cilı́ndrica de la superficie de la esfera es y 2 + r2 = a2 .
Un sistema de coordenadas esféricas en el espacio se introduce fijando un semiplano,

el semiplano esférico, delimitado por una recta numérica, eje esférico, cuyo origen
es el punto O. Las coordenadas esféricas (ρ, θ, φ) de un punto P , se obtienen de-
terminando la distancia ρ del origen O al punto P , el ángulo diedro θ (en sentido
dextrógiro) formado por el semiplano esférico y el semiplano delimitado por el eje
esférico y que contiene al punto P , y el ángulo φ que forma el segmento dirigido OP
con el semieje esférico positivo.
eje esférico
φ •P (ρ, θ, φ)
ρ
O
co
ri
θ
fé
es
no
la
ip
m
se
Figura 2.7
Al igual que en el caso anterior, se puede establecer una equivalencia entre las coor-
denadas cilı́ndricas y las coordenadas esférica de un punto en el espacio. Escogiendo,
por ejemplo, el eje W como eje esférico, el semiplano delimitado por el eje W que
contiene el eje polar como semiplano esférico y el origen del sistema como el polo O,
la relación entre las coordenadas (r, θ, φ) y las coordenadas (ρ, θ, φ) de un punto P
se obtiene de la siguiente manera. De acuerdo con la siguiente gráfica
Q r
•P (ρ, θ, φ)
w φ
ρ
O
Figura 2.8
25
el triángulo OQP es rectángulo y el ángulo QOP es φ por lo tanto w = |OQ| =

ρ cos φ, y r = |QP | = ρ sen φ. Obsérvese que en estas circunstancias el ángulo θ es
el mismo en los dos sistemas de coordenadas. En resumen,
√
r = ρ sen φ, ρ ∈ [0, ∞) ρ = r2 + w2
θ = θ, θ ∈ [0, 2π) y θ = θ (2.10)
r
w = ρ cos φ, φ ∈ [0, π] φ = arctan .
w
Teniendo en cuenta lo establecido antes entre los sistemas de coordenadas cilı́ndri-
cas y cartesianas obtenemos la siguiente relación entre las coordenadas esféricas y
cartesianas:
p
ρ = x2 + y 2 + z 2
x = ρ sen φ cos θ, y
y θ = arctan (2.11)
y = ρ sen φ sen θ, x
p
2
x +y 2
z = ρ cos φ, φ = arctan .
z
EJERCICIOS
1. Determine la ecuación de la superficie de la esfera con centro en el origen y

radio 2.
2. Determine la ecuación de la superficie de la esfera con centro en el punto

C(3, 2, 4) y radio 5.

C(−1, 2, −5) y que contiene al origen.

(−1, 3, −2) y tangente al plano XY .

(−1, 3, −2) y tangente al plano Y Z.

(1, 3, 2) y tangente a la superficie de la esfera de ecuación x2 − 2x + y 2 + 4y +
z 2 + 2z + 2 = 0.
7. En cada numeral determine si la ecuación corresponde o no a la superficie de

una esfera. Si lo hace, encuentre su centro y su radio:
a) x2 + y 2 + z 2 − 4x − 6y − 10z + 37 = 0
b) x2 + y 2 + z 2 + 2x − 2y − 8z + 19 = 0
c) 2x2 + 2y 2 + 2z 2 + 4x + 4y + 4z − 44 = 0
d ) x2 + y 2 − z 2 + 12x + 2y − 4z + 32 = 0
e) x2 + y 2 + z 2 − 4x + 6z = 3.
f ) 2x2 + 2y 2 + 2z 2 + 4x − y + 3z = 0.
8. En cada numeral escriba una ecuación del plano que contiene al punto Q y es
perpendicular al vector n:
a) Q(5, 1, −2), n = h4, −4, 3i b) Q(6, −2, 0), n = h2, 6, 4i
9. En cada numeral escriba una ecuación del plano que contiene los puntos dados:
27
a) (1, 0, 3), (1, 2, −1), (6, 1, 6) b) (−3, 1, −3), (4, −4, 3), (0, 0, 1)
10. Determine la ecuación del plano paralelo al plano de ecuación x − 2y + 7z = 4

y que contiene al punto (7, 3, −4).
11. ¿Es posible determinar la ecuación de un plano que contiene al punto (2, −1, −3),
es perpendicular al plano de ecuación x − y + 2z = 5 y al vector h−1, 8, 4i?
12. Determine la ecuación del plano que contiene al punto (2, −1, −3), es perpen-
dicular al plano de ecuación x − y + 2z = 5 y al vector h0, 8, 4i.
13. En cada numeral calcule la distancia del punto Q al plano P.
a) Q(4, 1, 2), P: 3x − y − 5z + 8 = 0
b) Q(0, 2, 0), P: −5x + 2y − 7z + 1 = 0
14. Determine la ecuación vectorial de la recta de intersección (si la hay) de los

planos dados:
a) x + 3y + 2z − 6 = 0, 2x − y + z + 2 = 0
b) 3x + y − 5z = 0, x + 2y + z + 4 = 0
15. Determine las coordenadas de los puntos de intersección (si los hay) de la recta
x−6
= y + 3 = z con el plano x + 3y + 2z − 6 = 0.
4
16. Las siguientes ecuaciones representan superficies cuadráticas, para cada una
dibuje las trazas y a partir de ellas haga un bosquejo de la gráfica. De los
siguientes nombres deduzca cuál le corresponde a cada ecuación: Elipsoide,
Paraboloide, Paraboloide Hiperbólico, Hiperboloide de una hoja, Hiperboloide
de dos hojas, Cono de doble hoja.
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
x2 y 2 z 2
f) 2 + 2 + 2 = 1
a b c
17. Determine las coordenadas de los puntos de intersección de la superficie de la

esfera (x − 3)2 + (y + 1)2 + (z − 3)2 = 9 y la recta x = −1 + 2t, y = −2 − 3t,
z = 3 + t.
18. Describa la intersección de las superficies de las esferas (x − 4)2 + (y + 2)2 +

(z − 4)2 = 9 y x2 + y 2 + z 2 = 9.
19. Describa la intersección de la superficie de la esfera x2 + y 2 + z 2 = 9 y el

cilindro x2 + y 2 = 4, represéntela mediante ecuaciones.
x2 y2 z2
20. Determine la traza del hiperboloide de una hoja a2
+ b2
− c2
= 1 en el plano
x = a, y la traza en el plano y = b.
x2 y2 z
21. Determine la traza del paraboloide hiperbólico a2
− b2
= c en el plano XY .
22. Dados tres puntos, no colineales, P , Q, R en un plano describa un procedi-

miento que le permita determinar el centro y el radio de la circunferencia que
los contiene.
23. Determine el centro y el radio de la circunferencia que contiene los puntos,

(con coordenadas rectangulares), P (2, 3), Q(1, 5) y R(3, 0).
24. Dados cuatro puntos, no coplanares, P , Q, R y S, describa un procedimiento

geométrico que le permita determinar el centro y el radio de la esfera que los
contiene.
25. Determine la ecuación de la superficie de la esfera que contiene los puntos

(con coordenadas cartesianas) P (1, −1, 2), Q(0, 0, 2), R(1, −4, 3) y S(0, −1, 3)
utilizando el procedimiento descrito en el numeral anterior.
26. Muestre que el hiperboloide de una sola hoja es una superficie doblemente
reglada, es decir, cada punto de la superficie está en dos rectas contenidas
completamente en la superficie.
27. Muestre que el paraboloide hiperbólico es una superficie doblemente reglada.
28. Sea S la superficie de la esfera con radio 1 y centro (0, 0, 1), y sea S ∗ , la su-
perficie de la esfera S sin el punto N (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 algún punto
29
Q(x, y, 0). Encuentre este punto Q en términos de a, b y c. El punto Q es la

proyección estereográfica de P . Esta correspondencia “identifica” S ∗ con R2 .
29. Cuando se está acostumbrado a trabajar exclusivamente en sistemas de coorde-

nadas cartesianas se dificulta pensar en otros sistemas independientemente de
los primeros y se tiende a pasar de forma mecánica, por ejemplo, de un sistema
de coordenadas polares a un sistema de coordenadas cartesianas y viceversa,
sin apreciar el potencial del sistema de coordenadas polares para resover pro-
blemas particulares. Se propone hacer los siguientes ejercicios sin cambiar a
otro sistema de coordenadas, para darle sentido al sistema de coordenadas
polares independientemente de los sistemas de coordenadas cartesianas.
a) Determine la distancia entre los puntos P (r1 , θ1 ) y Q(r2 , θ2 ), donde r1 6= 0

y r2 6= 0.
b) Dados los puntos P (r1 , θ1 ) y Q(r2 , θ2 ), determine las coordenadas polares
del punto R tal que el segmento dirigido OR tenga la misma magnitud,
la misma dirección y el mismo sentido que el segmento dirigido PQ.
c) ¿Qué condición deben satisfacer las coordenadas polares r y θ de un punto
P para que pertenezca a la recta perpendicular al eje polar que pasa por
el punto de coordenadas polares (3, 0)?
d ) ¿Qué condición deben satisfacer las coordenadas polares r y θ de un
punto P para que pertenezca a la recta paralela al eje polar que pasa por
el punto de coordenadas poplares (1, π/4)?
e) ¿Qué condición deben satisfacer las coordenadas polares r y θ de un punto
P para que pertenezca a la circunferencia de radio a cuyo centro tiene
coordenadas polares (a, 0)?
f ) ¿Qué condición deben satisfacer las coordenadas polares r y θ de un punto
P para que pertenezca a la circunferencia de radio a y cuyo centro tiene
coordenadas polares (a, π/2)?
g) ¿Qué condición deben satisfacer las coordenadas polares r y θ de un punto
P para que pertenezca a la circunferencia de radio b y cuyo centro tiene
coordenadas polares (a, α), con b < a?
30. Encuentre las coordenadas (a) cilı́ndricas y (b) esféricas de los puntos cuyas
coordenadas cartesianas se dan, explicando de qué manera hizo coincidir los
dos sistemas:
√
a) (2, 2 3, −1) b) (−5, 5, 6)
√ √ √
c) ( 21, − 7, 0) d ) (0, 2, 2)
31. Dado un sistema de coordenadas en el espacio donde las coordenadas de un

punto son (α, β, γ), las superficies de ecuación α = c1 , β = c2 o γ = c3
se denominan superficies coordenadas y la intersección entre dos superficies
coordenadas se denominan curvas coordenadas. Para cada uno de los sistemas
estudiados describa las superficies coordenadas y las curvas coordenadas y
haga un bosquejo de sus gráficas.
32. Describa el lugar geométrico representado por cada una de las siguientes ecua-
ciones, en un sistema de coordenadas rectangulares, y haga un bosquejo de la
gráfica:
a) x = 2
b) y = π/4
c) z = −1
d) x = y
ciones, en un sistema de coordenadas cilı́ndricas, y haga un bosquejo de la
gráfica:
a) r = 2
b) θ = π/4
c) w = −1
d) w = r
e) r = cos θ
f ) r = sen θ
ciones, en un sistema de coordenadas esféricas, y haga un bosquejo de la gráfica:
31
a) ρ = 2
b) φ = π/4
c) θ = 1
d ) ρ = cos φ
e) ρ = sen φ
f ) ρ = cos θ
g) ρ = sen θ
35. En cada uno de los siguientes numerales, describa el lugar geométrico de todos
los puntos Q que satisfacen la condición dada y, para cada lugar geométrico,
dé una ecuación utilizando un sistema de coordenadas tridimensional adecua-
do.
a) Para un punto P fijo, la magnitud de PQ es 2.

b) Para un punto P fijo y una recta fija l, la distancia de P a Q es igual a
la distancia de Q a l (en el espacio).
c) Para puntos P y A diferentes y fijos, el ángulo formado por los vectores
PA y PQ es π6 .
d ) Para puntos P , R y S fijos, PQ es paralelo a RS.
e) Para puntos P , R y S fijos, PQ es perpendicular a RS.
f ) Para puntos P y R fijos, P , Q y R son vértices de un triángulo equilátero.
g) Para puntos P y R fijos, P , Q y R son vértices de un triángulo isosceles.
h) Para puntos P y R fijos tales que la magnitud de PR es 5 unidades, P ,
Q y R son vértices de un triángulo de área 10 unidades cuadradas.
i ) Dada una recta fija L y un puto P de la recta. el ángulo formado por la
recta y PQ es π6 .
j ) Dado un plano P y una recta L perpendicular al plano, la distancia de Q
a la recta es igual a la distancia de Q al plano.
k ) Dado un plano P y una recta L perpendicular al plano, la distancia de Q
a la recta es el doble de la distancia de Q al plano.
l ) Dado un plano P y un punto P exterior al plano, la distancia de P a Q
es igual a la distancia de Q al plano.
m) Dados dos puntos fijos P y R y una constante K > kPRk la suma de las
distancias de P a Q y de Q a R es igual a K.
36. En un sistema de coordenadas cartesianas determine una ecuación para el

conjunto de todos los puntos P (x, y, z) del espacio tales que la suma de las
distancias a los puntos F1 (2, 0, 0) y F2 (−2, 0, 0) es 6.
37. En un sistema de coordenadas cartesianas determine una ecuación para el

conjunto de todos los puntos P (x, y, z) del espacio tales que la distancia de P
al punto F (0, 0, 3) es igual a la distancia de P al plano de ecuacin z = −3.
38. En un sistema de coordenadas cilı́ndricas determine una ecuación para el

conjunto de todos los puntos P (r, θ, w) tales que la distancia de P al pun-
to F (0, 0, 3) es igual a la distancia de P al plano w = 5.
39. En un sistema de coordenadas cilı́ndricas determine una ecuación para el con-

junto de todos los puntos P (r, θ, w) que pertenecen al plano perpendicular al
eje polar, que contiene al punto Q(5, 0, 0).

junto de todos los puntos P (r, θ, w) que pertenecen al plano paralelo tanto al
eje polar como al eje W y que contiene al punto Q(5, π2 , 0).

junto de todos los puntos P (r, θ, w) que pertenecen a la superficie de la esfera
con centro en el polo y radio 3.

junto de todos los puntos P (r, θ, w) que pertenecen a la superficie de la esfera
con centro en el punto C(4, π4 , 5) y radio 3.
43. En un sistema de coordenadas esféricas determine una ecuación para el con-

junto de todos los puntos P (ρ, θ, φ) que pertenecen a la superficie cónica que
tiene como eje de simetrı́a el semieje esférico y ángulo de la generatriz π6 .

junto de todos los puntos P (ρ, θ, φ) que pertenecen a la superficie de la esfera
con centro en el origen y radio 3.
33

junto de todos los puntos P (ρ, θ, φ) que pertenecen a la superficie de la esfera
con centro en el punto C(5, π2 , π3 ) y radio 3.

junto de todos los puntos P (ρ, θ, φ) que pertenecen al plano perpendicular al
eje esférico y que contiene al punto Q(4, 0, 0).
47. Escriba las ecuaciones dadas en coordenadas (a) cilı́ndricas, haciendo corres-
ponder el eje Z con el eje W y el eje polar con el semieje positivo X y (b)
esféricas, haciendo corresponder el eje Z con el eje esférico y el semiplano
esférico con el semiplano XZ en el cual x ≥ 0:
a) x2 + y 2 + z 2 = 25 b) x2 + y 2 = 2y c) x2 + y 2 + 9z 2 = 36
48. Describa la intersección de las superficies cuyas ecuaciones en coordenadas

esféricas son θ = π2 y φ = π4 .
49. Muestre que para a 6= 0, la ecuación ρ = 2a sin φ cos θ en coordenadas esféricas

describe la superficie de la esfera con centro (a, 0, 0) y radio |a|.
50. Sea P (a, θ, φ) un punto en coordenadas esféricas, con a > 0 y 0 < φ < π.
Entonces P está en la superficie de la esfera ρ = a. Como 0 < φ < π, el
segmento de recta desde el origen hasta P se puede extender hasta intersectar el
cilindro dado por r = a en coordenadas cilı́ndricas. Encuentre las coordenadas
cilı́ndricas del punto de intersección.
51. Describa el tetrahedro sólido de vértices (0, 0, 0), (2, 0, 0), (0, 3, 0), y (0, 0, 5)
usando coordenadas cartesianas.
52. Represente el sólido en el espacio, de puntos con coordenadas cilı́ndricas (r, θ, w)

tales que 1 ≤ r ≤ 2, π/6 ≤ θ ≤ π/4 y 1 ≤ w ≤ 2.
53. Los siguientes sólidos están encerrados por superficies en el espacio. Describalos
usando coordenadas cilı́ndricas.
a) Un cubo sólido de arista 1.

b) Un cilindro circular recto de altura 5 y radio 3.
c) La superficie de una esfera de radio 4.

d ) La superficie de una esfera de radio 4 perforada simétricamente por un
cilı́ndro de radio 2.
e) El sólido dentro de la superficie de una esfera de radio 4 y fuera de la
superficie de la esfera de radio 2 de igual centro que la primera.
54. Sean P1 y P2 puntos cuyas coordenadas esféricas son (ρ1 , θ1 , φ1 ) y (ρ2 , θ2 , φ2 ),

respectivamente. Sea v1 el vector desde el origen hasta P1 , y sea v2 el vector
desde el origen hasta P2 . Para el ángulo γ entre v1 y v2 , muestre que
cos γ = cos φ1 cos φ2 + sin φ1 sin φ2 cos( θ2 − θ1 ).
55. Muestre que la distancia d entre los puntos P1 y P2 con coordenadas cilı́ndricas
(r1 , θ1 , z1 ) y (r2 , θ2 , z2 ), respectivamente, es
q
d = r12 + r22 − 2r1 r2 cos( θ2 − θ1 ) + (z2 − z1 )2 .
56. Muestre que la distancia d entre los puntos P1 y P2 con coordenadas esféricas
(ρ1 , θ1 , φ1 ) y (ρ2 , θ2 , φ2 ), respectivamente, es
q
d = ρ21 + ρ22 − 2ρ1 ρ2 [sin φ1 sin φ2 cos( θ2 − θ1 ) + cos φ1 cos φ2 ] .
LECCIÓN 3
Funciones vectoriales
En esta lección estudiaremos funciones definidas en conjuntos especiales de Rn de-

nominados celdas.
Una 1-celda I es un intervalo cerrado [a, b] de números reales1 , es decir,
I = [a, b] = {x ∈ R : a ≤ x ≤ b}. (3.1)
Una 2-celda R es el producto cartesiano de dos 1-celdas [a, b] y [c, d], es decir
R = [a, b] × [c, d] = {(x, y) : a ≤ x ≤ b, c ≤ y ≤ d} . (3.2)
Una 3-celda P es el producto cartesiano de tres 1-celdas [a, b], [c, d] y [e, h], es decir,
P = {(x, y, z) : a ≤ x ≤ b, c ≤ y ≤ d, e ≤ z ≤ h} = [a, b] × [c, d] × [e, h] . (3.3)

1
Las 1-celdas serán los dominios de funciones vectoriales de una sola variable. Aunque, en los
libros de cálculo se consideran las funciones vectoriales definidas en dominios más generales, en
esta lección consideraremos los dominios de funciones vectoriales de una variable como 1-celdas.
En algunos casos permitiremos también que las 1-celdas sean de las siguientes formas: (a, b), (a, b],
[a, b), [a, b], [a, +∞), (a, +∞), (−∞, b], (−∞, b) o (−∞, +∞).
35
36 LECCIÓN 3. FUNCIONES VECTORIALES
Una función vectorial 2 es una función definida en alguna celda y que toma valores
en el espacio de vectores bidimensionales o tridimensionales 3 . La función vectorial
es de una variable si su dominio es una 1-celda, de dos variables si su dominio es
una 2-celda y de tres variables si su dominio es una 3-celda.
Una función vectorial tridimesional de una variable es una función definida en una
1-celda I que toma valores en R3 :
r : I −→ R3 , dada por r(t) = hx(t), y(t), z(t)i (3.4)

Por ser una función que toma valores en R3 , tiene tres funciones componentes, a
saber, x, y, z. Cada una de éstas es una función escalar definida en un intervalo
I y que toma valores en R. Debido a esto, los conceptos estudiados en los cursos
de cálculo diferencial e integral, para funciones reales de valor real se extienden
naturalmente a estas funciones vectoriales. Por ejemplo, si los lı́mites lı́mt→a x(t),
lı́mt→a y(t) y lı́mt→a z(t) existen, el lı́mite lı́mt→a r (t) existe y
D E
lı́m r (t) = lı́m x(t), lı́m y(t), lı́m z(t) . (3.5)
t→a t→a t→a t→a
Como consecuencia de lo anterior la función vectorial r =r (t) es diferenciable en I si

dx dy
cada una de sus componentes escalares lo es. Además, si las derivadas (a), (a)
dt dt
dz
y (a) existen, la derivada de la función vectorial r =r (t) en t = a es
dt

dr dx dy dz
(a) = (a), (a), (a) (3.6)
dt dt dt dt
y, si las funciones componentes son continuas en a, la función r es continua en a.
Por otro lado, si cada una de las componentes escalares de r son integrables sobre
[a, b], entonces r es integrable sobre [a, b] y
Z b Z b Z b Z b
r(t)dt = x(t)dt, y(t)dt, z(t)dt . (3.7)
a a a a
2
En realidad, una función vectorial es una función definida en algún conjunto y que toma valores
en algún espacio vectorial. Sin embargo, en estas lecciones trabajaremos con funciones vectoriales
definidas en celdas y con valores en R2 o en R3 .
3
En casos especiales podrı́amos considerar también espacios de vectores unidimensionales
37
A partir de dos funciones vectoriales r(t) = hx(t), y(t), z(t)i y s(t) = hex(t), ye(t), ze(t)i
definidas en un intervalo I se definen naturalmente en dicho intervalo las funciones
vectoriales
(r + s)(t) = hx(t) + x
e(t), y(t) + ye(t), z(t) + ze(t)i , (3.8)
(r − s)(t) = hx(t) − x
e(t), y(t) − ye(t), z(t) − ze(t)i , (3.9)
(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 , (3.10)
la función escalar
(r · s)(t) = x(t)e
x(t) + y(t)e
y (t) + z(t)e
z (t), (3.11)
si f es una función escalar definida en el mismo intervalo, la función vectorial
f r(t) = hf (t)x(t), f (t)y(t), f (t)z(t)i , (3.12)
si r es derivable en I, la función vectorial

r′ (t) = x′ (t), y ′ (t), z ′ (t) , (3.13)
y si existen las antiderivadas de las componentes, la familia de funciones vectoriales

Z Z Z Z
r(t)dt = x(t)dt, y(t)dt, z(t)dt . (3.14)
Supongamos ahora que r es diferenciable en su dominio. Para cada valor de t, el vec-

tor r(t) = hx(t), y(t), z(t)i se puede interpretar como el vector posición de un punto
P de coordenadas (x(t), y(t), z(t)). Ası́, cuando t recorre la 1-celda I, el punto P
describe una curva en el espacio, se dice en este caso que las ecuaciones paramétricas
de la curva son
x = x(t), y = y(t), z = z(t) (3.15)
o que r es una parametrización de la curva.

En otras palabras, el conjunto

P (x, y, z) ∈ R3 : OP = hx(t), y(t), z(t)i , t∈I , (3.16)
que es el conjunto de puntos que la función vectorial r señala a partir del origen
cuando t recorre la 1-celda I, representa una curva en el espacio.
Por ejemplo, consideremos el segmento de recta AB, que es el conjunto de puntos

P (x, y, z) de la recta que pasa por A(a0 , b0 , c0 ) y B(a1 , b1 , c1 ) y que están entre A y
B. Esta recta tiene vector director ha1 − a0 , b1 − b0 , c1 − c0 i y por consiguiente las
coordenadas de P deben satisfacer las ecuaciones paramétricas de la recta
x = a0 + t(a1 − a0 ), y = b0 + t(b1 − b0 ), z = c0 + t(c1 − c0 ),
o
x = (1 − t)a0 + ta1 , y = (1 − t)b0 + tb1 , z = (1 − t)c0 + tc1 ,
imponiendo la condición de que t esté en [0, 1] para garantizar que P esté entre A y B.
Podemos ası́ decir que la función vectorial
r(t) = h(1 − t)a0 + ta1 , (1 − t)b0 + tb1 , (1 − t)c0 + tc1 i

= (1 − t) ha0 , b0 , c0 i + t ha1 , b1 , c1 i (3.17)
= (1 − t)OA + tOB
t ∈ [0, 1], parametriza el segmento de recta AB.
Parametricemos ahora la curva de intersección de dos superficies. Considérese la

curva Γ de intersección de la superficie cilı́ndrica recta circular, de ecuación x2 +y 2 =
4, con el plano de ecuación z = x + 3. En otras palabras Γ = {(x, y, z) : x2 + y 2 =
4 ∧ z = x + 3}. Si para algún t ∈ [0, 2π], (x, y, z) satisface las ecuaciones x = 2 cos t,
y = 2 sen t, el punto P (x, y, z) pertenecerá a la superficie x2 + y 2 = 4, ya que, en
este caso, x2 + y 2 = 4 cos2 t + 4 sen2 t = 4. Si además, si z = (2 cos t) + 3, el punto
P (x, y, z) pertenecerá al plano z = x + 3. Ası́,
x = 2 cos t, y = 2 sen t, z = 2 cos t + 3,
con t ∈ [0, 2π], son ecuaciones paramétricas de Γ y
r(t) = h2 cos t, 2 sen t, 2 cos t + 3i
con t ∈ [0, 2π], parametriza a Γ.

Es de anotar que una curva no tiene una única parametrización; existen muchas
funciones vectoriales que recorren una curva determinada. Por ejemplo, las funcio-
nes vectoriales r1 (t) = cos ti + sen tj, con t ∈ [0, 2π], y r2 (t) = cos 2ti + sen 2tj, con
39
t ∈ [0, π], parametrizan la circunferencia con centro en (0, 0) y radio 1.
Una función vectorial tridimesional 4 de dos variables es una función definida en una
2-celda R que toma valores en R3 :
r : R −→ R3 , dada por r(u, v) = hx(u, v), y(u, v), z(u, v)i (3.18)
Si R = [a, b] × [c, d], para cada valor u0 ∈ [a, b] y cada valor v0 ∈ [c, d] se obtienen
funciones vectoriales de una sola variable
r1 : [c, d] −→ R3 , (3.19)
y
r2 : [a, b] −→ R3 , (3.20)
definidas por r1 (v) = r(u0 , v) y r2 (u) = r(u, v0 ), respectivamente. Si para todo

(u0 , v0 ) ∈ R, las funciones r1 y r2 son diferenciables, el conjunto

P (x, y, z) ∈ R3 : OP = hx(u, v), y(u, v), z(u, v)i , (u, v) ∈ R , (3.21)
que es el conjunto de puntos que la función vectorial r señala desde el origen cuan-
do (u, v) recorre R, representa una superficie en el espacio y se dice que r es una
parametrización de la superficie.
Para ilustrar lo dicho, parametricemos la superficie σ del plano limitada por las
gráficas de las funciones f (x) = x2 y g(x) = −x2 + 2x + 4 (haga una gráfica de σ).
En este caso como la superficie es un subconjunto de R2 , la función vectorial que
buscamos toma valores en R2 , y por lo tanto tendrá sólamente dos componentes.
Procedemos de la siguiente manera:
1. Los puntos de intersección de las gráficas son (−1, 1) y (2, 4). Todos los puntos
de σ tienen abscisa x entre -1 y 2, por lo que podemos tomar como primer
parámetro u = x con −1 ≤ u ≤ 2.
2. Para cada u ∈ [−1, 2] consideramos los puntos A(u, f (u)) en la gráfica de f y

B(u, g(u)) en la gráfica de g.
4
Pueden considerarse también funciones vectoriales bidimensionales de dos variables, es decir
definidas en 2-celdas y que toman valores en R2
b B
b b A
b b
−1 2
Figura 1.3
3. El segmento AB, que está contenido en la región, se puede parametrizar me-

diante la función vectorial (ver (3.17))
r(u, v) = (1 − v) hu, f (u)i + v hu, g(u)i

= (1 − v) u, u2 + v u, −u2 + 2u + 4

= u, u2 − 2vu2 + 2vu + 4v
con 0 ≤ v ≤ 1, lo que nos da una parametrización de σ limitada por las curvas

dadas.
4. Obsérvese que el dominio de r es [−1, 2] × [0, 1]
Parametricemos ahora la superficie σ que es la parte del paraboloide de ecuación

z = x2 +y 2 que se encuentra dentro del cilindro circular recto de ecuación x2 +y 2 = 4,

en otras palabras σ = (x, y, z) : x2 + y 2 ≤ 4 ∧ z = x2 + y 2 . Procedemos de la
siguiente manera:
Para cada u ∈ R+ , x = u cos v, y = u sen v con v variando entre 0 y 2π parametriza
una circunferencia con centro en el origen y radio u. Por lo tanto, si además de variar
v, hacemos variar a u entre 0 y 2, x = u cos v, y = u sen v parametriza el disco en el
plano XY con centro en el origen y radio 2. Ası́, para cada punto (x, y, z) de σ existe
un valor de u entre 0 y 2 y uno de v entre 0 y 2π tales que x = u cos v, y = u sen v y,
como z = x2 + y 2 , entonces z = u2 . En consecuencia, una parametrización de σ es
r : [0, 2] × [0, 2π] −→ R3

41
definida por

r(u, v) = u cos v, u sen v, u2 .
Al igual que en la parametrización de curvas, la parametrización de una superficie
no es única, ya que existen muchas funciónes vectorial de dos variables que recorren
una superficie determinada.
Una función vectorial tridimesional de tres variables es una función definida en una
3-celda P que toma valores en R3 :
r : P −→ R3 , dada por r(u, v, w) = hx(u, v, w), y(u, v, w), z(u, v, w)i (3.22)
Si P = [a, b] × [c, d] × [e, h], para cada valor u0 ∈ [a, b], cada valor v0 ∈ [c, d] y cada
valor w0 ∈ [e, h] se obtienen funciones vectoriales de una sola variable
r1 : [e, h] −→ R3 , (3.23)
r2 : [c, d] −→ R3 (3.24)
y
r3 : [a, b] −→ R3 , (3.25)
definidas por r1 (w) = r(u0 , v0 , w), r2 (v) = r(u0 , v, w0 ) y r3 (u) = r(u, v0 , w0 ), respec-
tivamente. Si para todo (u0 , v0 , w0 ) ∈ P las funciones r1 , r2 y r3 son diferenciables,
el conjunto

P (x, y, z) ∈ R3 : OP = hx(u, v, w), y(u, v, w), z(u, v, w)i , (u, v, w) ∈ P , (3.26)
que es el conjunto de puntos que la función vectorial r señala desde el origen, repre-
senta una sólido en el espacio.
En muchos casos es posible encontrar una función vectorial de tres variables que
recorra un sólido Λ determinado. Dicha función vectorial es una parametrización de
Λ. Por ejemplo, la función vectorial r : [0, 2] × [0, 2π] × [0, π] −→ R3 definida por
r(u, v, w) = hu cos v sen w, u sen v sen w, u cos wi , (3.27)

es una parametrización de la esfera sólida Λ = (x, y, z) : x2 + y 2 + z 2 ≤ 4 . Es claro
que
Λ = {(u cos v sen w, u sen v sen w, u cos w) : (u, v, w) ∈ [0, 2] × [0, 2π] × [0, π]} .
(3.28)
Las funciones vectoriales y su cálculo se usan frecuentemente en dos contextos que

estudiaremos a continuación.
Geometrı́a de curvas y superficies
Supongamos que la curva Γ está parametrizada por la función vectorial diferenciable
r : I −→ R3 , r(t) = hx(t), y(t), z(t)i . (3.29)
Diremos que Γ es regular si r′ (t) 6= 0, para todo t ∈ I y, en este caso, que el vector
r′ (t)
T(t) = , (3.30)
kr′ (t)k
es el vector unitario tangente 5 a Γ en el punto P (x(t), y(t), z(t)). Como T(t) y

r′ (t) son paralelos, entonces r′ (t) también es un vector tangente a Γ 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 fácilmente al calcular la derivada
de la función constante (T · T)(t) = 1. En efecto,
d
0= [T(t) · T(t)] = 2T(t) · T′ (t) (3.31)
dt
Fijemos un punto P sobre una curva Γ junto con una orientación determinada y
tomemos otro punto Q sobre ella. El procedimiento que consiste en escribir las
coordenadas de Q en términos de la longitud del arco de curva desde P hasta Q, se
conoce como parametrización por longitud de arco de la curva Γ. La longitud de este
arco es precisamente la distancia que hay que recorrer sobre la curva partiendo de P
y llegando a Q. Encontrar directamente la parametrización por longitud de arco de
una curva es muy difı́cil (si no imposible) por lo que se recurre a encontrar primero
una parametrización conveniente y utilizarla para calcular la longitud del arco que
parte de P y llega a Q. En cursos de cálculo anteriores se dedujo una fórmula para
5
En realidad, hay dos vectores tangentes unitarios. El otro se obtiene recorriendo Γ en sentido
contrario mediante una parametrización conveniente. Sin embargo, podrı́amos preguntarnos si con
otra parametrización no serı́a posible obtener otros vectores tangentes unitarios. ¿Podrı́a ser?
43
calcular la longitud de estos arcos para curvas en el plano. De manera similar se

puede deducir una fórmula para calcular longitudes de curvas en el espacio.
Recordemos que la longitud de arco entre P (x(a), y(a)) y Q(x(t0 ), y(t0 )) de una
curva Γ en el plano, parametrizada por
x = x(t), y = y(t), t ∈ I, (3.32)
es Z t0 p
x′ (τ )2 + y ′ (τ )2 dτ, (3.33)
a
y para el caso de una curva en el espacio, haciendo un razonamiento análogo,

al que se hace para obtener esta fórmula, se tiene que la longitud de arco entre
P (x(a), y(a), z(a)) y Q(x(t0 ), y(t0 ), z(t0 )) de una curva Γ en el espacio parametriza-
da por
x = x(t), y = y(t), z = z(t) t ∈ I. (3.34)
es Z t0 p
x′ (τ )2 + y ′ (τ )2 + z ′ (τ )2 dτ. (3.35)
a
Obsérvese que, siendo la función vectorial r (t) = hx(t), y(t), z(t)i una parametriza-
ción de Γ, la integral anterior se puede escribir en la forma
Z t0 ′
r (τ ) dτ. (3.36)
a
Evidentemente el valor de esta integral depende de la posición del punto Q sobre la

curva y por lo tanto es en realidad una función s = s(t):
Z t
s(t) = r′ (τ ) dτ. (3.37)
a
Esta función siempre es invertible(¿por qué?), por lo tanto podemos escribir las
coordenadas del punto Q en términos 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
parametrización por longitud de arco de Γ:
x=x
e(s) = x(t(s)), y = ye(s) = y(t(s)) z = ze(s) = z(t(s)) t ∈ I. (3.38)
La función vectorial que parametriza a Γ es ahora
r(s) = he
e x(s), ye(s), ze(s)i .
Tenemos entonces que r (t) = e

r(s(t)) y por lo tanto,
d de
r ds de
r
r′ (t) = r(s(t))] = (s(t)) (t) = (s(t)) r′ (t)
[e
dt ds dt ds
ya que la derivada de s es precisamente kr′ (t)k, por el Teorema Fundamental del
Cálculo. De aquı́ que
de
r r′ (t)
(s(t)) = ′ = T(t), (3.39)
ds kr (t)k
lo que muestra que el vector tangente que se obtiene mediante parametrizaciones
por longitud de arco es unitario.
La dirección de la curva Γ en cada punto Q(x(t), y(t), z(t)) está representada por
el vector T (t). La curvatura de Γ es una medida de la rapidez con que cambia
la dirección de la curva con respecto al desplazamiento sobre ella. Esta medida
está dada por
dT(t(s))
κ= ds
(3.40)
dT(t(s))
donde s es el parámetro de longitud de arco. Si desarrollamos la derivada
ds
usando la regla de la cadena obtenemos.
dT
dT(t(s)) dT dt dT 1 (t(s))
= (t(s)) (s) = (t(s)) = dt′ .
ds dt ds dt ds kr (t(s))k
(t(s))
dt
De aquı́ se obtiene la relación

dT ′

dt (t(s)) = κ r (t(s)) (3.41)
que necesitaremos más adelante.
La fórmula que define a κ es dispendiosa de calcular, porque involucra el cálculo

del vector tangente unitario y la reparametrización por longitud de arco. Hay otras
45
maneras de calcular la curvatura sin necesidad de hacer ésto. Como r′ = kr′ k T =

ds
T, entonces
dt
ds d2 s
r′′ = T′ + 2 T
dt dt
y por lo tanto,
ds 2
r′ × r′′ = ( )2 T × T′ = r′ T × T′ .
dt
Calculando la norma de este producto, teniendo en cuenta que T y T ′ son perpen-
diculares y que kT′ k = κ kr′ k se obtiene
′ ′ 2
r × r′′ = 2 3
r T × T′ = r′ T′ = r′ κ.
Despejando κ obtenemos una fórmula para calcular la curvatura de Γ a partir de r ,

r ′ y r ′′ :
kr′ × r′′ k
κ= . (3.42)
kr′ k3
T′ (t)
Consideremos ahora el vector N(t) = , que es un vector unitario perpen-
kT′ (t)k
dicular a T (t). Es decir, N es un vector unitario perpendicular a Γ en el punto
P (x(t), y(t), z(t)) que se conoce como vector normal unitario. Este vector indica la
dirección 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 Γ en el punto
P (x(t), y(t), z(t)). Un vector normal a este plano es
B = T × N, (3.43)
que también es unitario y que se conoce como vector binormal de Γ 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 Γ en P y tiene radio r = 1/κ se
denomina circunferencia osculadora y es la que mejor aproxima a Γ en el punto P ,
esto último quiere decir que la curva Γ y la circunferencia tienen la misma curvatura
en el punto P .
Sea σ una superficie en R3 parametrizada por una función vectorial de dos va-
riables r : [a, b] × [c, d] −→ R3 dada por r(u, v) = hx(u, v), y(u, v), z(u, v)i, y
P (x(u0 , v0 ), y(u0 , v0 ), z(u0 , v0 )) un punto de la superficie. Obsérvese que las fun-

ciones vectoriales
r1 : [c, d] −→ R3
y
r2 : [a, b] −→ R3
definidas por r1 (v) = r(u0 , v) y r2 (u) = r(u, v0 ), respectivamente, parametrizan dos
curvas Γu0 y Γv0 contenidas en la superficie σ y que se cortan en P . Al calcular las
derivadas de r1 y r2 en v0 y u0 , respectivamente, obtenemos dos vectores tangentes a
Γu0 y Γv0 , respectivamente, en el punto P . Explı́citamente, estos vectores tangentes
estan dados por
r′1 (v0 ) = rv (u0 , v0 ) = hxv (u0 , v0 ), yv (u0 , v0 ), zv (u0 , v0 )i (3.44)

y
r′2 (u0 ) = ru (u0 , v0 ) = hxu (u0 , v0 ), yu (u0 , v0 ), zu (u0 , v0 )i (3.45)

donde el subı́ndice v indica que se ha derivado con respecto a la variable v, y el
subı́ndice u indica que se ha derivado con respecto a la variable u. Como rv (u0 , v0 )
y ru (u0 , v0 ) son tangentes a Γu0 y Γv0 , respectivamente, en el punto P , también lo
son a la superficie σ en el punto P . Por consiguiente, el vector rv (u0 , v0 ) × ru (u0 , v0 )
es perpendicular a σ en el punto P , siempre y cuando rv (u0 , v0 ) × ru (u0 , v0 ) 6= 0.
Diremos que σ es regular en P si rv (u0 , v0 ) × ru (u0 , v0 ) 6= 0, y que σ es regular, si lo
es en todos sus puntos.
Si σ es una superficie regular parametrizada por una función vectorial de dos va-
riables r : [a, b] × [c, d] −→ R3 y P (x(u0 , v0 ), y(u0 , v0 ), z(u0 , v0 )) es un punto de σ,
entonces σ tiene plano tangente en P y cuya ecuación es
(rv (u0 , v0 ) × ru (u0 , v0 )) · hx − x(u0 , v0 ), y − y(u0 , v0 ), z − z(u0 , v0 )i = 0. (3.46)
Al igual que en el caso de curvas en el espacio, se tiene el concepto de curvatura

de una superficie. En realidad hay dos conceptos de curvatura que se conocen como
curvatura gaussiana K(P ) de σ en P y la curvatura media H(P ) de σ en P . En
el contexto de este curso no profundizaremos más en la geometrı́a de superficies,
47
pero si daremos las expresiones explı́citas de estas dos curvaturas en términos de la

parametrización de la superficie en los ejercicios.
Movimiento en el espacio
En el contexto de la cinemática la posición de un objeto que se mueve en el espacio se

puede modelar mediante una función vectorial r = r(t), t ∈ I, tomando el parámetro
t como el transcurrir del tiempo.
Ası́, la derivada r′ (t) = v(t) representa la velocidad del objeto que es la razón de
cambio instantáneo de la posición en el instante de tiempo t y contiene la informa-
ción de la dirección de su movimiento y la rapidez kv(t)k con que lo hace en cada
instante de tiempo t.
La derivada v′ (t) = a(t), representa la aceleración del objeto que es la razón de

cambio instantáneo de la velocidad en cada instante de tiempo t y contiene la infor-
mación de cómo cambian la dirección de movimiento y la rapidez con el transcurrir
del tiempo en cada punto P (x(t), y(t), z(t)) en el que se encuentra.
La razón de cambio instantáneo de la dirección del movimiento se denomina acele-

ración normal aN y la razón de cambio instantáneo de la rapidez del movimiento se
denomina aceleración tangencial aT .
Para encontrar la descomposición de la aceleración en términos de la aceleración

normal y la tangencial, se escribe la aceleración a como combinación lineal de T y
N. Es decir, debemos encontrar aT y aN tales que
a = aT T + aN N.
Resulta que aT y aN van a depender de la rapidez con que se mueve el objeto y de

la curvatura de la curva. Esta descomposición se deja como ejercicio.
EJERCICIOS
dr r(a + h) − r(a)
1. Sabiendo que (a) = lı́m , utilice las propiedades de la suma
dt h→0 h
de vectores, la multiplicación por escalar y (3.5) para demostrar (3.6).
2. Muestre que si las funciones vectoriales r y s son derivables en a entonces

d(r × s) dr ds
a) (a) = (a) × s(a) + r(a) × (a)
dt dt dt
d(r · s) dr ds
b) (a) = (a) · ×s(a) + r(a) · (a)
dt dt dt
d(f r) df dr
c) (a) = (a)r(a) + f (a) (a)
dt dt dt
3. El Teorema del Valor Medio no se satisface en general para funciones vecto-
riales. 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
4. Para un vector c 6= 0, la función r(t) = tc parametriza la recta que contiene
al origen y tiene vector director c.
a) ¿ Qué curva es parametrizada por p(t) = t3 c? Explique.

b) ¿ Qué curva es parametrizada por q(t) = et c? Explique.
c) Compare r ′ (0) y p ′ (0). ¿ Cómo explica la diferencia en las dos derivadas?

d dr d2 r
5. Muestre que r× = r × 2.
dt dt dt

d dr dp dq
6. Muestre que (r · (p × q)) = · (p × q) + r · × q +r· p × .
dt dt dt dt
7. Determine la derivada de cada una de las siguientes funciones en t = 3/2.

a) r(t) = t2 , ln t, t .
D√ E
b) r(t) = t2 − 4, t, 1/t .

8. Determine la integral de la función r(t) = t2 , 2 sen t, ln t en el intervalo [1, π]
9. Determine una función vectorial r tal que r′ (t) = hsen t, cos t, ti y r(π/3) =
h2, 5, −1i.
49
D √ E
10. Determine una función vectorial r tal que r′′ (t) = t3 + 2, t t2 + 9, t , r′ (3) =
h2, 3, 0i y r(0) = h1, −1, 4i.
11. Determine tres parametrizaciones diferentes del segmento RS si las coordena-

das de R y de S son (1, 1, 0) (−5, 5, 10), respectivamente. Sugerencia: use tres
parametrizaciones diferentes de la recta que pasa por los puntos R y S.
12. Determine tres parametrizaciones diferentes del segmento RS si las coordena-

das de R y de S son (1, −1, 7) (3, −1, 4), respectivamente.
13. Describa la curva de intersección de la superficie de la esfera de centro en el

origen y radio 3 con el plano z = 2 y determine una parametrización de la
misma.
14. Describa la curva de intersección de las superficies cilı́ndricas de ecuaciones

x2 + y 2 = 4 y y 2 + z 2 = 4 y determine una parametrización de la misma.
15. Determine una parametrización para la recta de intersección de los planos

x + y = 5 y z − y = 2.
16. Determine una parametrización para la recta de intersección de los planos

x + y − 3z = 5 y x − 3y + z = 2.
17. Realice la siguiente construcción: Trace dos circunferencias con centro en el

origen O de un plano cartesiano, una de radio a y otra de radio b con a > b,
ubique un punto B sobre la circunferencia de radio b y trace la semirrecta
OB. Nombre A el punto de intersección de la semirrecta con la circunferencia
de radio a y θ el ángulo formado por el semieje positivo X y la semirrecta
OB medido en sentido positivo y en radianes. Trace la recta paralela al eje X
que contiene al punto B y la recta paralela al eje Y que contiene al punto A,
nombre P al punto de intersección de las dos rectas. Dibuje el lugar geométrico
generado por el punto P al variar el ángulo θ. Describa este lugar geométrico y
parametrı́celo usando como parámetro el ángulo θ. Muestre que efectivamente
es el lugar geométrico que describió.
18. Realice la siguiente construcción: Trace una circunferencia con centro en el

punto (0, a) y radio a, los puntos de intersección entre la circunferencia y el
eje Y son el origen O y el punto A. Trace la recta L1 , paralela al eje X por el
punto A, trace una semirrecta con origen el punto O y que forme un ángulo θ
con el semieje positivo X, Sea C la intersección entre la semirrecta y la recta
L1 . Trace la recta L2 paralela al eje Y por el punto C. Sea B la intersección
de la semirrecta con la circunferencia, trace la recta L3 paralela al eje X por
el punto B. Sea P el punto de intersección de L2 y L3 . La curva generada por
el punto P al variar θ entre 0 y π se conoce como Bruja de Agnesi.
a) Parametrice la curva.
8a3
b) Muestre que la ecuación rectangular de la curva es y = .
x2 + 4a2
19. Dadas una recta y una circunferencia de radio a, coplanares, tales que el centro
de la circunferencia se encuentra a una distancia b de la recta.
a) Describa la superficie6 que genera la circunferencia al ser rotada tomando

como eje la recta si:
1) Si b > a
2) Si b = a
3) Si b < a
b) Si se introduce un sistema de coordenadas rectangulares de tal manera
que la recta coincida con el eje Z, parametrice la superficie generada para
cada uno de los casos del numeral anterior.
c) Parametrice el sólido interior a la superficie7 generada para el caso en que
b > a.
20. Parametrice la región del plano limitada por las curvas dadas, de manera que
el dominio de la parametrización 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
6
Esta superficie se denomina ((torus))
7
este sólido se denomina ((toroide))
51
√
f ) y = x2 , y = x
g) x + y = 3, x + y = 7, x − y = 4, x − y = 1.
21. Sea σ la superficie de la esfera con centro en el origen y radio 2 y σ ′ la superficie

p
cónica de ecuación z = x2 + y 2 , parametrice cada una de las siguientes
regiones del espacio de manera que el dominio de los parámetros sea una 1-
celda (intervalo), una 2-celda o una 3-celda, dependiendo de la naturaleza de
la región8 .
a) La superficie σ.
b) El sólido limitado por σ, es decir, la esfera.
c) La parte de la superficie σ que está dentro de la superficie cónica σ ′ .
d ) La parte de la superficie cónica σ ′ que se encuentra dentro de la superficie
σ.
e) La curva de intersección de las dos superficies, σ y σ ′ .
f ) El sólido que está dentro de σ y fuera de σ ′ .
22. Siendo σ la parte del plano de ecuación 2x + 4y + 3z = 6 que se encuentra en

el primer octante, parametrice cada una de las siguientes regiones del espacio
de manera que el dominio de la parametrización sea una 1-celda, una 2-celda
o una 3-celda, dependiendo de la naturaleza de la región.
a) La superficie σ.
b) El sólido limitado por la superficie σ y los planos coordenados.
c) El borde de σ.
23. Usando las parametrizaciones de las superficies del ejercicio 21 encuentre la

ecuación del plano tangente en los siguientes puntos:
√ √ √
a) El punto (2 3/3, 2 3/3, 2 3/3) de la superficie σ.
√ √
b) El punto ( 2/2, 2/2, 1) de la parte de la superficie σ ′ .
8
σ ′ no tiene dentro ni fuera porque no es una superficie cerrada. En este ejercicio, fuera se
p
refiere a los puntos de coordenadas (x, y, z) que satisfacen la condición z < x2 + y 2 y dentro hace
p
referencia a los puntos de coordenadas (x, y, z) que satisfacen la condición z > x2 + y 2
x2 y 2
24. Sea Γ la elipse de ecuación cartesiana + = 1.
4 9
a) Determine una parametrización de Γ.
b) Determine las ecuaciones de las circunferencias osculadoras de Γ en los
p
puntos (2, 0) y (1, 27/4)
c) Dibuje Γ y las circunferencias osculadoras del literal anterior.
25. Determine las ecuaciones paramétricas de la recta tangente a cada una de las
curvas Γ parametrizadas por r en el punto en el que t = 0.

D E
a) r(t) = t + 1, t2 + 1, t3 + 1 2
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
26. Calcule la velocidad v(t) y la aceleración a(t) de un objeto con el vector posi-
ción r(t) dado:
a) r(t) = ht, t − sin t, 1 − cos ti b) r(t) = h3 cos t, 2 sin t, 1i
27. Si v representa la velocidad de un objeto en movimiento, podemos escribir

v ds ds
v(t) = kvk = (t)T(t), para explicitar la rapidez y la dirección de
kvk dt dt
ds
movimiento T(t). Use la expresión v(t) = (t)T(t) para calcular la derivada
dt 2
d2 s ds
de v y muestre que a(t) = 2 (t)T(t) + (t) κN(t) donde κ = κ(t) es
dt dt
la curvatura de la curva que sigue el objeto en su camino y N(t) es el vector
normal de la curva en el instante de tiempo t. Esto muestra que aceleración
normal está dada por
aN = (s′ )2 κN,
y la aceleración tangencial está dada por
aT = s′′ T.

cos t sin t −at
28. Sea r(t) = √ ,√ ,√ , con a 6= 0.
1 + a 2 t2 1 + a 2 t2 1 + a 2 t2
a) Muestre que kr(t)k = 1 para todo t.
53
b) Muestre directamente que r ′ (t) · r(t) = 0 para todo t.
29. Si r′ (t) = 0 para todo t en algún intervalo (a, b), muestre que r(t) es constante
en (a, b).
30. Considere una partı́cula de masa constante m con vector posición r(t), ve-
locidad v(t), aceleración a(t) y momento p(t) en el instante t. El momento
angular L(t) de la partı́cula 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 partı́cula en
el instante t, se define el torque N(t) con respecto al origen, actuando sobre la
partı́cula por N(t) = r(t) × F(t). Muestre que L ′ (t) = N(t).
31. Sea r(t) el vector posición de una partı́cula que se mueve en R3 . Muestre que
d
(r × (v × r)) = krk2 a + (r · v)v − (kvk2 + r · a)r.
dt
32. Sea r(t) el vector posición en R3 de una partı́cula 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.
33. 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]

√
b) r(t) = (t2 + 1) cos t, (t2 + 1) sin t, 2 2t en [0, 1]

c) r(t) = 2 cos 3t, 2 sin 3t, 2t3/2 en [0, 1]
34. Determine la trayectoria r(t) de un punto P del borde de un disco de radio a

que rueda a lo largo de una recta L con una velocidad constante v.
35. Suponga que un objeto que sigue la trayectoria r(t) = ti + sin tj + cos tk se sale
por la tangente en t = 7π
3 . Determinar la posición de la partı́cula en t = 12.
36. Suponga que sobre el objeto, al que se hace referencia en el ejercicio anterior,
se ejerce una fuerza F = −10mi (donde m es la masa del objeto) a partir del
momento en que se sale por la tangente, determine la posición del objeto en
cada instante de tiempo, la altura máxima del objeto sobre el plano x = 0, el
punto de impacto con el plano x = 0.
37. Determine la posición, en cada instante de tiempo, de un objeto de masa m

que se mueve con rapidez constante v, siguiendo una trayectoria circular de
radio a.
38. Dado que la fuerza centrı́peta es el producto de la masa por la aceleración de un

objeto que semueve con rapidez constante siguiendo una trayectoria circular,
determine la fuerza centrı́peta de un objeto de 2 kg de masa que sigue una
trayectoria circular de radio 9 m con una frecuendia de 4 revoluciones por
segundo.
39. Si T, N y B son respectivamente los vectores tangente unitario, normal y

binormal, muestre que
dB
a) ·B=0
dt
dB
b) ·T=0
dt
dB
c) es múltiplo escalar de N
dt
40. Un objeto se mueve siguiendo la trayectoria r(t) = 2 sen t2 i + 2 cos t2 j, t ≥ 0.
a) Determine el vector velocidad y la rapidez de la partı́cula como funciones

del parámetro t.
b) ¿En qué punto de la circunferencia deberı́a liberarse el objeto con el fin
de alcanzar un blanco situado en el punto (4, 0, 0).
c) ¿En qué instante de tiempo deberá ser liberada la partı́cula?
d ) ¿Cuál es la velocidad y la rapidez en el instante de la liberación?
e) ¿En qué instante se alcanza el blanco?
41. 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).
42. 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
55
Ayuda: Para calcular dt d

(kr(t)k), calcule de dos formas la derivada de kr(t)k2
y luego use el ejercicio 33 de la Lección 1.
43. Sea r(t) una curva suave tal que r ′ (t) 6= 0 para todo t y T es el vector tangente
unitario. Muestre que
r ′ (t) × (r′′ (t) × r ′ (t))
T ′ (t) = .
kr ′ (t)k3
Ayuda: use el ejercicio anterior.
44. 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
r′ (t) × (r′′ (t) × r′ (t))
N(t) = .
kr′ (t)k kr ′′ (t) × r ′ (t)k
Ayuda: use el ejercicio anterior.
45. Si B es el vector binormal unitario, muestre que

r ′ (t) × r ′′ (t)
B(t) = .
kr ′ (t) × r ′′ (t)k
Ayuda: Use el ejercicio anterior y el ejercicio 33 de la Lección 1. Los vectores
T(t), N(t) y B(t) forman un sistema de vectores mutuamente perpendiculares

de vectores unitarios (vectores ortonormales) en cada punto de la curva r(t).
46. Encuentre T, N, B y κ en cada punto de la hélice r(t) = hcos t, sin t, ti.
47. Muestre que la longitud de arco L de una curva cuyas coordenadas esféricas
son ρ = ρ(t), θ = θ(t) y φ = φ(t) para t en un intervalo [a, b] es
Z bq
L= ρ ′ (t)2 + (ρ(t)2 sin2 φ(t)) θ ′ (t)2 + ρ(t)2 φ ′ (t)2 dt.
a
48. Un objeto de masa m se desliza sin rozamiento sobre una rampa circular
(cuarto de circunferencia de radio R). Introduzca un sistema de coordenadas
cartesianas conveniente y encuentre una funci’on vectorial que describa la po-
sición del objeto en función del tiempo transcurrido, si éste partió del reposo
desde la parte más alta de la rampa.
49. Repita el ejercicio anterior pero suponiendo que hay una fuerza de rozamiento
cuya magnitud es proporcional a la fuerza normal.
50. La curva Γ de Bézier pasa por cuatro puntos no colineales B0 , B1 , B2 y B3 en

R3 y está parametrizada mediante la función vectorial b30 : [0, 1] −→ R3 , a par-
tir de los vectores posición b0 , b1 , b2 , b3 de B0 , B1 , B2 y B3 , respectivamente,
que está definida mediante el siguiente procedimiento:
Se definen las funciones vectoriales:

1) b10 (t) = (1 − t)b0 + tb1
2) b11 (t) = (1 − t)b1 + tb2
3) b12 (t) = (1 − t)b2 + tb3
A partir de las anteriores se definen las funciones vectoriales
1) b20 (t) = (1 − t)b10 (t) + tb11 (t)
2) b21 (t) = (1 − t)b11 (t) + tb12 (t)
Por último, a partir de las anteriores se define la función vectorial b30 que
parametriza a Γ
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 expresión explı́cita de la parametrización de la curva de Bézier
para los puntos B0 = (0, 0, 0), B1 = (0, 1, 1), B2 = (2, 3, 0), B3 = (4, 5, 2),
y represéntela usando un programa graficador.
51. Sea σ una superficie regular en R3 parametrizada por una función vectorial de
dos variables r : [a, b]×[c, d] −→ R3 dada por r(u, v) = hx(u, v), y(u, v), z(u, v)i,
y P (x(u0 , v0 ), y(u0 , v0 ), z(u0 , v0 )) un punto de la superficie. Consideremos el
vector unitario normal a σ en P dado por
rv (u0 , v0 ) × ru (u0 , v0 )
N(u0 , v0 ) = .
krv (u0 , v0 ) × ru (u0 , v0 )k
Entonces, la curvatura gaussiana de σ en P está dada por
(N · ruu )(N · rvv ) − N · ruv

K(P ) = ,
krv × ru k2
57
todas las funciones calculadas en (u0 , v0 ). Por otro lado, la curvatura media
de σ en P está dada por
krv k2 (N · ruu ) + kru k2 (N · rvv ) − 2(ru · rv )(N · ruv )

H(P ) = ,
2krv × ru k2
todas las funciones calculadas en (u0 , v0 ).
a) Calcule las curvaturas gaussiana y media de un plano.

b) Calcule las curvaturas gaussiana y media de la superficie de una esfera.
c) Defina esfera osculadora de una superficie en un punto.
d ) Encuentre la esfera osculadora de un paraboloide en su vértice.
52. Determine la ecuación del plano tangente a la superficie parametrizada por
r(u, v) = h(2 + sin v) cos u, (2 + sin v) sin u, cos vi,
u ∈ [0, 2π], v ∈ [0, 2π], en el punto r(0, 0). ¿Qué superficie parametriza r?
53. Determine el vector normal unitario y la ecuación del plano tangente a la

superficie parametrizada por
r(u, v) = h3 cos u, 3 sin u, vi,
u ∈ [0, 2π], v ∈ [−4, 4], en el punto r(π, 2), ¿Qué superficie parametriza r?
Detemine todos los puntos de intersección del plano tangente con la superficie.
54. La gráfica de una función f de dos variables se puede parametrizar usando

la parametrización r(u, v) = hu, v, f (u, v)i tomando una 2-celda convenien-
te como dominio de la parametrización. Parametrice la gráfica de la función
f (x, y) = x2 − y 2 y encuentre la ecuación del plano tangente a la gráfica en el
punto (1, 2, −3).
55. Determine la ecuación del plano tangente a la superficie parametrizada por

r(u, v) en el punto r(p, q) para (p, q) dado.
a) r = hv sin u, v cos u, vi, (p, q) = (π/4, 2).

b) r = hcos u, sin u, vi, (p, q) = (π/4, 3)
c) r = hv sin u, v cos u, uvi, (p, q) = (π/3, 1)
d ) r = hsin v sin u, cos v sin u, cos ui, (p, q) = (π/4, π/4)

e) r = hev sin u, ev cos u, ev i, (p, q) = (π, 1)
f ) r = hv cos u, v sin u, vi, (p, q) = (π/2, 1)
g) r = hcos u, sin u, vi, (p, q) = (π/2, 1)
h) r = hv sin u, v 2 , v cos ui, (p, q) = (π/4, 1)
i ) r = hsin v sin u, cos v sin u, cos ui, (p, q) = (π/3, π/4)
j ) r = hsin u cosh v, sin u sinh v, cos ui, (p, q) = (π, ln2)
56. La gráfica de una función f es una superficie que tiene ecuación z = f (x, y). De-
termine la ecuación del plano tangente a la gráfica de f en el punto (p, q, f (p, q)),
parametrizándola convenientemente (considere como dominio de f una 2-celda
que contenga al punto (p, q)).
a) f (x, y) = x2 − y, (p, q) = (2, 1)

b) f (x, y) = 3x + 2y + 1, (p, q) = (0, 0)
c) f (x, y) = xexy , (p, q) = (2, 0)
d ) f (x, y) = x2 + y 2 , (p, q) = (1, 1)
e) f (x, y) = xy, (p, q) = (1, 1)
f ) f (x, y) = ln(x2 + y 2 ), (p, q) = (1, 1)
LECCIÓN 4
Funciones de varias variables
Las funciones de las que trata esta lección son funciones definidas en algún subcon-
junto D de Rn y que toman valores en R. Una tal función f es una función escalar
que a cada punto (x1 , x2 , . . . , xn ) de D le asigna un número real f (x1 , x2 , . . . , xn ).
En este texto, solamente estudiaremos funciones escalares de dos o tres variables, es
decir, funciones definidas en subconjuntos de R2 o de R3 que toman valores en R. Los
dominios de estas funciones pueden ser de diferente naturaleza. Sin embargo, hay
ciertos dominios que juegan un papel muy importante en el estudio de propiedades
de estas funciones. Estudiaremos algunas caracterı́sticas de dichos conjuntos en R2 .
(Todas las definiciones se extienden de manera natural a R3 ).
Dado un subconjunto D de R2 definiremos el interior de D, el exterior de D y la

frontera de D usando el concepto de vecindad de un punto en el plano.
Si P es un punto de R2 y α es un número real positivo, la vecindad de P de tamaño

α, es el conjunto

V (P, α) = Q ∈ R2 : d(P, Q) < α ,
cuya representación en el plano es un disco con centro en P y radio α, excluyendo

los puntos de la circunferencia. Un punto P es un punto interior de D si existe una
59
60 LECCIÓN 4. FUNCIONES DE VARIAS VARIABLES
vecindad de P que esté contenida en D. El conjunto de todos los puntos interio-

res de un conjunto D se denomina el interior de D y se escribe intD. Un punto
P es un punto exterior de D si existe una vecindad de P que no contiene puntos
de D. El conjunto de todos los puntos exteriores de un conjunto D se denomina el
exterior de D y se escribe extD. Un punto P es un punto frontera de D si toda ve-
cindad de P contiene puntos de D y puntos que no están en D. El conjunto de todos
los puntos frontera de un conjunto D se denomina la frontera de D y se escribe froD.
Todo subconjunto de R2 tiene interior, exterior y frontera, que son subconjuntos

disjuntos de R2 y cuya unión es todo R2 .
Un conjunto D es abierto si todos sus puntos son interiores, es decir, si D ⊂ intD. D

es cerrado si contiene todos sus puntos frontera, es decir si froD ⊂ D. D es acotado
si existe un punto P y una vecindad V (P, α) tal que D ⊂ V (P, α).
En la mayorı́a de los libros de texto las funciones se dan por su regla de asignación z =
f (x, y) y no se hace mención explı́cita de su dominio. En este caso, se supondrá que el
dominio de f es el conjunto de puntos (x, y) para los cuales el valor f (x, y) está bien
definido, es decir, para los cuales f (x, y) es un número real. Explı́citamente:
domf = {(x, y) : f (x, y) ∈ R} .
La imagen de f es el conjunto de todos valores que toma la función al evaluarla en

todos los puntos de su dominio. Es decir,
imf = {f (x, y) : (x, y) ∈ domf } .
La gárfica de f es el conjunto
grff = {(x, y, z) : (x, y) ∈ domf, z = f (x, y)} .
En general, para las funciones que se tratarán en el curso, la representación del

conjunto grff en el espacio es una superficie de ecuación z = f (x, y). Por otro lado,
si se tiene una ecuación de la forma z = f (x, y), su gráfica (que en general es una
superficie), es la gráfica de la función f .
61
Para hacer el bosquejo de la gráfica de una función se puede recurrir a la técnica de

las trazas, que se usó para graficar una superficie a partir de su ecuación. En este
caso, las trazas que se obtienen al intersectar la gráfica de la función f con planos
paralelos al plano XY (esto es planos con ecuación z = k) toman el nombre de
curvas de nivel1 de la función f :
Nk f = {(x, y) : f (x, y) = k} es la curva de nivel k de f,
que tiene ecuación f (x, y) = k. Nótese que el conjunto Nk f es diferente de vacı́o si

y sólo si k ∈ imf .
La representación de varias de las curvas de nivel de f en el plano XY es un mapa
de contorno de f y nos permite construir la gráfica de la función.
p
Veamos el siguiente ejemplo. Dada la función f (x, y) = 1 − x2 − y 2 , f (x, y) ∈ R
si y solo si 1 − x2 − y 2 ≥ 0, es decir si x2 + y 2 ≤ 1, ası́

domf = (x, y) : x2 + y 2 ≤ 1 ,
que es el disco cerrado con centro en (0, 0) y radio 1. Las curvas de nivel de f ,
p
Nk f = {(x, y) : 1 − x2 − y 2 = k} = {(x, y) : x2 + y 2 = 1 − k 2 },
√
si 0 ≤ k ≤ 1 son circunferencias con centro en (0, 0) y radio 1 − k 2 , de donde, un
mapa de contorno de f es:
Y
k=0
k = 2/3
k = 5/6
k = 13/14
1
Si la función f es una función de tres variables, Nk f {(x, y, z) : f (x, y, z) = k} es la superficie
de nivel k de f . La ecuación de esta superficie es f (x, y, z) = k
De aquı́, se puede reconstruir la gráfica de la función f , que no es otra cosa que una
semiesfera:
x y
Veamos ahora las caracterı́sticas del dominio de f . Se puede verificar que si D =domf :
1. Todos los puntos (x, y) tales que x2 + y 2 < 1 son puntos interiores de D , es
decir,
intD = {(x, y) ∈ R2 : x2 + y 2 < 1}.
2. Todos los puntos (x, y) tales que x2 + y 2 > 1 son puntos exteriores de D, es
decir,
extD = {(x, y) ∈ R2 : x2 + y 2 > 1}.
3. Todos los puntos (x, y) tales que x2 + y 2 = 1 son puntos frontera de D, es

decir,
froD = {(x, y) ∈ R2 : x2 + y 2 = 1}.
Como todos los puntos frontera de D pertenecen a D, se tiene que D es un conjunto

cerrado. Por otro lado, D está contenido en la vecindad con centro en (0, 0) y radio
2, por lo tanto D es un conjunto acotado. D no es un conjunto abierto ya que hay
puntos de D que no son interiores (¿cuáles?)
63
Lı́mites
Para extender el concepto de lı́mite, usaremos la idea de los lı́mites laterales de

funciones de una sola variable. En otras palabras, recordemos que el lı́mite de una
función de una variable en un punto existe si los lı́mites laterales (por la izquierda y
por la derecha) existen y son iguales. Es decir, para estudiar el lı́mite de la función
en un punto, nos acercamos al punto por dos caminos posibles: por la izquierda y
por la derecha. Esta idea es un poco más complicada para funciones de dos (o más)
variables, ya que hay muchı́simas formas de acercarse a un punto en el plano. Sea
f : D −→ R, con D ⊂ R2 y P (a, b) un punto en el plano, que no necesariamente
pertenece a D. Un acercamiento a P en D es una función vectorial aP : [−1, 0) −→
R2 con las siguientes dos propieedades:
i) aP (t) ∈ D − {P }, para todo t ∈ [−1, 0),
ii) lı́mt→0− aP (t) = OP.
Obsérvese que el item (i) dice que un acercamiento es un camino que conduce a P
usando puntos de D y que el item (ii) dice que lı́mt→0− x(t) = a y lı́mt→0− y(t) = b.
Ahora, si existe por lo menos un acercamiento a P podemos definir el lı́mite de f
cuando nos acercamos a P (a, b). Diremos que L es el lı́mite de f cuando (x, y) se
acerca a (a, b) y escribiremos lı́m(x,y)→(a,b) f (x, y) = L, si para todo acercamiento
aP se tiene que
lı́m f (aP (t)) = L.
t→0−
Si para algún acercamiento aP se tiene que lı́mt→0− f (aP (t)) = L, no podemos

asegurar que lı́m(x,y)→(a,b) f (x, y) exista, pero podemos asegurar que si existiese,
éste deberı́a ser L. Por ejemplo, consideremos la función f : R2 − {(0, 0)} −→ R
3x2 y
definida por f (x, y) = 2 y tomemos el acercamiento a P (0, 0), aP (t) = ht, ti.
x + y2
Entonces,
3t3 3t
lı́m f (aP (t)) = lı́m 2 = lı́m =0
t→0− t→0− 2t t→0− 2
Este cálculo no demuestra la existencia del lı́mite. Sin embargo, nos dice que si e-
xistiese el lı́mite deberı́a ser 0. Veamos que efectivamente existe y es 0. Sea aP (t) =
hx(t), y(t)i cualquier acercamiento a P , es decir, el acercamiento es tal que lı́mt→0− x(t) =
0 y lı́mt→0− y(t) = 0 . Entonces,

3x(t)2 y(t) x(t)2

0 ≤ |f (aP (t))| = =3 |y(t)| ≤ 3 |y(t)|
x(t)2 + y(t)2 x(t)2 + y(t)2
Usando el Teorema del Emparedado y el hecho de que lı́mt→0− y(t) = 0, se tiene que
lı́mt→0− f (aP (t)) = 0 y por consiguiente lı́m(x,y)→(0,0) f (x, y) = 0.
De acuerdo con la definición de lı́mite, lı́m(x,y)→(a,b) f (x, y) no existirá si encontramos

algún acercamiento aP para el que lı́mt→0− f (aP (t)) no exista, o si encontramos por
lo menos dos acercamientos a1P y a2P para los que
lı́m f (a1P (t)) 6= lı́m f (a2P (t))

t→0− t→0−
xy
Por ejemplo, si f : R2 − {(0, 0)} −→ R definida por f (x, y) = y tomamos
x2
+ y2
los acercamientos a P (0, 0), a1P (t) = ht, 0i y a2P (t) = ht, ti, se tiene que
t2 1
lı́m f (a1P (t)) = lı́m 0 = 0 6= lı́m f (a2P (t)) = lı́m 2
= ,
t→0− t→0− t→0− t→0− 2t 2
lo que muestra que el lı́mite no existe.
Una consecuencia inmediata de la definición de lı́mite, es que todas las propiedades

de los lı́mites que se estudiaron en cursos anteriores de cálculo son válidas: lı́mite
de una suma es la suma de los lı́mites, lı́mite de un producto es el producto de los
lı́mites, etc. Además, el concepto de continuidad se extiende a funciones de varias
variables de manera natural.
Sea f una función de D ⊂ R2 en R, P (a, b) ∈ D y suponga que exista por lo

menos un acercamiento a P en D . La función f es continua en P si y solamente
si lı́m(x,y)→(a,b) f (x, y) = f (a, b) (Si no existe ningún acercamiento a P aceptaremos
que f es continua en P ). Esta definición trae como consecuencia que la continuidad
de dos funciones en un punto implica la continuidad de la suma y del producto en
ese punto y las demás propiedades de las funciones continuas que se estudiaron en
cursos anteriores de cálculo.
Veamos algunos ejemplos.

65
1. Sea f : R2 −→ R, definida por f (x, y) = x. Veamos que
lı́m f (x, y) = a.
(x,y)→(a,b)
En efecto, sea aP (t) = hx(t), y(t)i un acercamiento a P (a, b). Es decir, lı́mt→0− x(t) =
a y lı́mt→0− y(t) = b. Entonces,
lı́m f (aP (t)) = lı́m f (x(t), y(t)) = lı́m x(t) = a

t→0− t→0− t→0−
Obsérvese que como lı́m(x,y)→(a,b) f (x, y) = a = f (a, b), la función f es continua

en (a, b).
2. Sea f : R2 −→ R, definida por f (x, y) = y. Vemos que
lı́m f (x, y) = b = f (a, b).

(x,y)→(a,b)
Es decir, f es continua en (a, b)
3. Sea f : R2 −→ R una función polinómica, es decir, f (x, y) es la suma de

términos de la forma cxn y m . f es continua en (a, b) pues ésta se obtiene de
productos y sumas de las funciones de los ejemplos 1 y 2 que son continuas en
(a, b).
4. Las funciones racionales, es decir, cocientes de funciones polinómicas, son con-

tinuas en todo punto de su dominio. Esto es una consecuencia de la continuidad
del cociente, cuando el numerador y el denominador lo son.
EJERCICIOS
1. Defina los conceptos de vecindad, interior, exterior, frontera, conjunto abierto,

conjunto cerrado y conjunto acotado en R3 .
2. Sea f (x, y) = ln(x + y − 1)
a) Evalúe f (1, 1) y f (e, 1).

b) Determine el dominio de f .
c) Determine, si existen, puntos del dominio cuya imagen sea cada uno de
los siguientes números reales: 0, 3, 10.000, -1.325.478, α.
d ) Determine la imagen de f .
e) Represente gráficamente el dominio D de f : ¿es abierto D? ¿es cerrado
D? ¿no es ni abierto ni cerrado D? ¿es acotado D? ¿Por qué?
2 −y
3. Sea f (x, y) = ex
a) Evalúe f (2, 4) y f (0, 1).

b) Determine el dominio de f .
los siguientes números reales: 0, 1, 10.000, -20, α.
d ) Determine la imagen de f .
e) Represente gráficamente el dominio D de f : ¿es abierto D? ¿es cerrado
p
4. Sea g(x, y) = 36 − 9x2 − 4y 2
a) Evalúe g(1, 2) y g(0, 1).

b) Determine el dominio de g.
√
los siguientes números reales: 0, 3, -1, 1235, -π, 2, α.
d ) Determine la imagen de g.
e) Represente gráficamente el dominio D de g: ¿es abierto D? ¿es cerrado
67
p
5. Sea g(x, y) = (xy 2 − 1)(x2 + y 2 − 4)
a) Evalúe g(1, 2) y g(−1, 1.5).

b) Determine el dominio de g.
√
los siguientes números reales: 0, 3, -1, 1235, π, 2, α.
d ) Determine la imagen de g.
e) Represente gráficamente el dominio D de g: ¿es abierto D? ¿es cerrado
6. Sea f (x, y, z) = x2 ln(x − y + z)
a) Evalúe f (3, 6, 4)
b) Determine el dominio D de f .
c) Determine la imagen de f .
d ) ¿Es abierto D? ¿es cerrado D? ¿no es ni abierto ni cerrado D? ¿es acotado
D? ¿Por qué?
1
7. Sea f (x, y, z) = p
x + y + z2 − 1
2 2
a) Evalúe f (1, 3, −4)

b) Determine el dominio D de f .
c) Determine la imagen de f .
d ) ¿Es abierto D? ¿es cerrado D? ¿no es ni abierto ni cerrado D? ¿es acotado
D? ¿Por qué?
8. Determine 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) = x+y
√ √
b) f (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
j ) f (x, y, z) = ln(16 − 4x2 − 4y 2 − z 2 )
1
k ) f (x, y) = x2 + y 2 − 1 2
x
1
l ) f (x, y) = 2
x + y2
p x2
m) f (x, y) = x2 + y 2 − 4
y
x2 + 1
n) f (x, y) =
y
ñ) f (x, y, z) = sin(xyz)
p
o) f (x, y, z) = (x − 1)(yz − 1)
9. Dé por lo menos dos acercamientos convenientes, calcule los lı́mites al restringir
la función a los acercamientos que dió, Determine el lı́mite de las funciones si
existen, o demuestre que el lı́mite no existe
a) lı́m(x,y)→(5,−2) (x5 + 4x3 y − 5xy 2 )

b) lı́m(x,y)→(6,3) xy cos(x − 2y)
x2
c) lı́m(x,y)→(0,0)
x2 + y 2
(x + y)2
d ) lı́m(x,y)→(0,0)
x2 + y 2
8x2 y 2
e) lı́m(x,y)→(0,0)
x4 + y 4
x3 + xy 2
f ) lı́m(x,y)→(0,0)
x2 + y 2
69
xy
g) lı́m(x,y)→(0,0) p
x2 + y 2
xy + 1
h) lı́m(x,y)→(0,0) 2
x + y2 + 1
2x2 y
i ) lı́m(x,y)→(0,0) 4
x + y2
x3 y 2
j ) lı́m(x,y)→(0,0) 2
x + y2
x2 + y 2
k ) lı́m(x,y)→(0,0) p
x2 + y 2 + 1 − 1
xy − 2y
l ) lı́m(x,y)→(0,0) 2
x + y 2 − 4x + 4
m) lı́m(x,y,z)→(3,0,1) e−xy sin(πz/2)
x2 + 2y 2 + 3z 2
n) lı́m(x,y,z)→(0,0,0)
x2 + y 2 + z 2
xy + yz + xz
ñ) lı́m(x,y,z)→(0,0,0)
x2 + y 2 + z 2
10. Determine el domino de la función, y diga si ésta es continua en su dominio:
1
a) F (x, y) =
x2
−y
x−y
b) F (x, y) = 2
x + y2 + 1
√
c) F (x, y) = arctan(x + y)
d ) F (x, y) = ln(2x + 3y)
p p
e) G(x, y) = (x + y) − (x − y)
f ) G(x, y) = sin−1 (x2 + y 2 )
xyz
g) f (x, y, z) = 2
x + y2 − z
√
h) f (x, y, z) = x + y + z

2 3
 x y

si (x, y) 6= (0, 0)
i ) f (x, y) = 2x2 + y 2

1 si (x, y) = (0, 0)
 xy
 si (x, y) 6= (0, 0)
2 2
j ) f (x, y) = x + xy + y
0 si (x, y) = (0, 0)
11. Suponga que lı́m f (x, y) y lı́m g(x, y) existen, y que k es un escalar.
(x,y)→(a,b) (x,y)→(a,b)
Demuestre que:

(a) lı́m
[f (x, y) + g(x, y)] = lı́m f (x, y) + lı́m g(x, y)
(x,y)→(a,b) (x,y)→(a,b) (x,y)→(a,b)

(b) lı́m k f (x, y) = k lı́m f (x, y)
(x,y)→(a,b) (x,y)→(a,b)

(c) lı́m [f (x, y)g(x, y)] = lı́m f (x, y) lı́m g(x, y)
(x,y)→(a,b) (x,y)→(a,b) (x,y)→(a,b)
lı́m f (x, y)
f (x, y) (x,y)→(a,b)
(d) lı́m = si lı́m g(x, y) 6= 0
(x,y)→(a,b) g(x, y) lı́m g(x, y) (x,y)→(a,b)
(x,y)→(a,b)
(e) Si |f (x, y) − L| ≤ g(x, y) para todo (x, y) y si lı́m g(x, y) = 0, enton-

(x,y)→(a,b)
ces lı́m f (x, y) = L.
(x,y)→(a,b)
12. Calcule los siguientes lı́mites:
a) lı́m cos(xy) b) lı́m exy

(x,y)→(0,0) (x,y)→(0,0)
x2 − y 2 xy 2
c) lı́m d) lı́m
(x,y)→(0,0) x2 + y 2 (x,y)→(0,0) x2 + y 4
x2 − 2xy + y 2 xy 2
e) lı́m f) lı́m
(x,y)→(1,−1) x−y (x,y)→(0,0) x2 + y 2
x2 − y 2 x2 − 2xy + y 2
g) lı́m h) lı́m
(x,y)→(1,1) x−y (x,y)→(0,0) x−y

y 4 sin(xy) 2 2 1
i) lı́m j) lı́m (x + y ) cos
(x,y)→(0,0) x2 + y 2 (x,y)→(0,0) xy
x
k) lı́m 1
(x,y)→(0,0) y l) lı́m cos
(x,y)→(0,0) xy
1 2 2 2
−(x +y )/2σ , para σ > 0, es constante en la circun-
13. Muestre que f (x, y) = 2πσ 2e
ferencia de radio r > 0 centrada en el origen. Esta función se llama Gaussian

blur, y se usa como filtro en los programas de procesamiento de imágenes para
producir un efecto de “desenfoque”.
71
14. 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 .
15. La advertencia de Thomas. (Ver George B. Thomas, Jr., Calculo en varias

variables, 12 Ed., Addison-Wesley, 2010, página 763) En algunas ocasiones se
recurre a escribir la función en coordenadas polares para obtener el lı́mite de f
en (0, 0), es decir, se escribe x = r cos θ, y y = r sen θ y se estudia el lı́mite de la
función cuando r tiende a 0. De acuerdo con la definición de lı́mite de funciones
de una sola variable, se trata de determinar si existe L tal que, para todo ǫ > 0
hay un δ > 0 tal que si 0 < |r| < δ, se tenga |f (x(r, θ), y(r, θ)) − L| < ǫ, para
todo θ. Por ejemplo,
xy 2 r3 cos θ sin2 θ
lı́m = lı́m = lı́m r cos θ sin2 θ = 0.
(x,y)→(0,0) x2 + y 2 r→0 r2 r→0
Sin embargo, debemos comprobar esta afirmación recurriendo a la definición.

xy 2
En este caso, f (x, y) = 2 = r cos θ sin2 θ, si r 6= 0. Es decir, necesitamos
x + y2
demostrar que, para todo ǫ > 0 hay un δ > 0 tal que, si 0 < |r| < δ, se tenga
|r cos θ sin2 θ − L| < ǫ, para todo θ. Se puede ver que esto es cierto, si tomamos
δ = ǫ ya que | cos θ sin2 θ| ≤ 1 para todo θ.
x2
a) Estudie el lı́mite de f (x, y) = en (0, 0) cambiando las variables a
x2
+ y2
coordenadas polares. ¿Qué se puede concluir en este caso?
b) Los ejemplos anteriores muestran que en algunas ocasiones cambiar a

coordenadas polares puede ser útil. Sin embargo, hay que proceder con
mucha cautela al usar este procedimiento ya que puede conducir a erro-
res si se usa de manera mecánica. Considere la función f definida por
2x2 y
f (x, y) = 4 y estudie su lı́mite en (0, 0) usando acercamientos y
x + y2
cambiando las variables a coordenadas polares. Compare sus respuestas
y explique por qué la diferencia en los resultados.
c) Cambie las variables a coordenadas polares para mostrar que
p
sin x2 + y 2
lı́m p = 1.
(x,y)→(0,0) x2 + y 2
LECCIÓN 5
Diferenciación Parcial
Una de las herramientas más poderosas del cálculo es la diferenciación, con la que se
estudia la forma en que los valores de una función varı́an con respecto a su variable,
mediante el cálculo de la razón de cambio instantánea en puntos de su dominio.
En la lección 3, recurrimos a la diferenciación para estudiar el comportamiento
de funciones vectoriales, lo que nos permitió estudiar el movimiento de objetos en
el espacio y la geometrĂa de las curvas y las superficies. Para poder extender la
diferenciación a funciones de varias variables recurriremos al estudio de la variación
de funciones en direcciones determinadas. Es ası́ que le damos sentido a la razón
de cambio instantanea de una función de varias variables en un punto de su dominio.
Con el propósito de precisar lo que significa la razón de cambio instantanea de una
función en un punto en una dirección determinada introduciremos a continuación el
concepto de paso por el punto P en la dirección de un vector v . Sean P (a, b) un
punto en el plano, ǫ un número real positivo y v = hv1 , v2 i un vector unitario. El
paso por el punto P en la dirección del vector v es la función vectorial
pvP : (−ǫ, ǫ) −→ R2
dada por
pvP (t) = ha + v1 t, b + v2 ti .
73
74 LECCIÓN 5. DIFERENCIACIÓN PARCIAL
Sea f : D −→ R una función definida en un subconjunto abierto D del plano, P (a, b)

un punto de D y v = hv1 , v2 i un vector unitario. La derivada direccional de f en el
df
punto P (a, b) en la dirección de v , que se escribe como (a, b) o como Dv f (a, b),
dv
se define por:
df f (a + v1 t, b + v2 t) − f (a, b)
(a, b) = lı́m .
dv t→0 t
df
simpre y cuando este lı́mite exista. Obsérvese que (a, b) no es otra cosa que la
v
dv
derivada de f ◦ pP en t = 0:

df d v d
(a, b) = f (pP (t)) = f (a + v1 t, b + v2 t).
dv dt t=0 dt t=0
Claramente (−ǫ, ǫ), dominio del paso, debe ser tal que (a + v1 t, b + v2 t) ∈ D para
todo t. Podemos definir la función derivada direccional de f en la dirección de v que
df
se nota como cuyo dominio es el conjunto de todos los puntos (x, y) de D para
dv
los cuales existe la derivada de f en (x, y) en la dirección de v.
Como para cada punto P (x, y) y dirección fija v la derivada direccional es la derivada
de una función de una sola variable, es natural que se cumplan todas las propiedades
algebraicas de las derivadas de funciones de una variable. Por ejemplo, se tiene que
d(f + g) df dg
(x, y) = (x, y) + (x, y)
dv dv dv
y
d(f g) df dg
(x, y) = g(x, y) + f (x, y).
dv dv dv
Las derivadas direccionales de f en las direcciones de los vectores i y j se conocen
como las derivadas parciales de f 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 más adelante todas las derivadas direccionales de funciones diferen-
ciables se van a poder expresar en términos 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 D√ √ E
la derivada (x, y) en los casos en que v = 22 , 22 , v =i y v =j . Tenemos que :
dv
75
D√ √ E D √ √ E
2 2 2 2
En el caso en que v = 2 , 2 , pvP (t) = x + 2 t, y + 2 t y

df d √
2 3
√
2 2
(x, y) = (x + 2 t) − 2(y + 2 t)
dv dt t=0
√ √ √ √
2 2 2 2 2
= 3 2 (x + 2 t) −4 2 (y + 2 t)t=0
√
2 2
√
= 3 2 x − 2 2y,
en el caso en que v =i ,

df d
(x, y) = (x + t)3 − 2y 2
di dt t=0

= 3(x + t)2 t=0
∂f
= 3x2 = (x, y),
∂x
y en el caso en que v =j ,

df d
(x, y) = x3 − 2(y + t)2
dj dt t=0
= −4(y + t)|t=0
∂f
= −4y = (x, y).
∂y
Obsérvese que las derivadas parciales se pueden calcular derivando la función f co-
mo si fuese de una sola variable, fijando como constante la otra. Ası́ por ejemplo, si
queremos calcular la derivada parcial con respecto a x, basta derivar con respecto a
∂f ∂f
x fijando y como si fuese constante. Las derivadas parciales y de f se suelen
∂x ∂y
escribir también como fx , fy , respectivamente. En lo que sigue supondremos que las
funciones fx y fy son continuas en D.
Considere un paso por un punto P (a, b) del dominio D de una función f en la

dirección de un vector unitario v . La curva de ecuaciones paramétricas
x = a + v1 t, y = b + v2 t, z = f (a + v1 t, b + v2 t), t ∈ (−ǫ, ǫ),
está contenida en la gráfica

de f . El vector tangente a esta curva en el punto
df
Q(a, b, f (a, b)) es v e = v1 , v2 , (a, b) , que también es tangente a la gráfica de
dv
la función f en el punto Q(a, b, f (a, b)). La pendiente de la recta tangente, es decir,
de la que tiene vector director v e y pasa por Q(a, b, f (a, b)), con respecto al vector
df
hv1 , v2 , 0i es precisamente (a, b).
dv
Hay dos de estos vectores tangentes que nos interesan. Éstos son ei = h1, 0, fx (a, b)i
y ej = h0, 1, fy (a, b)i que nos dan las pendientes de las rectas tangentes a la gráfica
de f en el punto Q(a, b, f (a, b)) con respecto a i y j , respectivamente. A partir de
estos vectores podemos encontrar la ecuación del plano tangente a la gárfica de f
en el punto Q(a, b, f (a, b)). En efecto, el vector ei × ej = h−fx (a, b), −fy (a, b), 1i es un
vector normal al plano tangente, por lo que dicha ecuación 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).
Esta última ecuación corresponde a la ecuación de la función afı́n

(a,b)
Lf (x, y) = f (a, b) + fx (a, b)(x − a) + fy (a, b)(y − b),
que aproxima a f en una vecindad del punto Q. Esta función se usa frecuentemente
para aproximar valores de f al hacer incrementos pequeños ∆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 valores 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 pequeño cuando los incrementos ∆x y ∆y son pequeños, pero se re-
quiere, para que la aproximación sea de utilidad, que éste sea mucho más pequeño
en el siguiente sentido. Diremos que E es bien pequeño si existe hǫ1 , ǫ2 i tal que
1. E = h∆x, ∆yi · hǫ1 , ǫ2 i.

77
2. Si h∆x, ∆yi → h0, 0i se tiene que hǫ1 , ǫ2 i → h0, 0i.
No todas las funciones se pueden aproximar mediante su función afı́n en (a, b) con
un error bien pequeño. Las que satisfacen esta condición se llaman funciones dife-
renciables en (a, b). En otras palabras, la función f es diferenciable en (a, b) si E
es bien pequeño. En general, no es fácil encontrar hǫ1 , ǫ2 i, por lo que verificar si
una función es diferenciable en (a, b) se vuelve una tarea dispendiosa. Sin embargo,
contamos con el siguiente resultado:
Si las derivadas parciales de f existen en una vecindad de (a, b) y son continuas en

(a, b) entonces f es diferenciable en (a, b).
Demostraremos que existe hǫ1 , ǫ2 i tal que E = hǫ1 , ǫ2 i · h∆x, ∆yi y tal que
lı́m hǫ1 , ǫ2 i = h0, 0i .

(∆x,∆y)→(0,0)
En efecto,
E = f (a + ∆x, b + ∆y) − f (a, b) − fx (a, b)∆x − fy (a, b)∆y
= f (a + ∆x, b + ∆y) − f (a + ∆x, b) + f (a + ∆x, b) − f (a, b)

−fx (a, b)∆x − fy (a, b)∆y
= fx (α, b)∆x + fy (a + ∆x, β)∆y − fx (a, b)∆x − fy (a, b)∆y
= (fx (α, b) − fx (a, b))∆x + (fy (a + ∆x, β) − fy (a, b))∆y
= hfx (α, b) − fx (a, b), fy (a + ∆x, β) − fy (a, b)i · h∆x, ∆yi
= hǫ1 , ǫ2 i · h∆x, ∆yi
En el tercer renglón hemos usado el Teorema del Valor Medio para derivadas aplica-
do a las derivadas parciales de f . Como α → a, y β → b cuando ∆x → 0 y ∆y → 0,
y fx y fy son continuas en (a, b), se tiene que lı́m(∆x,∆y)→(0,0) hǫ1 , ǫ2 i = h0, 0i, lo que
demuestra la diferenciabilidad de f en (a, b).
La continuidad de las derivadas parciales es mucho más fácil de verificar.
Nota. Una función f de dos variables que tiene derivadas parciales en el punto
(a, b) se dice que es diferenciable en (a, b) si es bien pequeño el error que se comete al
aproximar f (a+∆x, b+∆y) mediante su aproximación afı́n Lf(a,b) (a+∆x, b+∆y) en
una vecindad de (a, b). El error E = f (a+∆x, b+∆y)−Lf(a,b) (a+∆x, b+∆y) es bien
pequeño, si éste se puede escribir en la forma E = hε1 , ε2 i · h∆x, ∆yi de tal manera
que lı́m(∆x,∆y)→(0,0) ε1 = lı́m(∆x,∆y)→(0,0) ε2 = 0. Se harán varios ejercicios en los
que se mostrará que algunas funciones son diferenciables en puntos particulares, de
forma directa o verificando la continuidad de las derivadas parciales. Sin embargo,
mostrar que una función no es diferenciable en un punto no es una tarea sencilla
a partir de la definición. Esta nota sobre diferenciabilidad nos dá dos herramientas
muy útiles para mostrar que una función no es diferenciable en un punto en caso de
que no lo sea. El primer resultado que destacamos aquı́, es que si una función f es
diferenciable en un punto (a, b) entonces ésta debe ser continua en el mismo punto.
En efecto, si f es diferenciable en (a, b), podemos escribir
f (a + ∆x, b + ∆y) = f (a, b) + fx (a, b)∆x + fy (a, b)∆y + E,
Siendo E bien pequeño, se tiene que
E = hε1 , ε2 i · h∆x, ∆yi
tiende a cero cuando ∆x y ∆y tienden a cero. Por consiguiente,
lı́m f (a + ∆x, b + ∆y) = f (a, b),

(∆x,∆y)→(0,0)
lo que muestra la continuidad de f en (a, b).
Esto quiere decir, que si una función f no es continua en en un punto (a, b), ésta no
puede ser diferenciable en ese punto. Por otro lado, si f es diferenciable en (a, b),
entonces, como
|E|
0≤ ≤ khε1 , ε2 ik,
kh∆x, ∆yik
se tiene que
|E|
lı́m = 0.
(∆x,∆y)→(0,0) kh∆x, ∆yik
79
Ası́, si
|E|
lı́m 6= 0
(∆x,∆y)→(0,0) kh∆x, ∆yik
se concluye automáticamente que f no puede ser diferenciable en (a, b).

(a,b)
La parte lineal de Lf , es decir fx (a, b)∆x + fy (a, b)∆y se conoce como diferencial
de f , se escribe df y se usa para aproximar el incremento de f cuando se hacen
incrementos ∆x y ∆y de a y de b respectivamente. Claramente, df es una muy
buena aproximación 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 ilustración del uso de las diferenciales, estimemos el error máximo 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 medición de 0.1cm. El error cometido en la medición
del volumen V = πr2 h/3 es ∆V , que por lo que hemos visto es aproximadamente
igual a
dV = Vr ∆r + Vh ∆h = (2πrh/3)∆r + (πr2 /3)∆h.
Como |∆r| ≤ 0,1 y |∆h| ≤ 0,1. Entonces,
500π 100π
|dV | = (2πrh/3)∆r + (πr2 /3)∆h ≤ (0,1) + (0,1) = 20π
3 3
Es decir, el máximo error cometido al calcular el volumen a partir de las mediciones
es de aproximadamente 63 cm3 .
Ahora bien, las derivadas parciales de una función f de varias variables son nue-
vamente funciones de varias variables, por lo que son susceptibles de ser deriva-
bles parcialmente. En este caso obtendremos derivadas parciales de orden superior.
Ası́, las derivadas
parciales de segundo
orden 2de una función f de dos variables
∂ ∂f ∂2f ∂ ∂f ∂ f ∂ ∂f ∂2f
serán: = = f xx , = = f xy , = fyx y
∂x ∂x2 ∂x2 ∂y ∂x ∂y∂x ∂x ∂y ∂x∂y
∂ ∂f ∂ f
= = fyy . Este proceso de diferenciación parcial se podrı́a continuar
∂y ∂y ∂y 2
para obtener derivadas parciales de orden tres, orden cuatro, etc.
Otro resultado importante en esta teorı́a es el siguiente:
Cuando las derivadas parciales fx , fy , fyx y fxy de una función f están definidas
en un conjunto abierto D y son continuas, se tiene que las derivadas cruzadas son
iguales, es decir, se tiene que fxy = fyx .
La justificación de esta afirmación se puede hacer mediante la aplicación sucesiva

del Teorema del Valor Medio a ciertas funciones convenientes de una sola variable.
Sea (a, b) un punto de D y sean h y k incrementos en la primera y segunda variables
de modo que el rectángulo de vértices (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 común que denotaremos Γ.
Además, 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 función 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,

81
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 querı́amos mostrar.
EJERCICIOS
1. Calcule las derivadas direccionales de las siguientes funciones en el punto dado

y la dirección dada. No olvide que la dirección la dá un vector unitario.
a) f (x, y) = x2 + y 2 en el punto (1, 1) en la dirección del vector v= h1, 1i

b) f (x, y) = cos(x + y) en el punto (0, π) en la dirección del vector v= h2, 1i
p
c) f (x, y) = x2 + y + 4 en el punto (1, 0) en la dirección del vector v=
h−1, 1i
x+1
d ) f (x, y) = en el punto (1, 5) en la dirección del vector v= h−1, −2i
y+1
∂f ∂f
2. Calcule las derivadas parciales 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
p
e) f (x, y) = x2 + y 2 f ) f (x, y) = sin(x + y)
p
g) f (x, y) = 3
x2 + y + 4 xy + 1
h) f (x, y) =
x+y
2 +y 2 )
i ) f (x, y) = e−(x j ) f (x, y) = ln(xy)
k ) f (x, y) = sin(xy) l ) f (x, y) = tan(x + y)
∂f ∂f ∂f
3. Calcule las derivadas parciales , y de las siguientes funciones:
∂x ∂y ∂z
p
a) f (x, y, z) = x − y2 + z2
b) f (x, y, z) = tanh(x + y + 3z)
c) f (x, y, z) = arcsec (x + yz)
2 +y 3 )/z
d ) f (x, y, z) = xe(x
4. Determine la ecuación del plano tangente a la superficie de ecuación z = f (x, y)

en el punto P :
83
a) f (x, y) = x2 + y 3 , P (1, 1, 2) b) f (x, y) = xy, P (1, −1, −1)
c) f (x, y) = x2 y, P (−1, 1, 1) d ) f (x, y) = xey , P (1, 0, 1)

p
e) f (x, y) = x + 2y, P (2, 1, 4) f ) f (x, y) = x2 + y 2 , P (3, 4, 5)
5. Determine la función afin de cada una de las siguientes funciones en los puntos
dados:
a) f (x, y) = x2 + y 2 en (0, 0) y en (1, 1).

b) f (x, y) = (x + y + 2)2 en (0, 0) y en (1, 2).
c) f (x, y) = ex cos(2y) en (0, 0) y en (0, π/2).
d ) f (x, y, z) = xy + xz + yz en (1, 1, 1) y en (1, 0, 0).
sin(xy)
e) f (x, y, z) = en (π/2, 1, 1) y en (2, 0, 1).
z
f ) f (x, y, z) = arctan(xyz) en (1, 0, 0) y en (1, 1, 0).
6. Sea f (x, y) = −(1 − x2 − y 2 )1/2 si x2 + y 2 < 1, muestre que el plano tan-

gente a la gárfica de f en el punto (x0 , y0 , f (x0 , y0 )) es ortogonal al vector
hx0 , y0 , f (x0 , y0 )i.
7. Sea f la función definida por f (x, y) = x2 + y 2 .

(a,b)
a) Escriba explı́citamente la función afı́n Lf de f en el punto (a, b).
b) Evalúe la función afı́n al en (a + ∆x, b + ∆y).
c) Escriba explı́citamente el error E que se comete al aproximar f (a+∆x, b+
(a,b)
∆y) por Lf (a + ∆x, b + ∆y).
d ) Encuentre hǫ1 , ǫ2 i tal que E = h∆x, ∆yi · hǫ1 , ǫ2 i
e) Muestre que E es bien pequeño y concluya que f es diferenciable en (a, b).
8. Sea f la función definida por f (x, y) = x + 2y.

(a,b)
(a,b)
∆y) por Lf (a + ∆x, b + ∆y).

x+1
9. Sea f la función definida por f (x, y) = .
y+1
(a,b)
(a,b)
∆y) por Lf (a + ∆x, b + ∆y).
10. Si f (x, y) = ϕ((x + y)/(x − y)) para alguna función diferenciable ϕ : R −→ R

muestre que xfx + yfy = 0
11. Muestre que la función


2 3
 3x y − y

si (x, y) 6= (0, 0)
f (x, y) = x2 + y

0 si (x, y) = (0, 0)
tiene derivadas direccionales en todas las direcciones en (0, 0), pero no es di-
ferenciable en (0, 0).
12. Sea f definida por

 xy
 si (x, y) 6= (0, 0)
f (x, y) = x2 + y2
0 si (x, y) = (0, 0)
a) Muestre que las derivadas parciales de f existen en todo su dominio.

b) Encuentre los puntos en que f es diferenciable y los puntos en los que f
no lo es.
13. Considere la función f definida por


2 2
 sin(x y )

si (x, y) 6= (0, 0)
f (x, y) = x2 + y 2

0 si (x, y) = (0, 0).
85
a) Explique lo que significarı́a que f sea diferenciable en (0, 0).

b) Muestre que las derivadas parciales fx (0,0) y fy (0,0) existen y calcúlelas.
c) ¿Es f diferenciable en (0, 0)?
14. Considere la función f definida por


 sin(xyz)

si (x, y, z) 6= (0, 0, 0)
f (x, y, z) = x4 + y 4 + z 4

0 si (x, y, z) = (0, 0, 0).
a) Muestre que las derivadas parciales de f existen en todo el dominio.

b) Encuentre los puntos en que f es diferenciable y los puntos en los que f
no lo es.
15. Al efectuar incrementos ∆x y ∆y en las variables x y y respectivamente, a

partir del punto (x0 , y0 ), el cambio absoluto de una función diferenciable f
es ∆f = f (x0 + ∆x, y0 + ∆y) − f (x0 , y0 ), que se puede estimar usando el
diferencial df = fx (x0 , y0 )∆x + fy (x0 , y0 )∆y, ya que E = ∆f − df es un error
∆f
bien pequeño. Ası́, el cambio relativo, definido por , se puede estimar
f (x0 , y0 )
df ∆f
mediante y el cambio porcentual, definido por × 100, se
f (x0 , y0 ) f (x0 , y0 )
df
puede estimar mediante × 100.
f (x0 , y0 )
a) Suponga que el radio r0 = 1 y la altura h0 = 5 de un cilindro circular recto

experimentan variaciones ∆r = 0, 03 unidades y ∆h = −0, 1 unidades.
Estime los cambios absoluto, relativo y porcentual de su volumen.
b) El volumen V = πr2 h de un cilindro circular recto va a calcularse a partir
de los valores medidos de r y h. Si r se mide con un error de no más del
2 % y h con un error de no más de 0,5 %, estime el posible porcentaje de
error resultante en el cálculo de V .
c) El área de un lote triangular es (1/2)ab sin C, donde a y b son las lon-
gitudes de dos de sus lados y C es la medida del ángulo formado por
esos dos lados. Al efectuar medidas topográficas 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 cálculo del área del lote si la me-

dida a, b de sus lados tienen errores de 15 cm cada una y C un error de
2◦ .
16. Calcule todas las derivadas parciales de segundo orden de las siguientes fun-
ciones.
a) f (x, y) = sin(x2 − 3xy)

b) f (x, y) = x2 y 2 e2xy
c) f (x, y) = arctan xy
17. Una función u = f (x, y) es una función armónica si sus derivadas parciales
2 2
segundas existen, son continuas y satisfacen la ecuación de Laplace ∂∂xu2 + ∂∂yu2 =
0. Muestre que la función f (x, y) = x3 − 3xy 2 es armónica.
xy(x2 −y 2 )
18. Sea f (x, y) = x2 +y 2
si (x, y) 6= (0, 0) y f (0, 0) = 0. Calcule
∂f ∂f
a) ∂x y ∂y si (x, y) 6= (0, 0)
∂f ∂f
b) Verifique que ∂x (0, 0) = ∂y (0, 0) = 0 (usando la definición).
∂2f ∂2f
c) Verifique que ∂x∂y (0, 0) = 1 y que ∂y∂x (0, 0) = −1.
d ) ¿Qué pasa con las derivadas parciales cruzadas?
∂2f ∂2f ∂2f
19. Encuentre ,
∂x2 ∂y 2
y ∂y ∂x :
a) f (x, y) = x2 + y 2 b) f (x, y) = cos(x + y)

p
c) f (x, y) = x2 + y + 4 x+1
d ) f (x, y) =
y+1
e) f (x, y) = exy + xy f ) f (x, y) = x2 − y 2 + 6xy + 4x − 8y + 2
g) f (x, y) = x4 h) f (x, y) = x + 2y
i ) f (x, y) = ln(xy) j ) f (x, y) = sin(xy)
20. Muestre que la función f (x, y) = sin(x + y) + cos(x − y) satisface la ecuación

de onda
∂2f ∂2f
− = 0.
∂x2 ∂y 2
87
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 solución
de la ecuación general de onda unidimensional.
∂2f 1 ∂2f
− = 0.
∂x2 c2 ∂y 2
LECCIÓN 6
Regla de la cadena
En general, en cálculo se trabaja con funciones que se obtienen de componer, sumar

o multiplicar las funciones elementales. En el caso de funciones de una sola variable
se tiene un instrumento muy valioso para el cálculo de la derivada de una función
compuesta. Este instrumento es la Regla de la Cadena. En esta lección generaliza-
remos la regla de la cadena para el caso en que la función que se quiere derivar es
la compuesta de una función vectorial con una función escalar.
Supongamos que f : D −→ R es diferenciable en todos los puntos de su dominio

(que supondremos abierto) y que la función vectorial r : I −→ D, dada por r(t) =
hx(t), y(t)i, es derivable en I. La compuesta de estas dos funciones es una función
z : I −→ R, dada por z(t) = f (r(t)) y estamos interesados en estudiar la derivada
de z con respecto a t. Si hacemos un cambio ∆t en la variable t obtenemos cambios
∆x = x(t + ∆t) − x(t) y ∆y = y(t + ∆t) − y(t), ası́, usando el hecho de que f es
89
90 LECCIÓN 6. REGLA DE LA CADENA
diferenciable, tenemos
∆z z(t + ∆t) − z(t)
=
∆t ∆t
f (x(t + ∆t), y(t + ∆t)) − f (x(t), y(t))

=
∆t
fx (x(t), y(t))∆x + fy (x(t), y(t))∆y + E

=
∆t
∆x ∆y E
= fx (x(t), y(t)) + fy (x(t), y(t)) + ,
∆t ∆t ∆t
donde E = ∆f − df que es bien pequeño ya que f es diferenciable. Como

E ∆x ∆y dx dy
lı́m = lı́m hǫ1 , ǫ2 i · , = h0, 0i · , = 0,
∆t→0 ∆t ∆t→0 ∆t ∆t dt dt
se obtiene
dz dx dy
= fx + fy ,
dt dt dt
que es la Regla de la Cadena en esta situación 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, más explı́citamente,
dz
(t) = ∇f (r(t)) · r′ (t).
dt
Teniendo en cuenta la regla de la cadena, podemos reformular la derivada direccional
de una función diferenciable f en la dirección de un vector unitario v ası́:

df d
(x, y) = f (x + v1 t, y + v2 t)
dv dt t=0
= ∇f (x + v1 t, y + v2 t) · hv1 , v2 i|t=0
= ∇f (x, y) · v
91
A partir de esta última, se deducen varias propiedades del vector gradiente. En

efecto,
df
(x, y) = ∇f (x, y) · v
dv
= k∇f (x, y)k kvk cos θ
= k∇f (x, y)k cos θ,
siendo θ el ángulo entre ∇f (x, y) y v. De aquı́ que, la dirección en la que la derivada

df
direccional (x, y) es máxima es aquella en la que v tiene la misma dirección y
dv
df
sentido de ∇f (x, y). En este caso, (x, y) = k∇f (x, y)k es la razón de máximo
dv
crecimiento de la función en el punto (x, y) . De igual manera, en el sentido opuesto
al del gradiente se obtiene el máximo decrecimiento de la función en el punto y
df
(x, y) = − k∇f (x, y)k es la razón de máximo decrecimiento de la función en el
dv
punto (x, y).
Por otro lado, si r = r(t) es una parametrización 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 más variables. Por ejemplo, mediante un
razonamiento muy similar al último aquı́ presentado, se puede concluir que para una
función diferenciable f : D ⊆ 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 per-
mite determinar la ecuación del plano tangente a Nkf en el punto (x0 , y0 , z0 ) ∈ Nkf ,
ya que un vector normal al plano tangente en ese punto es ∇f (x0 , y0 , z0 ). Explı́cita-
mente, la ecuación de dicho plano es
fx (x0 , y0 , z0 )(x − x0 ) + fy (x0 , y0 , z0 )(y − y0 ) + fz (x0 , y0 , z0 )(z − z0 ) = 0.
Supóngase ahora que f es una función 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
incremento ∆s en la variable s obtenemos:

∆z = z(s + ∆s, t) − z(s, t)
= f (x(s + ∆s, t), y(s + ∆s, t)) − f (x(s, t), y(s, t))
∂f ∂f
= (x, y)(x(s + ∆s, t) − x(s, t)) + (x, y)(y(s + ∆s, t) − y(s, t)) + E
∂x ∂y
Por lo tanto
∆z ∂f (x(s + ∆s, t) − x(s, t)) ∂f (y(s + ∆s, t) − y(s, t)) E
= (x, y) + (x, y) +
∆s ∂x ∆s ∂y ∆s ∆s
∂f (x(s + ∆s, t) − x(s, t)) ∂f (y(s + ∆s, t) − y(s, t))

= (x, y) + (x, y)
∂x ∆s ∂y ∆s

(x(s + ∆s, t) − x(s, t)) (y(s + ∆s, t) − y(s, t))
+ , · hǫ1 , ǫ2 i
∆s ∆s
∂f ∂x ∂f ∂y
que tiende a (x, y) + (x, y) , cuando ∆s tiende a cero. Ası́,
∂x ∂s ∂y ∂s
∂z ∂f ∂x ∂f ∂y
= + .
∂s ∂x ∂s ∂y ∂s
De igual manera se puede verificar que:
∂z ∂f ∂x ∂f ∂y
= + .
∂t ∂x ∂t ∂y ∂t
Por último: uno de los usos más frecuentes de la regla de la cadena se dá en el cálculo
de derivadas de funciones definidas implı́citamente. Por ejemplo, si f es una función
diferenciable de dos variables y la ecuación 0 = f (x, y) define implı́citamente a y
como función diferenciable de x, y = g(x), podemos usar la regla de la cadena para
calcular la derivada de g. Sea F (t) = f (t, g(t)). Se tiene entonces que F (t) = 0 y por
lo tanto que su derivada es cero:
∂f dx ∂f dy
0 = F ′ (t) = +
∂x dt ∂y dt
∂f ∂f ′
= + g (t)
∂x ∂y
Por lo tanto, se tiene que
∂f
(t, g(t))
g ′ (t) = − ∂x
∂f
(t, g(t))
∂y
93
∂f
siempre y cuando (t, g(t)) 6= 0.
∂y
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 tem-

peratura es 0.
7π
b) Determine la razón de cambio de la temperatura del objeto en t = 3 .
c) Determine el punto o los puntos en la trayectoria del objeto en los cuales

la razón de cambio de la temperatura es máxima y los puntos en los cuales
es mı́nima. ¿Cuál es el valor de la temperatura en esos puntos?
2. En el instante t = 0, un objeto sale despedido de la superficie x2 +2y 2 +3z 2 = 6,

en el punto (1, 1, 1), en dirección normal a la superficie, con rapidez de 10
unidades por segundo.
a) ¿En qué instante de tiempo atraviesa el objeto cada uno de los planos
coordenados?
b) ¿En qué punto se encuentra el objeto en el instante t = 7?
c) ¿En qué punto atraviesa, el objeto, la superficie x2 + y 2 + z 2 = 103?
d ) Si a partir del momento en que el objeto sale despedido de la superficie,
se ejerce sobre él la fuerza F = 4mj − 10mk (donde m es la masa del
objeto), cuál es la altura máxima que alcanza el objeto sobre el plano
z = 0? y ¿cuál es el punto de impacto del objeto con el plano z = 0?
e) Si a partir del momento en que el objeto sale despedido de la superficie, se
ejerce sobre él la fuerza F = 4mj−10mk (donde m es la masa del objeto),
¿cuál es la altura máxima que alcanza el objeto sobre el plano tangente
a la superficie en el punto (1, 1, 1) ? y ¿cuál es el punto de impacto del
objeto con este plano?
3. El capitán Ralph tiene problemas cerca de la cara iluminada de Mercurio. La

temperatura en el casco de su nave cuando se encuentra en el punto (x, y, z)
2 2 2
es de T (x, y, z) = e−x −2y −3z , donde x, y, z se miden en metros. En este
momento está en el punto (1, −2, 1).
95
a) ¿En qué dirección debe moverse para que la temperatura baje lo más
rápidamente posible?
b) Si la nave se mueve con una rapidez de e6 metros por segundo, ¿a qué ve-
locidad bajará la temperatura cuando se desplace en esa dirección?
c) Desafortunadamente, el metal del casco se fracturar si la temperatura
√
cambia a con una rapidez mayor de 14e4 grados por segundo. Describir
el conjunto de direcciones según las cuales puede desplazarse para bajar
la temperatura a un ritmo inferior al lı́mite permitido.
2 2
4. Comprobar la regla de la cadena para la función f (x, y, z) = ln 1+2x +3z
1+y 2
y la
2
trayectoria r(t) = (1 − t)i + t j + cos tk.
5. En cada literal determine la ecuación 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) b) f (x, y) = xy, P (1, −1, −1)
c) f (x, y) = x2 y, P (−1, 1, 1) d ) f (x, y) = xey , P (1, 0, 1)

p
e) f (x, y) = x + 2y, P (2, 1, 4) f ) f (x, y) = x2 + y 2 , P (3, 4, 5)
6. En cada literal determine la ecuación del plano tangente a la superficie dada

en el punto P
√
x2 y2 z2 b) x2 + y 2 + z 2 = 9, P (0, 0, 3)
a) 4 + 9 + 16 = 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
1
a) f (x, y) = x2 + y 2 − 1 b) f (x, y) =
x2 + y2
p
c) f (x, y) = x2 + y 2 + 4 d ) f (x, y) = x2 ey
e) f (x, y) = ln(xy) f ) f (x, y) = 2x + 5y
g) f (x, y, z) = sin(xyz) h) f (x, y, z) = x2 eyz

p
i ) f (x, y, z) = x2 + y 2 + z 2 j ) f (x, y, z) = x2 + y 2 + z 2
8. Calcule
la
derivada direccional de f en el punto P en la dirección de v =
1
√ ,√1
2 2
.
1
a) f (x, y) = x2 + y 2 − 1, P (1, 1) b) f (x, y) = , P (1, 1)
x2 + y2
p
c) f (x, y) = x2 + y 2 + 4, P (1, 1) d ) f (x, y) = x2 ey , P (1, 1)
9.
Calcule la derivada
direccional de f en el punto P en la direccón de v =
1 1
√ ,√ ,√ 1
3 3 3
.
a) f (x, y, z) = sin(xyz), P (1, 1, 1) b) f (x, y, z) = x2 eyz , P (1, 1, 1)
10. ¿En qué dirección la función f (x, y) = xy 2 + x3 y crece más rápidamente desde
el punto (1, 2)? ¿En qué dirección decrece más rápidamente desde el mismo
punto?
11. La temperatura T en un sólido está dada por la función T (x, y, z) = e−x +

e−2y + e4z , donde x, y, z son las coordenadas espaciales relativas al centro del
sólido. ¿En qué dirección desde el punto (1, 1, 1) decrecerá más rápidamente
la temperatura?
12. Sean f y g funciones continuamente diferenciables, c una constante, y v un

vector unitario en R2 . Muestre que:
a) ∇(cf ) = c ∇f b) ∇(f + g) = ∇f + ∇g
g ∇f − f ∇g
c) ∇(f g) = f ∇g + g ∇f d ) ∇(f /g) = , g(x, y) 6= 0
g2
97
df df d(cf ) df
e) =− f) =c
d(−v) dv dv dv
d(f + g) df dg d(f g) dg df
g) = + h) =f +g
dv dv dv dv dv dv
p
13. La función r(x, y) = x2 + y 2 es la longitud del vector posición 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.
14. La Regla de la Cadena no se puede utilizar en los casos en que la función f no

es diferenciable. Veamos un ejemplo: sea f definida por

2
 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).

b) Si g(t) = hat, bti para constantes a y b, muestre que f ◦ g es diferenciable
y que (f ◦ g)′ (0) = ab/(a2 + b2 ).
c) Sin embargo, muestre que ∇f (0, 0) · g′ (0) = 0.
2 2
15. Para la función z(x, y) = f (u(x, y), v(x, y)), con f (u, v) = uu2 −v
+v
2 , u(x, y) =
e −x−y xy
y v(x, y) = e . Calcule cada una de las siguientes derivadas usando la
regla de la cadena.
∂2z
a) ∂z c) ∂x∂y
∂x
∂2z
b) ∂z d) ∂y 2
∂y
16. Para la función z(x, y) = f (u(x, y), v(x, y)). Determine cada una de las si-
guientes derivadas usando la regla de la cadena.
∂2z
a) ∂z c) ∂x∂y
∂x
∂2z
b) ∂z d) ∂y 2
∂y
17. Dada la función z(u, v) = f (x(u, v), y(u, v)), con f (x, y) = x2 + 2xy + 3y 3 ,
u ∂z ∂z
x(u, v) = y y(u, v) = ln(u − v). Calcule los valores de y cuando
v ∂u ∂v
a) u = 4 y v = 2.
b) x = 2 y y = 0.
x−y
18. Dada la función z(u, v) = f (x(u, v), y(u, v)), con f (x, y) = arctan( ), con
x+y
∂z ∂z
x(u, v) = uv y y(u, v) = eu/v . Calcule los valores de y cuando
∂u ∂v
a) u = 4 y v = 2.
b) x = 4 y y = e
19. En cada literal, la ecuación define implı́citamente a y como función de x.

Calcule dy/dx en el punto dado.
a) 2x3 − 3y 2 + xy = 0 en el punto (1, 1).

b) xy + y 2 − 3x − 3 = 0 en el punto (−1, 1).
c) yex + sin(xy) + x − 3 = 0 en el punto (3, 0).
20. Muestre que si f es una función diferenciable de tres variables y la ecuación

0 = f (x, y, z) define implı́citamente a z como función diferenciable de x y y,
z = g(x, y), use la regla de la cadena para calcular las derivadas parciales de
g y muestre que
∂f
∂g (s, t, g(s, t))
(s, t) = − ∂x
∂x ∂f
(s, t, g(s, t))
∂z
y que
∂f
(s, t, g(s, t))
∂g ∂y
(s, t) = −
∂y ∂f
(s, t, g(s, t))
∂z
∂f
siempre y cuando (s, t, g(s, t)) 6= 0.
∂z
21. Calcule ∂z/∂x y ∂z/∂y en los puntos dados, suponiendo que las ecuaciones
dadas definen implı́citamente a z como función de x y y.
a) y 3 − xz + yz + z 3 − 2 = 0 en el punto (1, 1, 1)
1 1 1
b) + + − 1 = 0 en el punto (2, 3, 6).
x y z
c) cos(x + y) + cos(y + z) + cos(x + z) + 3 = 0 en el punto (π/2, π/2, π/2).
99
22. Sea f una función diferenciable y w(r, θ) = f (r cos θ, r sin θ).
a) Calcule ∂w/∂r y ∂w/∂θ en términos de ∂f /∂x y ∂f /∂y.

b) Escriba ∂f /∂x y ∂f /∂y en términos de ∂w/∂r y ∂w/∂θ.
c) Use los items anteriores para demostrar que
2 2 2 2
∂f ∂f ∂w 1 ∂w
+ = + 2
∂x ∂y ∂r r ∂θ
LECCIÓN 7
Extremos de funciones
Uno de los propósitos de estudiar cálculo diferencial en una sola variable fue dispo-
ner de herramientas para encontrar y clasificar los puntos crı́ticos de una función.
Los conceptos que hemos estudiado hasta el momento permiten hacer lo mismo con
funciones de varias variables. Nuevamente en esta lección haremos este estudio en
funciones de dos variables que se puede extender a funciones de más variables.
Sea f : D −→ R, D ⊂ R2 .
1. f tiene un valor máximo local en (a, b) ∈ D, si existe una vecindad V de (a, b)

tal que f (x, y) ≤ f (a, b), para todo (x, y) ∈ D ∩ V .
2. f tiene un valor mı́nimo local en (a, b) ∈ D, si existe una vecindad V de (a, b)

tal que f (x, y) ≥ f (a, b), para todo (x, y) ∈ D ∩ V .
3. f tiene un valor máximo absoluto en (a, b) ∈ D, si f (x, y) ≤ f (a, b), para

todo (x, y) ∈ D.
4. f tiene un valor mı́nimo absoluto en (a, b) ∈ D, si f (x, y) ≥ f (a, b), para todo
(x, y) ∈ D.
101
102 LECCIÓN 7. EXTREMOS DE FUNCIONES
En algunos contextos se usan las palabras relativo en lugar de local, global en lu-
gar de absoluto y a un valor máximo o un valor mı́nimo se le denomina valor extremo.
No es difı́cil demostrar que si f tiene un valor máximo local en un punto P (a, b)

interior de D y f tiene derivada en el punto P (a, b) en la dirección de v, entonces
df
(a, b) = 0. En efecto, sea p : (−ǫ, ǫ) −→ D un paso por P (a, b) en la dirección de
dv
v, como f tiene un máximo local en P (a, b) entonces la función f ◦ p, que es de una
sola variable, tiene un máximo 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 dirección de v, para todo v y (a, b) = ∇f (a, b) · v, por lo tanto po-
dv
demos concluir que si f tiene un máximo 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 análisis completamente análogo al anterior se puede hacer para el caso en que f

tenga un mı́nimo local en (a, b).
Nótese que si f es una función de una sola variable que tiene un máximo local o
un mı́nimo local en un punto a y f es diferenciable en a, entonces f ′ (a) = 0. Por
lo tanto, la recta tangente a la gráfica en el punto (a, f (a)) es paralela al eje X y
tiene ecuación y = f (a). De igual manera, si f es una función de dos variables que
tiene un valor máximo local o un mı́nimo local en (a, b) y f es diferenciable en (a, b)
entonces ∇f (a, b) = h0, 0i. Por lo tanto, la ecuación del plano tangente a la gráfica
de f en el punto (a, b, f (a, b)) es z = f (a, b), que es la ecuación de un plano
paralelo al plano XY .
Un punto P (a, b) ∈ D es un punto crı́tico 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 función alcanza sus
valores máximos y mı́nimos relativos, se debe buscar entre los puntos crı́ticos. Una
vez identificados estos puntos, viene la tarea de su clasificación. Es decir, determinar
103
en cuáles de estos puntos f tiene un valor mı́nimo local, en cuáles f tiene un valor
máximo local, y en cuáles no tiene máximo ni mı́nimo local.
A continuación se hará la deducción del criterio de clasificación de puntos crı́ticos

para funciones de dos variables, herramienta muy útil en muchos casos en los que la
función f tiene segundas derivadas parciales continuas en el punto crı́tico (a, b).
Supóngase que (a, b) es un punto crı́tico de la función f en el que f es diferenciable

y tiene segundas derivadas parciales continuas. Sea p un paso por el punto (a, b) en
la dirección 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 además 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 crı́tico de g. Aplicando nuevamente la regla de la cadena se obtiene:
d
g ′′ (t) = {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 expresión y
completar el cuadrado, obteniendo

′′
{v1 fxx (p(t)) + v2 fxy (p(t))}2 + fyy (p(t))fxx (p(t)) − fxy
2 (p(t)) v 2
2
g (t) = ,
fxx (p(t))
y evaluando en t = 0, podemos escribir:

′′
{v1 fxx (a, b) + v2 fxy (a, b)}2 + fyy (a, b)fxx (a, b) − fxy
2 (a, b) v 2
2
g (0) = .
fxx (a, b)
Haciendo
δ(a, b) = fxx (a, b)fyy (a, b) − (fxy (a, b))2 ,
de la última expresión 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,
independientemente del vector v. Como g es una función de una variable, por el

criterio de la segunda derivada se tiene que en el primer caso g alcanza un valor
mı́nimo local en 0, y en el segundo caso g alcanza un valor máximo local en 0. Como
g(t) = f (p(t)), donde p(t) es un paso por (a, b), se puede concluir que en el primer
caso f alcanza un valor mı́nimo local en (a, b) en cualquier dirección v , y en el
segundo caso f alcanza un valor máximo local en (a, b) en cualquier dirección v . Si
δ < 0 se tiene que hay direcciones en las que f alcanza un mı́nimo local en (a, b), y
direcciones en las que f alcanza un máximo local en (a, b), en cuyo caso se dice que
f tiene un punto de silla o un punto de ensilladura en (a, b).
En conclusión se obtiene el siguiente criterio para clasificar los puntos crı́ticos, de la

función f en los que f es diferenciable y tiene segundas derivadas parciales continuas.
1. Si δ > 0 y fxx (a, b) > 0, entonces f toma un valor mı́nimo local en (a, b),
2. Si δ > 0 y fxx (a, b) < 0, entonces f toma un valor máximo local en (a, b),
3. Si δ < 0, entonces f no toma ni un valor máximo local, ni un valor mı́nimo

local en (a, b). En este caso (a, b) es un punto de ensilladura.
Hay resultados sobre valores extremos absolutos que se obtuvieron para funciones
de una sola variable definidas en intervalos cerrados y acotados y que se pueden
extender también a funciones de varias variables. Por ejemplo, si f es una función
continua, definida en un intervalo cerrado y acotado I, se tiene que existen puntos
x0 , x1 ∈ I en los que la función alcanza un valor mı́nimo absoluto f (x0 ) y un valor
máximo absoluto f (x1 ). Para funciones de dos variables se tiene que si f está de-
finida en un conjunto cerrado y acotado D y f es continua en D, entonces existen
puntos (x0 , y0 ), (x1 , y1 ) ∈ D en los que f alcanza un valor mı́nimo absoluto f (x0 , y0 )
y un valor máximo absoluto f (x1 , y1 ), respectivamente.
Las técnicas de demostración de estos resultados son un poco sofisticadas para el

nivel de estas lecciones por lo que no las incluiremos aquı́. En cambio, proporciona-
remos un protocolo para determinar los puntos en los que la función toma los valores
extremos absolutos:
Si f : D −→ R es una función continua en el conjunto cerrado y acotado D, para
determinar los puntos en los que f toma sus valores extremos absolutos.
105
1. Encuentre puntos destacados del dominio de la función (por ejemplo, si el

dominio es un polı́gono, considere sus vértices),
2. Encuentre los puntos crı́ticos de f que se encuentren en D,
3. Encuentre los puntos de la frontera de D en los que f tiene valores extremos

(hay varias formas de hacerlo, si se conoce una parametrización, o una ecuación
de la frontera, se puede restringir la función y buscar puntos crı́ticos de la
restricción, por ejemplo.),
4. 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 identificación de los puntos de la frontera en los que la función tiene

valores extremos no es tarea fácil. Sin embargo, contamos con la técnica de los
multiplicadores de Lagrange que en muchos casos es de gran ayuda, sobre todo
cuando la frontera tiene una ecuación de la forma g(x, y) = k, es decir cuando esta
frontera es una curva de nivel Nkg de una cierta función g.
En realidad, la técnica de los Multiplicadores de Lagrange se usa para resolver

problemas de optimización de una función en los que se cuenta con restricciones
representables mediante curvas o superficies de nivel de otras funciones. Más pre-
cisamente, para hallar los puntos de una superficie de nivel k de una función g en
los que una función f alcanza valores máximos (o mı́nimos), se usa la técnica de los
Multiplicadores de Lagrange.
Supongamos que la función f tiene un valor extremo en el punto P (x0 , y0 ) de la curva

de nivel k de la función g 1 . Si parametrizamos Nkg mediante una función vectorial
r : (−ǫ, ǫ) −→ Nkg de tal forma que r(0) = hx0 , y0 i, se tiene, por un lado, que en 0
la función z = f ◦ r tiene un valor extremo 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).
1
Suponemos que las funciones involucradas son diferenciables
Esto quiere decir que en los puntos (x0 , y0 ) de Nkg 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 Nkg , 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 ecuación
anterior.
La identificación de puntos crı́ticos en la frontera del dominio de una función es

un problema de optimización con restricciones. En otras palabras, si la frontera del
dominio de la función f es la curva de nivel de una función g, se puede restringir f
al conjunto Nkg = {(x, y); g(x, y) = k} y encontrar puntos crı́ticos de f en Nkg .
En el caso en que la función f de tres variables esté sometida a dos restricciones

dadas por superficies de nivel de otro par de funciones g1 y g2 , se tiene también una
técnica de Multiplicadores de Lagrange. Explı́citamente, supongamos que hemos
restringido la función f al conjunto de puntos
R = {(x, y, z) : g1 (x, y, z) = k1 y g2 (x, y, z) = k2 } ,
es decir, a la intersección de las superficies de nivel Nkg11 y Nkg22 . Supóngase además

que f tiene un valor máximo local en P (a, b, c) ∈ R (local quiere decir comparado
con puntos de R en una vecindad de P ) y que podemos parametrizar la parte de
R que contiene el punto P mediante una función vectorial r : (−ǫ, ǫ) −→ D, con
r(0) = ha, b, ci. Entonces, la función f ◦ r tiene un valor máximo relativo en 0 y por
consiguiente tiene derivada nula. Tenemos entonces, usando la regla de la cadena,
que
1. 0 = ∇f (a, b, c) · r′ (0)
2. 0 = ∇g1 (a, b, c) · r′ (0)
3. 0 = ∇g2 (a, b, c) · r′ (0)

107
Es decir, el gradiente de f en (a, b, c) está en el plano definido por los gradientes de

g1 y g2 en (a, b, c). Este hecho se puede escribir ası́:
∇f (a, b, c) = λ∇g1 (a, b, c) + µ∇g2 (a, b, c)
La conclusión es que, si estamos buscando puntos donde la función tiene valores

extremos sujeta a las dos restricciones en cuestión, debemos buscar entre aquellos
puntos en los que se satisface la ecuación anterior para valores convenientes de λ y
µ, conocidos como multiplicadores de Lagrange.
EJERCICIOS
1. Para cada una de las siguientes funciones determine los puntos crı́ticos y cla-
sifı́quelos
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
e) f (x, y) = 2x3 + 6xy + 3y 2 f ) f (x, y) = 2x3 − 6xy + y 2

p
g) f (x, y) = x2 + y 2 h) f (x, y) = −4x2 + 4xy − 2y 2 + 16x −
12y
i ) f (x, y) = x + 2y j ) f (x, y) = 4x2 − 4xy + 2y 2 + 10x − 6y
3x4 − 4x3 − 12x2 + 18 l ) f (x, y) = (2 + cos πx) sen πy

k ) f (x, y) =
12(1 + 4y 2 )
m) f (x, y) = y sen πx
2. Halle el punto del plano 2x−y +3z = 20 más próximo al origen, de tres formas
diferentes.
3. Determine los valores máximo y mı́nimo absolutos de la función f (x, y) =

cos x + sin y en el rectángulo 0 ≤ x ≤ 2π, −2π ≤ y ≤ 2π.
4. Determine los valores máximo y mı́nimo absolutos de la función f (x, y) =

x2 + y 2 − 3y − xy en el disco x2 + y 2 ≤ 9
5. Sea A una matriz simétrica 3 × 3 diferente de 0, considerando la función

f (x) = 21 (Ax) · x.
a) Determine ∇f .

b) Determine el máximo y el mı́nimo de f en B = (x, y, z) : x2 + y 2 + z 2 = 1
6. Para un sólido rectangular de 1000 metros cúbicos de volumen, encuentre las

dimensiones que maximizan el área de su superficie.
109
7. Demuestre que si (a, b) es un punto donde f alcanza un valor máximo local o

un valor mı́nimo local y f es diferenciable en (a, b), entonces el plano tangente
a la superficie z = f (x, y) en el punto (a, b, f (a, b)) es paralelo al plano XY .
8. Encuentre tres números positivos x, y, z cuya suma es 10 y tales que x2 y 2 z

sea máxima.
9. Encuentre los valores máximo y mı́nimo de f (x, y) = 2x+y dado que x2 +y 2 =

4.
10. Encuentre los valores máximo y mı́nimo de f (x, y) = xy dado que x2 +3y 2 = 6.
11. Encuentre los puntos de la circunferencia x2 +y 2 = 100 que están más cercanos
y más lejanos del punto (2, 3).
12. Encuentre los valores máximo y mı́nimo de f (x, y, z) = x + y 2 + 2z dado que

4x2 + 9y 2 − 36z 2 = 36.
13. Encuentre el volumen del paralelepı́pedo rectangular más grande que puede
inscribirse en el elipsoide.
x2 y 2 z 2
+ 2 + 2 =1.
a2 b c
14. Determine las dimensiones del cono de mayor área superficial que puede ins-
cribirse en una esfera de radio a.
15. La intersección entre la superficie cilı́ndrica de ecuación x2 + y 2 = 4 y el plano

de ecuación 2x + 3y + 6z = 24 es una elipse. Determine los puntos, de esta
elipse, mas cercanos y los más lejanos del origen del sistema de coordenadas.
16. Un alambre de longitud A se corta en dos pedazos. Uno se dobla en forma

de cuadrado y el otro en forma de cı́rculo, encuentre las dimensiones de los
dos pedazos de manera que la suma de las áreas del cuadrado y el cı́rculo sea
máxima y las dimensiones de los dos pedazos de manera que la suma de las
áreas sea mı́nima.
17. Resuelva el ejercicio anterior en el caso en que los alambres se doblan en forma
de triángulo equilátero y cuadrado.
18. Resuelva el ejercicio anterior en el caso en que los alambres se doblan en forma
de triángulo rectángulo isósceles y cuadrado.
19. Una caja rectangular sin tapa tiene base cuadrada, encontrar las dimensiones
de la caja de volumen constante y de mı́nima área superficial.
20. La longitud de generatriz de un cono circular recto es a, determine la medida

del ángulo formado por la generatriz y el eje de simetrı́a del cono de volumen
máximo.
21. Muestre que el producto de los senos de los ángulos de un triángulo es máximo
cuando el triángulo es equilátero.
22. Determine los máximos y mı́nimos absolutos de la función f (x, y) = 2x2 −4x+
y 2 − 4y + 1 en la región triangular de vértices (0, 0), (0, 2) y (1, 2)
23. Determine la distancia más corta del punto (a, 0), a 6= 0, a la curva de ecuación
x = y2 .
24. Determine el punto de la curva de intersección del plano 2y + 4z = 5 y el cono

z 2 = 4x2 + 4y 2 más cercano al origen.
25. Descomponer el número 97 en 4 sumandos de manera que la suma de sus

cuadrados sea mı́nima.
26. Descomponer el número positivo p en n factores positivos, de manera que la

suma de sus inversos sea mı́nima.
27. Determine las dimensiones del tetraedro, en el primer octante, de volumen

mı́nimo limitado por los planos coordenados y un plano tangente al elipsoide
2
x2 2
a2
+ yb2 + zc2 = 1.
28. Determine las dimensiones del recipiente cónico con tapa semiesférica con ca-
pacidad de 1 litro y área superficial máxima, y del de área superficial mı́nima.
29. Hay una versión del algoritmo de Newton para resolver sistemas de ecuaciones
f1 (x, y) = 0 and f2 (x, y) = 0 ,

111
donde f1 (x, y) y f2 (x, y) don funciones de valor real:

Tome un punto inicial (x0 , y0 ). Para n = 0, 1, 2, 3, . . . , defina:

f1 (xn , yn ) f2 (xn , yn )

∂f1
(xn , yn ) ∂f2 (xn , yn )
∂y ∂y
xn+1 = xn − ,
D(xn , yn )

f1 (xn , yn ) f2 (xn , yn )

∂f1 ∂f2

(xn , yn ) (xn , yn )
∂x ∂x
yn+1 = yn + ,
D(xn , yn )
donde D(xn , yn ) = ∂f ∂f2 ∂f1 ∂f2

∂x (xn , yn ) ∂y (xn , yn )− ∂y (xn , yn ) ∂x (xn , yn ) . Entonces
1
la sucesión de puntos (xn , yn )∞

n=1 converge a una solución. Escriba un programa
de computador que use este algoritmo para encontrar soluciones aproximadas
del sistema de ecuaciones:
sin(xy) − x − y = 0 y e2x − 2x + 3y = 0 .
Muestre que se obtienen soluciones diferentes cuando se usan (0, 0) y (1, 1)

como puntos iniciales (x0 , y0 ).
LECCIÓN 8
Integrales múltiples
En cursos de cálculo anteriores se estudió el concepto de integral de una función f

de una sola variable definida en un intervalo cerrado [a, b]. Los conceptos relaciona-
dos son los de partición de un intervalo, suma de Riemann y lı́mite de una suma
de Riemann, conceptos que vamos a extender a funciones de varias variables. Dos
resultados importantes que se usan frecuentemente son:
Una función continua definida en un intervalo cerrado es integrable, es decir

que el lı́mite de las sumas de Riemman de la función cuando la norma de la
partición tiende a cero existe.
El Teorema Fundamental del Cálculo: El valor de la integral de una función

continua sobre un intervalo cerrado es igual a la diferencia entre los valores de
una primitiva en los extremos del intervalo. (Este resultado permite darle la
vuelta al problema de calcular sumas de Riemann de la función, cuando sea
posible encontrar una primitiva.)
Sea f : [a, b] −→ R una función continua. Una partición de la 1-celda [a, b] es un

conjunto finito de puntos P = {x0 , x1 , x2 , . . . , xn } tales que x0 = a < x1 < x2 <
. . . < xn = b. Podemos suponer que estos puntos son equidistantes, es decir que
113
114 LECCIÓN 8. INTEGRALES MÚLTIPLES
xk+1 − xk = (b − a)/n = ∆x, en este caso se dice que la partición es uniforme y

depende del número entero positivo n. Una suma de Riemann R(f, n, ∗) está definida
por una partición y por la escogencia de un punto de muestra x∗k en cada uno de las
subceldas [xk , xk+1 ] de [a, b] definidos por la partición, ası́:
R(f, n, ∗) = f (x∗1 )(x1 − x0 ) + f (x∗2 )(x2 − x1 ) + · · · + f (x∗n )(xn − xn−1 )
= f (x∗1 )∆x + f (x∗2 )∆x + f (x∗3 )∆x + · · · + f (x∗n )∆x
= (f (x∗1 ) + f (x∗2 ) + f (x∗3 ) + · · · + f (x∗n ))∆x
n
X
= f (x∗k )∆x.
k=1
R Rb
La integral de f sobre la 1-celda I = [a, b], que se escribe I f (x)dx o a f (x)dx, es
el lı́mite de estas sumas de Riemann cuando n tiende a infinito,
Z b n
X
f (x)dx = lı́m R(f, n, ∗) = lı́m f (x∗k )∆x.
a n→∞ n→∞
k=1
Este lı́mite existe para funciones continuas independientemente de los puntos de

muestra escogidos.
Ahora extenderemos el concepto de integral a funciones de más de una variable de

manera natural.
Sea R la 2-celda [a, b] × [c, d]. Una partición uniforme de R, se obtiene escogiendo
particiones uniformes x0 = a < x1 < x2 < . . . < xn = b y y0 = c < y1 < y2 <
. . . < ym = d de las 1-celdas [a, b] y [c, d], respectivamente, tales que xk+1 − xk =
(b−a)/n = ∆x y que yj+1 −yj = (c−d)/m = ∆y. Ası́, la 2-celda R queda dividido en
mn 2-celdas de la forma [xk , xk+1 ]×[yj , yj+1 ], de área ∆A = (xk+1 −xk )(yj+1 −yj ) =
((b − a)/n)((c − d)/m) = ∆x∆y. Escogemos un punto de muestra (x∗kj , ykj ∗ ) en cada
una de las mn 2-celdas [xk , xk+1 ] × [yj , yj+1 ], y definimos la suma de Riemann de
una función continua f : R ⊂ R2 −→ R, correspondiente a esta partición y esta
escogencia de punto de muestra ası́:
115
R(f, n, m, ∗)
= f (x∗11 , y11
∗ )∆A + f (x∗ , y ∗ )∆A + · · · + f (x∗ , y ∗ )∆A
21 21 n1 n1
∗ )∆A + f (x∗ , y ∗ )∆A + · · · + f (x∗ , y ∗ )∆A

+ f (x∗12 , y12 22 22 n2 n2
...
+ f (x∗1m , y1m
∗ )∆A + f (x∗ , y ∗ )∆A + · · · + f (x∗ , y ∗ )∆A
2m 2m nm nm
n X
X m
= f (x∗kj , ykj
∗
)∆A
k=1 j=1
R
La integral de f sobre R que se escribe R f (x, y)dA es el lı́mite de estas sumas de
Riemann cuando n y m tienden a infinito:
Z
f (x, y)dA = lı́mn,m→∞ R(f, n, m, ∗)
R
n X
X m
= lı́mn,m→∞ f (x∗kj , ykj
∗
)∆A.
k=1 j=1
Este lı́mite existe para funciones continuas e independientemente de los puntos de

muestra escogidos.
Ahora, sea P la 3-celda [a, b] × [c, d] × [e, h]. Una partición uniforme de P se obtiene
escogiendo particiones uniformes
x0 = a < x1 < x2 < . . . < xn = b, y0 = c < y1 < y2 < . . . < ym = d y z0 = e < z1 <
z2 < . . . < zl = d de los intervalos [a, b], [c, d] y [e, h], respectivamente. Ası́, la 3-celda
P queda dividida en mnl 3-celdas de la forma [xk , xk+1 ] × [yj , yj+1 ] × [zi , zi+1 ], con
el mismo volumen ∆V = ∆x∆y∆z. Escogemos un punto de muestra (x∗kji , ykji ∗ , z∗ )
kji
en cada 3-celda [xk , xk+1 ] × [yj , yj+1 ] × [zi , zi+1 ] y definimos la suma de Riemann de
una función continua f : P −→ R, correspondiente a esta partición y esta escogencia
de punto de muestra, ası́:
X l
m X
n X
R(f, n, m, l, ∗) = ∗
f (x∗kji , ykji ∗
zkji )∆V.
k=1 j=1 i=1
R
La integral de f sobre P que se escribe P f (x, y, z)dV es el lı́mite de estas sumas
de Riemann cuando n, m y l tienden a infinito:
Z
f (x, y, z)dV =
P
= lı́mn,m,l→∞ R(f, n, m, l, ∗)
X l
m X
n X
= lı́mn,m,l→∞ f (x∗kji , ykji
∗ ∗
, zkji )∆V
k=1 j=1 i=1
Este lı́mite existe para funciones continuas e independientemente de los puntos de

muestra escogidos.
Evidentemente, el cálculo de integrales de funciones sobre 2-celdas o 3-celdas usan-

do sumas de Riemann, es una tarea muy dispendiosa. Sin embargo, las sumas de
Riemann son un instrumento de aproximación muy útil de estas integrales. Veremos
ahora un método para calcularlas sin necesidad de recurrir a las sumas de Riemann.
Este consiste en escribir la integral como una integral iterada, es decir una secuencia
de integrales parciales unidimensionales que se calcularán recurriendo al Teorema
Fundamental del Cálculo. Definamos primero lo que es la integral iterada de una fun-
ción de dos variables (el lector puede extender esta idea a funciones de tres variables).
Sea f : R −→ R, R = [a, b] × [c, d], una función continua. Para cada x ∈ [a, b]
Rd
definimos A(x) = c f (x, y)dy. Resulta ser que la función A : [a, b] −→ R es continua,
y por consiguiente se puede calcular su integral:
Z b Z b Z d
A(x)dx = f (x, y)dy dx.
a a c
Ésta se conoce como integral iterada. El Teorema de Fubini establece que la integral
de f sobre R es igual a ésta integral iterada. Explı́citamente, si f es continua en R
se tiene que
R R b R d
R f (x, y)dA = a c f (x, y)dy dx
R d R b
= c a f (x, y)dx dy
El Teorema de Fubini es una generalización del Principio de Cavalieri para calcular

volúmenes de sólidos (¿qué dice el principio de Cavalieri?). El Teorema de Fubini
117
para una función continua f definida en una 3-celda P = [a, b] × [c, d] × [e, h] dice
que:
R R b R d R h
P f (x, y, z)dV = a c e f (x, y, z)dz dy dx
R b R h R d
= a e c f (x, y, z)dy dz dx
R h R b R d
= e a c f (x, y, z)dy dx dz
= ...
los puntos suspensivos son para las otras tres integrales iteradas que faltan (¿cúales?).
Para dar cuenta de integrales de funciones definidas en regiones más generales que
las celdas, recurriremos a una estrategia de parametrización, concepto que no es muy
utilizado en los libros de cálculo pero que ’‘mata varios pájaros de un solo tiro”(ver
Lección 3).
Sean f : σ −→ R una función continua, definida en una superficie σ del plano y

r : [a, b] × [c, d] −→ R2 , r(u, v) = hx(u, v), y(u, v)i, una parametrización de σ. Su-
pondremos además que x = x(u, v) y y = y(u, v) tienen derivadas parciales continuas
en (a, b) × (c, d).
Una partición de la superficie σ en subregiones σij se obtiene a partir de una partición

de [a, b] × [c, d] en 2-celdas Rij = [ui−1 , ui ] × [vj−1 , vj ]:
σij = {(x(u, v), y(u, v)) : (u, v) ∈ [ui−1 , ui ] × [vj−1 , vj ]} .
Escogemos en cada subregión σij el punto Pij (x(ui−1 , vj−1 ), y(ui−1 , vj−1 )). Una suma
de Riemann de f sobre σ con respecto a la parametrización r y a una partición de
σ en mn subregiones, está definida por
P m Pn
R(f, m, n, r) = i=1 j=1 f (Pij ) k∆u ru × ∆v rv k ,
P m Pn
= i=1 j=1 f (Pij ) kru × rv k ∆u∆v,
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 aparición de k∆u ru × ∆v rv k en
la suma de Riemann se debe a que ésta es el área del paralelogramo generado por los
vectores ∆u ru (Pij ) y ∆v rv (Pij ) , que es una aproximación del área de la subregión
σij .
La integral de f sobre σ es el lı́mite de R(f, m, n, r) cuando m y n tienden a infinito:
ZZ m X
X n
f dA = lı́mm,n→∞ f (Pi,j ) kru × rv k ∆u∆v
σ i=1 j=1
ZZ
= f (x(u, v), y(u, v)) kru × rv k dudv
[ab]×[cd]
Obsérvese que si f : [a, b] × [c, d] −→ R es continua y r : [a, b] × [c, d] −→ R2

está definida por r(u, v) = hu, vi, entonces kru × rv k = 1, y por consiguiente
ZZ ZZ
f dA = f (u, v)dudv.
[a,b]×[c,d] [a,b]×[c,d]
Es decir, las integrales sobre 2-celdas son un caso particular de integrales sobre su-
perficies planas más generales.

Calculemos la integral de f (x, y) = x2 + y 2 sobre el disco σ = (x, y) : x2 + y 2 ≤ 4
con respecto a la parametrización de σ que dimos en la lección 3. Tenemos que
kru × rv k = khcos v, sen vi × h−u sen v, u cos vik = u,
y por lo tanto,
ZZ ZZ Z 2 Z 2π
2 3
f dA = u ududv = u dv du = 8π
D [0,2]×[0,2π] 0 0
Consideremos ahora una función continua f : Λ −→ R, donde Λ es una sóli-

do del espacio junto con una parametrización r : [a, b] × [c, d] × [e, h] −→ R3 ,
r(u, v, w) = hx(u, v, w), y(u, v, w), z(u, v, w)i. Supondremos que las funciones x =
x(u, v, w), y = y(u, v, w), z = z(u, v, w) tiene derivadas parciales continuas en
(a, b) × (c, d) × (e, h).
Una partición del sólido Λ en subregiones Λijk se obtiene a partir de una partición
de [a, b] × [c, d] × [e, h] en 3-celdas Rijk = [ui−1 , ui ] × [vj−1 , vj ] × [wk−1 , wk ]:
Λijk = {(x(u, v, w), y(u, v, w), z(u, v, w)) : (u, v, w) ∈ [ui−1 , ui ] × [vj−1 , vj ] × [wk−1 , wk ]} .
119
Escogemos en cada subregión Λijk el punto
Pijk (x(ui−1 , vj−1 , wk−1 ), y(ui−1 , vj−1 , wk−1 ), z(ui−1 , vj−1 , wk−1 )).
Una suma de Riemann de f sobre Λ con respecto a la parametrización r y a una

partición de Λ en mnl subregiones, está definida por
Pm P n Pl
R(f, m, n, l, r) = i=1 j=1 k=1 f (Pijk ) |(∆u ru × ∆v rv ) · ∆w rw |
Pm P n Pl
= i=1 j=1 k=1 f (Pijk ) |(ru × rv ) · rw | ∆u∆v∆w,
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
aparición de |(∆u ru × ∆v rv ) · ∆w rw | en la suma de Riemann se debe a que éste
es el volumen del paralelepı́pedo generado por ∆u ru , ∆v rv y ∆w rw , que es una
aproximación del volumen de la subregión Λijk .
La integral de f sobre Λ con respecto a r es el lı́mite de R(f, m, n, l, r) cuando m,

n y l tienden a infinito:
ZZZ
f dV =
Λ
Pm P n Pl
= lı́mm,n,l→∞ i=1 j=1 k=1 f (Pijk ) |(ru × rv ) · rw | ∆u∆v∆w
ZZZ
= f (x(u, v, w), y(u, v, w), z(u, v, w)) |(ru × rv ) · rw | dudvdw
[ab]×[cd]×[e,h]
Obsérvese que si f : [a, b] × [c, d] × [e, h] −→ R es continua y r : [a, b] × [c, d] ×

[e, h] −→ R3 está definida por r(u, v, w) = hu, v, wi, entonces |(ru × rv ) · rw | = 1, y
por consiguiente
ZZZ ZZZ
f dr = f (u, v, w)dudvdw.
[a,b]×[c,d]×[e,h] [a,b]×[c,d]×[e,h]
Es decir, las integrales sobre 3-celdas son un caso particular de integrales sobre sóli-
dos más generales.
Calculemos la integral de la función f (x, y, z) = x2 + y 2 + z 2 sobre la esfera Λ =

(x, y, z) : x2 + y 2 + z 2 ≤ 4 con respecto a la parametrización de Λ que dimos en
la Lección 3. Tenemos que
|(ru × rv ) · rw | =
= |(hcos v sen w, sen v sen w, cos wi × h−u 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,
y por lo tanto,
ZZZ
f dV =
Λ ZZZ
= u2 u2 sen wdudvdw
[0,2]×[0,2π]×[o,π]
Z 2 Z 2π Z π
4
= u sen wdw dv du
0 0 0
= 128π/5
Si el dominio de f es una superficie σ en el espacio podemos proceder exactamente
de la misma manera como lo hicimos en el caso de superficies planas. Sea f : σ −→ R

una función continua, definida en una superficie acotada σ del espacio y
r : [a, b] × [c, d] −→ R3 , r(u, v) = hx(u, v), y(u, v), z(u, v)i ,
una parametrización continua de σ 1 . Supondremos además que r tiene derivadas

parciales continuas en (a, b) × (c, d). La integral de f sobre σ con respecto a r
está dada por
ZZ
f dA =
σ ZZ
= f (x(u, v), y(u, v)) kru × rv k dudv.
[a,b]×[c,d]
1
Debemos suponer tambı́en que la restricción de r al rectángulo abierto (a, b)×(c, d) sea inyectiva
121
donde ru y rv son las derivadas parciales de r con respecto a u y a v.
Calculemos la integral de f (x, y, z) = x2 + y 2 + z 2 sobre el la superficie de la esfera

σ = (x, y, z) : x2 + y 2 + z 2 = 4 . Una parametrización de σ 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, π]} .
Tenemos que
kru × rv k =
= kh−2 sen u sen v, 2 cos u sen v, 0i × h2 cos u cos v, 2 sen u cos v, −2 sen vik
= 8 sen v,
y por lo tanto,
ZZ ZZ Z 2π Z π
f dA = 32 sen vdudv = 32 sen vdv du = 16π.
σ [0,π]×[0,2π] 0 0
Nos queda estudiar el caso de una función definida en una curva. Sea f : Γ −→ R
una función continua, definida en una curva acotada Γ del espacio (o del plano).
Una parametrización de Γ es una función vectorial
r : [a, b] −→ R3 , r(u) = hx(u), y(u), z(u)i ,
continua e inyectiva y tal que
Γ = {(x(u), y(u), z(u)) : u ∈ [a, b]} .
Supondremos además que r es derivable en (a, b).
Una partición de Γ en subarcos Γi se obtiene a partir de una partición de [a, b] en

subintervalos Ii = [ui−1 , ui ]:
Γi = {(x(u), y(u), z(u)) : u ∈ [ui−1 , ui ]} .
Escogemos en cada subarco Γi el punto Pi (x(ui−1 ), y(ui−1 ), z(ui−1 )). Una suma de
Riemann de f sobre Γ con respecto a la parametrización r y a una partición de Γ
en m subarcos, está definida por

m
X
R(f, m, r) = f (Pi ) ∆u r′ (ui−1 ) ,
i=1
m
!
X ′
= f (Pi ) r (ui−1 ) ∆u,
i=1
donde r′ (ui−1 ) es la derivada r en Pi . La aparición de k∆u r′ (ui−1 )k en la suma de

Riemann se debe a que ésta es la longitud del vector tangente ∆u r′ (ui−1 ) que es
una aproximación de la longitud del subarco Γi .
La integral de f sobre Γ con respecto a r es el lı́mite de R(f, m, r) cuando m tiende
a infinito: Z !
Xm
′
f dL = lı́mm→∞ f (Pi ) r (ui−1 ) ∆u
Γ i=1
Z

= f (x(u), y(u), z(u)) r′ (u) du
[a,b]
Calculemos la integral de f (x, y, z) = x2 + y 2 + z 2 sobre la circunferencia Γ =

(x, y, z) : x2 + y 2 = 4, z = 5 . Una parametrización de Γ está dada por la función
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π]} .
Tenemos que
′
r = kh−2 sen u, 2 cos u, 0ik = 2,
y por lo tanto,
Z Z Z 2π
f dL = 29 · 2 du = 58 du = 116π.
Γ [0,2π] 0
123
EJERCICIOS
1. Calcule la suma de Riemman R (f, n, m, ∗) para la función f (x, y) = x2 y + y,

la región 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
a) x2 y 2 z 2 dzdydx
0 2 1
Z 2Z 3
b) (x2 y + xy 3 )dydx
1 0
Z 1Z 4
x √3
c) ( + x2 )dydx
0 2 y
Z 1 Z π/4
d) (r cos θ − tan θ)dθdr
0 0
Z bZ dp
e) s2 + t2 dsdt
a c
3. Parametrizar la región del plano limitada por las curvas dadas, de manera que
el dominio de la parametrización 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 cónica

p
de ecuación z = x2 + y 2 , parametrice cada una de las regiones del espacio
de manera que el dominio de los parámetros sea una 1-celda (intervalo), una
2-celda o una 3-celda, dependiendo de la naturaleza de la región.
a) La superficie S.
b) El sólido limitado por S.

c) La parte de la superficie S que está dentro de la superficie cónica P .
d ) La parte de la superficie cónica que se encuentra dentro de la esfera.
e) La curva de intersección de las dos superficies.
f ) El sólido interior a la esfera y exterior al cono
5. Dado que S es la parte del plano de ecuación 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 parametrización sea una 1-celda, una 2-celda o
una 3-celda, dependiendo de la naturaleza de la región.
a) La superficie S.
b) El sólido 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 región.
R
a) C (x + y + z) dl donde C es el borde de la región 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 sólido descrito en el numeral 4b).
R
m) S (x + y + z)dV donde S es sólido descrito en el numeral 4f).
125
R
n) S (x + y + z)dV donde S es sólido descrito en el numeral 5b).
7. Teorema de Clairaut. Si f : D ⊆ R2 −→ R donde D es un disco que contiene

al punto (a, b), en el cual fxy y fyx son continuas entonces fxy (a, b) = fyx (a, b).
Teorema de Fubini. Sea f : R2 −→ R continua sobre la región R = [a, b]×[c, d].
RbRd RdRb RR
Entonces a c f (x, y)dydx = c a f (x, y)dxdy = R f (x, y)dA
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 = .
∂x∂y ∂y∂x
Use este ejemplo para establecer la relación que existe entre el teorema
de Clairaut y el de Fubini.
R 1 R 1 x2 − y 2
b) Determine el valor de cada una de las siguientes integrales 0 0 2 dxdy
(x + y 2 )
R 1 R 1 x2 − y 2
y 0 0 2 dydx ¿Contradice ésto el teorema de Fubini?
(x + y 2 )
8. Explique por qué el volumen de la región limitada por la superficie z = f (x, y),
el plano z = 0 cuando (x, y) ∈ R, siendo R una 2-celda, está dado por la
R
integral doble R f (x, y)dA.
9. Encuentre el volumen de la región limitada por la superficie z = f (x, y), el

plano z = 0 cuando (x, y) ∈ R:
a) f (x, y) = 4xy, R = [0, 1] × [0, 1]

b) f (x, y) = ex+y , R = [0, 1] × [−1, 1]
c) f (x, y) = x3 + y 2 , R = [0, 1] × [0, 1]
d ) f (x, y) = x4 + xy + y 3 , R = [1, 2] × [0, 2]
10. Evalúe la integral doble dada:

Z 1Z 2 Z 1Z 2
2
a) (1 − y)x dx dy b) x(x + y) dx dy
0 1 0 0
Z 2Z 1 Z 2 Z 1
c) (x + 2) dx dy d) x(xy + sin x) dx dy
0 0 −1 −1
Z π/2 Z 1 Z π Z π/2
e) xy cos(x2 y) dx dy f) sin x cos(y − π) dx dy
0 0 0 0
Z 2Z 4 Z 1 Z 2
g) xy dx dy h) 1 dx dy
0 1 −1 −1
RdRb
11. Si M es una constante muestre que c a M dx dy = M (d − c)(b − a).
12. Evalúe la integral doble dada:

Z 1Z 1
Z π Z y
2 b) sin x dx dy
a) √
24x y dy dx
0 x 0 0
Z 2 Z ln x Z 2 Z 2y
2
c) 4x dy dx d) ey dx dy
1 0 0 0
Z Z ∞Z
π/2 Z y ∞
2 +y 2 )
e) cos x sin y dx dy f) xye−(x dx dy
0 0 0 0
Z 2Z y Z 1 Z x2
g) 1 dx dy h) 2 dy dx
0 0 0 0
13. Encuentre el volumen V del sólido acotado por los tres planos coordenados y
el plano x + y + z = 1.
14. Encuentre el volumen V del sólido acotado por los tres planos coordenados y
el plano 3x + 2y + 5z = 6.
RR
15. Explique por qué la integral doble lı́mR 1 dA da el área de la región R.
16. Demuestre que el volumen de un tetraedro con lados adyacentes mutuamente

perpendiculares de longitudes a, b, y c es abc
6 .
17. Evalúe la integral triple:

Z 3Z 2Z 1 Z 1Z xZ y
a) xyz dx dy dz b) xyz dz dy dx
0 0 0 0 0 0
Z π Z x Z xy Z 1Z z Z y
2 2
c) x sin z dz dy dx d) zey dx dy dz
0 0 0 0 0 0
127
Z eZ y Z 1/y Z 2 Z y2 Z z2
e) x2 z dx dz dy f) yz dx dz dy
1 0 0 1 0 0
Z 2Z 4Z 3 Z 1 Z 1−x Z 1−x−y
g) 1 dx dy dz h) 1 dz dy dx
1 2 0 0 0 0
R z 2 R y2 R x 2
18. Si M es una constante. muestre que z1 y1 x1 M dx dy dz = M (z2 − z1 )(y2 −
y1 )(x2 − x1 ).
19. Encuentre el volumen V del sólido S acotado por los planos coordenados, por
encima por el plano x + y + z = 2, y por debajo por el plano z = x + y.
Z bZ zZ y Z b
(b−x)2
20. Muestre que f (x) dx dy dz = 2 f (x) dx.
a a a a
21. Escriba un programa que use el método de Monte Carlo para aproximar la
RR xy
integral doble e dA,donde R = [0, 1] × [0, 1]. Muestre el resultado del
R
programa para N = 10, 100, 1000, 10000, 100000 y 1000000 puntos aleatorios.
22. Escriba un programa que use el método de Monte Carlo para aproximar la
RRR xyz
integral triple e dV , donde S = [0, 1] × [0, 1] × [0, 1]. Muestre el resul-
S
tado del programa para N = 10, 100, 1000, 10000, 100000 y 1000000 puntos
aleatorios.
23. Encuentre el volumen V dentro del paraboloide z = x2 + y 2 para 0 ≤ z ≤ 4.

p
24. Encuentre el volumen V dentro del cono z = x2 + y 2 para 0 ≤ z ≤ 3.
25. Encuentre el volumen V del sólido dentro de x2 + y 2 + z 2 = 4 y x2 + y 2 = 1.
26. Encuentre el volumen V dentro de la esfera x2 + y 2 + z 2 = 1 y el cono z =

p
x2 + y 2 .
RR
27. Calcule sin x+y
2 cos x−y
2 dA, donde R es el triángulo con vértices (0, 0),
R
(2, 0) y (1, 1).
28. Encuentre el volumen del sólido acotado por z = x2 + y 2 y z 2 = 4(x2 + y 2 ).

x2 y2
29. Encuentre el volumen del sólido acotado por el cilindro elı́ptico a2
+ b2
=1
para 0 ≤ z ≤ 2.
2
x2 2
30. Muestre que el volumen del sólido delimitado por el elipsoide a2
+ yb2 + zc2 = 1
es 4πabc
3 .
31. Muestre que la función Beta definida por

Z 1
B(x, y) = tx−1 (1 − t)y−1 dt , for x > 0, y > 0,
0
satisface la relación B(y, x) = B(x, y) para x > 0, y > 0.
32. Use la sustitución t = u/(u + 1), para mostrar que la función Beta se puede
escribir como
Z ∞
ux−1
B(x, y) = du , for x > 0, y > 0.
0 (u + 1)x+y
33. Encuentre el centro de masa de la región R con función 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
p
d ) R = {(x, y) : y ≥ 0, x ≥ 0, 1 ≤ x2 + y 2 ≤ 4 }, δ(x, y) = x2 + y 2
e) R = {(x, y) : y ≥ 0, x2 + y 2 ≤ 1 }, δ(x, y) = y
34. Encuentre el centro de masa del sólido S con función de densidad dada δ(x, y, z).
a) S = {(x, y, z) : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1 }, δ(x, y, z) = xyz

b) S = {(x, y, z) : z ≥ 0, x2 + y 2 + z 2 ≤ a2 }, δ(x, y, z) = x2 + y 2 + z 2
c) S = {(x, y, z) : x ≥ 0, y ≥ 0, z ≥ 0, x2 + y 2 + z 2 ≤ a2 }, δ(x, y, z) = 1
d ) S = {(x, y, z) : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1 }, δ(x, y, z) = x2 +y 2 +z 2
e) S = {(x, y, z) : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1 − x − y}, δ(x, y, z) = 1
R∞ 2
35. Calcule la integral −∞ e−x dx.
R∞ 2 2
36. Para σ > 0 y µ > 0, evalue −∞ σ√12π e−(x−µ) /2σ dx.
LECCIÓN 9
Integrales múltiples
En este capı́tulo deduciremos una forma de calcular integrales múltiples de funciones

reales sobre algunos tipos especiales de dominios.
Estudiaremos dos tipos de dominios del plano denominados regiones del tipo I y re-
giones del tipo II. Las regiones del tipo I son aquellas superficies del plano limitadas
por rectas x = a y x = b y por las gráficas de dos funciones y = g1 (x) y y = g2 (x)
diferenciables en [a, b]. Las regiones del tipo II son aquellas superficies del plano
limitadas por las rectas y = c y y = d y por las gráficas de dos funciones x = h1 (y)
y x = h2 (y) diferenciables en [c, d].
Las regiones del tipo I son superficies planas de la forma:
σ = {(x, y) : a ≤ x ≤ b, g1 (x) ≤ y ≤ g2 (x)} .
129
y = g2 (x)
y = g1 (x)
X
a b
De acuerdo con lo estudiado en la Lección 3, una parametrización para esta superficie

σ es r : [a, b] × [0, 1] −→ R2 dada por
r(u, v) = h u , (1 − v)g1 (u) + vg2 (u) i .
Si f : σ ⊆ R2 −→ R es una función continua, definida en la superficie σ entonces

ZZ ZZ
f dA = f (x(u, v), y(u, v)) kru × rv k dudv.
σ [a,b]×[c,d]
Se tiene ası́ que

ru = h 1 , (1 − v)g1′ (u) + vg2′ (u) i ,
rv = h 0 , g2 (u) − g1 (u) i ,
y por consiguiente que
kru × rv k = g2 (u) − g1 (u).
Ası́,
ZZ ZZ
f dA = f (u, (1 − v)g1 (u) + vg2 (u)) (g2 (u) − g1 (u)) dudv
σ [a,b]×[0,1]
Z b Z 1
= f (u, (1 − v)g1 (u) + vg2 (u)) (g2 (u) − g1 (u))dv du.
a 0
131
Haciendo el cambio de variable y = (1−v)g1 (u)+vg2 (u) y tomando x = u se obtiene

ZZ Z Z !
b g2 (x)
f dA = f (x, y)dy dx.
σ a g1 (x)
Las regiones del plano del tipo II son superficies del plano de la forma:
σ = {(x, y) : h1 (y) ≤ x ≤ h2 (y), c ≤ y ≤ d} ,
x = h1 (y) x = h2 (y)
Una parametrización de esta superficie es r : [0, 1] × [c, d] −→ R2 dada por
r(u, v) = h (1 − u)h1 (v) + uh2 (v) , v i .
Si f : σ ⊆ R2 −→ R es una función continua definida en σ, procediendo de manera

análoga al caso anterior, se tiene que
ZZ Z Z !
d h2 (y)
f dA = f (x, y)dx dy
σ c h1 (y)
Para dominios del espacio, estudiaremos aquı́ regiones de del tipo I, del tipo II y
del tipo III. Las regiones de tipo I son sólidos determinados por una superficie σ del
plano XY y por las gráficas de dos funciones z = g1 (x, y) y z = g2 (x, y) diferenciables

en σ. Son de la forma
Λ = {(x, y, z) : (x, y) ∈ σ, g1 (x, y) ≤ z ≤ g2 (x, y)} .
Las regiones del tipo II son sólidos determinados por una superficie σ del plano XZ
y por las gráficas de dos funciones y = h1 (x, z) y y = h2 (x, z) diferenciables en la
superficie σ. Son de la forma
Λ = {(x, y, z) : (x, z) ∈ σ, h1 (x, z) ≤ y ≤ h2 (x, z)} .
Por último, las regiones del tipo III son sólidos determinados por una superficie σ del
plano Y Z y por las gráficas de dos funciones x = l1 (y, z) y x = l2 (y, z) diferenciables
en la superficie σ. Son de la forma
Λ = {(x, y, z) : (y, z) ∈ σ, l1 (y, z) ≤ x ≤ l2 (y, z)} .
Si ahora consideramos una función continua f : Λ −→ R donde Λ es uno de los

sólidos anteriores en el que σ es una región de uno de los tipos estudiados y para-
metrizamos de manera conveniente el sólido Λ podemos deducir que:
ZZZ ZZ Z !
g2 (x,y)
f dV = f (x, y, z)dz dA,
Λ σ g1 (x,y)
si Λ es del tipo I,
ZZZ ZZ Z !
h2 (x,z)
f dV = f (x, y, z)dy dA,
Λ σ h1 (x,z)
si Λ es del tipo II, y

ZZZ ZZ Z !
l2 (y,z)
f dV = f (x, y, z)dx dA,
Λ σ l1 (y,z)
si Λ es del tipo III. En cada caso, la integral doble sobre σ se calculará de acuerdo
a la naturaleza de σ.
133
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 I
4
3
2 II
III
1
−1
−1 1 2 3 4 5 6 7 8 9
a) Determine el centro de masa de una placa delgada con la forma de la

región I.
b) Determine el centro de masa de una placa delgada con la forma de la
región II.
c) Determine el centro de masa de una placa delgada con la forma de la
región III
d ) Determine el centro de masa de una placa delgada con la forma de la
región I unida con la región II.
e) Determine el centro de masa de una placa delgada con la forma de la
región I unida con la región II unida con la región III.
f ) Dé una forma general para determinar el centro de masa de una placa
delgada que es la unión no traslapada de n placas delgadas, en función
del centro de masa de cada una de las placas que se han unido.
2. Calcule cada una de las siguientes integrales.

Z 1Z 1
a) (x + y)2 dxdy
−1 |y|
Z 4Z 2
2
b) ex dxdy
y
0 2
Z 1 Z √9−y2
c) √ xdxdy
−3 − 9−y 2
Z 1Z π
4 x2 − y 2
d) dxdy
0 arctan y (x2 + y 2 )
3. Toda función continua f definida en una 2-celda R = [a, b] × [c, d] alcanza

su valor máximo M y su valor mı́nimo m en puntos de R. Además, si ϕ es
otra función continua y no negativa definida en R, se tiene que mϕ(x, y) ≤
G(x, y)ϕ(x, y) ≤ M ϕ(x, y) para todo (x, y) ∈ R. Por consiguiente,
ZZ ZZ ZZ
m ϕ(x, y) dA ≤ G(x, y)ϕ(x, y) dA ≤ M ϕ(x, y) dA,
R R R
desigualdad que se conoce como Desigualdad del Valor medio. Si ϕ no se anula

RR
en R, R ϕ(x, y) dA > 0 y ası́, dividiendo a ambos lados de esta desigualdad
RR
por R ϕ(x, y) dA, se tiene que
RR
G(x, y)ϕ(x, y) dA
RRR
m≤ ≤ M.
R ϕ(x, y) dA
Por el Teorema del Valor intermedio, existe (a, b) ∈ R tal que
RR
G(x, y)ϕ(x, y) dA
RRR
G(a, b) =
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

ZZ
1 1
≤ 2 f (x, y)dA ≤ e.
e 4π R
dA
b) Si f (x, y) = y R = [−1, 1] × [−1, 2] entonces
x2 + y 2 + 1
ZZ
dA
1≤ 2 2
≤6
R x +y +1
135
RRR R1RxRy
4. Si Λ f (x, y, z)dV = 0 0 0 f (x, y, z)dzdydx, dibuje el sólido Λ y plantee
las otras 5 integrales triples iteradas correspondientes.
RRR
5. Calcular Λ f (x, y, z)dV si
p 2 2 2
a) f (x, y, z) = x2 + y 2 + z 2 e−(x +y +z ) y Λ es el sólido limitado por las
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
x2 y 2 z 2
6. Sea Λ el elipsoide sólido + 2 + 2 ≤ 1 donde a, b, c son positivos.
a2 b c
a) Determine el volumen de Λ.
2
RRR x y2 z2
b) Calcule Λ a 2 + b2 + c 2 dV
7. Sea Ω la región en el plano XY interior a la gráfica de x2 + y 2 = 1 y exterior

a la gráfica de x2 + y 2 = 2y con x, y ≥ 0.
a) Dibujar la región Ω.
b) Si u = x2 + y 2 y v = x2 + y 2 − 2y, dibujar la región Ω∗ que corresponde
a Ω bajo esta transformación.
RR
c) Calcular Ω xey dA usando la transformación.
8. Determine la masa del sólido acotado por las gráficas de x2 + y 2 = 2x y

p
z 2 = x2 + y 2 si la densidad es δ(x, y) = x2 + y 2 .
9. Determine el valor medio (o promedio) de la función f (x, y, z) = sin2 πz cos2 πx

sobre la 3-celda [0, 2] × [0, 4] × [0, 6].
10. Determine el valor medio (o promedio) de la función f (x, y, z) = e−z sobre la

esfera de centro en el origen y radio a.
RARBRC
11. Calcular a b c f (x, y, z)dzdydx si f (x, y, z) = Fxyz (x, y, z) y a, b, c, A, B, C
son constantes.
RRR m n p
12. Calcular Λ x y z dV si Λ es la esfera sólida con centro en el origen y radio
1 y m, n, p son números 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
números positivos si Λ es el sólido 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́ razón 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 sólida con centro en el origen
y radio 1, determinando dos sólidos. Encuentre el valor de α para el cual la
razón entre los volúmenes de los dos sólidos es 3.
16. Dada la esfera de ecuación x2 + y 2 + z 2 = a2 y el cono de ecuación z =

p
k x2 + y 2 , determine el valor de la constante k para el cual la razón entre las
áreas 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 parábola en el plano XY de ecuación y = x2 desde el punto

(−1, 1) hasta el punto (2, 4).
b) La parte de la élice de ecuación r(t) = cos ti + sen tj + tk desde el punto
(1, 0, 0) hasta el punto (1, 0, 2π).
x2 y2
c) La elı́pse en el plano XY de ecuación a2
+ b2
= 1.
18. Determine el área 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 ecuación z = x2 + y 2 .
b) La parte de la superficie cilı́ndrica de ecuación z = 4 − y 2 cuyos puntos
satisfacen la desigualdad x2 + y 2 ≤ 4.
c) La parte del paraboloide de ecuación z = x2 + y 2 cuyos puntos satisfacen
las desigualdades |x| ≤ 1 y |y| ≤ 1
x2 y2
d ) La parte del plano XY acotada por la elı́pse de ecuación a2
+ b2
=1
137
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 sólido acotado por el cono de ecuación z = x2 + y 2 y por el plano
z = 1 si la densidad en cada punto está dada por δ = z 2 .(Thomas)
c) Una lámina 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 lámina delgada modelada por la superficie del cono z = x2 + y 2
cortada por el cilı́ndro circular x2 + y 2 = 2x. (Thomas)
x2 y2 z2
e) Un sólido acotado por el el elipsoide a2
+ b2
+ c2
= 1 si la densidad en
cada punto está dada por δ = z
LECCIÓN 10
Campos vectoriales
Los campos vectoriales son funciones que se usan frecuentemente para modelar
fenómenos de asignación vectorial, como por ejemplo el campo de velocidades de
un fluido, el campo eléctrico, el campo magnético y el campo gravitacional. Desde
el punto de vista matemático un campo vectorial es una función definida en algún
conjunto y que toma valores en un espacio vectorial. Nos ocuparemos aquı́ del estu-
dio de campos vectoriales bidimensionales y tridimensionales.
Un campo vectorial bidimensional es una función F : D −→ R2 donde D es un

subconjunto de R2 :
F(x, y) = hP (x, y), Q(x, y)i ; (x, y) ∈ D.
Un campo vectorial tridimensional es una función F : D −→ R3 donde D es un

subconjunto de R3 :
F(x, y, z) = hP (x, y, z), Q(x, y, z), R(x, y, z)i ; (x, y, z) ∈ D.
Las funciones P = P (x, y) y Q = Q(x, y) en el caso bidimensional o las funciones

P = P (x, y, z), Q = Q(x, y, z) y R = R(x, y, z) en el caso tridimensional, son
funciones escalares y se denominan las funciones componentes del campo. El estudio
139
140 LECCIÓN 10. CAMPOS VECTORIALES
de los campos vectoriales se hace de manera natural mediante el estudio de sus

funciones componentes:
El dominio D de un campo vectorial es la intersección de los dominios de sus funcio-
nes componentes. El rango es el conjunto de las imágenes de los puntos del dominio.
La gráfica es el conjunto de todas las parejas formadas por cada elemento del domi-
nio y su respectiva imagen.
El lı́mite de un campo vectorial en un punto de su dominio existe si existen los lı́mi-

tes de sus funciones componentes y las derivadas parciales de un campo vectorial
existen si existen las derivadas parciales de sus funciones componentes, los cálculos
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 gráficamente un campo vectorial bidimensional es

dibujando a partir de unos cuantos puntos (x, y) ∈ D los vectores F(x, y) corres-
pondientes. Para ilustrar la representación gráfica 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
X
141
La circulación y el flujo.
Los dos conceptos más importantes, desde el punto de vista fı́sico, relacionados con
campos vectoriales son el flujo y la circulación, y que se hacen evidentes al estudiar
la interacción de los campos con curvas y superficies. Estudiaremos primero estos
dos conceptos en campos vectoriales bidimensionales para los cuales nos interesa
estudiar su interacción con curvas contenidas en sus dominios.
Consideremos un campo vectorial bidimensional F : D −→ R2 definido en un sub-

conjunto D de R2 y una curva Γ en D parametrizada por una función vectorial
r : [a, b] −→ R2 . Definimos la densidad de circulación, δ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)).
F T
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
en ese punto. Obsérvese que las densidades de circulación δC y de flujo δF de un
campo son funciones escalares con dominio Γ, la curva parametrizada por la función
vectorial r. Por ser N el vector que se obtiene al rotar T 90◦ en sentido positivo, se
tiene que, como
* +
x′ (t) y ′ (t)
T= p ,p
x′ (t)2 + y ′ (t)2 x′ (t)2 + y ′ (t)2
entonces * +
y ′ (t) x′ (t)
N= −p ,p .
x′ (t)2 + y ′ (t)2 x′ (t)2 + y ′ (t)2
Ası́ como al integrar la densidad de masa se obtiene la masa, en este caso al integrar la
densidad de circulación del campo sobre la curva se obtiene la circulación del campo
a lo largo de la curva. Ası́ mismo, al integrar la densidad de flujo del campo sobre la
curva se obtiene el flujo del campo a traves de la curva. Es decir, la circulación de F
R
a lo largo de Γ es Γ δC dL, y ésta se puede escribir en términos de las componentes
escalares de F ası́:
Z Z b
δC dL = F(r(t)) · T(t) r′ (t) dt
Γ a
Z b
r′ (t)
r′ (t) dt
= F(r(t)) · ′
a kr (t)k
Z b
= F(r(t)) · r′ (t)dt
a
Z b
= (P (r(t))x′ (t) + Q(r(t))y ′ (t))dt.
a
Debido a la forma de la última integral, la circulación de F a lo largo de Γ se suele

R R
escribir como Γ P dx + Qdy o como Γ F · dL.
143
R
El flujo de F a través de Γ es Γ δf dL, y se puede escribir en términos de las com-
ponentes escalares de F ası́:
Z Z b
δF dL = F(r(t)) · N(t) r′ (t) dt
Γ a
Z * +
b
y ′ (t) x′ (t) ′
r (t) dt
= F(r(t)) · −p ,p
a x′ (t)2 + y ′ (t)2 x′ (t)2 + y ′ (t)2
Z b

= hP (r(t)), Q(r(t))i · −y ′ (t), x′ (t) dt
a
Z b
= (Q(r(t))x′ (t) − P (r(t))y ′ (t))dt.
a
Debido a la forma de la última integral, el flujo de F a través de Γ se suele escribir

R R
como Γ Qdx − P dy o como Γ F e · dL, donde F e = hQ, −P i se conoce como el campo
dual del campo F. Lo que se muestra en la cadena de igualdades arriba, es que el
flujo de un campo a través de una curva es igual a la circulación de su campo dual
a lo largo de la misma curva.
Densidad rotacional y densidad de expansión.
Los conceptos de densidad rotacional y de densidad de expansión de un campo vec-

torial 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
circulación y el flujo por unidad de área, 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 además que las componentes escalares del campo, P y Q, tienen
derivadas parciales continuas en una vecindad de (0, 0) y consideremos la circulación
1
Estas densidades se conocen tambı́en como microcirculación y microflujo del campo, respecti-
vamente.
de F a lo largo del borde Γ del rectángulo Ω de vértices V1 (∆x, ∆y), V2 (−∆x, ∆y),
V3 (−∆x, −∆y) y V4 (∆x, −∆y), orientado positivamente2 , donde V1 , V2 , V3 y V4
están en la vecindad . La circulación de F a lo largo del lado V1 V2 , que escribiremos
CV1 V2 F es, aproximadamente, el producto de la componente tangencial de F en
el punto medio de V1 V2 y de la longitud de V1 V2 , es decir CV1 V2 F ≈ {F(0, ∆y) ·
(−i)}{2∆x}. La circulación del campo a lo largo de los demás lados se puede expresar
de la misma forma, y ası́, se tiene que la circulación de F a lo largo de Γ, CΓ F es:
CΓ F = CV 1 V 2 F + CV 2 V 3 F + CV 3 V 4 F + CV 4 V 1 F
≈ {F(0, ∆y) · (−i)}{2∆x} + {F(−∆x, 0) · (−j)}{2∆y}
+{F(0, −∆y) · i}{2∆x} + {F(∆x, 0) · j}{2∆y}
= 2∆y{Q(∆x, 0) − Q(−∆x, 0)} − 2∆x{P (0, ∆y) − P (0, −∆y)}
= 4∆x∆yQx (α, 0) − 4∆x∆yPy (0, β)
= 4∆x∆y{Qx (α, 0) − Py (0, β)},

donde α ∈ (−∆x, ∆x) y donde β ∈ (−∆y, ∆y). En el penúltimo paso se hizo uso del
Teorema del Valor medio para derivadas. Si usamos la continuidad de las derivadas
parciales de P y Q obtenemos que la densidad rotacional de F en (0, 0), que es la
circulación por unidad de área de F en (0,0) y se denota por δR (0, 0), es:
Circulación de F a lo largo de Γ
δR (0, 0) = lı́m(∆x,∆y)→(0,0)
área de Ω
4∆x∆y{Qx (α, 0) − Py (0, β)}

= lı́m(∆x,∆y)→(0,0)
4∆x∆y
= Qx (0, 0) − Py (0, 0)
Calculemos ahora la densidad de expansión de F en el punto (0, 0), con las mismas
suposiciones hechas anteriormente. Para ésto, consideremos el flujo de F a través
del borde Γ del rectángulo Ω de vértices V1 (∆x, ∆y), V2 (−∆x, ∆y), V3 (−∆x, −∆y)
2
es decir, en el sentido contrario al del movimiento de las manecillas del reloj
145
y V4 (∆x, −∆y), orientado negativamente. El flujo de F a través del lado V1 V2 , que

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}{2∆x}. El flujo del campo a través de los demás lados se puede expresar de la
misma forma, y ası́, se tiene que el flujo de F a través 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}{2∆x} + {F(−∆x, 0) · (−i)}{2∆y}
+{F(0, −∆y) · (−j)}{2∆x} + {F(∆x, 0) · i}{2∆y}
= 2∆y{P (∆x, 0) − P (−∆x, 0)} + 2∆x{Q(0, ∆y) − Q(0, −∆y)}
= 4∆x∆yPx (α, 0) + 4∆x∆yQy (0, β)
= 4∆x∆y{Px (α, 0) + Qy (0, β)},
donde α ∈ (−∆x, ∆x) y donde β ∈ (−∆y, ∆y). En el penúltimo 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 expansión de F en (0, 0), que es el
flujo por unidad de área de F en (0, 0) y denotaremos por δE (0, 0), es:
Flujo de F a través de Γ
δE (0, 0) = lı́m(∆x,∆y)→(0,0)
área de Ω
4∆x∆y{Px (α, 0)) + Qy (0, β)}

= lı́m(∆x,∆y)→(0,0)
4∆x∆y
= Px (0, 0) + Qy (0, 0)
Las densidades rotacional y de expansión se pueden calcular en todos los puntos
del dominio del campo F de la misma forma en que se calculó en el punto (0, 0),
para encontrar que la densidad rotacional de F en (x, y), δR (x, y) y la densidad de
espansión de F en (x, y), δE (x, y), están dadas por:
δR (x, y) = Qx (x, y) − Py (x, y), y δE (x, y) = Px (x, y) + Qy (x, y),

respectivamente.
Debemos señalar que el concepto de circulación a lo largo de una curva se puede

extender de forma natural a campos vectoriales tridimensionales, pero el concepto
de flujo a través de una curva no (¿por qué?). En la siguiente lección discutiremos los
conceptos, flujo, densidad de expansión y densidad de rotación de campos vectoriales
tridimensionales.
Los campos gradiente
Un ejemplo de campo vectorial es el que proviene del gradiente de una función. El

argumento que presentaremos a continuación es válido para funciones de dos o de
tres variables, y demuestra que la circulación de un campo gradiente de una función
f depende exclusivamente de los valores de f en los puntos inicial y final de la
curva. Considere una función f de varias variables con dominio D y Γ una curva en
D parametrizada por una función vectorial r : [a, b] −→ D continua. Calculemos la
circulación del campo vectorial ∇f a lo largo de Γ:
Z Z b
∇f · dL = ∇f (r(t)) · r′ (t)dt
Γ a
Z b
d(f ◦ r)
= (t)dt
a dt
= f (r(b)) − f (r(a))
Este resultado significa que si F es un campo gradiente, la circulación de F a lo

largo de cualquier curva dentro del dominio sólamente depende de los puntos inicial
y final de la curva, recı́procamente, si la circulación de un campo F a lo largo de
cualquier curva dentro del dominio sólamente depende de los puntos inicial y final
de la curva, entonces el campo vectorial es un campo gradiente. En efecto, debido
R
a la suposición hecha, está bien definida la función 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 función f tal que F = ∇f , se llaman
147
campos conservativos y, en este caso, a la función f se le llama función potencial de f .
Para campos vectoriales tridimensionales se cuenta con criterios para decidir si un

campo vectorial es conservativo. En realidad contamos con la equivalencia de las si-
guientes condiciones, siempre y cuando F sea un campo vectorial cuyas componentes
tengan primeras derivadas parciales continuas, excepto posiblemente en un número
finito de puntos,
1. F es el gradiente de una función f definida en los puntos en los que las com-
ponentes de F son diferenciables.
R
2. Γ F · dL = 0, para cualquier curva cerrada simple Γ.
R R
3. Γ1 F · dL = Γ2 F · dL, para cualquer dos curvas simples Γ1 y Γ2 con el mismo
punto inicial y mismo punto final.
4. rotF = 0
La demostración de la equivalencia de estas condiciones se puede encontrar en el

libro de Cálculo Vectorial de Marsden y Tromba, Addison Wesley, 2004, quinta
edición. Sin embargo, debemos tener en cuenta que estas condiciones no son todas
equivalentes si consideramos campos vectoriales bidimensionales. El contraejemplo
−y x
lo dá el campo vectorial F = h 2 , 2 i que está definido en R2 − {(0, 0)} y
x + y x + y2
2
sus componentes tienen derivadas parciales de todos los órdenes en todos los puntos
de su dominio (es decir, todos los de R2 excepto un punto). Diga cuáles de las
condiciones anteriores no son equivalentes.
Para terminar esta sección 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.
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 expansión de F en (0, 0, 0). Ayuda: imite el cálculo de la densidad
de expansión de un campo bidimensional considerando una 3-celda conveniente
con centro en (0, 0, 0).
R
2. Muestre que si Γ F · dL depende únicamente de los puntos inicial y final de
Γ, entonces F = ∇f para alguna función 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 · dL = 0.
R R
b) Si F(r(t)) es paralelo a r′ (t), muestre que Γ F · dL = ± Γ kFk dL.
4. Si la curva Γ está parametrizada por r(t) = cos3 ti + sin3 (t)j, con 0 ≤ t ≤ 2π,
y F = y 1/3 i + x1/3 j calcule la circulación d e F a lo largo de Γ y el flujo de F
a través de Γ.
5. Calcule la circulación 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 vértices (±1, ±1, ±1).
2 2 2
6. Suponga que ∇f (x, y, z) = 2xyzex i + zex j + yex k. Si f (0, 0, 0) = 5, hallar
f (1, 1, 2)
R
7. Calcule Γ f dL para la función f y la curva Γ dadas:
a) f (x, y) = xy; Γ: x = cos t, y = sin t, 0 ≤ t ≤ π/2

x
b) f (x, y) = 2 ; Γ : x = t, y = 0, 0 ≤ t ≤ 1
x +1
c) f (x, y) = 2x + y; Γ: camino poligonal desde (0, 0) hasta (3, 0) hasta
(3, 2)
d ) f (x, y) = x + y 2 ; Γ: camino desde (2, 0) en el sentido antihorario a lo
largo de la circunferencia x2 + y 2 = 4 hasta el punto (−2, 0) y de regreso
al punto (2, 0) a lo largo del eje X.
149
8. Para cada uno de los siguientes campos vectoriales haga una representación
R
gráfica del campo y la curva, y calcule Γ F · dL.
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; Γ : x = cos t, y = sin t, 0 ≤ t ≤ 2π
e) F(x, y) = xy 2 i + xy 3 j; Γ : el camino poligonal desde (0, 0) hasta (1, 0)
hasta (0, 1) hasta (0, 0)
f ) F(x, y) = (x2 + y 2 ) i; Γ : x = 2 + cos t, y = sin t, 0 ≤ t ≤ 2π
R R
9. Demuestre que Γ f (x, y) dL = −Γ 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 · dL| ≤ M L.
Z
11. Calcule (x2 + y 2 ) dx + 2xy dy para Γ : x = cos t, y = sin t, 0 ≤ t ≤ 2π.
Γ
Z
12. Calcule (x2 + y 2 ) dx + 2xy dy para Γ : x = cos t, y = sin t, 0 ≤ t ≤ π.
Γ
13. En cada uno de los numerales determine si el campo F tiene una función
potencial y si la tiene determı́nela
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
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) · dL = a F · dL ± b G · dL .
Γ Γ Γ
R
15. Sea Γ una curva cuya longitud de arco es λ. Muestre que Γ 1 dL = λ.
16. Sean f y g funciones de valor real continuamente diferenciables en una región

Ω. Muestre que Z Z
f ∇g · dL = − g ∇f · dL
Γ Γ
para cualquier curva cerrada Γ en Ω.
17. Suponga que F(x, y) = x2−y +y 2

x
i + x2 +y 2 j para todo (x, y) 6= (0, 0) es un campo
de fuerzas, y sea Γ : x = cos t, y = sin t, 0 ≤ t ≤ 2π.
a) Determine el trabajo realizado por el campo al desplazar un objeto sobre

la curva Γ desde el punto (1, 0) hasta el puno (0, 1) en sentido positivo
(esto es sentido antihorario).
b) Determine el trabajo realizado por el campo al desplazar un objeto sobre
la curva Γ desde el punto (1, 0) hasta el puno (0, 1) en sentido negativo.
c) Determine el trabajo realizado por el campo al desplazar un objeto sobre
la curva Γ una vuelta completa desde algún punto de la misma.
d ) Determine el trabajo realizado por el campo al desplazar un objeto sobre
la curva Γ una vuelta completa y media desde algún punto.
e) Sea f (x, y) = tan−1 (y/x). Calcule el gradiente de f y discuta la relación
que éste tiene con el campo vectorial F.
18. Sean g y h funciones diferenciables de una sola variable, y sea F(x, y) = h(y) i+
g(x) j. Si F tiene un potencial f , encuéntrelo.
19. Dado el campo vectorial F(x, y) = (2x sen y + x3 )i + (x2 cos y + tan y)j, dé tres
curvas que tengan como punto inicial el punto A(2, 0) y como punto final el
√ √
punto B(− 2, − 2), y calcule la circulación del campo sobre cada una de las
curvas. ¿Qué puede afirmar acerca de de la circulación de este campo sobre
cualquier curva que tenga como punto inicial A y como punto final B?
20. Dado el campo vectorial F(x, y) = P (x, y)i + Q(x, y)j ¿Qué condiciones deben
satisfacer las funciones P y Q para que el campo sea un campo gradiente y
cómo se puede saber si si lo es?
LECCIÓN 11
Divergencia y rotacional
En esta sección presentaremos la generalización de los conceptos de densidad de ex-

pansión y densidad rotacional a campos vectoriales tridimensionales y enunciaremos
los Teoremas de la Divergencia y del Rotacional.
Densidad de expansión para campos tridimensionales
Ası́ como la densidad de expansión para un campo vectorial bidimensional es el flujo

del campo por unidad de área, la densidad de expansión de un campo tridimensional
es el flujo del campo por unidad de volumen. Para hacer esta generalización debemos
precisar el concepto de flujo de un campo a través 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 tridimen-
sional y σ una superficie contenida en el dominio del campo. Supongamos que σ es
una superficie orientable, esto es, se puede parametrizar por una función 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 través de la superficie σ es la integral de δf
151
152 LECCIÓN 11. DIVERGENCIA Y ROTACIONAL
sobre la superficie:
R RR
σ δf dA = σ F · NdA
RbRd ru × rv
= a c F(r(u, v)) · kru × rv k dvdu
kru × rv k
RbRd
= a c F(r(u, v)) · (ru × rv )dvdu
Con este concepto de flujo de un campo a través de una superficie podemos ge-
neralizar el concepto de densidad de expansión 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 pequeños para que la 3-celda esté contenida en el domi-
nio del campo. Procediendo de manera similar a como se hizo en la lección anterior,
podemos hacer una aproximación del flujo del campo F a través de la superficie de
la 3-celda, dividir esta cantidad por el volumen de la 3-celda y hacer tender las lon-
gitudes de las aristas a cero. Se Obtiene ası́ la densidad de expansión δ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 cálculos 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 expansión 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)
Densidad rotacional para campos tridimensionales

153
La circulación de un campo tiene sentido sólamente a lo largo de curvas, debido a que

para cada punto en una curva se puede determinar un único vector tangente unita-
rio, 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 circulación por unidad de área. 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 área determinaremos
la densidad de circulación referida a un plano que contenga al punto y haremos los
cálculos 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 ecuación ax + by − z = 0, que se puede
parametrizar con la función r(x, y) = hx, y, ax + byi, tomemos el cuadrilátero de
vértices V1 (∆x, ∆y, a∆x+b∆y), V2 (−∆x, ∆y, −a∆x+b∆y), V3 (−∆x, −∆y, −a∆x−
b∆y) y V4 (∆x, −∆y, a∆x − b∆y). Sea Γ la curva cerrada formada por los segmentos
V1 V2 , V2 V3 , V3 V4 y V4 V1 , recorrida en el sentido descrito y Ω la región del plano
limitada por Γ. Calculemos ahora una aproximación de la circulación del campo F
a lo largo Γ, de manera idéntica a como se hizo en el caso bidimensional. Si CΓ F es
la circulación de F a lo largo de Γ:
√ 1
CΓ F ≈ 2∆y 1 + b2 F(∆x, 0, a∆x) · √ h0, 1, bi
1 + b2
√ 1
+2∆x 1 + a2 F(0, ∆y, b∆y) · √ h−1, 0, −ai
1 + a2
√ 1
+2∆y 1 + b2 F(−∆x, 0, −a − ∆x) · √ h0, −1, −bi
1 + b2
√ 1
+2∆x 1 + a2 F(0, −∆y, −b∆) · √ h1, 0, ai
1 + a2
= 2∆y(Q(∆x, 0, a∆x)) − Q(−∆x, 0, −a∆x) + bR(∆x, 0, a∆x)) − bR(−∆x, 0, −a∆x))
= 4∆x∆y(Qx (α2 , 0, a∆x) + aQz (−∆x, 0, β2 ) + bRx (α3 , 0, a∆x) + abRz (−∆x, 0, β3 ))
−4∆x∆y(Py (0, α1 , b∆y) + bPz (0, −∆y, β1 ) + aRy (0, α4 , b∆y) + abRz (0, −∆y, β4 )),
donde α2 , β2 , α3 , β3 ∈ (−∆x, ∆x) y donde α1 , β1 , α4 , β4 ∈ (−∆y, ∆y). En el último paso se

hizo uso del Teorema del Valor medio para derivadas. Como el área del paralelogramo es
√
A(Ω) = 4∆x∆y 1 + a2 + b2 , se tiene que la densidad rotacional δR (0, 0, 0) del campo F,
referida al plano dado, en el punto (0, 0, 0) es:
CΓ F
δR (0, 0, 0) = lı́m(∆x,∆y)→(0,0)
A(Ω)
= lı́m(∆x,∆y)→(0,0)

Qx (α2 , 0, a∆x) + aQz (−∆x, 0, β2 ) + bRx (α3 , 0, a∆x) + abRz (−∆x, 0, β3 )
√
1 + a 2 + b2

Py (0, α1 , b∆y) + bPz (0, −∆y, β1 ) + aRy (0, α4 , b∆y) + abPz (0, −∆y, β4 )
− √
1 + a 2 + b2
1
= √ (a(Qz − Ry ) + b(Rx − Pz ) + (Qx − Py ))
1 + a 2 + b2
1
= √ h−a, −b, 1i · hRy − Qz , −(Rx − Pz ), Qx − Py i |(0,0,0)
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 ecuación
ax + by − z=0, es:
1
δR (x, y, z) = √ h−a, −b, 1i · hRy − Qz , −(Rx − Pz ), Qx − Py i |(x,y,z)
1 + a 2 + b2
Obsérvese que la densidad rotacional en un punto depende de dos vectores, el vector
1
√ h−a, −b, 1i
1 + a 2 + b2
que es un vector unitario normal al plano z = ax + by y el vector
hRy (x, y, z) − Qz (x, y, z), −(Rx (x, y, z) − Pz (x, y, z)), Qx (x, y, z) − Py (x, y, z)i ,
que depende únicamente del campo. A este último vector se le denomina el rotacional del
campo F en el punto (x, y, z) y se nota como rotF(x, y, z), esto es
rotF(x, y, z) = hRy (x, y, z) − Qz (x, y, z), −(Rx (x, y, z) − Pz (x, y, z)), Qx (x, y, z) − Py (x, y, z)i .
De ésta última observación podemos concluir que la densidad rotacional en un punto (x, y, z)
de un campo vectorial F referida a un plano con vector normal unitario N es
δR (x, y, z) = rotF(x, y, z) · N.
155
Para recordar las componentes del vector rotF se usa la siguiente fórmula mnemotécnica

i j k

rotF = ∂x ∂y ∂z

P Q R
Podemos ahora calcular la densidad rotacional de un campo F referida a una superficie orien-
table σ, 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 rotaciónal
de F. Explı́citamente 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 nota-
ción. Primero que todo, recordemos que para definir el concepto de integral múltiple sobre
curvas, superficies y sólidos, introdujimos el concepto de parametrización mediante funcio-
nes vectoriales r de una, dos y tres variables respectivamente. Para poder diferenciar la
dimensión de la integral introducimos aquı́ los siguientes elementos de integración:
1. Elemento de longitud de arco:
dL = r′ (t)dt y dL = kr′ (t)k dt
2. Elemento de área de superficie:
dA = (ru × rv )dudv y dA = kru × rv k dudv
3. Elemento de volumen de sólido:
dV = |(ru · (rv × rw )| dudvdw
Con esta terminologı́a podemos escribir:
1. Si Γ es una curva parametrizada por la función vectorial derivable r = r(t), a ≤ t ≤ b

Z Z b
F · dL = F(r(t)) · r′ (t)dt
Γ a
y
Z Z b
f dL = f (r(t)) kr′ (t)k dt.
Γ a
2. Si σ es una superficie parametrizada por la función vectorial diferenciable r = r(u, v),

a ≤ u ≤ b, c ≤ v ≤ d
ZZ Z bZ d
F · dA = F(r(u, v)) · (ru × rv )dvdu
σ a c
y
ZZ Z b Z d
f dA = f (r(u, v)) kru × rv k dvdu.
σ a c
3. Si Λ es una sólido parametrizado por la función vectorial diferenciable r = r(u, v, w),

a ≤ u ≤ b, c ≤ v ≤ d, e ≤ w ≤ h
ZZZ Z bZ dZ h
f dV = f (r(u, v, w)) |ru · (rv × rw )| dwdvdu.
Λ a c e
Usando esta notación podemos escribir la rotación total de un campo F sobre la superficie
σ: RR RR
δ dA =
σ R σ
rotF · NdA
RbRd ru × rv
= a c
rotF(r(u, v)) · kru × rv k dvdu
kru × rv k
RbRd
= a c
rotF(r(u, v)) · ru × rv dvdu
RR
= σ
rotF · dA
La expansión total de un campo F sobre el sólido Λ:
RRR RRR
δ dV =
Λ E Λ
divFdV
RbRdRh
= a c e
divF(r(u, v, w)) |(ru × rv ) · rw | dwdvdu,
y la circulación de un campo F a lo largo de la curva Γ:

R R
δ dL = Γ F · TdL
Γ c
Rb r′ (t)
= a
F(r(t)) · kr′ (t)k dt
kr′ (t)k
Rb
= a
F(r(t)) · r′ (t)dt
R
= Γ
F · dL
Estamos ahora en la capacidad de enunciar los teoremas fundamentales: el de la Divergencia

y el de Stokes.
157
Teorema de la divergencia: Suponga que F= hP, Q, Ri es un campo vectorial y σ una

superficie cerrada que encierra un sólido Λ contenidos en el dominio del campo y que Λ
se puede parametrizar mediante una función vectorial diferenciable r : P −→ R3 . Si las
componentes escalares de F tienen derivadas parciales continuas, entonces, la expansión
total del campo en Λ es igual al flujo de F a través de la frontera σ de Λ. Más precisamente
se tiene que: ZZZ ZZ
divF dV = ± F · dA
Λ σ
El signo depende de la parametrización.
Teorema del rotacional(Stokes): Suponga que F= hP, Q, Ri es un campo vectorial y que

σ es una superficie con borde Γ contenidos en el dominio del campo y que σ se puede
parametrizar mediante una función vectorial diferenciable r : R −→ R3 . Si las componentes
escalares de F tienen derivadas parciales continuas, entonces, la rotación total del campo
sobre σ es igual a la circulación de F a lo largo del borde Γ de σ. Más precisamente se tiene
que: ZZ Z
rotF · dA = F · dL
σ Γ
La orientacón de Γ la dá la parametrización de σ.
Nota: En el caso en que Γ y σ estén contenidas en un plano y F sea un campo bidimensional,

el Teorema de Stokes se conoce también como el Teorema de Green y en ciertos contextos,
el Teorema de la Divergencia se conoce como el Teorema de Gauss.
EJERCICIOS
1. Use el Teorema del Rotacional1 para calcular las siguientes integrales. Se considera Γ
orientada positivamente.
Z
a) (x2 − y 2 ) dx + 2xy dy; Γ es la frontera de Ω = { (x, y) : 0 ≤ x ≤ 1, 2x2 ≤ y ≤
Γ
2x }
Z
b) x2 y dx + 2xy dy; Γ es la frontera de Ω = { (x, y) : 0 ≤ x ≤ 1, x2 ≤ y ≤ x }
Γ
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 triángulo de vértices (0, 0),
Γ
(4, 0) y (0, 4).
Z
e) Muestre que para constantes a, b y cualquier curva cerrada simple Γ, a dx +
Γ
b dy = 0.
2. Use el Teorema del Rotacional para demostrar que si Γ es una curva plana, cerrada,
simple (que no se cruza ella misma), orientada positivamente y Ω es la región del
plano acotada por Γ entonces el área de Ω, A(Ω), está dada por
Z Z Z
1
A(Ω) = − y dx = x dy = xdy − ydx.
Γ Γ 2 Γ
Z
3. Calcule ex sin y dx + (y 3 + ex cos y) dy, donde Γ es la frontera del rectángulo con
Γ
vértices (1, −1), (1, 1), (−1, 1) y (−1, −1), orientada positivamente.
4. Calcule el flujo del campo F a través de la superficie σ
a) F(x, y, z) = xi + yj + zk, σ : borde del cubo sólido Λ = { (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.
d ) F(x, y, z) = x2 yi + j + yzk, σ es la parte de la superficie cilı́ndrica x2 + y 2 = a2
limitada por los planos z = 0 y z = 3.
1
Para superficies planas el Teorema del Rotacional se conoce como Teorema de Green
159
5. En cada numeral use el Teorema de la Divergencia para calcular el flujo del campo F
a través de las superficies σ dada:
a) F(x, y, z) = xi + 2yj + 3zk, σ : x2 + y 2 + z 2 = 9

b) F(x, y, z) = x3 i + y 3 j + z 3 k, σ : x2 + y 2 + z 2 = 1
c) F(x, y, z) = 2i + 3j + 5k, σ : x2 + y 2 + z 2 = 1
6. Muestre que el flujo de un campo vectorial constante a través de cualquier superficie

cerrada es cero.
7. Encuentre un campo vectorial tal que su flujo a través de la superficie de una esfera
de radio r sea el área de esa esfera.
2
x2 2
8. El elipsoide a2 + yb2 + zc2 = 1 se puede parametrizar usando las coordenadas elipsoidales
x = a sin φ cos θ , y = b sin φ sin θ , z = c cos φ , for 0 ≤ θ ≤ 2π y 0 ≤ φ ≤ π.
Use esta parametrización para plantear una integral que permita calcular el área del
elipsoide. (Nota: La integral que usted va a plantear no puede evaluarse por me-
dios elementales. Para valores dados de las constantes se puede calcular por métodos
numéricos.)
R
9. Calculate Γ f (x, y, z) dL para la función f y la curva Γ dadas.
a) f (x, y, z) = z; Γ : x = cos t, y = sin t, z = t, 0 ≤ t ≤ 2π

x
b) f (x, y, z) = + y + 2yz; Γ : x = t2 , y = t, z = 1, 1 ≤ t ≤ 2
y
√
2 2 3/2
c) f (x, y, z) = z 2 ; Γ : x = t sin t, y = t cos t, z = 3 t , 0≤t≤1
10. En cada numeral calcule la circulación 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; Γ : x = cos t, y = sin t, z = t, 0 ≤ t ≤ 2π
c) F(x, y, z) = x i + y j + z k; Γ : x = cos t, y = sin t, z = 2, 0 ≤ t ≤ 2π
d ) F(x, y, z) = (y − 2z) i + xy j + (2xz + y) k; Γ : x = t, y = 2t, z = t2 − 1,
0≤t≤1
e) F(x, y, z) = yz i + xz j + xy k; Γ : el camino poligonal desde (0, 0, 0) a (1, 0, 0)
y hasta (1, 2, 0).
f ) F(x, y, z) = xy i + (z − x) j + 2yz k; Γ : el camino poligonal desde (0, 0, 0) a
(1, 0, 0) a (1, 2, 0) y hasta (1, 2, −2).
11. 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; σ : x2 + y 2 + z 2 = 1, z ≥ 0
b) F(x, y, z) = xy i + xz j + yz k; σ : z = x2 + y 2 , z ≤ 1
12. Sea σ una superficie cerrada y F un campo vectorial. Muestre que

ZZ
(rot F ) · N dA = 0.
σ
13. Para F(r, θ, z) = r er + z sin θ eθ + rz ez en coordenadas cilı́ndicas, encuentre div F y

rot F.
14. Para F(ρ, θ, φ) = eρ + ρ cos θ eθ + ρ eφ en coordenadas esféricas, encuentre div F y
rot F.
15. Demuestre las siguientes fórmulas siendo r = krk la norma del campo vectorial de
posición 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
16. 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)
17. Muestre que si un campo vectorial es conservativo, su rotacional es el vector cero.

18. Muestre que si F es un campo vectorial tridimensional, con segundas derivadas par-
ciales continuas entonces la divergencia del rotacional de F es cero.
19. Para cada uno de los siguientes campos F muestre que no existe ningún campo G
cuyo rotacional sea F.
161
a) F(x, y, z) = (2x + y)i + (y − 2z)3/2 j + y 3 k

3/2
b) F(x, y, z) = cos(x − 2z)i + sen( z)
x j + y3 k
20. Si F es un campo vectorial tridimensional, con derivadas parciales continuas, deter-

mine la densidad rotacional del campo para:
a) Un punto arbitrario del plano XY

b) Un punto arbitrario del plano XZ
c) Un punto arbitrario del plano Y Z
∂f ∂f
21. Deduzca la expresión del gradiente en coordenadas cilı́ndricas: ∇f = ∂r er + 1r ∂θ eθ +
∂f
∂z ez
22. Ponga F = u ∇v en el Teorema de la Divergencia para demostrar:

RRR RR
a) Primera identidad de Green (u ∆ v + (∇u) · (∇v)) dV = (u ∇v) · dA
Λ σ
RRR RR
b) Segunda identidad de Green (u ∆ v − v ∆ u) dV = (u ∇v − v ∇u) · dA
Λ σ
−y x
23. Sea F el campo definido en R3 −{(x, y, z) ∈ R3 : x2 +y 2 = 0} por F = h , , 0i.
x2 + y 2 x2 + y 2
Demostrar que rotF = 0 pero que la circulación de F a lo largo de cualquier circun-
ferencia paralela al plano XY , con centro en el eje Z no es nula, y por consiguiente
no puede existir una función f tal que F = ∇f .
24. Si f : [a, b] −→ R es diferenciable, no negativa y su gráfica de en el plano XY se gira

alrededor del eje X en R3 , se genera una superficie σ . Demostrar que el área de σ es
Z b p
2π f (t) 1 + (f ′ )2 (t)dt.
a
25. Sean, Ω ⊂ R2 abierto, f : Ω −→ R con derivadas parciales contiuas y σ = {(x, y, z) ∈

R3 : z = f (x, y)}. Parametrice σ, encuentre el elemento de área correspondiente y
demuestre que ZZ ZZ q
dA
= 1 + fx2 + fy2 dA
Ω cos θ Ω
siendo θ = θ(x, y, z) el ángulo que un vector normal a σ, en el punto (x, y, z), con
tercera componente positiva, forma con el eje Z.
p 1
26. Sean r = x2 + y 2 + z 2 , E = R3 − {0} y f : E −→ R dada por f (x, y, z) = .
r
a) Demostrar que f es armónica en E, es decir que △f = 0.
b) Calcular Z
df
dA,
σR dN
donde σR es la superficie de la esfera de radio R centrada en el origen.
c) Demostrar que no puede existir un campo vectorial F tal que ∇f = rotF.
d ) Calcular la circulación de ∇f a lo largo de las curvas: (a) x2 + y 2 + z 2 = R2 ,
√
z = a, 0 < a < R y (b) x2 + y 2 + z 2 = R2 , y = 3x
27. Sean, a, b ∈ R, a < b < 4a, Λ = {(ρ cos θ, y, ρ sin θ) : a ≤ ρ ≤ b, z − 4a < y <
4a − z, −π/6 < θ < π + π/6} y F el campo vectorial tridimensional definido por
F = (2a − z)k.
a) Calcular el flujo de F a través de ∂Λ.

b) Calcular el flujo de F a través de σ1 = ∂Λ ∩ {(ρ, y, θ) : y = −ρ sin θ + 4a}.
c) Calcular el flujo de F a través de σ2 = ∂Λ ∩ {(ρ, y, θ) : ρ = a}.
d ) Calcular la circulación de F a lo largo de σ1 ∩ σ2 .

Doce Lecciones de Calculo Vectorial Libr PDF

Cargado por

Copyright:

Formatos disponibles

Doce Lecciones de Calculo Vectorial Libr PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Doce Lecciones de Calculo Vectorial Libr PDF

Cargado por

Copyright:

Formatos disponibles

DOCE LECCIONES DE CÁLCULO VECTORIAL

Ernesto Acosta Gempeler

Escuela Colombiana de Ingenierı́a

2. Ecuaciones de superficies en el espacio 17

4. Funciones de varias variables 59

7. Extremos de funciones 101

8. Integrales múltiples 113

9. Integrales múltiples 129

10.Campos vectoriales 139

Hacemos un recuento muy rápido de los conceptos fundamentales de coordenadas de

Para localizar puntos en el plano, se introducen dos ejes numéricos perpendiculares

El mismo procedimiento se puede emplear para introducir un sistema de coordenadas

Dos puntos diferentes P y Q definen los segmentos de recta dirigidos PQ y QP.

Si introducimos un sistema de coordenadas en el espacio y (x1 , y1 , z1 ), (x2 , y2 , z2 )

La igualdad de las componentes escalares correspondientes entre segmentos dirigidos

este número representa también la longitud |RS| del segmento RS y la distancia

y la multiplicación del escalar a por el vector u es

au = hau1 , au2 , au3 i . (1.5)

Dividir un vector por un escalar no nulo consiste en multiplicar el vector por el

u · v = kuk kvk cos θ (1.6)

Al proyectar ortogonalmente el vector u sobre el vector v se obtiene un vector de-

al escalar kuk cos θ se le denomina la componente escalar de u en la dirección de v

ku × vk = kuk kvk sen θ, (1.9)

±vol(u, v, w) = ku × vk kwk cos γ

siendo γ el ángulo entre u × v y w .

ecuaciones o inecuaciones en dos o tres variables.

Por ejemplo, la recta que contiene un punto dado P y es paralela a un vector v es el

Si P tiene coordenadas (a, b, c) y v = hv1 , v2 , v3 i, un punto Q de coordenadas (x, y, z)

hx, y, zi = ha, b, ci + t hv1 , v2 , v3 i (1.12)

x = a + tv1 , y = b + tv2 , z = c + tv3 . (1.13)

A partir de las ecuaciones paramétricas se pueden escribir las ecuaciones simétricas

1. Para cada par de puntos P y Q determine las componentes escalares de los

a) P (1, 2, 1) y Q(3, 1, 2).

a) El segmento dirigido PQ.

3. En cada numeral determine si los segmentos PQ y RS son equivalentes, jus-

5. Dados los puntos P (p1 , p2 , p3 ), Q(q1 , q2 , q3 ) y R(r1 , r2 , r3 ), determine todos los

6. Muestre que la suma de vectores es conmutativa.

7. Muestre que el producto escalar de dos vectores es conmutativo.

8. Muestre que el producto vectorial no es conmutativo ni asociativo. ¿Qué rela-

9. Calcule la magnitud de cada uno de los siguientes vectores:

10. Sean v = h−1, 5, −2i y w = h3, 1, 1i. Calcule

El producto vectorial de dos vectores también se define como: Un vector que

22. Sea P el paralelepı́pedo en R3 definido por

a) Dibuje el paralelepı́pedo en un sistema de coordenadas cartesianas.

26. Sean v = h−1, 5, −2i y w = h3, 1, 1i.

a) Determine el vector u tal que u + v + w = i.

c) ¿Existe un escalar m tal que m(v + 2 w) = k?

27. Calcule v · w, v × w, la componente escalar de v en la dirección de w y la

a) v = h5, 1, −2i, w = h4, −4, 3i . b) v = − i+2 j+k, w = −3 i+6 j+3 k.

c) v = h1, 3, 2i, w = h7, 2, −10i . d ) v = i, w = 3 i + 2 j + 4k.

a) v = h5, 1, −2i, w = h4, −4, 3i . b) v = h7, 2, −10i, w = h2, 6, 4i .

c) v = h2, 1, 4i, w = h1, −2, 0i . d ) v = i, w = 3 i + 2 j + 4k

29. En cada numeral determine si v y w son o no perpendiculares. Justifique su

a) v = h8, 4, 3i y w = h−2, 1, 4i.

30. Para cada una de las siguientes afirmaciones, demuéstrela si es verdadera o

31. Sean α, β y γ los ángulos entre un vector no nulo v en R3 y los vectores i, j,

32. En cada numeral calcule u · (v × w) y u × (v × w)

a) u = h1, 1, 1i, v = h3, 0, 2i, w = h2, 2, 2i.

34. Si v y w son vectores unitarios en R3 , ¿bajo que condiciones se tendrı́a que