Cal Vectc9
Cal Vectc9
Cal Vectc9
Para aplicar el calculo a muchas situaciones reales y a las matematicas avanzadas, nece-
sitamos una descripcion matematica del espacio tridimensional.
Coordenadas rectangulares
1
a2 y a3 se denominan componentes de a en las direcciones x, y y z, respectivamente.
Considere un punto P(x, y, z) en el espacio. El vector r que parte del origen O hacia el
punto P se llama vector de posicion (o radio vector). Entonces, podemos escribir
p
r = xi + yj + zk, |r| = x2 + y 2 + z 2
1. x2 + y 2 + z 2 < 4
2. x2 + y 2 + z 2 4
3. x2 + y 2 + z 2 = 4, z 0
4. x2 + y 2 + z 2 + 3x 4z + 1 = 0
P0 X = tv (x x0 )i + (y y0 )j + (z z0 )k = t(v1 i + v2 j + v3 k)
Cuales son las ecuaciones parametricas de una recta L?
2
Planos Sea un punto P Rn y dos vectores c, Rn diferentes de cero y no paralelos.
Diremos que el conjunto de puntos X que determinan vectores P X que son combinacion
lineal de los vectores c y , es el plano P que pasa por el punto P y tiene como vectores
directores a c y .
2
EJEM. Encontremos una ecuacion vectorial del plano que contiene los puntos P = 5 ,
3
0 2
Q = 2 y R = 0
1 3
Cuando dos planos son paralelos?, Cuando dos planos son iguales? Cuando una recta L
con vector director Rn es paralella a un plano P con vectores directores c1 , 1 Rn ?.
Cuando una recta L con vector director Rn y un plano P con vectores directores
c1 , 1 Rn son ortogonales.
3
Coordenadas Cilindricas
Las coordenadas cilndricas (r, , z) de un punto P en el espacio son un hbrido natural de sus
coordenadas polares y rectangulares. Se usan las coordenadas polares (r, ) de un punto en el
plano con coordenadas rectangulares (x, y) y se usa la misma coordenada z como en las coordenadas
rectangulares.
De manera que las relaciones entre las coordenadas rectangulares (x, y, z) y las coordenadas polares
(r, , z) de un punto P en el espacio estan dadas por
Observese que la ecuacion r = r0 (r0 > O) es, en el sistema cilndrico (r0 , , z), la ecuacion
de un cilindro (circular recto) cuyo eje es el eje z. La proyeccion en el plano xy siempre dista
del origen una cantidad constante igual a r0 .
4
Ejemplo 3. Convertir el punto r = 4, 5
6 , 3 a coordenadas rectangulares.
Ejemplo 4. Convertir el punto (x, y, z)r = 1, 3, 2 a coordenadas cilndricas
Los planos verticales que contienen el eje z y los planos horizontales tambien tienen
ecuaciones simples de coordenadas cilndricas, como se muestra en la figura
5
Ejemplo 5. Hallar una ecuacion en coordenadas cilndricas para la superficie representada por
cada ecuacion rectangular, a) x2 + y 2 = 4z 2 , y b) y 2 = x
SOL: a) la grafica de la superficie x2 + y 2 = 4z 2 es un cono de dos hojas, con su eje a lo largo del
eje z, como se muestra en la figura. Luego,
SOL b) La grafica de la superficie es un cilindro parabolico con rectas generatrices paralelas al eje
z. Observe que,
Hay que observar que esta ecuacion comprende un punto en el que por lo cual nada se pierde al
dividir cada lado entre el factor r
Ejemplo 6. Hallar una ecuacion en coordenadas rectangulares de la superficie representada por
la ecuacion cilndrica, r2 cos 2 + z 2 + 1 = 0
SOL:
6
Coordenadas Esfericas
En el sistema de coordenadas esfericas, cada punto se representa por una terna ordenada:
la primera coordenada es una distancia, la segunda y la tercera coordenadas son angulos.
Un punto P en el espacio se representa por medio de una terna ordenada (, , ).
El nombre coordenadas esfericas se utiliza porque la ecuacion = const es una esfera con
mas precision, una superficie esferica de radio c centrada en el origen. La ecuacion = Const
describe (una parte de) un cono si 0 < c < 2 o si 2 < c < . La ecuacion esferica del plano
xy es = 2 .
Por tanto un punto P en el espacio esta determinado por la interseccion de un cono =
const, un plano = const y una esfera = const; de ah surge el nombre de coordenadas
esfericas. La relacion entre coordenadas rectangulares y esfericas es dada por:
Esfericas a rectangulares:
Rectangulares a esfericas:
y z
2 = x2 + y 2 + z 2 , tan = , arc cos p
x x2 + y 2 + z 2
Para cambiar entre los sistemas de coordenadas cilndricas y esfericas, usar lo siguiente.
Esfericas a rectangulares r 0:
r2 = 2 sin2 , = , z = cos
Cilindricas a esfericas (r 0)
p z
= r2 + z 2 , = , arc cos
r2 + z 2
Ejemplo 7. Hallar una ecuacion en coordenadas esfericas para la superficie representada por cada
una de las ecuaciones rectangulares a) Cono x2 + y 2 = z 2 b) Esfera x2 + y 2 + z 2 4z = 0
2 4 cos = 0, ( 4 cos ) = 0
= 4 cos
Hay que observar que el conjunto solucion de esta ecuacion comprende un punto en el cual = 0
de manera que no se pierde nada al eliminar el factor .
7
EJERCICIO I: Convierta las coordenadas esfericas (6, /4, /3)en coordenadas:
(a)
Coordenadas rectangulares. (b) Coordenadas cilindricas. SOL : (a) ( 23 2, 32 6, 3 2) (b)
(3 2, /3, 3 2)
EJERCICIO II: a) Encuentre las coordenadas rectangulares del punto P que tiene coor-
denadas esfericas (8, 5/6, /3)
b) Encuentre las coordenadas esfericas aproximadas del punto Q con coordenadas rectan-
gulares (3, 4, 12).
c) Encuentre la ecuacion en coordenadas esfericas del paraboloide con ecuacion en coor-
denadas rectangulares z = x2 + y 2 .
d) Determine la grafica de la ecuacion en coordenadas esfericas = 2 cos .
e) Determine la grafica de la ecuacion en coordenadas esfericas = sin sen .
F (x, y, z) = 0
es el conjunto de todos los puntos cuyas coordenadas (x, y, z) satisfacen esta ecuacion y recibe
el nombre de superficie. Por ejemplo,
Esf eras : (x x0 )2 + (y y0 )2 + (z z0 )2 = r2
P lanos : ax + by + cz + d = 0
Un tercer tipo de superficies que estudiaremos en el espacio son las llamadas superficies
cilndricas, superficies cuadradicas y superficies de revolucion.
Planos y trazas
Al graficar manualmente una superficie es util analizar las curvas que se forman al cortar
la superficie con planos paralelos a los planos coordenados. Estas curvas se conocen como
trazas.
Para bosquejar una superficie S, en ocasiones es util examinar sus intersecciones con varios
planos. La traza de la superficie S en el plano P es la interseccion de P y S. Por ejemplo, si
S es una esfera, se puede ver que la traza de S con un plano P es una circunferencia.
8
Cuando queremos visualizar una superficie especfica en el espacio, suele ser suficiente
analizar sus trazas en los planos coordenados y posiblemente unos cuantos planos paralelos
a ellos
SOL: Para hallar la traza en el eje xy hacemos z = 0. La ecuacion se reduce a la ecuacion 3x+2y = 6
de una recta en el plano xy. De manera similar, con y = 0 obtenemos la recta 3x + 2z = 6 como
la traza del plano dado en el plano xz. Para encontrar la traza en el plano yz, hacemos x = 0 y
obtenemos la recta y + z = 3. La figura muestra las porciones de estas tres trazas que se encuentran
en el primer octante. Todas juntas dan una buena idea de como esta situado el plano 3x+2y+2z = 6
en el espacio.
Superfices Cilindricas
{(x, y, z) : x2 + y 2 = 1, z arbitraria}
De modo similar, la grafica de una ecuacion tal como y + 2z = 2 es una recta en el espacio
bidimensional (el plano yz), pero en el espacio tridimensional la grafica del conjunto
{(x, y, z) : y + 2z = 2, x arbitraria}
Las superficies de este tipo reciben un nombre especial cilindros. Usamos el termino cilin-
dro en un sentido mas general que el de un cilindro circular recto. Especficamente, si C es
una curva en un plano y L es una recta no paralela al plano, entonces el conjunto de todos
los puntos (x, y, z) generado al mover una lnea que recorra a C paralela a L se denomina
cilindro. En palabras mas coloniales, un cilindro se genera al deslizar la curva C en la misma
direccon de la recta L, donde la recta L es representada por la variable que falta en su
ecuacion. La curva C recibe el nombre de generatriz del cilindro.
9
Como sugiere las graficas anteriores, cualquier curva
definen un cilindro paralelo al eje z, eje y y eje x, respectivamente cuya ecuacion tambien
es f (x, y) = c1 , g(x, z) = c2 y h(y, z) = c3 . Por tanto, conluimos que una curva en un plano,
cuando se consideran tres dimensiones, es un cilindro perpendicular a ese plano.
Ejemplo 9. Cilindros
Superficies Cuadricas
10
Hay seis tipos basicos de superficies cuadricas: elipsoide, hiperboloide de una hoja, hiper-
boloide de dos hojas, cono elptico, paraboloide elptico y paraboloide hiperbolico.
Elipsoide: La grafica de cualquier ecuacion de la forma
x2 y2 z2
2
+ 2 + 2 = 1,
a b c
corta a los ejes coordenados en (a, 0, 0), (0, b, 0) y (0, 0, c) los numeros reales positivos
a, b y c se llaman semiejes del elipsoide. La superficie es simetrica con respecto a cada uno
de los planos coordenados, ya que en la ecuacion que la define, cada variable esta elevada al
cuadrado.
x2 y2 z2
+ = ,
a2 b2 c2
recibe el nombre de cono elptico (o circular si el cono a = b)
x2 y2 z
+ = ,
a2 b2 c
es simetrico con respecto a los planos x = 0 y y = 0. La unica interseccion con los ejes es
el origen. Excepto por este punto, la superficie esta completamente por arriba (si c > 0) o
completamente debajo (si c < 0 ) del plano xy, segun el signo de c.
11
Hiperboloide de una hoja La grafica de una ecuacion de la forma
x2 y2 z2
+ = 1,
a2 b2 c2
se llama hiperboloide de una hoja. En este caso, un plano paralelo al plano xy, corta la
superficie en secciones transversales elpticas (o circulares si a = b) y es simetrico con respecto
a cada uno de los tres planos coordenados.
z2 x2 y2
= 1,
c2 a2 b2
se llama apropiadamente hiperboloide de dos hojas. Es simetrico con respecto a los tres planos
coordenados. El plano z = 0 no corta a la superficie; de hecho, para que un plano horizontal
corte a la superficie, debemos |z| c. tener
z y2 x2
= 2 2,
c b a
12
se conoce como paraboloide hiperbolico y tiene forma de una silla de montar.Su traza en el
plano horizontal z = z0 es una hiperbola (o dos rectas que se cruzan si z0 = 0).
Para clasificar una superficie cuadrica, se empieza por escribir la superficie en la forma
canonica o estandar. Despues, se determinan varias trazas en los planos coordenados o en
planos paralelos a los planos coordenados.
4x2 3y 2 + 12z 2 + 12 = 0
y2 x2 z2
=1
4 3 1
se puede concluir que la superficie es un hiperboloide de dos hojas con el eje y como su eje.
Para esbozar la grafica de esta superficie, conviene hallar las trazas en los planos coordenados.
y2 x2
T raza xy (z = 0) =1 Hiperbola
4 3
x2 z2
T raza xz (y = 0) + = 1 N o hay traza
3 1
y2 z2
T raza yz (x = 0) =1 Hiperbola
4 1
La superficie es un .
13
Ejemplo 12. Clasificar y dibujar la superficie dada por z = 4 x2 y 2 .
EJERCICIO:
Clasificar y dibujar las superficies dadas por
a) x2 + 2y 2 + z 2 4x + 4y 2z + 3 = 0.
b) 2x2 4y 2 + z 2 = 0
c) 2x2 + 4y 2 + z 2 = 36
SOL:
14
Superficies de Revolucion:
En Calculo integral vimos que una superficie S podra generarse rotando una curva plana C
alrededor de un eje. En la discusion que sigue se encontraran ecuaciones de superficies de revolucion
cuando C es una curva en un plano de coordenadas y el eje de revolucion es un eje de coordenadas
Ahora veamos un procedimiento para hallar su ecuacion. Consideremos la grafica de la funcion
radio
y = r(z), Curva generadora o directriz.
en el plano yz. Si esta grafica se gira sobre el eje z, forma una superficie de revolucion, como se
muestra en la figura. La traza de la superficie en el plano z = z0 es un crculo cuyo radio es r(z0 )
y cuya ecuacion es
x2 + y 2 = [r(z0 )]2
Sustituyendo z0 por z se obtiene una ecuacion que es valida para todos los valores de z. De
manera similar, se pueden obtener ecuaciones de superficies de revolucion para los otros dos ejes,
y los resultados se resumen como sigue.
Superficie de Revolucion: Si la grafica de una funcion radio r se gira sobre uno de los
ejes coordenados, la ecuacion de la superficie de revolucion resultante tiene una de las formas
siguientes.
1.) Girada sobre el eje x: y 2 + z 2 = [r(x)]2
2.) Girada sobre el eje y: x2 + z 2 = [r(y)]2
3.) Girada sobre el eje z: x2 + y 2 = [r(z)]2
Ejemplo 13. Hallar una ecuacion para la superficie de revolucion generada al girar la grafica
1
de y = en torno al eje z.
z
SOL: Como el giro es en el eje z entonces la ecuacion de la superficie es de la forma
r(z) = x2 + y 2 y como la funcion radio es r(z) = z1 . entonces la ecuacion de la superficie es:
1 2
x2 + y 2 =
z
Ejemplo 14. Hallar una ecuacion para la superficie de revolucion generada al girar la grafica
de 9x2 = y 3 en torno al eje y.
eje y.
15
Ejemplo 15. Hallar una curva generadora y el eje de revolucion de la superficie de revolucion
Hallar una directriz y el eje de revolucion de la superficie dada por x2 + 3y 2 + z 2 = 9
SOL: Sabemos que una superficie de revolucion tiene una de las siguientes ecuaciones
Como los coeficientes de x2 y z 2 son iguales, se debe elegir la tercera forma y escribir x2 +z 2 =
9 3y 2 . Entonces el eje y es el eje de revolucion. Se puede elegir una directriz de las trazas
siguientes
De manera que una ecuacion de una superficie generada al rotar una curva en un plano de coor-
denadas alrededor de
EJERCICIOS
1) Escriba una ecuacion del elipsoide de revolucion que se obtiene al girar la elipse 4y 2 +z 2 = 4
alrededor del eje z
2) Determine la grafica de la ecuacion z 2 = x2 + y 2 .
3) La grafica 4x2 + y 2 = 16 de se rota en torno al eje x. Encuentre una ecuacion de la
superficie de revolucion. 4) La grafica z = y, y 0 de se rota en torno al eje z. Encuentre
una ecuacion de la superficie de revolucion.
16
1.2. Funciones de varias variables
Hasta ahora en este texto, solo se han visto funciones de una sola variable (independiente).
Sin embargo, muchos problemas comunes son funciones de dos o mas variables. Por ejemplo,
1. Area de un rectangulo A = xy
4. Permetro de un rectangulo P = 2x + 2y
Nota: Puesto que se requieren cuatro dimensiones, no es posible graficar una funcion
de tres variables.
OBS:
Dominios e imagenes
Para definir una funcion de mas de una variable, seguimos la practica usual de excluir las
entradas que conducen a numeros complejos o a la division entre cero.
17
Ejemplo 16.
p
1. Dado que f (x, y) = 4 + x2 y 2 ) encuentre f (1, 0), f (5, 3) y f (4, 2).
SOL: (a)
f (1, 0) =
f (5, 3) =
f (4, 2) =
SOL: (b) Las coordenadas (x, y) debe cumplir x2 y 2 0, el dominio consiste en la region
sombreada de la figura.
Ejemplo 17. Calcule el dominio para las siguientes funiciones y si es posible indique su
grafica
p x
1. f (x, y) = 9 x2 y 2 4. g(x, y, z) =
9x2 y 2 z 2
2. f (x, y) = ln xy
2 2
3. f (x, y) =
x +y 9 5. g(x, y, z) = x+y+z
x 2x +y +z 2
2
SOL: (2) Dominio: las coordenadas (x, y) deben satisfacer . Es decir, el dominio de f
es
SOL: (3) Dominio: las coordenadas (x, y) de esta funcion deben satisfa-
cer . Es decir, el dominio de f son los puntos que estan
. OJO con el eje y
SOL: (5) Dominio: las coordenadas (x, y, z) de esta funcion deben satisfacer
. Es decir, el dominio
SOL: Para que f (x, y) este definida, xy 2 > 0 es decir, y 2 < x. Este dominio esta sombreado
en la figura. Observe que la parabola en aparece en lnea punteada para indicar que NO
esta incluida en el dominio de f .
Por otra parte f (x, y) = 1 implica que y 2 = 1, es decir, x = 2y 2 . Por lo tanto
xy
f (x, y) = 1 se obtiene en cada punto de la parabola x = 2y 2 [OJO (x, y) 6= (0, 0)].
18
Graficas
Como en el caso de las funciones de una sola variable, se puede saber mucho acerca del
comportamiento de una funcion de dos variables dibujando su grafica. La grafica de una
funcion f de dos variables es el conjunto de todos los puntos (x, y, z) para los que z = f (x, y)
y (x, y) esta en el dominio de f .
p
Ejemplo 19. Cual es el rango de f (x, y) = 16 4x2 y 2 . Describir la grafica de f .
El rango de f esta formado por todos los valores z = f (x, y) tales que .
Un punto (x, y, z) esta en la grafica de f si y solo si
Por lo tanto, la grafica de f es .
Curvas de nivel
En general, si una funcion de dos variables esta dada por z = f (x, y) entonces el conjunto
de puntos (x, y) en el dominio de f donde f (x, y) = c tiene un valor constante es una curva
de nivel de f . El conjunto de todos los puntos (x, y, f (x, y)) en el espacio, para (x, y) en
el dominio de f , se llama la grafica de f o tambien se les llama superficie z = f (x, y).
19
llama curva de contorno para distinguirla de la curva de nivel f (x, y) = c en el dominio
de f . Sin embargo, no todo mundo hace esta distincion; por ejemplo en la mayora de los
mapas, a las curvas que representan una elevacion (altura sobre el nivel del mar) se les llama
contornos, no curvas de nivel.
Ejemplo 21. Grafique f (x, y) = 100 x2 y 2 y trace las curvas de nivel f (x, y) = 0,
f (x, y) = 51 y f (x, y) = 75 en el dominio de f en el plano.
Ejemplo 22. Halle las curvas de nivel de una funcion polinomial f (x, y) = y 2 x2
SOL: las curvas de nivel son dadas por . Por tanto las curvas
de nivel son .
20
Superficies de nivel
Para una funcion de tres variables, w = f (x, y, z) las superficies definidas por f (x, y, z) = c
donde c es una constante, se llaman superficies de nivel de la funcion f . Aqu NO graficamos
la funcion pues los puntos (x, y, z, f (x, y, z)) que estan en R4 ; Sin embargo, podemos ver
como se comporta la funcion, estudiando sus superficies de nivel en su dominio, ya que ellas
muestran en que forma cambian los valores de la funcion al movernos por el dominio.
Ejemplo 24.
21
1.3. Lmites y continuidad
Para poder abordar adecuadamente el estudio de la diferenciabilidad de funciones de varias
variables es necesario tener algunos conceptos sobre limites y continuidad de estas funciones.
Aunque en estas notas de clase no voy hacer un acercamiento riguroso y exhaustivo sobre
limites y continuidad, debido a que a que se requiere un trabajo previo con la topologa del
espacio Rn (en el curso de analisis, exijan esta rigurosidad)
En el estudio de los lmites de las funciones de varias variables, se ponen al descubierto las
grandes dificultades de pasar del calculo de una variable al de varias variables: para funciones
de una variable sus dominios son pedazos de la recta, muchas veces intervalos. Para una
funcion de n variables, su dominio es un pedazo de Rn , y... aqu empiezan los problemas.
Como son los subconjuntos de Rn .equivalentes.a los subconjuntos de R? La respuesta la en-
contramos en la matematica llamada topologa. Nosotros solamente estudiaremos aquellos
conceptos topologicosque nos hagan la vida mas eficiente para comprender los conceptos
de lmite y continuidad para una funcion f : U Rn R.
En el caso de funciones de una variable, en muchos casos es factible hacer un juicio acerca
de la existencia de lm f (x) a partir de la grafica de y = f (x). Tambien se aprovecha que
xa
lm f (x) existe si y solo si lm f (x) y lm+ f (x) existe y son iguales al mismo numero L,
xa xa xa
en cuyo caso lm f (x) = L. En esta seccion veremos que la situacion es mas DIFICIL en la
xa
consideracion de lmites de funciones de dos variables.
El estudio del lmite de una funcion de dos variables inicia definiendo el analogo bidi-
mensional de un intervalo en la recta real. Utilizando la formula para la distancia entre dos
puntos (x, y) y (x0 , y0 ) en el plano, se puede definir el entorno de (x0 , y0 ) como el disco con
radio > 0 centrado en (x0 , y0 )
n p o
D = (x, y) : (x x0 )2 + (y y0 )2 <
Cuando esta formula contiene el signo de desigualdad menor que, <, al disco se le llama
abierto, y cuando contiene el signo de desigualdad menor o igual que, al disco se le llama
cerrado.
Definition 2. Un conjunto V Rn se dice que es abierto si para cada x V
existe r > 0 tal que B(x, r) V . Un conjunto F Rn se dice que es cerrado si su
complemento F c = Rn F es un conjunto abierto.
5. Una region que contiene algunos pero no todos sus puntos frontera no es ni abierta ni
cerrada.
22
Las definiciones de interior, frontera, abierta, cerrada, acotada y no acotada para las
regiones en el espacio, son similares a las de las regiones en el plano. Para considerar la
dimension extra, ahora usamos esferas solidas de radio positivo en lugar de discos.
lm f (x, y) = L
(x,y)(a,b)
Para tener un poco mas de precision, los puntos en el espacio (x, y, f (x, y)) pueden hacerse
arbitrariamente cercanos a (a, b, L) siempre que (x, y) sea suficientemente cercano a (a, b)
La nocion de (x, y) aproximandose a un punto (a, b) no es tan simple como para funciones
de una variable donde x a significa que x puede acercarse a a solo desde la izquierda y
desde la derecha.
23
En el plano xy hay un numero infinito de maneras de aproximarse al punto (a, b) para
que lm f (x, y) exista, requerimos ahora que f se aproxime al mismo numero L a lo
(x,y)(a,b)
largo de cualquier trayectoria o curva posible que pase por (a, b).
OBS:
1. Si f (x, y) no se aproxima al mismo numero L por dos trayectorias diferentes a (a, b),
entonces lm f (x, y) no existe.
(x,y)(a,b)
x2 3y 2 3y 2 3
lm = lm =
(x,y)(0,0) x2 + 2y 2 (0,y)(0,0) 2y 2 2
x=0
xy 0
lm = lm =0
(x,y)(0,0) x2 + y 2 (0,y)(0,0) y 2
x=0
xy
Sin embargo, esto NOOOOO significa que lm
exista, ya que no se ha exa-
+ y2 (x,y)(0,0) x2
minado toda trayectoria a (0, 0). Ahora, usemos todas las rectas que pasan por el origen
y = mx.
xy mx2 m
lm = lm =
(x,y)(0,0) x2 + y 2 2 2
(x,mx)(0,0) x + m x 2 1 + m2
y=mx
24
Ahora el limite depende de la pendiente m de la recta, concluimos que el lmite no existe.
Por ejemplo, tomando las rectas y = x y en y = 2x tenemos,
xy x2 1
lm = lm =
(x,y)(0,0) x2 + y 2 (x,x)(0,0) x2 + x2 2
y=x
xy 2x2 2
lm = lm =
(x,y)(0,0) x2 + y 2 (x,2x)(0,0) x2 + 4x2 5
y=2x
x3 y
Ejemplo 28. Demuestre que lm no existe.
(x,y)(0,0) x6 + y 2
x3 y
lm =
(x,y)(0,0) x6 + y 2
y=0
x3 y
lm =
(x,y)(0,0) x6 + y2
x=0
x3 y
lm =
(x,y)(0,0)
y=mx
x6 + y2
x3 y
lm =
(x,y)(0,0) x6 + y 2
y=x 3
x 2 y 2 2
lm =
(x,y)(0,0) x2 + y 2
x=0
x 2 y 2 2
lm = =
(x,y)(0,0)
y=x
x2 + y 2
25
SOL: Usando las coordenadas polares x = r cos y y = r sen tenemos
10xy 2 10r3 cos sen2
= = 10r cos sin2
x2 + y 2 r2
Por tanto,
10xy 2
lm = lm 10r cos sin2 = 0
(x,y)(0,0) x2 + y 2 r0
Observe, que el limite es independiente del valor de . De ahi que el limite exista.
lm f (x, y) = L
(x,y)(a,b)
significa que para toda > 0, existe un numero > 0 tal que
p
|f (x, y) L| < , siempre que 0 < (x a)2 + (y b)2 <
10xy 2
Ejemplo 31. Demuestre que lm = 0.
(x,y)(0,0) x2 + y 2
Observe que
10xy 2 y2
0 = 10|x|
2
x + y2 x2 + y 2
y2 p
10|x| 2 = 10|x| = 10 x2 10 x2 + y 2
y
< 10 .
2
De modo que si se elige = 10 , logramos tener x10xy
2 +y 2 0 < 10 = . Esto demuestra que
10xy 2
lm =0
(x,y)(0,0) x2 + y 2
5x2 y
Ejemplo 32. Demuestre que lm = 0.
(x,y)(0,0) x2 + y2
p
SOL: Dado > 0 fijo, debemos encontar un > 0 tal que, si |(x, y)(0, 0)| = x2 + y 2 <
entonces |f (x, y) 0| < Observe que
5x2 y x2
|f (x, y) f (0, 0)| = | | = 5|y|
x2 + y 2 x2 + y 2
2
x p p
5|y| 2 = 5|y| = 5 y 2 5 x2 + y 2
x
= 5 .
26
Por tanto, se puede elegir = 5 y concluir que si |(x, y)(0, 0)| <
5 entonces |f (x, y)0| <
5x2 y
es decir, lm = 0.
(x,y)(0,0) x2 + y 2
SOL: Debemos
p demostrar que para cualquier > 0 existe un > 0 tal que si k(x, y)
(1, 2)k = (x 1)2 + (y 2)2 < entonces |(3x2 + y) 5| <
Primero observe que
|(3x2 + y) 5| = |3x2 3 + y 2| 3|x2 1| + |y 2| = 3|x 1||x + 1| + |y 2|
Ahora el objetivo es usar para controlar los terminos |x 1|, |x + 1| y |y 2|. No es dificil
ver que
p p
|x 1| = (x 1)2 (x 1)2 + (y 2)2 <
p p
|y 2| = (y 1)2 (x 1)2 + (y 2)2 <
Para controlar el termino |x+1| debemos imponer una restriccion sobre el . Dicha restriccion
consiste en elegir un valor para por ejemplo = 1. De manera que el que estamos buscando
debera cumplir dos condiciones
p
(x 1)2 + (y 2)2 < , y 1
Usando la segunda condicion podemos decir que el termino
|x 1| < = 1, 1 < x 1 < 1, 1 < x + 1 < 3, |x + 1| < 3
Por tanto
|(3x2 + y) 5| 3|x 1||x + 1| + |y 2| < (3)(3) + = 10
Como nuestro objetivo final es que |(3x2 +y)5| < entonces el que buscamos debe cumplir,
1
10 < , es decir 10 . Esto significa que hemos impuesto dos restricciones sobre el ;
1 1
1 y 10 . Para que ambas restricciones sean validas debemos tomar = mn{1, 10 }
1
con es tenemos demostrado que para cualquier > 0 elegimos a = mn{1, 10 } y entonces
si p
(x 1)2 + (y 2)2 < entonces |(3x2 + y) 5| < 10 = .
Esto demuestra que lm 3x2 + y = 5.
(x,y)(1,2)
sen(x2 + y 2 )
Ejemplo 34. Demuestre que lm = 1.
(x,y)(0,0) x2 + y 2
sen h
Observe que si h = x2 + y 2 entonces lm = 1, esto implica que dado > 0 podemos
h0 h
hallar un 1 > 0 con 0 < 1 < 1 tal que
sen(h)
1 < , siempre que 0 < |h 0| < 1
h
2
Observe
que, si 0 < |t| < 1 entonces 0 < 2t < 21 < 1 , sacando raiz a
ambos lados,
0 < t < 1 . Deshaciendo el reemplazo h = x + y , y tomando como = 1 logramos
concluir que
p sen(x2 + y 2 )
Si 0 < x2 + y 2 < entonces 1<
x2 + y 2
sen(x2 + y 2 )
Esto demuestra que el lm = 1.
(x,y)(0,0) x2 + y 2
27
1.4.1. Continuidad
lm f (x, y) = f (a, b)
(x,y)(a,b)
x4 y 4
Ejemplo 35. Demuestre que la funcion f (x, y) = es discontinua en (0, 0), pero
x2 + y 2
4 4
x y
(x, y) 6= (0, 0)
F (x, y) = x2 + y 2
0 (x, y) = (0, 0)
SOL: Claramente la funicion f (x, y) es discontinua en (0, 0), ya que f (0, 0) no esta definida.
Por otra parte, F (0, 0) = 0 y
x4 y 4 (x2 y 2 )(x2 + y 2 )
lm = lm =0
(x,y)(0,0) x2 + y 2 (x,y)(0,0) x2 + y 2
SOL Claramente la funcion f (x, y) es continua en todos puntos (x, y) que verifican x2 +
y 6= 1. En efecto, considere los puntos (x0 , y0 ) que verifican x20 + y02 6= 1
2
lm f (x, y) = lm 0 = 0 = f (x0 , y0 )
(x,y)(x0 ,y0 ) (x,y)(x0 ,y0 )
x20 +y02 >1
28
Ahora estudiemos el caso en que (x0 , y0 ) que verifican x20 +y02 = 1 y veamos si lm f (x, y)
(x,y)(x0 ,y0 )
existe Sea S1 = (x, y) : x2 + y 2 1 y S2 = (x, y) : x2 + y 2 > 1
lm f (x, y) = lm 0=0
(x,y)(x0 ,y0 ) (x,y)(x0 ,y0 )
P S2
Como los limites no coinciden entonces f es discontinua en todos los puntos (x, y) que verifican
x2 +y 2 = 1. De manera que la region de continuidad es todo R2 menos la circuferencia unitaria
x2 + y 2 = 1.
SOL: Sea x (a, b) y consideremos la distancia L = mn x a, b x > 0 entonces
B(x, L) = (x L, x + L) (a, b)
esto es x es un punto interior de (a, b) y como x es arbitrario, entonces todo punto de (a, b)
es interior. Por tanto (a, b) es un conjunto abierto de R
Ejemplo 38. La bola abierta con centro en x0 y radio R > 0, es el conjunto abierto:
n o
B(x0 , R) = x Rn : kx x0 k < R
SOL: Sea x B(x0 , R), luego por definicion kx x0 k < R. De aqu podemos decir que
r := R kx x0 k > 0 y as considerar la bola B(x, r).
Nuestro objetivo es demostrar que B(x, r) B(x0 , R), para ello tomemos un y B(x, r)
el cual verifica, (ky xk < r) Ahora observe que
ky x0 k = ky x + x x0 k ky xk + kx x0 k r + kx x0 k = R
SOL: Sea (x, y) A, entonces x > 0. Por tanto, si tomaremos como radio r = x. Si
(x1 , y1 ) B((x, y), r) tenemos
p p
|x1 x| = (x1 x)2 = (x1 x)2 + (y1 y)2 < r = x
y as x1 x < x, es decir x1 > 0 por tanto (x1 , y1 ) A. entonces B((x, y), r) A, y por
tanto A es abierto.
SOL: Sea (a1 , a2 ) A1 A2 tenemos que demostrar que existe un 2-bola B(a1 , a2 ; r)
A1 A2 . Puesto que A1 , y A2 son abiertos en R existe una 1-bola B(a1 , r) en A1 , y una
29
1-bola B(a2 , r2 ) en A2 . Sea r = mn{r1 , r2 }. Podemos facilmente demostrar que la 2-bola
(B(a1 , a2 ; r) A1 A2 . En efecto, si (x1 , x2 ) es un punto cualquiera de B((a1 , a2 ), r) entonces
es un conjunto abierto en X.
esto es
B(x, ) f 1 B(f (x), ) f 1 (V )
2. A = (x, y) R2 : xy 2
7. E = (x, y) R2 : y < x2 , y < x1 , x > 0
8. F = (x, y) R2 : xy y + 1
2 3
4. A = (x, y) R : y x
5. C = (x, y) R2 : |x| < 1, |y| 2 9. G = (x, y) R2 : (x 1)2 + y 2 1, x 1
SOL (1) El conjunto representa al disco de centro C = (0, 0) y radio menor o igual a 2.
Si consideramos la funcion f : R2 R definida por f (x, y) = x2 + y 2 es continua y positiva,
30
el conjunto A se puede escribir como
= f 1 [0, 2]
= (x, y) R2 : g(x, y) (, 2]
= g 1 (, 2]
A = (x, y) R2 : 0 < k(x, y) (1, 3)k < 2 = (x, y) R2 : 0 < f (x, y) < 2
= f 1 (0, 2)
f (x, y) = x3 y 0
= f 1 [0, )
Sol (5): La representacion grafica del conjunto A se puede ver el la figura. Los puntos P y
Q estan en la frontera de C. Como P / C, vemos que C no es cerrado y como Q C,
vemos que C no es abierto. Vemos que (C) C 6= , por lo que el conjunto no es abierto.
Ademas, C 6= C por lo que el conjunto no es cerrado.
31
SOL (6) Para reescribir este conjunto D consideremos las funciones definidas de R2 en R
Como los intervalos (, 0), (0, ) R son abiertos y cada fi es continua entonces afirma-
mos que, el conjunto A es abierto. Observe que D es acotado porque esta contenido en la
bola de centro (0, 0) y radio 1, y como D D = , el conjunto es abierto.
SOL (7) Para reescribir este conjunto E consideremos las funciones definidas de R2 en R
1
f1 (x, y) = y x2 f2 (x, y) = y f3 (x, y) = x
x
Claramente cada fi es continuas y el conjunto E se puede escribir como
1
E = (x, y) R2 : y < x2 , y < , x > 0
x
= (x, y) R2 : f1 (x, y) < 0, f2 (x, y) < 0, f3 (x, y) > 0,
Como los intervalos (, 0), (0, ) R son abiertos y cada fi es continua entonces afir-
mamos que, el conjunto E es abierto. Observe que E no es acotado y como E E = , el
conjunto es abierto.
Sol (8) Para este conjunto podemos definir la funcion g : R2 R dada por g(x, y) = xy y
claramente es continua, el conjunto F se puede escribir como
= (x, y) R2 : g(x, y) (, 1]
= g 1 (, 1]
32
SOL (9) Para reescribir este conjunto G consideremos las funciones definidas de R2 en R
G = (x, y) R2 : (x 1)2 + y 2 1, x 1
= f11 , 1] f21 , 1]
33
Calculo Vectorial
Diferenciacion Vectorial
2. Las funciones importantes a estudiar, bajo la optica del calculo, son las funciones dife-
renciables.
5. El signo de f (x0 ) nos habla del crecimiento y/o decrecimiento de la funcion alrededor
del punto, etc. Este es,
6. En realidad, uno de los objetivos principales del calculo: obtener informacion de una
funcion diferenciable a partir de su derivada.
f f (x + h, y) f (x, y)
(x, y) = lm
x h0 h
y la derivada parcial con respecto a y es
f f (x, y + h) f (x, y)
(x, y) = lm
y h0 h
siempre que exista el lmite.
f
Observe que el la definicion (x, y) la variable y no cambia en el proceso del lmite,
y
en otras palabras, y se mantiene fija.
f
Obverse que el la definicion de (x, y) la variable x no cambia en el proceso del lmite,
y
en otras palabras, y se mantiene fija.
34
En un nivel practico
f
Para calcular , hay que tomar a y como constante y derivar respecto a x.
x
f
Para calcular , hay que tomar a x como constante y derivar respecto a y.
y
Las derivadas parciales son conceptos puntuales, es decir, se habla de la derivada
parcial de una funcion en un punto dado (de su dominio).
f f
En general, se debe hacer explcito el punto (x0 , y0 ) donde estan evaluados y
x y
f f
escribiendo (x0 , y0 ) y (x0 , y0 ).
x y
Muchas veces calculamos las derivadas parciales en un punto cualquiera (x, y) de
f f
su dominio. En tal caso, basta escribir y (fx y fy , o D1 f y D2 f , respectiva-
x y
mente).
Ejemplo 43. Sea f : U R2 R, f (x, y) = x2 y 3 . Entonces
f f (x + h, y) f (x, y)
= lm
x h0 h
=
f f (x, y + h) f (x, y)
= lm
y h0 h
=
En este ejemplo se observa que, como era de esperarse, las derivadas parciales de una
funcion z = f (x, y) se obtienen derivando parcialmenteccada una de las variables, y conser-
vando la otra como constante (es decir, pensando en la funcion f como dependiente solo de
x o de y)
Ejemplo 44. Calcular las derivadas parciales, de las siguientes funciones
1. f (x, y) = x2 + 2xy 2 y 3 .
2. z = (x2 + y 2 )exy .
3. z = 4x3 y 24x2 + y 6 + 1,
NOTA Smbolos como x y y se denominan operadores de diferenciacion parcial
y denotan la operacion de tomar una derivada parcial, en este caso con respecto a x y y. Por
ejemplo,
3 2 3 2
(x y sin y) = (x y ) (sin y)
x x x
x4 y 5 x4 y 5
(e sin xy) = (e ) (sin xy)
y x y
Ejemplo 45. Encuentre fx y fy (o zx y zy ) de las siguientes funciones
1. f (x, y) = x4 y 6 cos(x3 y 2 y 3 ).
2. z = (ln(x2 y) + tan(xy 2 )
35
p
3. z = sin( 4x3 y 24x2 y 6 + 1),
No podemos negar que esta definicionde diferenciabilidad para funciones de varias va-
riables goza de una gran simplicidad. Sera bueno que esta fuera en verdad la definicion
que buscamos. Sin embargo, este anhelo pronto se viene abajo pues es facil convencerse
que la existencia de las derivadas parciales de una funcion, (siendo una condicion necesaria)
esta muy lejos de ser una condicion suficiente para que la funcion sea diferenciable en el
sentido buscado.
xy
2 + y2
(x, y) 6= (0, 0)
2
Ejemplo 46. Considere la funcion f : R R definida por f (x, y) = x
0 (x, y) = (0, 0)
f f
Calculemos x (0, 0) y y (0, 0). Segun la definicion de derivadas parciales, tenemos
f
(0, 0) =
x
f
(0, 0) =
y
de modo que la funcion tiene ambas derivadas parciales iguales a en el origen. Sin
embargo, no es difcil demostrar que esta funcion NOOOOO es continua en (0, 0). (Ver
Ejemplo 27)
36
SOL: Para x 6= 0, se tiene
f
fy (x, 0) = (x, 0) =
y
f f (x + h, b) f (x, b)
(x, b) = lm
x h0 h
Interpretacion geometrica de fx = f
dx = zx . . .
El valor fx (x, b) es la pendiente de la recta tangente en P (a, b, c) de la curva x que
se encuentra sobre la superficie z = f (x, y).
37
Fig 4 Fig 5 Fig 6
As, se observa que el significado geometrico de fy es el siguiente:
Interpretacion geometrica de fy = f
dy = zy . . .
El valor fy (x, b) es la pendiente de la recta tangente en P (a, b, c) de la curva y que
se encuentra sobre la superficie z = f (x, y).
Observe que ahora tenemos dos rectas tangentes asociadas a la superficie z = f (x, y) en
el punto P (x0 , y0 , f (x0 , y0 )) (figura 7). El plano determinado por estas rectas tangentes es
tangente a la superficie en P?
Fig 7
Veremos que s, pero tenemos que aprender mas sobre las derivadas parciales antes de
saber por que.
38
Ejemplo 48. Considere un plano x = 1 el cual corta al paraboloide z = x2 + y 2 en una
parabola. Determine la pendiente de la tangente a la parabola en (1, 2, 5)
z
SOL: La pendiente es el valor de la derivada parcial en (1, 2):
y
z
=
y
(1,2)
Ejemplo 49. Hallar las pendientes en las direcciones de x y de y de la superficie dada por
2
f (x, y) = x2 y 2 + 25 1
8 en el punto ( 2 , 1, 2).
f
=
y
( 12 ,1)
Ejemplo 50. Hallar las pendientes de la superficie dada por f (x, y) = 1 (x 1)2 (y 2)2
en el punto (1, 2, 1), en las direcciones de x y de y.
f
=
y
(1,2)
Sin importar cuantas variables haya, las derivadas parciales se pueden interpretar como
tasas, velocidades o razones de cambio.
Ejemplo 51. El area de un paralelogramo con lados adyacentes a y b entre los que se forma
un angulo esta dada por A = ab sen ,
39
SOL (1): Para hallar la tasa o la razon de cambio del area respecto de a, se mantienen b
y constantes y se deriva respecto de a para obtener:
A
=
a
(10,20, 6 )
SOL (2): Para hallar la tasa o la razon de cambio del area respecto de , se mantienen a y b
constantes y se deriva respecto de a para obtener:
A
=
(10,20, 6 )
f f (x + h, y, z) f (x, y, z)
= lm aqu y, z son constantes
x h0 h
f f (x, y + h, z) f (x, y, z)
= lm aqu x, z son constantes
y h0 h
f f (x, y, z + h) f (x, y, z)
= lm aqu x, y son constantes
z h0 h
Las derivadas parciales de las funciones con aun mas variables se definen de manera analoga.
Una funcion f (x1 , x2 , . . . , xn ) de n variables tiene n derivadas parciales, una respecto a cada
una de sus variables independientes. Por ejemplo,
f f (x1 , . . . , xi1 , xi + h, xi+1 , . . . , xn ) f (x1 , . . . , xn )
= lm aqu xi son constantes
xi h0 h
f f (x1 , . . . , xi , . . . , xn ) + (0, . . . , h, . . . , 0) f (x1 , . . . , xn )
= lm
xi h0 h
f f (x 1 , . . . , x n ) + h(0, . . . , 1, . . . , 0) f (x1 , . . . , xn )
= lm
xi h0 h
f f x + hei f (x)
= lm donde x = (x1 . . . , xn )
xi h0 h
40
Esta ultima ecuacion es conocida como la Derivada Direccional de la funcion f en la direccion
del vector ei .
Ejemplo 53. Si f (x, y, t) = e3t cos(4x) sen(6y), entonces las derivadas parciales con res-
pecto a x, y y t son:,
SOL:
f
=
x
f
=
y
f
=
t
Ejemplo 54. Si g(x, y, u, v) = eux cos(vy), entonces las derivadas parciales con respecto a
x, y, u y v son:,
SOL:
f
=
x
f
=
y
f
=
u
f
=
v
Si se escribe z = f (x, y), entonces se puede reemplazar cada ocurrencia del smbolo f por z.
41
4. El orden de los smbolos en los subndices de las parciales mixtas es justamente lo
opuesto al orden de los smbolos cuando se usa la notacion de operador de diferenciacion
parcial:
fx f 2f
(fx )y = fxy = = =
y y x yx
fy f 2f
(fy )x = fyx = = =
x x y xy
Ejemplo 55. Hallar las derivadas parciales de segundo orden de f (x, y) = 2xy 2 3x + 3x2 y 2
y determinar el valor de fxy (1, 2)
fxy = fyx
2w 2w ey
Ejemplo 57. Encuentre y si w = xy+ 2 . Que opinas de estos dos calculos?
xy yx y +1
SOL:
42
numero de veces que se pueda derivar una funcion, siempre que las derivadas implicadas
existan. As, podemos obtener derivadas de tercer y cuarto orden, denotadas por smbolos
como
3w 4w
= f yyx , = fyyxx
xy 2 x2 y 2
etcetera. Como con las derivadas de segundo orden, el orden de derivacion no importa, mien-
tras todas las derivadas hasta el orden en cuestion sean continuas.
p
Ejemplo 58. Si f (x, y, z) = x2 + y 4 + z 6 determine fyzz
SOL:
Ejemplo 59. Mostrar que fxz = fzx y fxzz = fzxz = fzzx para la funcion dada por
f (x, y, z) = yex + x ln z
SOL:
43
Diferenciabilidad de z = f (x, y)
En el caso de una funcion de una sola variable se asumio que y = f (x) era diferenciable en x0 , esto
es,
f (x + x) f (x)
f (x0 ) = lm (2.1)
x0 x
existe. Recuerde tambien que si f es diferenciable en x0 , tambien es continua en ese numero.
Al repetir la suposicion en (2.1), deseamos que una funcion de dos variables z = f (x, y) sea
diferenciable en un punto (x0 , y0 ). Aunque hemos considerado lo que significa que z = f (x, y)
posea derivadas parciales en un punto, aun no formulamos una definicion de diferenciabilidad
de una funcion de dos variables f en un punto.
z = f (x + x, y + y) f (x, y) (2.2)
Ejemplo 60. Encuentre z para la funcion polinomial z = x2 xy Cual es el cambio en
la funcion de (1, 1) a (1.2, 0.7)
Solucion.
z = f (x + x, y + y) f (x, y)
= (x + x)2 (x + x)(y + y) (x2 xy)
= (2x y)x xy + (x)2 (x)(y) (2.3)
Con x = 1, y = 1, x = 0.2 y y = 0.3,
z = 1(0.2) 1(0.3) + (0.2)2 (0.2)(0.3) = 0.6
Una breve reinspeccion del incremento z en (2.3) muestra que en los primeros dos
terminos los coeficientes de x y y son fx (x, y) y fy (x, y) respectivamente. El importante
teorema que sigue muestra que esto no es un accidente.
Teorema 8. Considere que z = f (x, y) tiene derivadas parciales continuas fx (x, y) y
fy (x, y) en una region rectangular abierta que esta definida por a < x < b, c < y < d.
Si (x, y) es cualquier punto en esta region, entonces existe una funcion r(x, y), tal que
donde
lm r(x, y) = 0
(x,y)(0,0)
Demostracion.
44
Una simple maniputacion de z encontramos
donde,
r(x, y) = fx (x0 , y + y) fx (x, y) x + fy (x, y0 ) fy (x, y) y,
Observe que,
donde
r(x, y)
lm =0
(x,y)(0,0) kx, yk
El siguiente teorema es una implicacion directa del Teorema 8 el cual proporciona una
condicion suficiente bajo la cual la existencia de las derivadas parciales implica diferenciabi-
lidad.
Teorema 10 (Condicion suficiente para la diferenciabilidad).
Si las primeras derivadas parciales fx y fy son continuas en un punto en una region abierta
R, entonces z = f (x, y) es diferenciable sobre R.
Veamos ahora que esta definicion SIIII garantiza la continuidad de la funcion en el punto
(x0 , y0 ) en donde es diferenciable (como acontece con funciones de una variable).
Tornando el lmite aqu cuando (x, y) (0, 0), y observando que de la condicion estable-
45
cida en la definicion para el residuo r(x, y), se deduce que
r(x, y)
lm =0 lm r(x, y) = 0
(x,y)(0,0) k(x, y)k (x,y)(0,0)
Por tanto,
lm [f (x0 +x, y0 +y) = lm f (x0 , y0 )+fx (x0 , y0 )x+fy (x0 , y0 )y+r(x, y) = f (x0 , y0 )
(x,y)(0,0) (x,y)(0,0)
SOL: Del Ejemplo 46 sabemos que f no es continua en (0, 0) por tanto NOOOOO es
diferenciable en (0, 0). Otra forma de ver la no diferenciabilidad (0, 0) es usando la propia
definicion de ser diferenciable, primero que todo las derivadas parciales SI existen fx (0, 0) =
fy (0, 0) = 0 (Ejemplo 46). Ahora veamos incremento en z en puntos cercanos de (0, 0)
de donde
xy
r(x, y) =
(x)2 + (y)2
y
r(x, y) xy
lm = lm
(x,y)(0,0) k(x, y)k (x,y)(0,0) (x)2 + (y)2 3/2
(x)2
= lm = (NO EXISTE)
(x,y)(0,0) 23/2 |x|3
x=y
Ejemplo 62. Verifiquemos que la funcion f (x, y) = xy 2 es diferenciable en punto (0, 0).
46
OBSERVACIONES
1. A medida que (x, y) tiende a cero, r(x, y) tambien tiende a cero (r sera llamado
error o residuo, pensando en que es lo que le falta al plano tangente para describir a
la superficie z = f (x, y) en los alrededores de (x0 , y0 )).
2. OJO la condicion de
lm r(x, y) = 0,
(x,y)(0,0)
r(x, y)
lm =0
(x,y)(0,0) kx, yk
Intuitivamente, podemos pensar que la unica manera como se puede dar esta condicion
del residuo, es cuando el plano tangente se embarra con la superficie en los alrededores
de (x0 , y0 ). En otras palabras, la superficie tiene que ser suave.en (x0 , y0 ), para que
se pueda ver localmente como un plano tangente.
5. Si z = f (x, y) la funcion es diferenciable en cada punto en una region R del plano xy,
entonces se dice que f es diferenciable en R.
6. Si f es diferenciable sobre plano xy completo, entonces f es diferenciable.
Ejemplo 63. Mostrar que la funcion dada por f (x, y) = x2 + 3y es diferenciable en todo
punto del plano.
z = f (x + x, y + y) f (x, y)
= (x2 + 2xx + x2 ) + 3(y + y) (x2 + 3y)
= 2x(x) + 3(y) + x(x) + 0(y)
= fx (x, y)x + fy (x, y) + r(x, y)
r(x, y)
donde r(x, y) = x(x) + 0(y). Obviamente, lm = 0. Por tanto
(x,y)(0,0) kx, yk
f es diferenciable en todo punto en el plano.
47
Linealizacion de una funcion z = f (x, y)
Empezamos el analisis de la derivada con el problema de encontrar la recta tangente a la
grafica de una funcion y = f (x) en un punto (a, f (a)). Intuitivamente, es de esperar que una
recta tangente este muy proxima a la grafica de f siempre que x este cerca del numero a. En
otras palabras, cuando x esta en una pequena vecindad de a, los valores de la funcion f (x)
estan muy proximos al valor de las coordenadas y de la recta tangente. As, al encontrar
una ecuacion de la recta tangente en (a, f (a)) podemos usarla para aproximar a f (x).
Sabemos que la ecuacion de la recta tangente esta dada por
Esta ultima ecuacion se escribira como L(x) = f (a) + f (a)(x a). Esta funcion lineal recibe
un nombre especial.
Definition 12 (Linealizacion de y = f (x)).
Si una funcion y = f (x) es diferenciable en un punto a entonces decimos que la funcion
f (x) L(x)
Empleando x = x0 + x y y = y0 + y obtenemos
48
Definition 13 (Linealizacion de z = f (x, y)).
Si una funcion z = f (x, y) es diferenciable en un punto (x0 , y0 ) entonces la funcion
se dice que es una linealizacion de f en (x0 , y0 ). Para un punto (x, y) cercano a (x0 , y0 )
la aproximacion
f (x, y) L(x, y)
se denomina una aproximacion lineal local de f en (x0 , y0 ).
OBSERVACIONES:
Esto demuestra que la formula de linealizacion z = L(x, y) es una ecuacion del plano
tangente a la grafica de z = f (x, y) en el punto (x0 , y0 , f (x0 , y0 )) y tiene como vector
normal a fx (x0 , y0 ), fy (x0 , y0 ), 1
Cuando fx y fy son continuas y cuando (x, y) (0, 0), entonces dz es una apro-
ximacion de z esto es, dz = z.
p
Ejemplo 65. Encuentre una linealizacion de f (x, y) = x2 + y 2 en (4, 3)
Suponga que y = f (x) es una funcion diferenciable en el punto a, suponemos que x cambia
de su valor original al nuevo valor a + x. El cambio en el valor de y es el incremento y,
calculado restando el nuevo valor de y menos su valor anterior, esto es,
y = f (a + x) f (a) (2.9)
y f (a + x) f (a)
=
x x
49
Pero para valores x de que estan proximos a 0, (esto es, x cerca de a), este cociente es una
aproximacion del valor de la derivada de f en a:
y
f (a) o bien y f (a)x
x
Las cantidades x y f (a)x se denominan diferenciales y se denotan por los smbolos dx
y dy, respectivamente. Es decir,
x = dx dy = f (a)dx
dy f (a)dx
= = f (a)
dx dx
y coincide totalmente con la notacion que se ha usado. De hecho, Leibniz origino la notacion
diferencial al visualizar incrementos infinitesimalesdx y dy, donde su razon dy/dx es la
pendiente de la recta tangente. La clave del descubrimiento independiente de Leibniz del
calculo diferencial en la decada de 1670 fue su entendimiento de que si dx y dy son suficien-
temente pequenos, entonces el segmento de la curva y = f (x) y el segmento de recta que
une (x, y) con (x + dx, y + dy) son virtualmente indistinguibles.
Esta concepto se ilustra con las amplificaciones sucesivas en las siguiente grafica de la
curva y = x2 cerca del punto (1, 1).
L = L(a + x) L(a)
= [f (a) + f (a)(a + x a)] + [f (a) + f (a)(a a)]
= f (a)dx
= dy
En el caso de una funcion f de dos variables tenemos naturalmente tres diferenciales. Los
cambios en las variables independientes x y y son dx y dy; los cambios en la linealizacion
50
L(x, y) se denotan por medio de dz. En el punto (x0 , y0 ) el cambio en la linealizacion es
L = L(x0 + x0 , y0 + y) L(x0 , y0 )
= [f (x0 , y0 ) + fx (x0 , y0 )(x0 + x x0 ) + fy (x0 , y0 )(y0 + y y0 )]
f (x0 , y0 ) + fx (x0 , y0 )(x0 x0 ) + fy (x0 , y0 )(y0 y0 )]
= fx (x0 , y0 )x + fy (x0 , y0 )y
por lo que en este ejemplo la diferencial parece ser una buena aproximacion al incremento.
Las definiciones, as como los teoremas visto para funciones de una variable y = f (x), y dos
variables z = f (x, y) se generalizan de la manera esperada a funciones de tres o mas variables.
A continuacion mencionare algunos puntos importantes para una funcon w = f (x, y, z).
w = f (x + x, y + y, z + z) f (x, y, z).
w = fx x + fy y + fz z + r(x, y, z)
donde
r(x, y, z)
lm = 0.
x,y,z)(0,0,0) k(x, y, z)k
51
3. Si f es diferenciable en (x0 , y0 , z0 ) entonces la linealizacion de f se define como
L(x, y, z) = f (x0 , y0 , z0 )+fx (x0 , y0 , z0 )(xx0 )+fy (x0 , y0 , z0 )(yy0 )+fz (x0 , y0 , z0 )(zz0 )
4. La diferencial total de f es
f f f
dw = dx + + dz
x y z
donde
r(h)
lm =0
h0 khk
f f
(aqu h := (x1 , . . . , xn ), f := ( x 1
, . . . x n
)).
x = g1 (u, v) y = g2 (u, v)
52
Ejemplo 68. Considere la funcion f (x, y, z) = 2x3 y 2 z + cos2 (x + y + z) y sean x = et ,
v = cos t, y z = t2
53
Ahora vamos a extender la regla de la cadena a funciones de varias variables. Si z = f (x, y) y
x y y son funciones de una sola variable t, entonces el siguiente teorema indica como calcular
dz
la derivada ordinaria.
dt
dz z dx z dy
= +
dt x dt y dt
dz
Ejemplo 72. Si z = x3 y y 4 y x = 2t2 , y = 5t2 6t calcule en en t = 1.
dt
Solucion. No es difcil ver que
dz z dx z dy
= + = (3x2 y)(4t) + (x3 4y 3 )(10t 6)
dt x dt y dt
dz
En este caso, en t = 1, x(1) = 2 y y(1) = 1, por lo que = 0.
dt t=1
NOTA: Aunque NO HAY necesidad de hacerlo de esa manera, tambien podemos encontrar
la derivada dz 2 2 3 4
dt al sustituir las funciones x = 2t , y = 5t 6t en z = x y y y despues
6 2 2 4
derivar la funcion resultante de una sola variable z = 8t (5t 6t) (5t 6t) con respecto
a t.
Ejemplo 73. Sea f : U R2 R la funcion dada por f (x, y) = x2 + 3y 2 , y sea g : U
R R2 la funcion g(t) = (et , cos t). Calcule sus derivadas parciales de f g
Solucion.
Regla de la cadena para derivadas parciales Para una funcion compuesta de dos
variables z = f (x, y) donde x= g(u, v) y y = h(u, v) se esperaran naturalmente dos formulas,
z z
ya que z = f g(u, v), h(u, v) y por ello pueden calcularse tanto u como v . La regla de la
cadena para funciones de dos variables se resume en el siguiente teorema.
54
Ahora bien, si
x = g(u, v + v) g(u, v) y y = h(u, v + v) h(u, v),
entonces como x = g(u, v) y y = h(u, v) tenemos que g(u, v + v) = x + x y h(u, v + v) =
y + y al reemplazar en z obtenemos que
z = f (x + x, y + y) f (x, y)
Puesto que f es diferenciable, entonces su incremento z puede escribirse como
z z
z = x + y + 1 x + 2 y
x y
donde 1 , 2 , son funciones de x y y con la propiedad de que 1 , 2 0 cuando (u, v)
(0, 0). Por tanto,
z z x z y x y
= + + 1 + 2
v x v y v v v
Ahora, tomando el lmite de la ultima lnea cuando v 0 obtenemos
z z x z y x y z x z y
= + +0 +0 = +
v x v y v v v x v y v
puesto que x y y se aproximan a cero cuando v 0.
2 3 2u3v 2 2 z z
Ejemplo 74. Si z = x y y x = e , y = sen(u v ), determine u y v
Solucion.
z x
1. = 4. =
x v
z y
2. = 5. =
y u
x y
3. = 6. =
u v
Por tanto,
z
=
x
z
=
y
2 2 2 2 2 3 3
Ejemplo 75. Sea u = 3x y z +6 sen(xyz), donde x = 4vw y = 5v +10w z = v . Calcule
u u
v , w
Solucion.
55
z z
Ejemplo 76. Consideremos la funcion z = f (x2 + y 3 , 2x2 y 2 ) y hallemos x y y
z f u f v
= + =
x u x v x
z z u f v
= + =
y u y v y
Ejemplo 77. Consideremos la funcion z = f ( xy ) y hallemos z
x y z
y
Solucion.
Generalizaciones regla de la cadena: Si z = f (x1 , x2 . . . , xn ) es diferenciable en
(x1 , x2 , . . . , xn ) y si xi , i = 1, . . . , n son funciones diferenciables de una sola variable t, enton-
ces
dz z dx1 z dx2 z dxn
= + + +
dt x1 dt x2 dt xn dt
De manera similar, si z = f (x1 , x2 . . . , xn ) y si cada xi , i = 1, . . . , n en una funcion de k
variables u1 , u2 , . . . , uk entonces
dz z dx1 z dx2 z dxn
= + + +
dui x1 dui x2 dui xn dui
donde i = 1, 2, . . . , k
56
r
Ejemplo 78. Si r = x2 + y 5 z 3 y x = uve2s , y = u2 v 2 s z = sen(uvs2 ), encuentre y
u
r
s
Solucion. En este caso r es una funcion de tres variables x, y y z, y cada una es en s misma
una funcion de tres variables u, y y s.
r
=
u
r
=
s
dz
Ejemplo 79. Si z = u2 v 3 w4 y u = t2 , v = 5t 8 y w = t3 + t, determine .
zt
Solucion. Siguiendo el diagrama de arbol tenemos
z
=
t
Diferenciacion implcita
57
y = f (x). Esta misma observacion la habamos hecho con superficies z = f (x, y), las cuales
pueden ser vistas como el nivel cero de la funcion w = F (x, y, z) = z f (x, y).
El punto que llama nuestro interes ante estas consideraciones es el planteamiento recproco
de ellas: por ejemplo, dada la funcion z = F (x, y), es su nivel cero una curva que se pueda
ver como la grafica de una funcion y = f (x)?, o de otra manera, de la expresion F (x, y) = 0
(que da el nivel cero de F ) podemos despejar a y en terminos de x y dejar as establecida
la funcion y = f (x)? Ahora nuestra pregunta es, de hecho, mas especfica: se quiere saber si
el nivel cero de z = F (x.y) es la grafica de una funcion y = f (x) Al pensar en la simple
funcion F (x, y) = x2 + y 2 1, vemos que la respuesta a la pregunta anterior es NO. De la
expresion x2 +l1 = 0 no podemos establecer una funcion y = f (x) despejando y enterminos
de x. Es bien sabido
que x2 + y 2 = 1 define dos funciones, a saber, y1 = f1 (x) = 1 x2 y
y = f2 (x) = 1 x2 .
Si la ecuacion F (x, y) = 0 define a una funcion y = f (x) de manera implcita, F (x, f (x)) =
0 entonces para toda x en el dominio de f . El siguiente teorema demuestra que la derivada
dy
dx tambien puede determinarse de la regla de la cadena.
dy Fx (x, y)
= , donde Fy (x, y) 6= 0
dx Fy (x, y)
z Fx (x, y, z) z Fy (x, y, z)
= =
x Fz (x, y, z) y Fz (x, y, z)
donde Fz (x, y, z) 6= 0
58
dy z
a) Encuentre si x2 4xy 3y 2 = 10 b) Encuentre si x2 y 5xy 2 = 2yz 4z 3 .
dx y
dy Fx (x, y) 2x 4y x 2y
De aqui deducimos que = = =
dx Fy (x, y) 4x 6y 2x + 3y
(b) Consideremos la funcion G(x, y, z) = x2 y 5xy 2 2yz + 4z 3 . Entonces definimos z
como una funcion de x y y mediante el uso de F (x, y, z) = 0. Es decir F depende solo de
las variables x, y. Luego podemos derivar respecto a x o y. En nuestro caso vamos a derivar
respecto a y. Usando la regla de la cadena encontramos que
dF dy dF dx dF dz
+ + =0
dy dy dx dy dz dy
y x z Fy (x, y, z) x2 10xy 2z
Como y = 1, y = 0 entonces = =
y Fz (x, y, z) 2y + 12z 2
Ejemplo 81.
Suponga que la ecuacion z 2 = x2 + xyz 2 define a z implcitamente como una funcion de x y
z z
y. Encuentre x y x .
SOL:
59
Derivada Direccional
en todo lugar en que existan estos lmites. Aqu x = (x, y), e1 = (1, 0) y e2 = (0, 1), como
es habitual. As, fx y fy representan tasas de cambio de z respecto a la distancia en las
direcciones de los vectores unitarios e1 y e2 . No hay razon para restringir nuestra atencion
solo a dos direcciones; podemos encontrar la tasa de cambio de una funcion diferencial en
cualquier direccion.
60
La derivada direccional de f en el punto (x0 , y0 ) D en la direccion del vector unitario
u = (u1 , u2 ), esto es, Du f (x0 , y0 ) mide la razon (o velocidad) de cambio instantaneo
del valor de la variable dependiente z = f (x, y) con respecto a la distancia en el plano
XY , medida en la direccion del vector unitario u.
La funcion f en la ecuacion (2.12) puede ser una funcion de dos, tres o mas variables.
las derivadas parciales de una funcion de dos variables x y y se escriben como
u = (cos , sin ), 0 2
f (x + h) f (x) f (x) h
lm =0
h0 |h|
En decir,
Du f (x) = f (x) u = kf kkuk cos = kf k cos
61
Propiedades de la Derivada Direccional:
f se lee como nabla. Otra notacion para el gradiente es gradf (x, y). OOJOO para cada
(x, y) el gradiente f (x, y) es un vector en el plano (no un vector en el espacio).
d
Nota: El smbolo no tiene ningun valor. Es un operador de la misma manera que dx es
un operador. Cuando opera sobre f (x, y) produce el vector f (x, y)
Ejemplo 82. Sea f (x, y) = exy y 2 . Halle la derivada direceional de f en cualquier punto
(x, y) Df , en la direccion del vector unitario u = ( 12 , 12 )
f (x, y) =
Por tanto, la derivada direccional de f en el punto (x, y) en la direccion u esta dada por
Du f (x, y) = f (x, y) u =
P1 Q 1
u=
kP1 Q1 k
62
f (x, y, z) =
f (2, 2, 4) =
Du f (2, 2, 4) = f (2, 2, 4) u =
Ejemplo 84. Determine las direcciones en que f (x, y) = (x2 /2) + (y 2 /2), crece mas rapida-
mente, decrece mas rapidamente en el punto (1, 1) y Cuales son las direcciones de cambio
nulo de f en (1, 1)?
Solucion. Como sabemos cualquier funcion f (x, y) crece mas rapidamente en la direccion de
f (x, y)
su f , es decir u = en nuestro caso, f (x, y) = y por tanto,
|f (x, y)|
f (1, 1)
f (1, 1) = u= =
|f (1, 1)|
Adicionalmente la funcion decrece mas rapidamente en la direccion de en f en (1, 1), es
decir u =
Por otra parte, las direcciones de cambio nulo en (1, 1) son las direcciones ortogonales a f .
es son:
n= n =
Ejemplo
p 85. Cual es el valor del angulo para el cual la derivada direccional de f (x, y) =
25 x y 2 en el punto (1, 2) es mnimo y cual es este valor mnimo?
2
Solucion. Aqu tenemos que usar usar el vector unitario u = (cos , sen ), y el gradiente de
f esta dado por
f (x, y) = f (1, 2) =
g() = Du f (1, 2) u =
g () =
g () =
g ( )=
63
Du f (1, 2) =
Teorema 22.
Si f (x, y) = c es la curva de nivel de una funcion diferenciable de dos variables z = f (x, y)
que pasa por un punto especificado (x0 , y0 ) entonces f (x0 , y0 ) es normal (ortogonal) a la
curva de nivel que pasa por (x0 , y0 ).
Demostracion. Suponga que la curva de nivel (t) se parametriza mediante las funciones
diferenciables
entonces por la regla de la cadena, la derivada de f ((t)) = f (g(t), h(t)) = c con respecto a
t esta dada por
f dg f dh f f dg dh
+ =0 , , = 0, f (t) = 0
x dt y dt x y dt dt
Hasta ahora las superficies en el espacio se han representado principalmente por medio
de ecuaciones de la forma
z = f (x, y), Ecuacion de una superficie S
Sin embargo, en el desarrollo que sigue, es conveniente utilizar la representacion mas general
F (x, y, z) Una superficie S dada por z = f (x, y) se puede convertir a la forma general
definiendo F como
F (x, y, z) = f (x, y) z
Puesto que f (x, y) z = 0, se puede considerar S como la superficie de nivel de F dada por
F (x, y, z) = 0 Ecuacion alternativa de la superficie S
Teorema 23.
Si F (x, y, z) = c es la superficie de nivel de una funcion diferenciable de tres variables
w = F (x, y, z) que pasa por un punto especificado (x0 , y0 , z0 ) entonces F (x0 , y0 , z0 ) es
normal (ortogonal) a la superficie de nivel que pasa por (x0 , y0 , z0 ).
Demostracion. Consideremos una curva C(t) derivable en esta superficie, con ecuaciones
parametricas x = x(t), y = y(t), y z = z(t) para las cuales C(t0 ) = (x0 , y0 , z0 ) y C (t) 6= 0,
entonces por la regla de la cadena, la derivada de F (C(t)) = F (x(t), y(t), z(t)) = c con
respecto a t esta dada por
F dx F dy F dz F F F dx dy dz
+ + =0 , , , , = 0,
x dt y dt z dt x y z dt dt dt
En particular, en t = t0 ,
F (x0 , y0 , z0 ) C (t0 ) = 0.
Puesto que este argumento se cumple para cualquier curva C diferenciable que pasa por
(x0 , y0 , z0 ) sobre la superficie, concluimos que: f es normal a la superficie de nivel en
(x0 , y0 , z0 ).
64
NOTA: Todas las rectas tangentes en S se encuentran en un plano que es normal a
F (x0 , y0 , z0 ) y contiene a como se muestra en la figura
Demostracion.
Sea Q(x, y, z) un punto arbitrario en el plano tangente. Entonces el vector P Q = (x
x0 , y y0 , z z0 ) se encuentra en el plano tangente. Como f (x0 , y0 , z0 ) es normal al plano
tangente en (x0 , y0 , z0 ) debe ser ortogonal a todo vector en el plano tangente. En particular
f (x0 , y0 , z0 ) (x x0 , y y0 , z z0 ) = 0
En otra palabras
Ejemplo 86. Dibujar la curva de nivel que corresponde a c = 0 para la funcion dada por
f (x, y) = y sen x y hallar un vector normal a varios puntos de la curva.
1. f (, 0) = 3. f (0, 0) =
2. f (/2, 1) = 4. f (/2, 1) =
65
Ejemplo 87. Un rastreador termico se encuentra en el punto (2, 3) sobre una placa metali-
ca cuya temperatura en (x, y) es T (x, y) = 20 4x2 y 2 . Hallar la trayectoria del rastreador,
si este se mueve continuamente en direccion de maximo incremento de temperatura.
Solucion. Denotamos por r(t) = (x(t), y(t)) la posicion del rastreador. Como el rastreador
busca el maximo incremento de temperatura, este se encuentra en la direccion del f (x, y).
Por tanto, las direcciones de r (t) = y f (x, y) =
son iguales en todo punto de la trayectoria.
2
Solucion. La elipse es una curva de nivel de la funcion f (x, y) = x4 + y 2 (para c = 2) Luego
el gradiente f (x, y) = es un vector normal a la recta tangente, luego la
ecuacion de la recta tangente en elpunto (2 1) (dada como hiperplano o ecuacion general)
es:
2 2 2
Ejemplo 89. Encuentre la superficie de nivel de F (x, y, z) = x +y +z y una ecuacion del
plano tangente que pasa por (1, 1, 1). Grafique el gradiente y el plano tangente en el punto.
Solucion. Puesto que F (1, 1, 1) = 3 la superficie de nivel que pasa por (1, 1, 1) es la esfera
x2 + y 2 + z 2 = 3. El gradiente de la funcion es f = y por ello
F (1, 1, 1) =
66
Ejemplo 90. Encuentre una ecuacion de un plano tangente a la grafica del paraboloide
z = 12 x2 + 12 y 2 + 4 en (1, 1, 5). Adicionalmente encuentre las ecuaciones parametricas para
la recta normal a la superficie en este mismo punto.
(x 1) (y + 1) (z 5) = 0 x+y+z =3
Ejemplo 91. Hallar una ecuacion del plano tangente al hiperboloide z 2 2x2 2y 2 = 12 en
el punto (1, 1, 4)
F (x, y, z) = F (1, 1, 4) =
Ejemplo 92. Hallar un conjunto de ecuaciones simetricas para la recta normal a la superficie
dada por xyz = 12 en el punto (2, 2, 3)
67
Recta tangente en interseccion de superficies
Si F y G son funciones de tres variables derivables continuamente, entonces la interseccion
de las superficies
por lo general sera alguna clase de curva C en el espacio, pues las ecuaciones en (2.14)
pueden resolverse de manera implicita para dos de las variables en terminos de la tercera.
(2 ecuaciones 3 incognitas). En cualquier caso, C es una curva suave que pasa por P . Como
esta curva se localiza en ambas superficies, su vector tangente en P es perpendicular a ambos
vectores normales F (P ) y G(P ). Si estos vectores NOOOO son colineales, entonces el
vector
T = F (P ) G(P )
es tangente a P en la curva C de la interseccion de las dos superficies F (x, y, z) = 0 y
G(x, y, z) = 0
Ejemplo 93. Describir la recta tangente a la curva de interseccion de las superficies x2 +
2y 2 + 2z 2 = 20 y x2 + y 2 + z = 4 en el punto (0, 1, 3)
Solucion.
Los gradientes de ambas superficies en el punto (0, 1, 3).
Por tanto, la recta tangente a la curva de interseccion de las dos superficies en el punto
(0, 1, 3) es
Elipsoide Paraboloide
F (x, y, z) = G(x, y, z) =
F (1, 1, 2) = G(1, 1, 2) =
68
El producto vectorial de estos dos gradientes es un vector tangente a ambas superficies en
el punto (1, 1, 2).
i j k
F (1, 1, 2) G(1, 1, 2) =
Por tanto, la recta tangente a la curva de interseccion de las dos superficies en el punto
(1, 1, 2) es
Perspectiva general
r(h)
f (x0 + h) = f (x0 ) + Df (x0 )(h) + r(h), donde lm =0
h0 khk
NOTA: Siendo r una funcion que toma valores en Rm , la interpretacion del lmite anterior
sera de que todas sus funciones componentes tenderan a cero cuando, al dividirlas por khk,
h tiende a cero.
En esta definicion aparece ya el termino transformacion lineal de Rn en Rm , que
se estudia en los cursos de Algebra Lineal. Sabemos que las transformaciones lineales tienen
representaciones concretas por medio de matrices. De hecho, la transformacion lineal T :
Rn Rm tiene asociada una matriz de orden m n que la representa, la cual esta construida
de la forma siguiente:
h i
AT = [T (e1 )]B m [T (e )]
2 B Rm [T (e )]
n B Rm
R
69
Sea f : U Rn Rm una funcion diferenciable en x0 U (Aqui, f = (f1 , f2 , . . . , fm ).
Luego existe una transformacion lineal Df (x0 ) : Rn Rm , tal que
r(h)
f (x0 + h) = f (x0 ) + Df (x0 )(h) + r(h), donde lm =0
h0 khk
Nuestro objetivo aqu es caracterizar la matriz ADf (x0 ) asociada a la transformacion
Df (x0 ),
es decir, vamos a calcular la columna j-esima de ADf (x0 ) (esto es, Df (x0 )(ej ) B m ).
R
Por tanto
f1
(x0 )
xj
f2
(x0 )
Df (x0 )(ej )
=
xj
B Rm ..
.
fm
, (x0 )
xj
f1 f1 f1
(x ) (x0 ) , (x0 )
x1 0 x2 xn
f2 f2 f2
x (x0 )
(x0 ) , (x0 )
1 x2 xn
= = Jf (x0 )
.. .. ..
. . .
f fm fm
m
(x0 ) (x0 ) , (x0 )
x1 x2 xn
A esta matriz m n, se le llama matriz Jacobianade la funcion f en x0 y se denota Jf (x0 ).
Esta es entonces la Derivada de la funcion diferenciable f en x0 (identificando, como siempre,
a la transformacion lineal Df (x0 ) : Rn Rm con la matriz que la representa Jf (x0 )).
Veamos, por ejemplo,
70
!
n n f f f
m = 1. f :U R R Df (x0 ) : R R Jf (x0 ) = (x0 ), (x0 ), , (x0 ) := f (x0 )
x1 x2 xn
f1 f1 f1
x1 (x0 ) (x0 ) , (x0 )
x2 xn
m = 2. f : U Rn R2 Df (x0 ) : Rn R2 Jf (x0 ) =
f2 f2 f2
(x0 ) (x0 ) , (x0 )
x1 x2 xn
Ejemplo 95. Calcule la matriz Jacobiana de la funcion f : R2 R2 dada por f (x, y) =
(x2 + 3y 2 , 5x3 + 2y 6 )
Solucion.
Ejemplo 96. Calcule la matriz Jacobiana en el punto (0, 0) de la funcion f : R2 R3 dada
por f (x, y) = (sin(x + y), xex+y , x + y)
Solucion.
71
o bien, en terminos matriciales
J(f g)(x0 ) = Jf g(x0 ) Jg(x0 )
el lado derecho de esta ultima expresion es una multiplicacion de matrices. Observese que
Jg(x0 ) es una matriz n mn y Jf (g(x0 )) una matriz p n, de modo que el producto
Jf g(x0 ) Jg(x0 ) esta bien determinado y da por resultado una matriz p m, como debe ser
la derivada de J(f g)(x0 )
Solucion.
72
Maximos y mnimos de z = f (x, y)
Primero consideraremos una funcion z = f (x, y) continua de dos variables. Suponga que
nos interesan los valores mas grande y mas pequeno que alcanza f (x, y) sobre una region
plana, acotada y cerrada R. Recuerdese que una region en el plano es cerrada si contiene
todos sus puntos frontera. Mientras que una region en el plano se le llama acotada si es
una subregion de un disco cerrado en el plano.
Definition 27. Sea z = f (x, y) una funcion definida en una region plana R
1. La funcion f tiene un mnimo local en (x0 , y0 ) si existe, una bola abierta B (x0 , y0 )
R, ( > 0), tal que
f (x0 , y0 ) f (x, y), (x, y) B (x0 , y0 )
2. La funcion f tiene un maximo local en (x0 , y0 ) si existe, una bola abierta B (x0 , y0 )
R, ( > 0), tal que
f (x, y) f (x0 , y0 ) (x, y) B (x0 , y0 )
1. Decir que f tiene un maximo relativo en (x0 , y0 ) significa que el punto (x0 , y0 , f (x0 , y0 ))
es por lo menos tan alto como todos los puntos cercanos en la grafica de z = f (x, y).
Como en el caso de las funciones de una sola variable, la clave para identificar los extremos
locales es un criterio de la primera derivada. Recordemos que para determinar los valores
extremos locales de una funcion f (x) de una variable, buscamos los puntos donde la grafica
y = f (x) tiene una recta tangente horizontal . En tales puntos, buscamos los maximos
locales, los mnimos locales y los puntos de inflexion.
Para una funcion f (x, y) de dos variables, buscamos los puntos donde la superficie z =
f (x, y) tiene un plano tangente horizontal . En tales puntos, buscamos los maximos locales,
los mnimos locales y los puntos silla (de estos ultimos daremos mas detalles en breve).
fx (a, b) = 0 fy (a, b) = 0
73
Demostracion. Si f tiene un extremo local en (a, b), entonces la funcion g(x) = f (x, b)
tiene un extremo local en x = a . Por tanto, g (a) = 0. Ahora, g (a) = fx (a, b) de modo que
fx (a, b) = 0 Un argumento similar con la funcion h(y) = f (a, y) muestra que fy (a, b) = 0.
z = f (a, b) = constante.
Por tanto, el teorema anterior nos dice que la superficie tiene un plano tangente horizontal
en un extremo local, si dicho plano existe.
Puntos crticos
Recuerde que en calculo I definimos un valor crtico c de una funcion f de una sola
variable x como un numero en su dominio para el cual f (c) = 0 o f (c) no existe. En la
definicion que sigue definimos un punto crtico de una funcion f de dos variables x y y.
Definition 29 (Puntos Crticos).
Un punto crtico de una funcion z = f (x, y) es un punto (a, b) en el dominio de f para el
cual fx (a, b) = 0 y fy (a, b) = 0, o si una de sus derivadas parciales no existe en el punto.
1. Los puntos crticos corresponden a puntos donde f podra posiblemente tener un ex-
tremo relativo.
2. En algunos libros los puntos crticos tambien reciben el nombre de puntos estacio-
narios.
3. En el caso en que las primeras derivadas parciales EXISTAN, notamos que un punto
crtico se encuentra al resolver simultaneamente las ecuaciones
fx (x, y) = 0 fy (x, y) = 0
4. El Teorema 28 dice que los unicos puntos donde una funcion f (x, y) puede asumir
valores extremos son los puntos crticos y los puntos frontera.
5. Al igual que en las funciones diferenciables de una sola variable, no todo punto crtico
da lugar a un extremo local.
6. Una funcion diferenciable de una variable podra tener un punto de inflexion. Una
funcion diferenciable de dos variables puede tener un punto silla.
74
Definition 30 (Punto Silla).
Un punto crtico (a, b) de una funcion z = f (x, y) es un punto silla si para toda bola B (a, b)
existen puntos (x, y) B tal que f (x, y) > f (a, b) y puntos (x, y) B tal que f (x, y) <
f (a, b).
Solucion. El dominio de f es todo el plano R2 (no hay puntos frontera). Ahora hallemos los
puntos crticos de f estudiando fx y fy ,
fx (x, y) = fy (x, y) =
estan definidas en todo R2 . Por tanto, los puntos extremos locales pueden ocurrir solamente
cuando , es decir,
Solucion. El dominio de f es todo el plano R2 (no hay puntos frontera). Ahora hallemos los
puntos crticos de f estudiando las fx y fy
fx (x, y) = fy (x, y)
estan definidas en todo R2 salvo en (0, 0). Observe que aqu las derivadas fx = fy = 0 no
pueden ser ambas 0. As que, buscamos el punto crtico en los (x, y) donde fx o fy NO
existan. Es por esto que (0, 0) es el unico punto crtico de f . No es dificil ver que f (0, 0) = 1
y que
f (x, y) = 1 (x2 + y 2 )1/3 < 1 = f (0, 0)
Por tanto, f tiene un maximo relativo en (0, 0)
Nota: Recuerdese que una de las derivadas parciales debe no existir o las dos deben ser 0
para tener un punto crtico.
75
Ejemplo 101. Hallar los extremos locales (si existen) de f (x, y) = y 2 x2
Solucion. El dominio de f es todo el plano R2 (no hay puntos frontera). Ahora hallemos los
puntos crticos de f estudiando las fx y fy ,
fx (x, y) = fy (x, y) =
claramente existen en todas partes. Por tanto, los puntos extremos locales pueden ocurrir
solamente cuando , es decir, solo cuando x = yy=
Luego el punto crtico es . y su valor crtico es . Ahora observe
que,
f (x, 0) = x2 < 0, f (0, y) = y 2 > 0
Es decir, f (x, y) es positiva y es negativa en puntos cercanos de , por eso
concluimos que f tiene un punto silla en . Observe que f la funcion NO tiene
valores extremos locales.
Ejemplo 102. Hallar los puntos crticos de f (x, y) = xy(3 x y) definida en el conjunto
D = (x, y) R2 : x > 0, y > 0, x + y < 3
Solucion. El dominio de f es la region de la figura, (no hay puntos frontera). Ahora hallemos
los puntos crticos de f estudiando fx y fy ,
fx (x, y) = fy (x, y) =
estan definidas en todo R2 . Por tanto, los puntos extremos locales pueden ocurrir solamente
cuando , es decir,
En las funciones de los ejemplos anteriores, fue relativamente facil determinar los extremos
locales, porque cada una de las funciones estaba dada, o se poda expresar, en forma de
cuadrado perfecto. Con funciones mas complicadas, los argumentos algebraicos son menos
adecuados y es mejor emplear los medios analticos presentados en el siguiente criterio de las
segundas derivadas parciales. Es el analogo, para funciones de dos variables, del criterio
de las segundas derivadas para las funciones de una variable.
Pero antes, definamos la matriz hessiana de una funcion de varias variables.
Definition 31 (Matriz Hessiana).
Sea f : D Rn R una funcion de n variables con dominio el conjunto abierto D tal que
fxi (p) y fxi xj (p) existen i, i = 1, 2, . . . , n y p D Se denomina matriz hessiana de la
funcion f en el punto p D, a la matriz dada por
fx1 x1 (p) fx1 x2 (p) fx1 xn (p)
fx2 x1 (p) fx2 x2 (p) fx2 xn (p)
Hess(f (p)) =
.. .. .. ..
. . . .
fxn x1 (p) fxn x2 (p) fxn xn (p)
76
fx (x, y) = fxx (x, y) = fxy (x, y) =
Hess(g(x, y, z)) =
OBS: Las menores principales de la matriz hessiana Hess(f (p)) son los determinantes
dados por
f
x1 x1 (p) fx1 x2 (p) fx1 x3 (p)
fx1 x1 (p) fx1 x2 (p)
H1 (f (p)) = fx1 x1 (p) H2 (f (p)) =
,
H3 (f (p)) = fx2 x1 (p) fx2 x2 (p) fx2 x3 (p)
fx2 x1 (p) fx2 x2 (p)
fx3 x1 (p) fx3 x2 (p) fx3 x3 (p)
fx x (p)
1 1 fx1 x2 (p) fx1 xn (p)
fx2 x1 (p) fx2 x2 (p) fx2 xn (p)
Hess(f (p)) =
.. .. .. ..
. . . .
fx x (p)
n 1
fxn x2 (p) fx x (p)
n n
77
Teorema 32. Sea f : D Rn R una funcion de n variables con dominio el conjunto
abierto D Rn , tal que fxi (p) y fxi xj (p) existen i, j = 1, 2, . . . , n y p D. Sea x0 D
un punto crtico de la grafica de f , entonces se tiene:
1. Si la matriz hessiana Hess(f (x0 )) es definida positiva, esto es, las menores
(es decir, los menores principales de orden impar son negativas y los de orden par
positivos), entonces f (x0 ) es un valor maximo local de f .
entonces (1, 3) corresponde a un mnimo local. El valor mnimo local es f (1, 3)) = 14
78
Para el punto (5, 1), se tiene
12 0 H (f (5, 1)) = 12 < 0
1
Hess(f (5, 1)) = , y como ,
0 12 H2 (f (5, 1)) = 144 > 0
entonces (5, 1) corresponde a un maximo local. El valor maximo local es f (5, 1)) = 50
Para el punto (5, 3), se tiene
12 0 H (f (5, 3)) = 12 < 0
1
Hess(f (5, 3)) = , y como ,
0 12 H2 (f (1, 3)) = 144 < 0
OBSERVACION:
1. Con el criterio de las segundas derivadas parciales pueden no hallarse los extremos
locales si alguna de las primeras derivadas parciales no existe, no se puede aplicar el
criterio.
2. Cuado el metodo no funciona se pueden tratar de hallar los extremos mediante la grafica
o mediante algun otro metodo.
79
Al igualar a cero estas derivadas parciales, se obtiene los puntos crticos
Ejemplo 107. Supongase que x e y representan las cantidades (en kilos) de los complementos
diferentes mezclados en un alimento balanceado para pollos. Los pesos resultantes (en kilos)
para vender los gallos y las gallinas se estiman por p = 6x 3y y q = 3y 9 respectivamente.
La utilidad (en miles pesos) obtenido de un lote de pollos se modela por la funcion
(se supone que el numero de gallos y gallinas en cada lote no vara). Determine las cantidades
x e y de cada complemento alimentario que maximiza la utilidad.
80
Solucion. Al aplicar la regla de la cadena, las derivadas parciales de la funcion utilidad con
respecto a x e y son
U
=
x
U
=
y
Al igualar a cero estas derivadas parciales, se obtiene un unico punto crtico La
matriz hessiana de f en un punto generico (x, y) es
fxx (x, y) fxy (x, y)
Hess(f (p)) = =
fyx (x, y) fyy (x, y)
Solucion.
Sean a, b y c las longitudes de las aristas de la caja, tal que P (a, b, c) es el vertice (opuesto al
origen) de la caja situado en el plano Q (ver figura). Como P (a, b, c) Q, entonces Luego,
V = abc
Va (a, b) =
Vb (a, b) =
Al igualar a cero estas derivadas parciales, se obtiene un unico punto crtico
La matriz hessiana de f en un punto generico (x, y) es
Vaa (a, b) Vab (a, b)
Hess(V (p)) = =
Vba (a, b) Vbb (a, b)
81
Extremos en conjuntos acotados cerrados
Recuerde que si una funcion f de una variable x es continua en un intervalo cerrado [a, b]
entonces f posee siempre un maximo global o absoluto y un mnimo global o absoluto en el
intervalo [a, b]. Tambien vimos que estos extremos globales sobre [a, b] ocurren en un punto
frontera x = a o x = b del o en un numero crtico c (a, b).
A continuacion se presenta el teorema del valor extremo para una funcion f de dos
variables x y y que es continua sobre un conjunto R cerrado y acotado en el plano xy.
Teorema 33. Una funcion f de dos variables x y y que es continua sobre un conjunto R
cerrado y acotado tiene siempre un maximo absoluto y un mnimo absoluto sobre R.
OBS: Los valores maximo o mnimo globales o absolutos de f (x, y) sobre un conjunto
R cerrado y acotado se encuentran como sigue:
Ejemplo 109. Determine los valores maximos y mnimos globales o absolutos de f (x, y) =
2 + 2x + 2y x2 y 2 en la region triangular del primer cuadrante acotada por las rectas x = 0
y y = 0, y = 9 x
Solucion. Como f es diferenciable, los unicos lugares donde f puede asumir estos valores son
los puntos interiores del triangulo donde fx = fy = 0 y los puntos en la frontera.
fx (x, y) = 2 2x fy (x, y) = 2 2y
Al igualar a cero estas derivadas parciales, se obtiene un unico punto crtico (1, 1) y el
valor de f ah es f (1, 1) = 4.
f (1, 0) = g(1) = 3.
82
b) Sobre el segmento OB, x = 0. La funcion
f (x, y) = 2 + 2y y 2 := h(y)
x=0
Ejemplo 110. Determine los valores maximos y mnimos globales o absolutos de f (x, y) =
6x2 8x + 2y 2 5 en la region cerrada R definido por x2 + y 2 1,
Solucion. Como f es diferenciable, los unicos lugares donde f puede asumir estos valores son
los puntos interiores del circulo x2 + y 2 < 1 donde fx = fy = 0 y los puntos en la frontera
x2 + y 2 = 1.
Al igualar a cero estas derivadas parciales, se obtiene un unico punto crtico ( 32 , 0). Como
2
2
3 + 02 < 1, el punto esta en el interior de R y el valor de f ah es f ( 23 , 0) = 23 3 .
f (1, 0) = g() = 9.
83
frontera de g: los unicos valores frontera de g son
t=0 donde = g(0) = f (1, 0) = 7 En
t = 2 donde = g(2) = f (1, 0) = 7
Multiplicadores de Lagrange
En ocasiones es necesario determinar los valores extremos de una funcion cuyo dominio
esta restringido a cierto subconjunto particular del plano (por ejemplo, un disco, una region
triangular cerrada, o a lo largo de una curva). Tales restricciones tienden a complicar los
problemas de optimizacion porque la solucion optima puede presentarse en un punto frontera
del dominio. En esta seccion exploraremos un poderoso metodo para determinar los valores
extremos de funciones restringidas: el metodo de multiplicadores de Lagrange.
Teorema 34 (Teorema de Lagrange).
Sean f y g funciones con primeras derivadas parciales continuas, y tales que f tiene un
extremo en un punto (x0 , y0 ) sobre la curva suave de restriccion g(x, y) = 0 Si g(x0 , y0 ) 6= 0
entonces existe un numero real tal que
f (x0 , y0 ) = g(x0 , y0 ).
dh f dx f dy
= + = f (x0 , y0 ) C (t0 ) = 0
dt x dt y dt
As, f (x0 , y0 ) es ortogonal a C (t0 ). Pero sabemos tambien que g, tambien es ortogonal
C (t). Por consiguiente, los gradientes f (x0 , y0 ) y g(x0 , y0 ) y son paralelos y debe existir
un escalar tal que
f (x0 , y0 ) = g(x0 , y0 ).
OBSERVACION:
1. Una interpretacion geometrica interesante del teorema 35. En las figuras se observa la
84
curva de restriccion g(x, y) = 0 junto con curvas de nivel comunes de la funcion f (x, y).
Como los vectores gradiente f y g son normales a las curvas de nivel de las funciones
f y g, respectivamente, se sigue que las curvas f (x, y) = c y g(x, y) = 0 son tangentes
una con otra en puntos donde los dos vectores gradiente son colineales y f toma su
valor maximo (o mnimo), c.
2. El criterio de los multiplicadores de Lagrange sirve para seleccionar, de entre las curvas
de nivel de f , aquella que es tangente a la curva de restriccion.
3. El metodo de los multiplicadores de Lagrange NO tiene un indicador integrado que
senale MAX o MIN cuando se encuentra un extremo.
4. El numero real para el cual f = g recibe el nombre de multiplicador de
Lagrange.
5. Al aplicar el metodo de los multiplicadores de Lagrange, en realidad no estamos intere-
sados en determinar los valores de que satisfacen el sistema f = g.
6. La ecuacion f = g es equivalente a
fx (x, y) = gx (x, y) fy (x, y) = gy (x, y)
Entre las soluciones (x, y, ) del sistema (2.15) estaran los puntos donde f tiene un
extremo. Cuando f tiene un maximo (mnimo), este sera el numero mas grande (o mas
pequeno) en la lista de los valores de la funcion f (xi , yi ).
Ejemplo 111. Encuentre los puntos de la hiperbola rectangular xy = 1 que sean los mas
cercanos al origen (0, 0).
p
Solucion. Se necesita minimizar la distancia d = x2 + y 2 que hay del origen a un punto
P (x, y) sobre la curva xy = 1. Pero el algebra es mas sencilla si en lugar de eso se minimiza
el cuadrado de esta distancia, es decir,
f (x, y) = x2 + y 2 sujeta a la restriccion g(x, y) = xy 1
Luego, f (x, y) = (2x, 2y) y g(x, y) = (y, x). Por tanto, la ecuacion f = g nos lleva
a
(i) 2x = x 2x2 = xy
(ii) 2y = y 2y 2 = xy
(iii) xy 1 = 0
Pero como xy = 1 > 0 implica que x y y tienen el mismo signo. Entonces, el hecho de que
x2 = y 2 implica que x = y. Al sustituir xy = 1 se obtiene x2 = 1, por lo que se sigue
finalmente que x = y = 1 o x = y = 1. Los dos resultados posibles (1, 1) y (1, 1) estan
indicados en la figura.
Fig. 1
85
Ejemplo 112. Determine los extremos f (x, y) = y 2 4x sujetos a x2 + y 2 = 9
Aqu la restriccion esta dada por g(x, y) = x2 + y 2 9 = 0. Luego, f (x, y) = (4, 2y) y
g(x, y) = (2x, 2y). Por tanto, la ecuacion f = g nos lleva a
(i) 4 = 2x
(ii) 2y = 2y y(1 ) = 0
2 2
(iii) x +y 9=0
Si y = 0, de (iii) tenemos x2 = 9 o x = 3. Por consiguiente, (3, 0) y (3, 0) son soluciones
del sistema.
2
Si =1, entonces de (i) tenemos x = 2. Al sustituir este valor en
(iii) obtenemos
y =5
o y = 5. luego tenemos dos o soluciones mas del sistema (2, 5) y (2, 5). De la
lista de valores de la funcion
f (3, 0) = 12 f (3, 0) = 12, f (2, 5) = 13, f (2, 5) = 13
De manera
que f tiene un mnimo con restricciones en (3, 0) y un maximo con restricciones
en (2, 5)
Los cuatro puntos que encontramos yacen en el plano xy sobre el crculo de radio 3; los
tres extremos con restricciones corresponden a los puntos (3, 0, 12), (2, 5, 13) en el
espacio tridimensional sobre la curva de interseccion de la superficie del cilindro circular.
Aqu la grafica muestra tres curvas de nivel de y 2 4x = c. Dos de las curvas de nivel son
tangentes al crculo x2 + y 2 = 9.
Ejemplo 113. Determine los extremos f (x, y, z) = x2 + y 2 + z 2 sujetos a 2x + 2y z = 5
86
Solucion. El tronco esta limitado por la elipse (x/2)2 +y 2 = 1 es decir, x2 +4y 2 = 4. Queremos
maximizar el area de la seccion transversal
A = f (x, y) =
de la viga sujeta a la restriccion . Luego, f (x, y) =
y g(x, y) = . Por tanto, la ecuacion f = g nos lleva a
Esta claro que ni x 6= 0 ni y 6= 0 dan el area maxima, por lo que podemos resuelver estas
dos ecuaciones de la siguiente manera
Ejemplo 115. Determine los puntos mas cercanos al origen sobre el cilindro hiperbolico
x2 z 2 1 = 0
87
Multiplicadores de Lagrange con dos restricciones
Muchos problemas nos piden determinar los valores extremos de una funcion diferenciable
f (x, y, z) cuyas variables estan sujetas a dos restricciones. Si las restricciones son
g1 (x, y, z) = 0 g2 (x, y, z) = 0
g1 (x, y, z) = 0 g2 (x, y, z) = 0
f = 1 g1 + 2 g2 , para algunas 1 y 2 .
Como los puntos que buscamos estan en ambas superficies, sus coordenadas tambien deben
satisfacer las ecuaciones g1 (x, y, z) = 0 y g(x, y, z) = 0 que son los requisitos restantes de las
ecuaciones (2.17).
Ejemplo 116. El plano x + y + z = 1 corta al cilindro x2 + y 2 = 1 en una elipse. Determine
los puntos, sobre la elipse, mas cercanos y mas lejanos del origen.
88
Luego,
f (x, y) = (2x, 2y, 2z), g1 (x, y, z) = (2x, 2y, 1), g2 (x, y, z) = (1, 1, 0).
2x = 21 x + 2 , 2y = 21 y + 2 , 2z = 2
2x = 21 x + 2z (1 1 )x = z
2y = 21 y + 2z (1 1 )y = z
x2 + x2 1 = 0 x+x+z1=0
2
x= z =1 2
2
Luego los puntos correspondientes sobre la elipse son P1 ( 22 , 22 , 1 2) y P2 ( 22 , 22 , 1+
2). Sin embargo, debemos tener CUIDADO. Aunque P1 y P2 dan maximos locales de f
sobre la elipse, P2 esta mas alejado del origen que P1 .
Los puntos sobre la elipse mas cercanos al origen son (1, 0, 0) y (0, 1, 0). El punto sobre la
elipse mas lejano del origen es P2 .
Solucion. Basandonos en las grafica, sugiere que existen dos de tales puntos P1 y P2 con
coordenadas z no negativas. La funcion f para la cual deseamos encontrar un maximo y un
mnimo es simplemente la distancia desde cada uno de estos puntos al plano xy, esto es,
f (x, y, z) = z. (distancia del plano xy a la curva intersecion) sujeta a las restricciones
89
Luego,
f (x, y) = (0, 0, 1), g1 (x, y, z) = (2x, 2y, 2z), g2 (x, y, z) = (1, 1, 3).
a) 0 = 21 x+2 , b) 0 = 21 y2 , c) 1 = 21 z+32
x2 + x2 + z 2 9 = 0 x + x + 3z 6 = 0
2 2
2x + z = 9 2x + 3z = 6 (2.20)
90
Calculo Vectorial
Integrales Multiples
Los pasos preliminares analogos que conducen al concepto de integral definida bidimensio-
nal, conocidos simplemente como integral doble de una funcion f de dos variables, se dan a
continuacion.
Sea z = f (x, y) una funcion definida en una region cerrada y acotada R del plano xy.
Considere los siguientes cuatro pasos:
91
Por medio de una retcula de lneas verticales y horizontales paralelas a los ejes de
coordenadas, forme una particion P de R en n subregiones rectangulares Rk de areas
Ak que esten por completo sobre R. Son los rectangulos que se muestran en rojo claro
en la figura.
Sea kP k la norma de la particion o la longitud de la diagonal mas grande de las n
subregiones rectangulares Rk .
Elija un punto muestra (xk , yk ) en cada subregion Rk
n
X
Forme la suma f (xk , yk )Ak
k=1
92
3. Cuando kP k 0 los rectangulos se vuelven mas angostos y mas cortos, y su
numero n aumenta, de modo que tambien podemos escribir
n
X n
X
lm f (xk , yk )Ak = lm f (xk , yk )Ak
kP k0 n
k=1 k=1
Volumen neto No toda integral doble produce volumen. Para la superficie que se muestra
en la figura, f (x, y)dA es un numero real pero no es el volumen puesto que f es no negativa
R
sobre R. Podemos interpretar la integral doble como la suma del volumen acotado entre la
grafica de f y la region R siempre que f (x, y) 0 y el negativo
del volumen entre la grafica
de f y la region R siempre que f (x, y) 0. En otras palabras, f (x, y)dA representa un
R
volumen neto entre la grafica de f y el plano xy sobre la region R.
93
Teorema 37 (Propiedades). Sean f y g funciones de dos variables que son integrables sobre
una region cerrada R del plano xy. Entonces
1. kf (x, y)dA = k f (x, y)dA
R R
2. [f (x, y) g(x, y)]dA = f (x, y)dA g(x, y)dA
R R R
3. f (x, y) = f (x, y)dA + f (x, y)dA donde R1 y R2 son subregiones que no
R R1 R2
se traslapan y R = R1 R2
4. Si f (x, y) g(x, y) sobre R entonces f (x, y)dA g(x, y)dA
R R
Ejemplo 118. Considere la region de integracion R en el primer cuadrante acotado por las
graficas de x + y = 2, y = 0 y x = 0. Aproxime la integral doble (5 x 2y)dA utilizando
R
una suma de Riemann, las Rk que se muestran en la Figura y los puntos muestra (xk , yk )
en el centro geometrico de cada Rk .
94
En otras palabras, para evaluar la integral parcial f (x, y)dy mantenemos x fija (como
si fuera una constante), en tanto que en f (x, y)dx mantenemos y fija.
Al evaluar una integral definida podemos prescindir de las funciones C1 (x) y C(y) en (3.2)
y (3.3).
Definition 39. Integracion parcial definida:
Si F (x, y) es una funcion tal que Fy = f (x, y) entonces la integral parcial definida con
respecto a y se define como
g2 (x) g2 (x) g2 (x)
f (x, y)dy = Fy (x, y)dy = F (x, y) = F (x, g1 (x) F (x, g2 (x))) (3.4)
g1 (x) g1 (x) g1 (x)
Si F (x, y) es una funcion tal que Fy = f (x, y) entonces la integral parcial definida con
respecto a x se define como
h2 (y) h2 (y) h2 (x)
f (x, y)dx = Fx (x, y)dx = F (x, y) = F (h1 (y), y) F (h2 (y), y)) (3.5)
h1 (y) h1 (y) h1 (x)
Las funciones g1 (x) y g2 (x) en (3.4) y las funciones h1 (y) y h2 (y) en (3.5) se denominan
los lmites de integracion. Desde luego los resultados en (3.4) y (3.5) se cumplen
cuando los lmites de integracion son constantes.
donde las funciones frontera g1 y g2 son continuas, se denomina region tipo I . Mientras que
la region
R:cyd h1 (y) x h2 (y)
donde h1 y h2 son continuas, se denomina region tipo II
95
Definition 40. Integral iterada:
Si f es continua sobre una region R de tipo I, definimos una integral iterada de f sobre
la region mediante
b g2 (x) " #
b g2 (x)
f (x, y)dydx = f (x, y)dy dx (3.6)
a g1 (x) a g1 (x)
Si f es continua sobre una region R de tipo II, definimos una integral iterada de f sobre
la region mediante
d h2 (y) " #
d h2 (y)
f (x, y)dxdy = f (x, y)dx dy = (3.7)
c h1 (y) c h1 (y)
b g2 (x)
1. Sobre una region de tipo I la integral iterada f (x, y)dydx es primero una
a g1 (x)
sumatoria en la direccion de y.
3. El dy situado antes del dx significa que los volumenesf (x, y)dydx de los paraleleppe-
dos construidos sobre los rectangulos se suman verticalmente con respecto a y desde la
curva frontera inferior y = g1 (x) hasta la curva frontera superior y = g2 (x).
Ejemplo 119. Evalue la integral iterada de f (x, y) = 2xy sobre la region cerrada por y =
x2 + 1, y = x x = 1 y x = 2.
y
Ejemplo 120. Evalue la integral iterada de f (x, y) = 8x + e sobre la region cerrada por
x = 2y, y = x x = 0 y x = 4.
96
Solucion. La region es de tipo la integral iterada es mas facil si usamos
8x + ey dA =
Ejemplo 121. Evalue la integral doble ex+3y sobre la region acotada por las graficas de
R
y = 1, y = 2 y = x y y = x + 5.
Solucion. Al dibujar la region dada concluimos que la region es de tipo II la integral iterada
es mas facil si usamos dxdy
ex+3y dA =
Ejemplo 122. Emplee la integral doble para determinar el area de la region acotada por las
graficas de y = x2 y y = 8 x2
Solucion. Primero que todo f (x, y) = 1. Ahora al dibujar la region dada concluimos que la
region es de tipo I, luego la doble integral usaremos
A= dA =
Ejemplo 123. Emplee la integral doble para determinar el volumen V del solido en el primer
octante que es que esta acotado por los planos de coordenadas y las graficas del plano z =
3 x y y el cilindro x2 + y 2 = 1.
97
Solucion. De las figuras vemos que el volumen esta dado por V = zdA y la region de
R
integracion R es tipo I, luego,
1 1x2
V = zdA = (3 x y)
0 0
Ejemplo 124. Evalue (x + y)dA sobre la region acotada por las graficas de x = y 2 y
R
y = 12 x 3
2
S
Solucion. Al dibujar la region vemos que puede escribirse como la union R = R1 R2 de
las dos regiones tipo I. o interpretarla de tipo II, calculemos la doble integral usando la
region tipo I y tipo II
I. (x + y)dA = (x + y)dA + (x + y)dA
R R1 R2
II. (x + y)dA =
R
Pregunta: Porque la respuesta en este ejemplo 4 NOOOO representa el volumen del
solido sobre R y debajo del plano z = x + y?
2
Ejemplo 125. Evalue xey dA sobre la region R en el primer cuadrante acotado por las
R
graficas de y = x2 , x = 0, y = 4.
98
Solucion: Al dibujar la region vemos que puede considerarse de tipo I o II
I. (x + y)dA =
R
II. (x + y)dA =
R
Ejemplo 126. Hallar el volumen de la region solida R acotada superiormente por el para-
boloide z = 1 x2 y 2 e inferiormente por el plano z = 1 y ver figura
El Volumen
es la diferencia entre el volumen bajo el parabolide y el volumen bajo el plano,
V = zsup zinf dA =
R
Ejemplo 127. Encuentre el volumen V del solido T limitado por los planos z = 6 y z = 2y,
y por los cilindros parabolicos y = x2 y y = 2 x2 . Este solido se ilustra en la figura
99
cuales se intersecan en los puntos (1, 1) y (1, 1).
V = zsup zinf dAdA =
R
Un cambio de variables para esta integral esta determinado por una transformacion T deri-
vable continuamente, del plano uv al plano xy -es decir, una funcion T que asocia al punto
(u, v) un punto (x, y) = T (u, v) dado por ecuaciones de la forma
El punto (x, y) se denomina la imagen del punto (u, v) con la transformacion T . Se dice
que una transformacion T es uno a uno si cada punto en R es la imagen bajo T del punto
unico en S. Dicho de otro modo, ningun par de puntos en S tiene la misma imagen en R.
En este caso T tiene una transformacion inversa T 1 tal que
T 1 (x0 , y0 ) = (u0 , v0 )
100
T convierte al
rectangulo S en la
figura curvilnea R
T (u, v) = (x, y)
x = g(u, v)
y = h(u, v)
Sea f (x, y) continua sobre R y sea {S1 , S2 , . . . , Sn } una particion interior de S en rectangu-
los cada uno de los cuales tiene dimensiones u por v. Cada rectangulo Si es transformado
por T en una figura curvilnea Ri en el plano xy. De manera que {R1 , R2 , . . . , Rn } constitu-
yen entonces una particion interna de la region R (aunque en figuras curvilneas en lugar de
rectangulos).
Sea (ui , vi ) el punto en la esquina inferior izquierda de Si , y sea
(xi , yi ) = T (ui , vi ) = g(ui , vi ), h(ui , vi )
101
Los vectores velocidad de estas curvas son
u (t) = Tu (t, vi ) = gu (t, vi ), hu (t, vi ) v (t) = Tu (ui , t) = gv (ui , t), hv (ui , t)
Podemos aproximar la imagen de la region R = T (S) por medio del paralelogramo formado
por los vectores secante,
El area del paralelogramo curvilneo Ri es entonces aproximadamente igual al area del para-
lelogramo Pi generado por los vectores a y b, es decir , Ai area(Pi ). Esta aproximacion
sera tanto mejor en cuanto u y v sean pequenos. Entonces
pero,
i j k
gu hu g gv
k = u
Tu Tv = gu hv 0 = k
gv hv hu hv
gv hv 0
El determinante de 2 2 en el lado derecho de la ultima ecuacion se denomina Jacobiano de
2
la transformacion T : Ruv R2xy , en honor del matematico aleman Carl Jacobi (1804-1851),
quien fue el primero que investigo los cambios generales de variables en las integrales dobles.
2
Definition 41. El Jacobiano de la transformacion derivable continuamente T : Ruv R2xy
2
es la funcion (de variable real) JT : Ruv R definida por
xu (u, v) xv (u, v)
J(u, v) = := (x, y)
yu (u, v) yv (u, v)
(u, v)
Por tanto, cuando planteamos sumas de Riemann para aproximar integrales dobles, encon-
tramos
n
X
f (x, y)dxdy f (xi , yi )Ai
R i=1
n
X
f g(ui , vi ), h(ui , vi ) |JT (ui , vi )|uv
i=1
f g(u, v), h(u, v)|JT (u, v)|dudv.
S
102
Este analisis es, en realidad, un bosquejo de una demostracion del siguiente teorema
general del cambio de variables. Para asegurar la existencia de las integrales dobles
indicadas, aceptamos que las fronteras de ambas regiones R y S consisten en un numero
finito de elementos curvos suaves.
2
Teorema 42. Suponga que la transformacion derivable continuamente T : Ruv R2xy toma
la region limitada S en el plano uv en la region limitada R en el plano xy, y es uno a uno
del interior de S al interior de R. Si F (x, y) es continua sobre R, entonces
f (x, y)dxdy = f T (u, v) |JT (u, v)|dudv
R S
OBSERVACION: Dada una integral doble particular f (x, y)dxdy, como encontra-
R
mos un cambio provechoso de variables? . Un enfoque estandar consiste en escoger una
transformacion T tal que la frontera de R consiste en curvas u y curvas v. En caso de que sea
mas conveniente expresar u y v en terminos de x y y, primero se calcula en forma explcita
J(u, v) y despues se encuentra el Jacobiano necesario J(x, y) con la formula
J(u, v)J(x, y) = 1
Esta ecuacion es una consecuencia de aplicar la regla de la cadena a T (T 1 (x, y)) = (x, y).
Ejemplo 128. Halle (x2 +y 2 )dxdy donde R es la region plana limitada por las hiperbolas
R
xy = 1, xy = 3, x2 y 2 , x2 y 2 = 4.
T (u, v) = (x(u, v), y(u, v))
u = xy
v = x2 + y 2
Sin embargo podemos calcular sin problema el Jacobiano J(u, v) para lograr una buena
sustitucion,
(u, v)
J(u, v) = =
(x, y)
(x, y)
J(x, y) = =.
(u, v)
Observe que,
(x2 + y 2 )2 = (x4 + 2x2 y 2 + y 4 ) = 4x2 y 2 + (x2 y 2 )2 = 4u2 + v 2 .
Ahora estamos listos para aplicar el teorema del cambio de variables, con las regiones S y R
4 3
(x2 + y 2 )dxdy =
R 1 1
103
Ejemplo 129. Halle 3xydA donde R es la region plana limitada por las x 2y = 0,
R
x 2y = 4, x + y = 4, x + y = 1.
(x, y)
J(u, v) = =
(u, v)
4 y/2+1
2x y
Ejemplo 130. Halle dxdy.
0 y/2 2
Solucion. Trazamos la region de integracion R en el plano xy para ello observe que los limites
de x son x = y2 y x = y2 + 1 y los lmites de y son y = 0 y y = 4. Luego las rectas que
limitan a R son curvas u = x y2 y v = y (aqu podemos tomar tambien v = y2 ), Observe
que aqu es facil despejar x, y
x=u+v y y = 2v
(x, y)
J(u, v) = =
(u, v)
1 1x
Ejemplo 131. Halle x + y(y 2x)2 dydx.
0 0
104
Solucion. Trazamos la region de integracion R en el plano xy para ello observe que los limites
de x son x = 0 y x = 1 y los lmites de y son y = 0 y y = 1 x esto implica que la frontera
esta limitada por las rectas x = 0, y = 0 y x + y = 0. El integrando sugiere la transformacion
u = x + y y v = y 2x. Observe que aqu es facil despejar x, y
u v 2u v
x= y y= + (3.9)
3 3 3 3
Por tanto, el Jacobiano J(u, u) esta dado por,
(x, y)
J(u, v) = =
(u, v)
Ahora hallemos la NUEVA region donde vamos a integrar, para ello usaremos la frontera
original x = 0, y = 0 y x + y = 0 y las ecuaciones (3.9), si x = 0 entonces v = u, si y = 0
entonces v = 2u y si x + y = 1 entonces u = 1. Una vez dibujada nos damos cuenta que la
Nueva region es tipo I.
Ahora estamos listos para aplicar el teorema del cambio de variables,
1 1x
x + y(y 2x)2 dydx =
0 0
1 1x2
Ejemplo 132. Halle (x2 + y 2 )dydx.
0 0
y ahora a . Esto nos motiva a buscar una buena sustitucion para mejorar el
integrando, nos es difcil que es una excelente opcion,
x= y=
(x, y)
J(u, v) = =
(u, v)
La NUEVA region esta descrita mediante desigualdades sencillas 0 r 1 y 0 2.
Ahora estamos listos para aplicar el teorema del cambio de variables,
1 1x2
(x2 + y 2 )dydx =
0 0
NOTA: La transformacion a coordenadas polares es muy eficaz en este caso pues x2 + y 2
se se simplifica como r2 y otra razon es que los lmites de integracion se vuelven cosntantes.
2 2
Ejemplo 133. Halle ex +y dydx donde R es la region semicircular acotada por el eje x
R
105
y la curva y = 1 x2 .
Solucion. Dado que la integrando tiene el factor x2 + y 2 , entonces usaremos las coordenadas
polares
x = r cos y = r sin
Por tanto, el Jacobiano J(u, u) esta dado por,
(x, y)
J(u, v) = =
(u, v)
Ejemplo
p 134. Halle el volumen de la region solida limitada superiormente por el hemisferio
z = 16 x2 y 2 e inferiormente por la region circular R dada por x2 + y 2 4.
Solucion. No es difcil ver que el volumen del solido esta dado por zdA donde R esta
R
2 2 2 2
dado por x + y 4. Dado que la base contiene el factor x + y , entonces usaremos las
coordenadas polares
x = r cos y = r sin
Por tanto, el Jacobiano J(u, u) esta dado por,
(x, y)
J(u, v) = =
(u, v)
EJER:
yx
1. Calcule e y+x dxdy, donde D es el triangulo limitado por la recta x + y = 2 y los
D
ejes coordenados.
106
Integrales en Coordenadas polares
Como hemos visto algunas integrales dobles son mucho mas faciles de evaluar en forma polar
que en forma rectangular. Esto es as especialmente cuando se trata de regiones circula-
res, cardioides y petalos de una curva rosa, y de integrandos que contienen el termino x2 +y 2 .
donde R es el rectangulo polar. Es decir, se quiere obtener el volumen del solido cuya base es
R y que se encuentra debajo de la superficie z = f (x, y). Anteriormente se definio la integral
doble como el lmite de sumas de Riemann con particiones que consisten en rectangulos
ordinarios. Tambien es posible definir la integral doble en terminos de particiones polares,
formadas por rectangulos polares.
Para definir una integral doble de una funcion continua en coordenadas polares,
1. Sea R una region limitada o acotada por las graficas de r = g1 () y r = g2 () y las rectas
= y = .
2. En lugar de hacer una particion de R en rectangulos pequenos, se utiliza una particion en
sectores polares pequenos.
3. A R se le superpone una red o cuadrcula polar formada por rayos y arcos circulares, como
se muestra en la figura.
4. Los sectores polares Ri que se encuentran completamente dentro de R forman una particion
polar interna cuya norma es la longitud de la diagonal mas larga en los n sectores polares.
5. Considerar un sector polar especfico Ri . Se puede mostrar que el area de Ri es Ai =
ri ri i donde ri = r2 r1 y i = 2 1 .
107
Ejemplo 135. Utilizar una integral doble para hallar el area encerrada por la grafica de
r = 3 cos(3)
Solucion. No es difcil ver que la region es r-simple y los lmites de un petalo de la curva son
6 6 y 0 r 3 cos(3). Por tanto, el area de un petalo es
1
A= dA =
3 R
Ejemplo 136. Hallar el area de la region acotada superiormente por la espiral r = 3 e
inferiormente por el eje polar, entre r = 1 y r = 2
Solucion. No es difcil ver que la region es -simple y los lmites polares son 0 3r y
1 r 2. Por tanto, el area de la region es
A= dA =
R
Solucion. No es difcil ver que la region es r-simple y los lmites polares son 0 2 y
0 r 1. Por tanto, el area de la region es
A= dA =
R
108
Ejemplo 138. Encuentre el volumen de la region solida que es interior tanto a la esfera
x2 + y 2 + z 2 = 4 de radio 2, como al cilindro (x 1)2 + y 2 = 1. Halle la region R limitada en
su interior por la circunferencia r = 1 y en el exterior por la curva conocida como caracol
de Pascalr = 2 + cos .
Solucion. No es difcil ver que el volumen del solido esta dado por
p
V = (zsup zinf )dA = 2 zsup = 4 x2 y 2 dA.
R R R
donde R esta dado por (x 1)2 + y 2 1. Pero esta integral sera muy difcil de evaluar en
coordenadas rectangulares, por lo que la cambiamos a coordenadas polares. x = r cos y
y = r sin
En este caso polar la region es r-simple y los lmites polares son 0 r 2 cos con
/2 /2. Por tanto, el volumen de la region es
p
V =2 4 x2 y 2 dA =
R
Masa
Si la lamina que corresponde a la region R, que se muestra en la figura, tiene una densidad
constante entonces la masa de la lamina esta dada por
Masa = A = dA = dA Densidad Constante
R R
Si no se especifica otra cosa, se supone que una lamina tiene densidad constante. Las inte-
grales dobles pueden usarse para calcular la masa de una lamina de densidad variable, donde
la densidad en (x, y) esta dada por la funcion de densidad
Definition 43. Si es una funcion de densidad continua sobre la lamina que corresponde
a una region plana R, entonces la masa m de la lamina esta dada por
m= (x, y)dA Densidad Variable
R
NOTA: La densidad se expresa normalmente como masa por unidad de volumen. Sin embar-
go, en una lamina plana la densidad es masa por unidad de area de superficie.
Ejemplo 139. Hallar la masa de la lamina triangular con vertices (0, 0), (0, 3) y (2, 3), dado
que la densidad en (x, y) es (x, y) = 2x + y
m= dA =
R
109
Ejemplo 140. Hallar la masa de la lamina correspondiente a la porcion en el primer
cuadrante del crculo x2 + y 2 = 4 donde la densidad en el punto (x, y) es proporcional a
la distancia entre el punto y el origen.
p
Solucion. En todo punto (x, y), la densidad de la lamina es = k (x 0)2 + (y 0)2 .
Ademas, la region es de tipo . Luego la masa de la lamina es
m= dA =
R
NOTA: Notese que la lamina plana esta sombreada; el sombreado mas oscuro corresponde
a la parte mas densa.
Momento y centro de masa
En laminas de densidad variable, los momentos de masa (M omento = (masa)(distancia)) se
definen de manera similar a la empleada en el caso de densidad uniforme. Dada una particion
P de una lamina, correspondiente a una region plana R, considerar el rectangulo i-esimo
Ri de area Ai como se muestra en la figura. Suponer que la masa de Ri se concentra en
uno de sus puntos interiores (xi , yi ). El momento de masa de Ri respecto al eje x puede
aproximarse por medio de
De manera similar, el momento de masa con respecto al eje y puede aproximarse por medio
de
M om(Ri )y = (masa)xi [(xi , yi )Ai ](xi )
Formando la suma de Riemann de todos estos M omentos(Ri ) y tomando lmites cuando
kP k 0, se obtienen las definiciones de momentos de masa con respecto a los ejes x y y.
Definition 44. Sea una funcion de densidad continua sobre la lamina plana R. Los
momentos de masa con respecto a los ejes x y y son
Mx = y(x, y)dA y My = x(x, y)dA
R R
110
Ejemplo 141. Hallar el centro de masa de la lamina que corresponde a la region parabolica
0 y 4 x2 donde la densidad en el punto (x, y) es proporcional a la distancia entre
(x, y) y el eje x,
Solucion. En todo punto (x, y), la densidad de la lamina es (x, y) = ky. Ademas, la region
es de tipo . Como la lamina es simetrica con respecto al eje y, x = 0. Para hallar y
primero debemos calcular la masa de la lamina.
m= dA =
R
Mx
As y = m = . y el centro de masa es (0, ).
Ejemplo 142. Una lamina ocupa la region acotada por la recta y = x + 2 y la parabola
y = x2 . La densidad de la lamina en el punto P (x, y) es proporcional al cuadrado de la
distancia de P a partir del eje y. Encuentre la masa y centroide de la lamina.
Solucion. En todo punto (x, y), la densidad de la lamina es (x, y) = ky. Ademas, la region
es de tipo , donde sus limites estan dados por
m= (x, y)dA =
R
111
My = x(x, y)dA =
R
63k
As, la lamina tiene una masa de 20 , y su centroide se localiza en el punto ( 87 , 118
49 )
Ejemplo 143. Una lamina tiene la forma de la region R acotada por la grafica de la elipse
x2 y2
4 + 16 = 1, 0 y 4 y y = 0. Encuentre su centro de masa si la densidad es (x, y) = |x|y.
Solucion. vemos que la region es simetrica con respecto al eje y, Ademas, puesto que
(x, y) = (x, y), la densidad es simetrica alrededor del eje y. Por tanto x = 0. Pa-
ra hallar y primero hallemos la masa integrando en la region tipo II
m= (x, y)dA =
R
As, la lamina tiene una masa de 16, y su centroide se localiza en el punto (0, 32
15 )
Observacion: Para que el centro de masa este sobre un eje de simetra de una lamina, no
region debe ser simetrica y la funcion de densidad tambien debe ser simetrica con respecto a
ese eje.
Ejemplo 144. Encuentre el centro de masa de la lamina que corresponde a la region acotada
por la curva llamada petalo de rosa r = 2 sin(2) en el primer cuadrante si la densidad en el
punto P en la lamina es directamente proporcional a la distancia desde el origen polar.
Solucion. La distancia desde el origen polar es d(0, P ) = |r|. Por tanto, la densidad de la
lamina es (r, ) = k|r|.
Ahora hallemos, la masa de la lamina
m= (x, y)dA = kkr|dA =
R R
112
Ahora, el momento con respecto al eje y.
My = x(r, )dA =
R
16
As, la lamina tiene una masa de 9 k, y su centroide se localiza en el punto ( 32 32
35 , 35 )
Momentos de inercia
I = md2 = (masa)(distancia)2
113
Igual que ocurre con los momentos de masa, se puede generalizar este concepto para obtener
los momentos de inercia de una lamina de densidad variable respecto de los ejes x y y. Estos
segundos momentos se denotan por Ix e Iy y en cada caso el momento es el producto de
una masa por el cuadrado de una distancia.
Solucion. En todo punto (x, y), la densidad de la lamina es (x, y) = ky, la region es de tipo
I, y de acuerdo con la definicion de momento de inercia, se tiene
Ix = (y 2 )dA =
R
Ejemplo 146. Calcule Ix para una lamina de densidad constante = 1 que ocupa la region
limitada por las curvas x = y 4 , 1 y 1. Adicionalmente halle el area de la region
Solucion. En todo punto (x, y), la densidad de la lamina es (x, y) = 1, la region es de tipo
, y de acuerdo con la definicion de momento de inercia, se tiene
Ix = (y 2 )(x, y)dA =
R
Resistencia de columna: Se sabe que la rigidez, de una columna horizontal es proporcional
al momento de inercia de su seccion transversal respecto a un eje horizontal a traves del
centroide de la seccion transversal de la columna.
114
Supongamos que la region del ejemplo anterior es una seccion trasversal de una columna y
comparemos su resistencia con otra columna rectangular con igual altura, 2, y area similar
4
5 . Al hallar el momento de inercia de su seccion transversal
1 1/5
4
Ix = (y 2 )dxdy =
1 1/5 15
Como la razon de 47 / 15
4
= 15
7 , se observa que la columna del ejemplo es mas del doble de
fuerte que una columna rectangular con la misma superficie de seccion transversal. A esta
resistencia se debe que sea comun que las columnas del ejemplo se usen en la construccion.
Ejemplo 147. Encuentre el momento polar de inercia de una lamina circular R de radio a
y densidad constante concentrada en el origen.
Solucion. En todo punto (x, y), la densidad de la lamina es (x, y) = 1, la region de la lamina
R ocupa la region plana x2 + y 2 a2 y de acuerdo con la definicion de momento de inercia,
se tiene
I0 = (x2 + y 2 )(x, y)dA =
R
Area de superficie
En esta seccion se vera como hallar el area de la superficie superior del solido. Para
empezar, considerar una superficie S dada por z = f (x, y) definida sobre una region R.
Suponer que R es cerrada y acotada y que f tiene primeras derivadas parciales continuas.
Para hallar el area de la superficie,
115
5. El area de la superficie de S esta dada por
n
X n
X
A(S) = Si Ti
i=1 i=1
Definition 45. Sea f una funcion para la cual las primeras derivadas parciales fx y fy son
continuas sobre una region cerrada R. Entonces el area de la superficie z = f (x, y) sobre
R esta dada por
s
z z
q
A(S) = 1 + [fx (x, y)]2 + [fy (x, y)]2 dA = 1 + [ ]2 + [ ]2 dA
R R x y
1. Si fuera mas comodo proyectar la superficie S sobre el plano xz, es decir, y = f (x, z)
se usa la formula p
A(S) = 1 + [fx ]2 + [fz ]2 dA
R
donde R es la region cerrada en el plano xz.
2. Si la proyeccion de la superficie S esta sobre el plano yz, es decir, x = f (y, z) se usa la
formula q
A(S) = 1 + [fy ]2 + [fz ]2 dA
R
donde R es la region cerrada en el plano yz
116
Ejemplo 149. Determine el area de la superficie de la porcion de la esfera x2 + y 2 + z 2 = 4
que estan dentro del cilindro (x 1)2 + y 2 = 1.
Ejemplo 150. Hallar el area de la superficie de la porcion del plano z = 2 x y que se
encuentra sobre el crculo x2 + y 2 1 en el primer cuadrante.
117
q
A(S) = 1 + [fx (x, y)]2 + [fy (x, y)]2 dA =
R
Ejemplo 152. Hallar el area de la superficie del paraboloide z = 1+x2 +y 2 que se encuentra
sobre el crculo unidad.
Integrales triples
Los pasos que conducen a la definicion de la integral triple, son bastante similares a los
de la integral doble.
118
1. Sea w = f (x, y, z) definida sobre una region cerrada y acotada D del espacio tridimen-
sional.
n
X
6. La suma f (xk , yk , zk )Vk donde (xk , yk , zk ) es un punto arbitrario dentro de cada
k=1
Dk y Vk denota el volumen de cada Dk recibe el nombre de suma de Riemann.
Definition 46. Sea f una funcion de tres variables definida sobre una region cerrada D del
espacio tridimensional. Entonces la integral triple de f sobre D, se define como
n
X
= lm f (xk , yk , zk )Vk
R kP k0
k=1
Limites de integrales triples Si la region D esta acotada por arriba por la grafica de
z = g2 (x, y) y acotada por abajo por la grafica de z = g1 (x, y) entonces
h g2 (x,y) i
f (x, y, z)dV = f (x, y, z) dA
D R g1 (x,y)
En una integral doble hay solo dos posibles ordenes de integracion: dydx y dxdy. Las
integrales triples tiene seis posible formas ordenes de integracion:
Las dos ultimas diferenciales nos indican el plano de coordenadas en el cual se localiza la
region R.
119
b h2 (x) h g2 (x,y) i
dzdydx Rxy tipo I f (x, y, z)dV = f (x, y, z) dzdydx
D a h1 (x) g1 (x,y)
d h2 (y) h g2 (x,y) i
dzdxdy Rxy tipo II f (x, y, z)dV = f (x, y, z) dzdxdy
D c h1 (y) g1 (x,y)
b h2 (x) h g2 (x,z) i
dydxdz Rxz tipo I f (x, y, z)dV = f (x, y, z) dydzdx
D a h1 (x) g1 (x,z)
k h2 (z) h g2 (x,z) i
dydzdx Rxz tipo II f (x, y, z)dV = f (x, y, z) dydxdz
D e h1 (z) g1 (x,z)
d h2 (y) h g2 (y,z) i
dxdydz Ryz tipo I f (x, y, z)dV = f (x, y, z) dxdzdy
D c h1 (y) g1 (y,z)
k h2 (z) h g2 (y,z) i
dxdzdy Ryz tipo II f (x, y, z)dV = f (x, y, z) dxdydz
D e h1 (z) g1 (y,z)
APLICACIONES
VOLUMEN: Si f (x, y, z) = 1 entonces el volumen del solido D es
V = dV
D
MASA: Si (x, y, z) es la densidad, entonces la masa del solido D esta dada por
m= (x, y, z)dV
D
CENTRO DE MASA Las coordenadas del centro de masa de D estan dadas por
Myz Mxz Mxy
x= , y= , z=
m m m
SEGUNDOS MOMENTOS: Los segundos momentos, o momentos de inercia,
de D alrededor de los ejes de coordenadas indicados por los subndices estan dados por
Ix = (y 2 +z 2 )(x, y, z)dV, Iy = (x2 +z 2 )(x, y, z)dV, Iz = (x2 +y 2 )(x, y, z)dV,
D D D
Ejemplo 153. Encuentre el volumen del solido en el primer octante acotado por las graficas
de z = 1 y 2 , y = 2x y x = 3.
Solucion. No es difcil ver que la proyeccion del solido sobre el plano xy es una region tipo
II. Luego, usaremos el orden dzdxdy para calcular el volumen de D
V = dV =
D
120
Ejemplo 154. Un solido tiene la forma determinada por las graficas del cilindro |x|+|y| = 1
y los planos z = 2 y z = 4.
1. Encuentre su centro de masa si la densidad esta dada por (x, y, z) = kz, con k una
constante.
Solucion. (a) No es difcil ver que el dominio del solido esta contenido en el plano xy.
Observese que la proyeccion del solido sobre el plano xy es la region dada por |x| + |y| = 1.
Esta region es equivalente a 4 rectas,
1 1x 4
Mxy = z(x, y, z)dzdydx = 4 kz 2 dzdydz
D 0 0 2
Mxy 28
Por consiguiente, z = m = 9 . Las coordenadas del centro de masa son entonces
(0, 0, 28
9 )
(b) Sabemos que el momento de inercia
1 1x 4
Iz = z(x2 + y 2 )(x, y, z)dzdydx = 4k (x2 + y 2 )zdzdydz
D 0 0 2
121
Ejemplo 155. Cambie el orden de integracion en
6 42x/3 3x/23y/4
f (x, y, z)dzdydx
0 0 0
a dydxdz
f (x, y, z)dzdydx =
D
Ejemplo 156. Evaluar
/2 /2 3
sin(y 2 )dzdydx
0 x 1
Solucion. Una integracion directa nos lleva a 2 sin(y 2 )dydx la cual no es elemental. As que
debemos
p cambiar el orden
p de integracion. La region esta acotada por las rectas x = 0,
x = /2, y = x, y = /2, z = 1 y z = 3
/2 /2 3
sin(y 2 )dzdydx =
0 x 1
122
Las coordenadas rectangulares de un punto se obtienen de las coordenadas cilndricas
mediante las ecuaciones.
Entonces, si f (r, , z) es una funcion continua sobre la region D, la integral triple de f sobre
D esta dada por
" h2 (r,) # g2 () h2 (r,)
f (r, , z)dV = f (r, , z)dz dA = f (r, , z)rdzdrd
S R h1 (r,) g1 () h1 (r,)
Ejemplo 157. Un solido en p el primer octante tiene la forma determinada por la grafica del
cono de un solo manto z = x2 + y 2 y los planos z = 1, x = 0 y Determine el centro de
masa si la densidad esta dada por (r, , z) = r.
Solucion. No es difcil ver que el cono tiene como ecuacion z = r. (Coordenadas cilndricas).
Por tanto, la masa esta dada por
m= dV =
D
123
Mxy = zdV =
D
Mxz = ydV =
D
Myz = xdV =
D
M Mxy
De manera que x = myz = 0.38 y = Mmxz = 0.38 z = m = 0.8. El centro de masa tiene las
coordenadas aproximadas (0.38, 0.38, 0.8).
Ejemplo 158. Calcule la integral triple quepproduce el volumen del solido que tiene la forma
determinada por el cono de un manto x = y 2 + z 2 y el paraboloide x = 6 y 2 z 2 .
Solucion. Primero debemos hallar la interseccion de las dos superficies para luego proyectarla
a un plano coordenado Observe que al sustituir y 2 + z 2 = x2 en y 2 + z 2 = 6 x encontramos
que x2 = 6 x o (x + 3)(x 2) = 0. As, las dos superficies se intersecan en el plano x = 2.
La proyeccion del solido sobre el plano yz de la curva de interseccion es y 2 + z 2 = 4. Luego
esta regon es tipo I o II
2
2 4y2 2
6y z
V = dxdzdy = 4 dxdzdy
D 0 0 y 2 +z 2
Ejemplo 159. Use coordenadas cilndricas para hallar el volumen de la region solida Q que
corta en la esfera x2 + y 2 + z 2 = 4 el cilindro r = 2 sen .
Solucion. Como x2 +y 2 +z 2 = r2 +z 2 = 4 los lmites z son z = 4 r2 Sea R la proyeccion
circular del solido sobre el plano r. Los lmites de R son 0 r 2 sin y 0 . Por
lo tanto,
V = dV = rdzdrd
Q
124
Coordenadas Esfericas
De la Figura vemos que el volumen de una cuna esferica esta dado por la aproximacion
V 2 sin
De tal modo, en una integral triple de una funcion continua en coordenadas esfericas f (, , )
la diferencial de volumen dV es
V = 2 sin ddd
Por consiguiente, una integral triple comun en coordenadas esfericas tiene la forma
" h2 (,) # g2 () h2 (,)
f (, , )dV = f (, , )d dA = f (, , )ddd
S R h1 (,) g1 () h1 (,)
Ejemplo 160. Un solido enpel primer octante tiene la forma determinada por la grafica del
cono de un solo manto z = x2 + y 2 y los planos z = 1, x = 0. Determine en coordenadas
esfericas para determinar el volumen del solido.
Solucion. Usando las ecuaciones. x = sin cos , y = sin sin , z = cos . Encontramos
que
125
Ejemplo 161. Hallar el volumen del elipsoide dado por 4x2 + 4y 2 + z 2 = 16
Ejemplo 162. Determine el momento de inercia en torno al eje z del solido homogeneo
acotado entre las esferas x2 + y 2 + z 2 = a2 y x2 + y 2 + z 2 = b2 , a < b.
Ejemplo 163. Hallar el volumen en coordenadas esfericas de la region solida D limitada
inferiormente por la hoja superior del cono z 2 = x2 + y 2 y superiormente por la esfera
x2 + y 2 + z 2 = 9.
126
como una transformacion T uno a uno de la region E en el espacio uvw a la region D en el
espacio xyz. Adicionalmente, el Jacobiano de T esta dado por
x
u xv xw
(x, y, v)
JT (u, v, w) = = yu yv yw
(u, v, w)
zu zv zw
Teorema 47. Si x = x(u, v, w), y = y(u, v, w), z = z(u, v, w) es una transformacion que
mapea una region E en el espacio uvw hacia una region D en el espacio xyz y f es una
funcion continua sobre E, entonces
f (x, y, z)dV = f x(u, v, w), y(u, v, w), z(u, v, w) |J(u, v, w)|dV
D E
127
Cuatro
Calculo Vectorial
Hasta este punto en nuestro estudio del calculo, hemos encontrado tres tipos de integrales:
la integral definida, la integral doble y la integral triple. En este captulo se presentan dos
nuevos tipos de integrales: las integrales de lnea y las integrales de superficie.
Suponga que C es una curva parametrizada por x = x(y), y = y(t) a t b, y que
A = (x(a), y(a)) y B(x(a), y(a)) son los puntos inicial y terminal,respectivamente. Afirmamos
que:
1. C es una curva suave si x (t) y y (t) son continuas sobre el intervalo cerrado [a, b] y
no son simultaneamente cero sobre el intervalo abierto (a, b)
2. C es una curva suave por partes si consiste en un numero finito de curvas suaves
c1 , C2 , . . . , Cn unidas extremo por extremo; esto es, C = C1 C2 Cn .
3. C es una curva cerrada si A = B.
4. C es una curva simple si no se cruza a s misma entre A y B.
5. C es una curva cerrada simple si A = B y la curva no se cruza a s misma.
6. Si C no es una curva cerrada, entonces la orientacion impuesta sobre C es la direccion
que corresponde a los valores crecientes de t.
7. Esta misma terminologa lleva de manera natural a las curvas en espacio tridimensional
C(t) = (x(y), y(t), z(t))
A = P0 , P1 , P2 , . . . , Pn = B
128
5. Sea kP k la longitud del subarco mas largo.
6. Escoja un punto muestra (xk , yk ) sobre cada subarco sk como se ilustra en la FIG.
9. Tomamos el lmite de estas tres sumas cuando kP k 0. Las integrales que resultan se
resumen a continuacion.
Definition 48. Sea f una funcion de dos variables x y y definida en una region del plano
que contiene una curva suave C.
Interpretacion Geometrica
Como Evaluar integrales de linea La idea basica es convertir la integral de lnea en una
integral definida de una sola variable. Para ello nos enfocaremos en como esta definida la
curva C (parametricamente o ecuacion explcita).
129
p
Ademas, recordando el hecho de que ds = dx2 + dy 2 y que C esta parametrizada, encon-
tramos que p p
ds = (x (t)dt)2 + (y (t)dt)2 = x (t)2 + y (t)2 dt
Por consiguiente, podemos concluir que
b p b
f (x, y)ds = f (x(t), y(t)) x (t)2 + y (t)2 dt = f (C(t))kC (t)kdt,
C a a
Curva Ingeniero Si la curva C esta definida por una funcion explcita y = g(x),
a x b (o x = h(y), c y d) es posible utilizar x o (y) como un parametro. Luego las
integrales de lnea son mas faciles
f (x, y)dx = , f (x, y)dy = , f (x, y)ds =
Cx Cx Cx
f (x, y)dx = , f (x, y)dy = , f (x, y)ds =
Cy Cy Cy
Una integral de lnea a lo largo de una curva C suave por partes se define como la suma
de las integrales sobre las distintas curvas suaves cuya union compone a C.
f (x, y)ds = f (x, y)ds + f (x, y)ds
C C1 C2
NOTA:
3. El smbolo C denota la curva que tiene los mismos puntos pero la orientacion opuesta
de C.
P dx + Qdy = P dx + Qdy
C C
130
Ejemplo 164. Sea C una curva que da vuelta alrededor de la circunferencia x2 + y 2 = 1 en
sentido contrario a las manecillas del reloj. Calcule 1 + xds
C
Ejemplo 165. Calcule xy 2 zds donde C es la curva interseccion de las superficies x2 +
C
y 2 + z 2 = 16 y x2 + y 2 = 4 en el primer octante.
xy 2 zds =
C
Un campo vectorial sobre una region T en el espacio es una funcion de variable vectorial
F que asocia con cada punto (x, y, z) de T un vector
Observe que las componentes P , Q y R de una funcion vectorial son funciones escalares
(de variable real). Un campo vectorial en el plano es similar, excepto que no se involucran
componentes z ni coordenadas z. As, un campo vectorial en la region plana R es una funcion
F de variable vectorial que asocia un vector a cada punto (x, y) de R
Es util visualizar un campo vectorial F dado. Una forma comun de hacerlo es dibujando
un conjunto de vectores comunes F (x, y), cada uno de los cuales esta representado por una
flecha de longitud kF (x, y)k y se localiza con (x, y) como su punto inicial.
131
Ejemplo 166. Describa el campo vectorial F(x, y) = xi + yj.
Solucion. Para cada punto (x, y) en el plano coordenado, F(x, y) = (x, y) es decir, es sim-
plemente su vector de posicion r. Apunta directamente hacia afuera del origen y tiene una
longitud de p
kF (x, y)k = kxi + yjk = x2 + y 2 = r,
igual a la distancia desde el origen (x, y).
OBSERVACION: Entre los campos vectoriales mas importantes debido a sus aplicaciones,
se hallan los campos vectoriales de velocidad. Imagine el flujo estable de un fluido, como el
agua de un ro o el viento solar. Por flujo estable queremos decir que el vector de velocidad
v(x, y, z) del fluido que pasa a traves de cada punto (x, y, z) es independiente del tiempo
(aunque no necesariamente independiente de x, y y z), por lo que el patron del flujo permanece
constante. Por lo tanto, v(x, y, z) es el campo vectorial de la velocidad del movimiento del
fluido.
Ejemplo 167. Suponga que el plano horizontal xy esta cubierto por una lamina delgada
de agua en movimiento giratorio (como en un remolino) alrededor del origen con velocidad
angular constante de radianes por segundo en direccion contraria a la del movimiento de
las manecillas del reloj. Describa el campo vectorial de la velocidad asociado.
Solucion. En este caso se tiene un campo vectorial de dos dimensiones v(x, y). El agua se
mueve en cadap punto (x, y) con velocidad v = r y en forma tangencial a la circunferencia
de radio r = x2 + y 2 . El campo vectorial
Ejemplo 168. Suponga que una masa M esta fija en el origen en el espacio. Cuando una
partcula de masa unitaria se coloca en el punto (x, y, z) distinto del origen, esta sujeta a
una fuerza F(x, y, z) de atraccion gravitacional dirigida hacia la masa M en el origen. De
2
acuerdo con pla ley de Newton del inverso del cuadrado, la magnitud de F es F = M G/r ,
donde r = x2 + y 2 + z 2 es la longitud del vector de posicion r = xi + y + zk. Se sigue, de
inmediato, que
kr
F(x, y, z) = 3
r
donde k = M G. Observe que F(x, y, z) no esta definido en el origen y que |F | + cuando
r 0+ .
132
El campo vectorial gradiente
Ejemplo 170. Dadas las funciones derivables f (x, y, z) y g(x, y, z), demuestre que
(f g) = f g + gf
Solucion. Aplicamos la definicion del gradiente y la regla del producto para la derivacion
parcial. encontramos que,
(f g) =
Ahora suponga que F = P i+Qj+Rk es un campo de fuerza definido sobre una region que
contiene a la curva C del punto A al B. Tambien suponga que C tiene una parametrizacion
Queremos aproximar el trabajo W realizado por el campo de fuerzas F cuando una partcula
se mueve a lo largo de la curva C, de A a B.
133
1. Subdivida C.
v dx dy dz
Tds = vdt = vdt = , , dt = dxi + dyj + dzk = dr
v dt dt dt
Por lo que,
Tds = dr
Con esta notacion, tenemos que
W = F T ds = F dr
C C
134
Por lo tanto
W = P dx + Qdy + Rdz
C
Este calculo revela una relacion importante entre los dos tipos de integrales de lnea que se
han definido aqu.
Teorema 49. Suponga que el campo vectorial F = P i+Qj+Rk tiene funciones componentes
continuas y que T es el vector unitario tangente a la curva suave C. Por lo cual
W = F Tds = P dx + Qdy + Rdz = F dr (4.3)
C C C
Solucion. No es difcil ver que dx = dt, dy = 2tdt, y dz = 3t2 dt, la sustitucion en terminos
de t produce
1
W = F dr =F Tds = ydx + zdy + xdz =
C C C 0
Ejemplo 172. Evalue xydx + x2 dy, donde C esta dada por y = x3 . 1 x 2.
C
Solucion. Como la curva C es una curva funcion en x entonces usaremos a x como el parame-
tro. Luego
xydx + x2 dy =
C
Ejemplo 173. Evalue y 2 dx x2 dy, sobre la curva cerrada C que se muestra en la fig 1.
C
Solucion. Como la curva C es suave por partes, simbolicamente, escribimos C = C1 + C2 + C3
donde C1 , C2 y C3 son las curvas, y = 0, x = 2 y y = x2 respectivamente.
135
Ejemplo 174. Evalue ydx + xdy + zdz, donde C es la helice x = 2 cos t, y = 2 sen t, z = t
C
0 t 2.
F(r(t) = r (t)dt =
Luego, F dr = F(r(t)) r (t)dt =
C C
Rotacional y Divergencia
Hemos visto que si un campo vectorial de fuerza F es conservativo, entonces puede escri-
birse como el gradiente de una funcion potencial :
F = = i+ j+ k.
x y x
El operador diferencial vectorial, u operador nabla
=i +j +k
x x x
combinarse con un campo vectorial
de dos modos diferentes: en un caso produciendo otro campo vectorial llamado Rotacional
de F y en el otro caso produciendo una funcion escalar llamada la Divergencia de F.
Definition 50 (Rotacional de un campo vectorial).
El rotacional de un campo vectorial F = P i + Qj + Rk es el campo vectorial
R Q P R Q P
rot F = i+ j+ k
y z z x x y
OBSERVACION: NOOO es necesario memorizar los complicados componentes del campo ro-
tacional de F, basta con interpretar el rotacional F como el producto cruz de y el vector
136
F: Es decir,
i j k
rot F = F =
x y z
P Q R
Solucion.
Divergencia Si F(x, y, z) = P (x, y, z)i + Q(x, y, z)j + R(x, y, z)k representa el campo de
velocidades de un fluido, entonces sabemos que el flujo F a traves del area S, es aproxima-
damente
137
Calculemos el flujo total de F a traves de sus seis lados de un paraleleppedo rectangular
en la direccion hacia fuera,
Calculamos primero el flujo total hacia el exterior de dos caras paralelas.
El area de la cara F1 es xz, y la normal unitaria hacia fuera es j, luego el flujo de F a
traves de F1 es
(F (j))xz = Q(x, y, z)xz
El flujo hacia fuera de la cara F2 , cuya normal hacia fuera es j, esta dado por
De manera analoga, vemos que las contribuciones al flujo total hacia fuera del para-
leleppedo desde las dos caras paralelas al plano yz, y desde las dos caras paralelas al
plano xy, originan
P R
xyz xyz
x z
Al sumar estos resultados, vemos que el flujo total de F hacia fuera del paraleleppedo
es aproximadamente
P Q R
+ + xyz
x y z
Dividiendo la ultima expresion entre xyz obtenemos el flujo hacia fuera de F por
unidad de volumen:
P Q R
+ +
x y z
Esta combinacion de derivadas parciales es una funcion escalar y recibe el nombre
especial de divergencia de F.
Definition 52 (Divergencia).
La divergencia de un campo vectorial F = P i + Qj + Rk es la funcion escalar
P Q R
divF = + +
x y z
OBS: La funcion escalar divF tambien puede escribirse en terminos del operador como
un producto punto:
divF = F = P (x, y, z) + Q(x, y, z) + R(x, y, z)
x y z
138
la divergencia de un campo de velocidades F cerca de un punto P (x, y, z) es el flujo por
unidad de volumen.
Si divF > 0 se dice que P es una fuente para F, ya que hay un flujo neto hacia fuera del
fluido cerca de P .
Si divF > 0, se afirma entonces que P es un sumidero para F, puesto que hay un flujo
neto hacia dentro del fluido cerca de P .
Si divF = 0 no hay fuentes o sumideros cerca de P .
El producto punto de consigo mismo obtenemos un importante operador diferencial escalar
de segundo orden llamado Laplaciano:
2 2 2
2 = = + +
x2 y 2 z 2
Cuando este operador actua sobre una funcion escalar f (x, y, z) el resultado se denomina
laplaciano tridimensional,
2f 2f 2f
2 f = + +
x2 y 2 z 2
y aparece en matematicas aplicadas en muchas ecuaciones diferenciales parciales. Una
de las ecuaciones diferenciales parciales mas famosas,
2f 2f 2f
+ + =0
x2 y 2 z 2
recibe el nombre de ecuacion de Laplace en tres dimensiones.
Solucion.
div F = F =
=
APENDICE
Funciones Vectoriales
Definition 53. Una funcion f : I R Rn cuya regla de correspondencia es
f (t) = f1 (t), f2 (t), . . . , fn (t)
139
4. Si esta regla de correspondencia la escribimos en la forma
x1 = f1 (t), x2 = f2 (t), xn = fn (t), t I. (4.4)
Los puntos (x1 , . . . , xn ) = (f1 (t), . . . , fn (t)), t I forman la curva C parametrizada en
el espacio Rn , y describen normalmente la trayectoria de una particula,
5. A las ecuaciones xi = fi (t) se llaman ecuaciones parametricas de una curva C. Si en las
ecuaciones parametricas xi = f (t) de la curva C se elimina el parametro t, logramos
encontrar las ecuaciones cartesianas de la curva (similar a las ecuacion simetricas
de una recta)
6. Si f : I Rn es una funcion vectorial tal que f (t) = (f1 (t), . . . , fn (t)) entonces f (t) es
el vector de posicion del punto P (f1 (t); f2 (t), . . . , fn (t)) en la curva C. El extremo del
vector de posicion f (t) traza la trayectoria de la curva C e indica su orientacion.
Las funciones f , g y h son las funciones componentes del vector de posicion r. Consideremos
la trayectoria de la partcula como la curva descrita por r durante el intervalo de tiempo I.
la curva.
140
Una funcion de la forma
r(t) = f (t)i + g(t)j
representa un curva el plano. Mientras que un funcion vectorial de la forma
representa una curva en el espacio, donde las funciones componentes f , g y h son funciones
del parametro t. Vease figura
Ejemplo 179. Grafique la funcion vectorial r(t) = cos ti + sin tj + tk
Como x2 + y 2 = 1 entonces afirmamos que la curva esta sobre un cilindro circular, pero la
curva sube cuando el componente en k, z = k aumenta. Adicionalmente, Cada vez que t
aumenta en 2, la curva completa una vuelta, solo en en el plano xy. De manera que la
trayectoria de esta curva es una Helice
La curva en el ejemplo anterior es una de tipos de curvas espaciales conocidas como curvas
helicoidales. En general, una funcion vectorial de la forma
describe una helice circular . El numero 2c/k recibe el nombre de horquilla de una
helice, (es la separacion vertical de los lazos de la helice). Una helice circular es solo un caso
especial de la funcion vectorial
141
se denomina helice conica. Por ultimo, una curva dada por
r(t) = a sin(kt) cos(t)i + a sin(kt) sin(t)j + a cos(kt)k
se llama helice esferica.En estas ecuaciones se supone que a, b, c y k son constantes posi-
tivas.
Ejemplo 180. Encuentre una funcion vectorial para el segmento de recta del punto P0 (3, 2, 1)
al punto P1 (1, 4, 5).
SOL: No es dificil ver que los vectores de posicion correspondientes a los puntos dados
son r0 = (3, 2 1) y r1 = (1, 4, 5). Entonces, una funcion vectorial para el segmento de recta
es:
Ejemplo 181. Grafique la curva trazada por la funcion vectorial r(t) = 2 cos ti+ 2 sen tj + 3k
SOL: Las ecuaciones parametricas de la curva son las componentes de la funcion vectorial
son
x = 2 cos t, y = 2 sen t, z=3
.
Ejemplo 182. Determine la funcion vectorial que describe la curva C de interseccion del
plano y = 2x y el paraboloide z = 9 x2 y 2
x = t, y = 2t, z = 9 5t2 ,
y por tanto una funcion vectorial que describe el trazo del paraboloide en el plano y = 2x
esta dada por
r(t) = ti + 2tj + (9 5t2 )k.
Ejemplo 183. Encuentre la funcion vectorial que describe la curva C de interseccion de los
cilindros y = x2 y z = x3
142
x2
Ejemplo 184. Dibujar la grafica C representada por la interseccion del semielipsoide 12 +
y2 z2 2
24 + = 1, z 0 y el cilindro parabolico y = x . Despues, hallar una funcion vectorial
4
que represente la grafica.
z2 t2 t4 (6 + t2 )(4 t2 )
=1 + = ,
4 12 24 24
Como la curva se encuentra sobre el plano xy, hay que elegir para z la raz cuadrada positiva.
Por tanto, la funcion vectorial resultante es:
r
2 (6 + t2 )(4 t2 )
r(t) = ti + t j + k, 2 t 2
6
Muchas de las tecnicas y definiciones utilizadas en el calculo de funciones reales se pueden
aplicar a funciones vectoriales. Por ejemplo, las funciones vectoriales se pueden sumar y
restar, multiplicar por un escalar, tomar su lmite, derivarlas, integrarlas y as sucesivamente.
La estrategia basica consiste en aprovechar la linealidad de las operaciones vectoriales y
extender las definiciones, componente por componente.
Esta extension, componente por componente, de las operaciones con funciones reales a
funciones vectoriales se ilustra mas ampliamente en la definicion siguiente del lmite de una
funcion vectorial.
143
Si r(t) tiende al vector L cuando t a, la longitud del vector r(t) L tiende a 0. Es
decir,
kr(t) Lk 0 Cuando ta
Ejemplo 185. Calcule lm f (t) (en caso exista) de las siguientes funciones vectoriales
tt0
!
1 t+1 t
1. f (t) = , , 2 , para t0 = 0 R/ (0, 0, 2)
t+2 t+1
!
et e ln t sen(t 1)
2. f (t) = , , , para t0 = 1 R/ (e, 1, 1)
t1 1t t1
!
1 cos(sen t) cos t cos(sen t) 1
3. f (t) = , , , para t0 = 0 R/ ( 12 , 0, 1 )
sen2 t t2 t
!
tan( sen( 5 t 1) t1
4. f (t) = (2 t) 2 t) , , , para t0 = 1 R/ (e2/ , 1, 12 )
tan( 5 t 1) t 1
144
!
sen t ln(1 + t) cos t 1
2. Dada la funcion vectorial f (t) = , , Determine si la funcion
t 1t t
vectorial es continua en t = 0. R/: NO
f (t + h) f (t)
f (t) = lm
ta h
para todo t para el cual existe el lmite.
NOTA: Ademas de la notacion f (t) otras notaciones para la derivada de una funcion
vectorial son
d df
Dt [f (t)], [f (t)],
dt dt
Velocidad y aceleracon
Si una partcula se mueve a lo largo de una curva C en el espacio Rn , de modo que su
vector posicion en el tiempo t es f (t) = (f1 (t), f2 (t), . . . , fn (t)) entonces, el vector velocidad
v(t) y el vector aceleracion a(t) de la partcula en el instante t son dadas por
v(t) = f (t) = f1 (t), f2 (t), . . . , fn (t)
a(t) = v (t) = f (t) = f1 (t), f2 (t), . . . , fn (t)
El vector velocidad v(t) tiene la direccion del vector tangente a la curva C en el punto f (t)
y el vector aceleracion a(t) apunta hacia el lado concavo de la curva C (lado hacia donde se
doble la curva).
145
Reglas de Derivacion de funciones vectoriales: Sean f , g : I Rn funciones
vectoriales derivables de t, c una constante real y : I R una funcion real derivable
de t. Entonces se tiene:
d dv
0= (v v) = 2v = 2v a Entonces, v(t) a(t) = 0 para todo t.
dt dt
Ejemplo 187. Si f (t) = (t, t2 ; 3 + t), g(t) = (cos t, sen t, ln(t + 1)) y (t) = e4t , calcule
146
Ejemplo 188. Considere la curva C dada por r(t) = cos(2t)i + sin(t)j, /2 t /2.
Encuentre la derivada r (t) y grafique los vectores r (0) y r (/6)
en concecuencia
1
r (0) = j, r (/6) = 3i + 3j
2
Como graficamos la curva? R/:
Ejemplo 189. Encuentre las ecuaciones parametricas de la recta tangente a la curva C
cuyas ecuaciones parametricas son x = t2 , y = t2 t, z = 7t en el punto correspondiente
a t = 3.
Definition 56.
Si f : [a, b] Rn es una funcion vectorial continua en el intervalo [a, b] tal que
f (t) = (f1 (t), f2 (t), . . . , fn (t)) entonces la integral indefinida de f es
!
f (t)dt = f1 (t)dt, f2 (t)dt, . . . , fn (t)dt
y la integral definida de f es
!
b b b b
f (t)dt = f1 (t)dt, f2 (t)dt, . . . , fn (t)dt
a a a a
1
Ejemplo 190. 1. Halle la integral indefinida de la funcion vectorial f (t) = cos t, , tet
1+t
1
1
2. Calcule la integral f (t)dt, donde f (t) = (2t, , tet )
0 1+t
147
Teorema 57 (Teorema Fund. del Calculo). Sea f : [a, b] Rn es una funcion
vectorial continua en [a, b] entonces
t
La funcion F definida por F(t) = f (u)du, a t b es derivable y F (t) =
a
f (t), t [a, b]
b
f (u)du = F(b) F(a),
a
/4
Ejemplo 191. Calcule f (t)dt h(0) , donde
0
r !
4 3 21 3 2 t
f (t) = tan t sec t, sen (2t) cos t sen (2t) sen t,
y !
1 1 1
5 5 10
t2 1 2 3
h(t) = e dt, (t t)dt, t dt , R:( , , )
1 0 0 288 2 63
Curvas regulares
Definition 58.
Se dice que una curva C Rn es una curva parametrizada, si existe una fun-
cion vectorial
: [a, b] Rn tal que a ([a, b]) = C A la funcion vectorial
(t) = 1 (t), 2 (t), . . . , n (t) se llama parametrizacion de la curva C.
148
1. Se dice que C es una curva con puntos dobles si (t1 ) = (t2 ), con t1 6= t2
2. Se dice que C es una curva simple si no tiene puntos dobles.
3. Se dice que C es una curva cerrada si (a) = (b)
4. Se dice que C es una curva regular, si la funcion vectorial (t) tiene derivada continua
y a (t) 6= 0 , t [a, b]
donde : [c; d] [a, b] es una funcion real derivable y sobreyectiva tal que (u) 6= 0,
u [c, d].
es una reparametrizacion de la curva (t). Com (u) = 2 > 0, entonces la curva (u)
mantiene la misma orientacion de la curva (t).
2. Consideremos la funcion : [0, 2] [0, 2] definida por (u) = 2 u entonces
149
es una reparametrizacion de (t). Como (u) = 1 < 0, entonces la curva (u) invierte
la orientacion de la curva (t).
1 1
2. (t) = (t, 1, t3 + ), desde t = 1 hasta t = 3. R/ 14
3
6 2t
150
Vectores unitarios: Tangente, Normal principal y Binormal
(t)
T(t) =
k (t)k
OBS: Como kTk = 1 entonces T(t) T(t) = 1, luego al derivar esta expresion, se tiene
T (t)
N(t) = siempre que kT (t)k 6= 0
kT (t)k
(t)
Luego, de la expresion del vector tangente unitario T(t) = k (t)k , se tiene
Esta ecuacion indica que la direccion del vector velocidad (t) es igual a la del vector
tangente unitario T(t) y la velocidad escalar o rapidez es dada por
l (t) = k (t)k
Por tanto, si un objeto se mueve a lo largo de una curva C, el vector tangente unitario
T(t) apunta en la direccion del movimiento, mientras que el vector normal principal N(t)
es ortogonal a T(t) y senala la direccion hacia donde gira el objeto (lado concavo de la curva
C). Ademas kN (t)k = 1, t [a, b].
OBS: Sea : [a, b] Rn una curva regular, tal que C = ([a, b]). Si (t) es derivable en
[a, b], entonces al derivar la expresion (4.7) resulta
Luego, el vector aceleracion a (t) es combinacion lineal de los vectores tangente unitario
T(t) y normal principal N(t).
Definition 63 (Vector Binormal). Sea : [a, b] Rn una curva regular tal que
(t) 6= 0, t [a, b]. El vector unitario dado por B(t) = T(t) N(t) se denomina
vector binormal a la curva : [a, b] Rn en el punto (t)
151
Los tres vectores unitarios T, N y B forman un conjunto de mano derecha de vectores
mutuamente ortogonales denominado triedro movil.
[ (t)]
T= , B= , N=
k k k k [ (t)]
y = (x 2)2 z = (x 2)2
Halle la ecuacion cartesiana del plano osculador a la curva C en los puntos (2, 0, 0), (3, 1, 1)
y (x0 , y0 , z0 ) R/: PO : y z = 0, B(t) es const
152
Curvatura y torsion de una curva
Reparametrizacion por Longitud de arco
Teorema 64. Sea : [a, b] Rn una curva regular tal que ([a, b]) = C y la longitud
b
de arco de la curva desde t = a hasta t = b es L = k (t)kdt Entonces, la funcion
a
longitud de arco l : [a, b] [0, L] dada por
t
s = l(t) = k (t)kdt
a
Ejemplo 198. Sea : [0, +) R3 una curva regular definida por (t) = (cos t, sen t, t).
Halle:
153
Ejemplo 199. Dada la curva C : (t) = t, ln(sec t), ln(sec t + tan t) , t [0, 4 ] Halle la
longitud de arco de C y reparametrice
! C en terminos de longitud de arco. R: a) 2 ln( 2+1),
2s 1
e 2
parte del b) (s) = arc sen 2s +1
e 2
Luego, el vector curvatura K(t) tiene la misma direccion que el vector unitario normal
principal N(t) y es ortogonal al vector tangente unitario.
154
Definition 67 (Curvatura).
La funcion escalar que multiplica a N(t) en (4.8) se denomina curvatura de la curva
C en el punto (t) y se denota por
kT (t)k
k(t) =
k (t)k
La curvatura k(t) es un numero real que nos indica que tanto se tuerce (o se dobla)
la curvatura C en el punto (t).
T(t) = T (t) =
kT (t)k
k(t) = =
kr (t)k
Ejemplo 201. En el espacio tridimensional la posicion de una partcula en movimiento
esta dada por la funcion vectorial r(t) = 2 cos t + 2sentj + 3tk. Encuentre los vectores T(t),
N(t) y B(t). Determine la curvatura (t).
Ahora el binormal es
T (t)k
Por ultimo encontramos la curvatura k(t) = =
kr (t)k
OBSERVACIONES:
155
4. Si C : : I R3 es una curva regular, entonces la curvatura de la curva C en punto
(t) en terminos de sus derivadas es dado por
k (t) (t)k
k(t) = , solo en R3
k (t)k3
Ejemplo 202. Sea C la curva de interseccion del cilindro x2 + y 2 + 2(y x) 2 = 0 con el
plano x y 2z 2 = 0 Determine la curvatura de C en el punto (3, 1, 1).
De donde resulta
156
3. Como los vectores T y N estan en el plano osculador, entonces la circunferencia de
curvatura se encuentra tambien sobre el plano osculador.
4. La circunferencia de curvatura (C1 ) esta en el lado concavo o interior de la curva C y
tiene la mism a curvatura que C en (t0 ).
5. El centro de curvatura de la curva C en el punto (t0 ) es ei centro de ia circunferencia
de curvatura (C1 ) y es dado por
Solucion. Como la curva interseca al plano dado, entonces las componentes de la funcion
vectorial (t) satisface la ecuacion del plano, esto es,
t t
1 2t sen u 1
+ e du + t = , esen u du = 0, t = 2
2 2 2 2
Luego, la curva C corta al plano dado en el punto (2) = P0 14 2 , 0, 2) Por otro lado, se
tiene
La ecuacion del plano osculador de la curva C que pasa por el punto (2) es
" #
1 4 1
P0 : (x, y, z) , 0, 2) (1, 0, 1) = 0, P0 : x + z =
2 2
Como las coordenadas del centro de curvatura satisfacen la ecuacion del plano osculador,
entonces la circunferencia de curvatura se encuentra sobre el plano osculador.
Definition 69 (Torsion).
Sea : [a, b] R3 una curva regular parametrizada por longitud de arco. Existe una
funcion real (t) que satisface
dB
= (t)N(t)
ds
Al numero real (t) se llama torsion de la curva C en el punto (t).
157
Supongamos que
t
(t) = (1 (t), 2 (t), 3 (t)) y s = l(t) = k (u)kdu, t [a, b]
0
dB
determina el grado de torsion de la curva C en el punto (t). Para los vectores unitarios,
ds
se tiene
Luego, los vectores N(t) y B (t) son paralelos. Al derivar el vector binormal con respecto al
param etro longitud de arco, se obtiene
dB dB(t) dt 1 dB(t) 1
= =
=
B (t)
ds dt ds k (t)k dt k (t)k
OBSERVACIONES:
158
Ejemplo 206. Dada la curva C : (t) = t sen t, 1 cos t, 4 sen 2t , halle la curvatura y
la torsion de la curva C en el punto donde el plano normal principal a la curva es paralelo
al plano z = 1. R: (0) = (0, 0, 0), k(0) = 14 , (0) = 21
donde T es el vector tangente unitario y l es la funcion longitud de arco. Por otra parte
sabemos que, el vector aceleracion es dado por
Definition 70 (Componentes aT y aN ).
La aceleracion de un particula esta dada por
kv(t)k = l (t)
159
2. La componente normal de la aceleracion es siempre positiva.
4. Esto explica por que un automovil que toma una curva cerrada a velocidad moderada o a
una curva suave a gran velocidad exige, en ambos casos, una fuerza normal (rozamiento
de los neumaticos) de gran magnitud para que el vehculo no se salga de la carretera.
Halle sus vectores velocidad y aceleracion, su velocidad escalar, los vectores unitarios T y N.
y los componentes normal y tangencial del vector aceleracion, todo para t = /3.
Solucion.
v(t) = (t) = (1, sec t, tan t), v( ) = (1, 2, 3)
3
2
a(t) = (t) = (0, sec t tan t, sec t), a( ) = (0, 2 3, 4)
3
l (t) = k (t)k = 2 sec t, l( )=2 2
3
l (t) = 2 sec t tan t, l ( )=2 6
3
La velocidad escalar, el vector tangente unitario y la curvatura en t = /3 son
v = k( )k = l ( ) = 2 2
3 3
( 3 ) 1
T( ) = = (1, 2, 3)
3 k ( 3 )k 2 2
(3) (3) 1
k( ) = =
3 k ( 3 )k3 4
2
Como a = l ( 3 )T + k l ( 3 ) N entonces N = ( 23 , 0, 12 ) y
2
aT ( ) = l ( ) = 2 6, aN ( ) = k l ( ) = 2
3 3 3 3
160
Teorema fundamental de las integrales de lnea
El estudio realizado anteriormente indica que el trabajo W realizado por un campo F sobre
un objeto que se mueve entre dos puntos generalmente depende de la trayectoria seguida por
el objeto. Ahora vamos a estudiar el trabajo W realizado por ciertos campos F , pero con
resultados sorprendentes, mas precisamente estudiaremos el teorema fundamental de las
integrales de lnea.
Para empezar, definamos lo que es un campo vectorial conservativo F
Definition 71. Un campo vectorial F se dice que es conservativo si existe una funcion
diferenciable tal que
F = .
La funcion recibe el nombre de funcion potencial de F.
Solucion. No es difcil ver que que la funcion potencial de F esta dada por (x, y) = xy, ya
que verifica que = F. Asi que F es un campo conservativo.
Ejemplo 209. La integral ydx + xdy tiene el mismo valor en cada trayectoria C entre
C
(0, 0) y (1, 1) que se muestra en la siguientes figura
Ejemplo 210. Hallar el trabajo realizado por el campo de fuerzas F(x, y) = 12 xyi + 14 x2 j
sobre una partcula que se mueve de (0, 0) a (1, 1) a lo largo de cada una de las trayectorias
1. C1 : y = x 2. C2 : x = y 2 3. C3 : y = x3
161
Solucion. a) No es dificil parametrizar C1 por el vector posicion r(t) =
para 0 t 1 Por lo tanto dr = y F(x, y) = . De
manera que el trabajo realizado es:
W = F dr =
C1
b) Una parametrizacion para el vector posicion de C2 es r(t) = ti + tj para 0 t 1.
Por lo tanto dr = y F(x, y) = . De manera que el trabajo
realizado es:
W = F dr =
C1
Por tanto, el trabajo realizado por un campo vectorial conservativo es el mismo para todas
las trayectorias, 14
OBSERVACIONES:
Sin embargo, hay excepciones muy importantes, como lo son los campos conservativos.
162
donde f (x) es una antiderivada de f (x). En el siguiente teorema, conocido como teorema
fundamental para integrales de lnea,
Teorema 72. Suponga que C es la trayectoria en una region abierta R del plano xy dada
por r(t) = x(t)i + y(t)j, a t b. Si F(x, y) = P (x, y)i + Q(x, y)j es un campo vectorial
conservativo en R y es una funcion potencial para F, entonces
F dr = dr = (B) (A)
C C
F = =
dx dy
Usando despues r (t) = i+ j es posible escribir la integral de lnea de F a lo largo
dt dt
de la trayectoria C como
F dr = F r (t)dt =
C C
b
= regla de la cadena
a
= (B) (A)
NOTA: Para curvas suaves por partes, la prueba anterior debe modificarse considerando
cada arco suave de la curva C.
1. Si el valor de una integral de lnea es el mismo para cada trayectoria en una region
que conecta el punto inicial A y el punto terminal B, entonces se dice que la integral
sera independiente de la trayectoria.
2. Una integral de lnea Fdr es independiente de la trayectoria si Fdr = Fdr
C C1 C2
para cualesquiera dos trayectorias C1 y C2 entre A y B.
3.
Si F es un campo vectorial conservativo en una region abierta en el espacio, entonces
F dr depende solo de los puntos inicial y terminal A y B, y no de la trayectoria C.
C
Ejemplo 211. Hallar el trabajo realizado por el campo de fuerzas F(x, y) = 12 xyi + 14 x2 j
sobre una partcula que se mueve de (0, 0) a (1, 1) a lo largo de cada una de las trayectorias
163
1. C1 : y = x 2. C2 : x = y 2 3. C3 : y = x3
(1,1) (1,1)
F dr = F dr = dr =
C (0,0) (0,0)
Ejemplo 212. Evalue ydx + xdy donde C es la trayectoria con punto inicial (0, 0) y
C
punto terminal (1, 1).
(1,1) (1,1)
ydx + xdy = F dr = dr
C (0,0) (0,0)
Ejemplo 213. Evalue F dr donde C es una curva suave a trozos desde (1, 4) hasta
C
(1, 2) y F(x, y) = 2xyi + (x2 y)j.
Solucion. No es difcil ver que F = 2xyi + (x2 y)j es un campo vectorial conservativo
definido en el plano xy y que (x, y) = es una funcion potencial
para F. Por tanto tenemos
(1,4) (1,4)
F dr = F dr = dr
C (1,2) (1,2)
Ejemplo 214. Evalue F dr donde C es una curva suave a trozos desde (1, 1, 0) hasta
C
(0, 2, 3) y F(x, y, z) = 2xyi + (x2 + z 2 )j + 2yzk.
(1,1,0) (1,1,0)
F dr = F dr = dr
C (0,2,3) (0,2,3)
164
Conceptos equivalentes
3. Si C es cualquier curva cerrada simple que yace por completo en R, entonces la region
en el interior de C tambien yace por completo en R.
5. Si la curva representativa C rodea a uno de los hoyos, entonces NO puede contraerse No conexa
hasta un punto sin dejar la region. Esta ultima region se dice que es multiplemente
conexa.
Teorema 73 (Conceptos Equivalentes).
En una region conexa abierta R, F dr es independiente de la trayectoria C si y solo si
C
el campo vectorial F es conservativo en R.
Multiplemente conexa
Demostracion () Si F es conservativo en R, entonces ya hemos visto que F dr =
C
(A) (B), es decir, la integral es independiente de la trayectoria C. ver Teorema 72
() Suponga que F dr es independiente de la trayectoria en la region R y que (x0 , y0 )
C
y (x1 , y1 ) son puntos arbitrarios en la region R. Sea la funcion definida como
(x,y)
(x, y) = F dr
(x0 ,y0 )
En este caso,
(x,y) (x,y)
=0+ P dx + Qdy = P dx = P (x, y)
x x (x1 ,y) x (x1 ,y)
puesto que la primera integral no depende de x. Sin embargo, sobre el segmento de recta
entre (x1 , y) y (x, y) y es constante de manera que dy = 0. De manera similar es posible
demostrar que = Q(x, y). En consecuencia, de
y
= i+ j = P i + Qj = F(x, y)
x y
concluimos que F es conservativo.
165
Integrales alrededor de trayectorias cerradas Recordemos que una trayectoria, o
curva, C se dice que es cerrada cuando su punto inicial A es el mismo que el punto terminal
B. Si C es una curva parametrica definida por la funcion vectorial r(t), a t b, entonces
C es cerrada cuando A = B, esto es, r(a) = r(b). El siguiente teorema es una consecuencia
inmediata
Teorema 74 (Conceptos Equivalentes).
En una region conexa abierta R, F dr es independiente de la trayectoria C si y solo si
C
donde C2 es
ahora una trayectoria de A a B. Debido a la independencia de la trayectoria,
F dr = F dr. De tal modo, que
C1 C2
F dr = F dr F dr = 0
C C1 C1
RESUMEN: Suponga que F es un campo vectorial conservativo definido sobre una region
conexa abierta y C es una trayectoria cerrada que yace por completo en la region.
F conservativo independencia de la trayectoria F dr = 0
C
Estas equivalencias nos lleva a conluir que: Si la integral de lnea F dr no es inde-
C
pendiente de la trayectoria, entonces el campo vectorial F no es conservativo. Ahora vamos
a determinar de una manera SENCILLA si un campo vectorial F es conservativo.
166
Teorema 75 (Prueba para un campo conservativo).
Suponga que F(x, y) = P (x, y)i + Q(x, y)j es un campo vectorial conservativo en una region
abierta R y que P y Q son continuas y tienen primeras derivadas parciales continuas en R.
Entonces
P Q
= (4.9)
y x
para todo (x, y) en R. Inversamente, si se cumple la igualdad (4.9) para todo (x, y) en una
region R simplemente conexa, entonces F = P i + Qj es conservativo en R.
Demostracion. Probamos la primera mitad del teorema. Suponemos que las funciones P y
Q del campo conservativo F = P i + Qj son continuas y tienen primeras derivadas parciales
continuas en una region abierta R. Puesto que F es conservativo, existe una funcion potencial
tal que = F,
i+ j = = F = P i + Qj
x y
As, P = yQ= . En este caso
x y
P 2 Q 2
= = y = =
y y x yx x x y xy
Ahora como P y Q tienen primeras derivadas parciales continuas en una region abierta R.
P Q
Entonces las derivadas parciales mixtas de segundo orden son iguales y por ello = .
y x
Ejemplo 215. Demuestre que los siguientes campos son conservativos
1. F(x, y) = yi + xj
2. F(x, y) = (x2 2y 3 )i + (x + 5y)j
3. F(x, y) = yexy i xexy j
167
y una curva espacial suave por partes r(t) = x(t)i + y(t)j + z(t)k, a t b, tenemos
que
F dr = dr = (B) (A)
C C
Si C es una curva en el espacio tridimensional, una integral de lnea F dr es inde-
C
pendiente de la trayectoria siempre que el campo vectorial tridimensional
sea conservativo.
Inversamente, si (4.10) se cumple en toda una region apropiada del espacio tridimen-
sional, entonces F es conservativo.
Ejemplo 216.
a) Demuestre que la integral de lnea (y + yz)dx + (x + 3z 3 + xz)dy + (9yz 2 + xy 1)dz
C
es independiente de la trayectoria entre (1, 1, 1) y (2, 1, 4).
(2,1,4)
b) Evalue F dr
(1,1,1)
= P, = Q, =R
168
Por tanto,
=
Por tanto si C es cualquier trayectoria con puntos inicial (1, 1, 1) y terminal (2, 1, 4) encon-
tramos que
(2,1,4)
F dr = (y + yz)dx + (x + 3z 3 + xz)dy + (9yz 2 + xy 1)dz
C (1,1,1)
Ejemplo 217. Calcule (x + z, y z, x y) dr , siendo C la curva de interseccion entre
C
la esfera x2 + y 2 + z 2 = 16 y el cilindro x2 + y 2 = 4x
Conclusiones:
El
trabajo W realizado por la fuerza F a lo largo de una trayectoria C cerrada es cero,
F dr = 0
C
Por esta razon, un campo de fuerzas F de este tipo se dice que es conservativo.
Teorema de Green
Vamos a examinar uno de los teoremas mas importantes del calculo integral vectorial.
Este teorema relaciona una integral de lnea alrededor de una curva cerrada simple suave por
partes con una integral doble sobre la region acotada por la curva.
Suponga que C es una curva cerrada simple suave por partes que forma la frontera de
una region simplemente conexa R. Decimos que la orientacion positiva alrededor de C es
la direccion en la que un punto sobre la curva debe moverse, o la direccion en la que una
persona debe caminar, para completar un solo recorrido de C manteniendo la region R a la
izquierda.
169
Las orientaciones positiva y negativa corresponden a recorridos de C en sentido contra-
rio al de las manecillas del reloj y en el sentido de las manecillas del reloj, respectivamente.
El siguiente teorema recibe el nombre de teorema de Green.
Teorema 76 (Teorema de Green).
Suponga que C es una curva cerrada simple suave por partes con una orientacion positiva
P Q
que limita una region simplemente conexa R. Si P , Q, y son continuas sobre R,
y x
entonces
Q P
F dr = P (x, y)dx + Q(x, y)dy = dA (4.11)
C C R x y
Demostracion. Solo haremos una pruba para una region R que es simultaneamente de tipo
I y de tipo II:
= P (x, y)dx
C
170
Si bien la prueba anterior NO es valida para regiones mas complicadas, el teorema se aplica
a esas regiones, tal como la que se ilustra en la FIGURA. La demostracion consiste en
descomponer R en un numero finito de subregiones para las cuales (4.11) puede aplicarse y
despues se suman los resultados.
El pequeno crculo sobrepuesto sobre el signo integral en el primer termino subraya el hecho
de que la integracion es a lo largo de una curva cerrada; la flecha en el crculo en el segundo
termino refuerza la nocion de que la integracion es a lo largo de una curva cerrada C con
orientacion positiva.
Ejemplo 218. Evalue (x2 y 2 )dx + (2y x)dy donde C consiste en la frontera de la
region en el primer cuadrante que esta acotada por las graficas de y = x2 y y = x3 .
P
Solucion: No es difcil ver que P (x, y) = y Q(x, y) = y que =
y
Q
y = son continuas. Por tanto del TG concluimos que
x
(x2 y 2 )dx + (2y x)dy =
R
3
Ejemplo 219. Evalue (x5 +3y)dx+(2xey )dy donde C es el crculo (x1)2 +(y5)2 = 4
C
P
Solucion. No es difcil ver que P (x, y) = y Q(x, y) = y que =
y
Q
y = son continuas. Por tanto del TG concluimos que
x
3
(x5 + 3y)dx + (2x ey )dy =
C R
171
Ejemplo 220. Determine el trabajo realizado por el campo de fuerza F = (16y +sen x2 )i+
(4ey + 3x2 )j que actua a lo largo de una curva cerrada simple C que se muestra en la Figura
P
No es difcil ver que P (x, y) = y Q(x, y) = y que =
y
Q
y = son continuas. Por tanto del TG concluimos que
x
2 y 2
W = (16y + sen x )dx + (4e + 3x )dy =
C R
W =
P Q
ya que P , Q, y y no son continuas en el origen.
y x
172
Teorema de Green para regiones multiplemente conexas:
El teorema de Green tambien puede extenderse a una region R con hoyos, esto es, una
region que no es simplemente conexa. Recuerde que una region con hoyos se dice que
es multiplemente conexa.
Observe que en la figura 1 se muestra una region R acotada por una curva C que
consiste en dos curvas cerradas simples C1 y C2 esto es C = C1 C2 .
El ultimo resultado sigue del hecho de que las integrales de lnea sobre los cortes
cruzados (trayectorias con orientaciones opuestas) se cancelaran entre s.
y x
Ejemplo 221. Evalue 2 + y2
dx + 2 dy donde C = C1 C2 es la frontera de la
C x x + y2
region sombreada R de la figura.
OBSERVACION IMPORTANTE: Como consecuencia del ejemplo anterior, podemos es-
tablecer un resultado para integrales de lnea que nos permite, en ciertas circunstancias,
sustituir una trayectoria cerrada complicada por una trayectoria que es mas simple.
173
Suponga, como se muestra en la Figura que C1 y C2 son dos trayectorias cerradas
simples suaves por partes que no se intersecan y que tienen la misma orientacion
positiva o contraria a la de las manecillas del reloj.
Suponga ademas que P y Q tienen primeras derivadas parciales continuas tales que
P Q
= en la region R acotada entre C1 y C2 . (F conservativo en esta region
y x
acotada)
Entonces del TG y como = se tiene
C C
0= F dr = P dx + Qdy + P dx + Qdy
C1 C2 C1 C2
= P dx + Qdy P dx + Qdy
C1 C2
Por tanto
P dx + Qdy = P dx + Qdy
C1 C2
y x
Ejemplo 222. Evalue 2 + y2
dx + 2 dy donde C la curva poligonal cerrada con-
C x x + y2
sistente en los cuatro segmentos de recta C1 , C2 , C3 y C4 que se muestran en la figura.
P Q
= =
y x
son continuas sobre la region acotada por C y C . Por consiguiente podemos acambiar de
curva de integracion
y x y x
2 2
dx + 2 dy = dx + 2 dy
C x +y x + y2 C x2 +y 2 x + y2
Una curva ingenieroC descrita por una funcion continua y = f (x) puede parametri-
174
zarse tomando como parametro a x = t de manera que las ecuaciones parametricas son
x = t y y = f (t).
Para parametrizar una superficie ingenieroS : z = g(x, y) necesitaremos dos variables
parametricas. Si x = u y y = v entonces las ecuaciones parametricas para S estan dadas
por
x=u y=v z = g(u, v)
p
Ejemplo 223. Encuentre ecuaciones parametricas para el cono de un manto z = x2 + y 2
Solucion.
p No es difcil darse cuenta que cono un manto es una funcion de dos variables
z = x2 + y 2 . Por tanto
el cono podemos parametrizarlos por las siguientes ecuaciones
x = u, y = v y z = u2 + v 2 . Alternativamente, el cono se describe mediante la funcion
vectorial p
r(u, y) = ui + yj + u2 + v 2 k
o tambien podemos parametrizarlo usando coordenadas polares, la parametrizacion
2 2
Ejemplo 224. Encuentre ecuaciones parametricas del cilindro circular y + z = 1 para
3x8
175
Solucion.
OBSERVACION: Resulta practicamente imposible identificar superficies cuando se dan
en forma parametrica o vectorial. Sin embargo, en algunos casos una superficie puede iden-
tificarse al eliminar los parametros.
Ejemplo 225. Identifique la superficie con la funcion vectorial r(u, v) = (2u v)i + (u +
v + 1)j + uk
x = 2u v y =u+v+1 z=u
Plano tangente a una superficie parametrica Recordemos que una superficie para-
metrica esta dada por
r(u, v) = x(u, v)i + y(u, v)j + z(u, v)k
Para generar un plano tangente en el punto (x0 , y0 , z0 ) de S, necesitamos o bien conocer
dos vectores directores o bien, un vector n normal. Los vectores directores del plano que
buscamos son los vectores velocidades de las llamadas curvas u y curvas v,
Por este motivo las derivadas parciales de r con respecto a u y v nos definen los vectores que
usaremos como vectores directores para generar cualquier plano en la superficie S
r x y z
= i+ j+ k Vector velocidad de la Curva u
u u u u
r x y z
= i+ j+ k Vector velocidad de la Curva v
v v v v
Observe que el plano tangente a S existe en (x0 , y0 , z0 ), si
r r
(x0 , y0 , z0 ) 6= 0 (x0 , y0 , z0 ) 6= 0.
u v
176
En muchas ocasiones es mejor encontrar el vector normal n a la superficie S y describir
el plano tangente por medio de la ecuacion normal. Para ello calculamos el producto cruz
r r
u v .
i j k
r r x y z
=
u v u u u
x y z
v v v
r r
La condicion de que el vector 6= 0 en (x0 , y0 , z0 ) asegura la existencia de un plano
u v
tangente en el punto (x0 , y0 , z0 ). Una vez conocido el vector normal podemos describir la
ecuacion del plano
n (x p) = 0
Suponga que S es una superficie parametrica cuya ecuacion vectorial r(u, v) tiene pri-
meras derivadas parciales continuas sobre una region R del plano uv.
r r
Se dice que la superficie S es suave si los vectores tangentes y de las curvas u
u v
r r
y v satisfacen 6= 0 en (u0 , v0 ).
u v
r r
Se afirma que la superficie S es suave sobre R si 6= 0 en todo R.
u v
En terminos generales, una superficie suave no tiene esquinas, puntos afilados o inte-
rrupciones.
Una superficie S suave por partes es una que puede escribirse como S = S1 S2
Sn donde las superficies Si son suaves.
Ejemplo 226. Encuentre una ecuacion del plano tangente a la superficie parametrica defi-
nida por x = u2 + v y = v, z = u + v 2 en el punto correspondiente a u = 3 y v = 0.
Solucion. No es difcil ver la superficie esta descrita por la funcion vectorial r(u, v) = u2 +
vi + vj + (u + v 2 )k y que en u = 3 y v = 0 el punto sobre la superficie es .
Entonces los vectores tangentes estan dados por
177
Area de una superficie parametrica.
El area de la superficie
n n
X X r r
A(S) = area(Sk )
u v uv
i=1 k=1
178
OBSERVACION: Si una superficie esta descrita por una funcion explcita z = g(x, y)
puede parametrizarse mediante las ecuaciones x = u, y = v y z = g(u, v). Para esta parame-
trizacion, nuestra ecuacion de area de superficie se reduce a
p
A(S) = 1 + (gu (u, v))2 + (gv (u, v))2 dudv
R
Ejemplo 227. Encuentre el area del cono r = (u cos v)i + (u sen v)j + uk, donde 0 u 1
y 0 u 2.
r r
Solucion. Lo primero es hallar a
u v
OBSERVACIONES
179
Si z = g(x, y) es la ecuacion de una superficie S, entonces
q
f (x, y, z)dS = f (x, y, g(x, y) 1 + [gx (x, y)]2 + [gy (x, y)]2 dA
S S
Suponga que (x, y, z) representa la densidad de una superficie S en el punto (x, y, z),
o la masa por unidad de area de superficie. Entonces la masa m de la superficie es
m= (x, y, z)dS
S
Solucion. No es dificil ver que (x, y, z) = y que la superficie en cuestion tiene una
ecuacion de la forma z = g(x, y) = 1 + x2 + y 2 Por tanto debemos calcular zx = y
zy = . Por tanto,
m= (x, y, z)dS =
S
180
Ejemplo 229. Evalue xz 2 dS donde S es la porcion del cilindro y = 2x2 +1 en el primer
S
octante acotado por x = 0, x = 2, z = 4 y z = 8.
Solucion. Observe que aqu la superficie S esta dada por y = g(x, z) = 2x2 + 1. Por tanto
necesitamos yx = y yz = . Por tanto,
xz 2 dS =
S
Definition 80 (Integral de Superficie parametrica).
Sea S una superficie parametrica suave definida por la funcion vectorial
donde (u, v) vara sobre la region R del parametro del plano uv y sea f (x, y, z) es continua
sobre S. Entonces la integral de superficie parametrica de f sobre S es
r r
f (x, y, z)dS = f (x(u, v), y(u, v), z(u, v))
dA
S R u v
p
Ejemplo 230. Evalue la integral de superficie 1 + x2 + y 2 dS, donde S es la superficie
S
parametrica dada por r(u, y) = u cos vi + u sen vj + vk, donde 0 u 2 y 0 v 4
r r
Solucion. Dado que S es parametrica debemos calcular . Las ecuaciones parametri-
u v
cas de S son:
x, y= z=
Por tanto, sustituyendo obtenemos
p
1 + x2 + y 2 =
r r
=
u v
Helicoide circular.
181
Orientacion de una Superfice:
Una superficie suave S es una superficie orientada si existe una funcion normal unitaria
continua n definida en cada punto sobre la superficie.
La cinta de Mobius NOOOO es una superficie orientada porque si una normal unitaria n
empieza en P sobre la superficie y se mueve una vez alrededor de la cinta sobre la curva
C, termina en el lado opuesto de la cinta en P y por ello apunta en la direccion opuesta.
Una superficie S definida por z = g(x, y) tiene una orientacion hacia arriba cuando
las normales unitarias estan dirigidas hacia arriba y tiene una orientacion hacia abajo
cuando las normales unitarias estan dirigidas hacia abajo.
Si una superficie suave S esta definida implcitamente por h(x, y, z) = 0 (es decir S es
la superficie de nivel 0 de h) entonces un campo normal de S es
h 1
n= = (hx , hy , hz )
khk khk
Si S esta definida por una funcion explcita z = g(x, y) entonces podemos usar h(x, y, z) =
z g(x, y) o h(x, y, z) = g(x, y) z dependiendo de la orientacion de S.
Una superficie cerrada se define como la frontera de un solido finito tal como la
superficie de una esfera.
182
En el caso de una superficie cerrada S, si n es la normal exterior (interior), entonces
(4.14) produce el volumen del fluido que fluye hacia fuera (hacia dentro) a traves de S
por unidad de tiempo.
Para evaluar una integral de flujo de una superficie dada por z = f (x, y) nos apoyamos
en una funcion G(x, y, z) = z g(x, y) para la cual la superficie S es una curva de nivel
cero de G. Esto implica que G = (fx , fy , 1) es una campo normal a S. Entonces,
G(x, y, z)
n dS = dS
kG(x, y, z)k
G(x, y, z)
q
=p (fx )2 + (fy )2 + 1 dA
(fx )2 + (fy )2 + 1
= G(x, y, z) dA
definida sobre una region R del plano uv, se puede definir la integral de flujo de F a
traves de S como
ru rv
Flujo = F n dS = F kru rv k dA (4.15)
S R kru rv k
Ejemplo 231. Considere que F(x, y, z) = zj+zk representa el flujo de un lquido. Determine
el flujo de F a traves de la superficie S dada por la parte del plano z = 6 3x 2y en el
primer octante orientado hacia arriba.
Solucion. Aqu la superficie (el plano) esta dado por por la superficie de nivel 0 de la funcion
h(x, y, z) = 3x + 2y + z 6. Por tanto el campo normas a la superficie (plano) es
h
n= =
khk
Luego
Flujo = (F.n)dS =
S
183
Ejemplo 232. Sea S la porcion del paraboloide z = 4 x2 y 2 que se encuentra sobre el
plano xy, orientado por medio de un vector unitario normal dirigido hacia arriba. Un fluido
de densidad constante fluye a traves de la superficie S de acuerdo con el campo vectorial
F(x, y, z) = xi + yj + zk. Hallar la tasa o ritmo de flujo de masa a traves de S.
Solucion. Aqu la superficie esta dado por la superficie de nivel 0 de la funcion G(x, y, z) =
. Por tanto el campo normal a la superficie es
n = G(x, y, z) =
Ejemplo 233. Hallar el flujo sobre la esfera S de radia a donde F es un campo cuadratico
inverso dado por
r
F(x, y, z) =
krk3
Supongase que S esta orientada hacia afuera.
donde 0 u y 0 v 2.
Por tanto el campo normal a la superficie es
n = ru rv =
Flujo = (F.n)dS = F (ru rv ) dA = 4
S R
184
Teorema de Stokes
Veamos que el teorema de Green podemos reescribirlo en forma vectorial tanto en R2
como en R3 .
Teorema de Green en R2 : Si F (x, y) = P (x, y)i + Q(x, y)j es un campo vectorial
bidimensional, entonces
i j k
Q P
rot F = F = = k
x y z x y
P Q 0
Ahora vamos a ver que el Teorema de Green en R3 relaciona una integral de lnea alrededor
de una curva C cerrada simple suave que forma la frontera de una superficie S con una integral
de superficie sobre S. es decir,
=
C S
Suponga que z = f (x, y) es una funcion continua cuya grafica es una superficie orien-
tada suave por partes sobre una region R en el plano xy
Por ultimo, sea T un vector tangente unitario a C que apunta en la direccion positiva.
La forma tridimensional del teorema de Green, la cual se presenta a continuacion,
recibe el nombre de Teorema de Stokes.
185
Teorema 82 (Teorema de Stokes).
Sea S una superficie orientada suave por partes acotada por una curva C cerrada simple
suave por partes. Sea
un campo vectorial para el cual P , Q, y R son continuas y tienen primeras derivadas parciales
continuas en una region abierta del espacio tridimensional que contiene a S. Si C se recorre
en la direccion positiva, entonces
F dr = (F T )ds = (rot F) n dS. (4.18)
C C S
Demostracion. Suponga que la superficie S se orienta hacia arriba y esta definida por una
funcion z = f (x, y) que tiene segundas derivadas parciales continuas. De la definicion
R Q P R Q P
rot F = i+ j+ k
y z z x x y
Por consiguiente,
" R
#
Q g P R g Q P
(rot F) n dS = + (4.19)
S R y z x z x y x y
Nuestro objetivo ahora es demostrar que la integral de lnea F dr se reduce a (4.19). Si
C
Cxy es la proyeccion de C sobre el plano xy y tiene las ecuaciones parametricas x = x(t),
y = y(t), a t b entonces ecuaciones parametricas para C son x = x(t), y = y(t),
186
z = g(x(t), y(t)),a t b. De tal manera
b" #
dx dy dz
F dr = P +Q +R dt
C a dt dt dt
b" #
dx dy g dx g dy
= P +Q +R + dt Regla de la cadena
a dt dt x dt y dt
g g
= P +R dx + Q + R dy
Cxy x y
" #
g g
= Q+R P +R dA T. Green (4.20)
R dx y dy x
2g 2g
Restando (4.22) de (4.21) y utilizando el hecho de que = vemos que, despues de
xy yx
rearreglar, (4.20) se convierte en
" #
R Q g P R g Q P
+
R y z x z x y x y
La ultima expresion es la misma que el lado derecho de (4.19), que era lo que se quera
demostrar
187
En este caso, si h(x, y, z) = z x2 1 define el cilindro para el nivel 0, la normal es
h 2xi + k
n= =
khk 4x2 + 1
Por tanto,
1 2
2xy x
(rot F n)dS = dS = (2xy x)dA = (2xy x)dydx = 2
S S 4x2 + 1 R 0 2
11 19
En consecuencia xydx + yzdy + xzdz = 0 +0 = 2. lo cual concuerda perfec-
C 15 15
tamente.
188
Ejemplo 235. Evalue zdx + xdy + ydz donde C es la traza del cilindro x2 + y 2 = 1 en
C
el plano y + z = 2. Oriente C en el sentido contrario al de las manecillas del reloj cuando
se observe desde arriba.
Solucion. Si F = zi + xj + yk entonces
i j k
rot F = =i+j+k
x y z
xy yz xz
En este caso definimos a h(x, y, z) = y + z 2 define el plano para el nivel 0, con orientacion
arriba y vector normal
h 1 1
n= = j+ k
khk 2 2
Por tanto,
" 1
#
1
F dr = (rot F n)dS = (i + j + k) j + k dS = 2 dS = 2 2ds = 2
C S S 2 2 S S
Interpretacion fsica del rotacional: Anteriormente vimos que si F es un campo de
velocidades de un fluido, entonces la circulacion de C F dr de F alrededor de C es una
medida de la cantidad por medio de la cual el fluido tiende a girar la curva C circulando
alrededor de ella. La circulacion de F se relaciona estrechamente a rot F. Para ver esto,
suponga que P0 (x0 , y0 , z0 ) es cualquier punto en el fluido y Cr es un pequeno crculo de
radio r centrado en P0 . Entonces por el teorema de Stokes,
F dr = (rot F) ndS (4.23)
Cr S
Ahora, en todos los puntos P (x, y, z) dentro del crculo pequeno Cr , si tomamos rot F(P )
rotF(P0 ), entonces (4.23) produce la aproximacion
F dr (rot F(P0 )) n(P0 )dS
Cr Sr
= (rot F(P0 )) n(P0 )) dS
Sr
= (rot F(P0 )) n(P0 )Ar (4.24)
1
(rot F(P0 )) n(P0 ) = lm F dr
r0 Ar C
r
189
De tal modo, vemos que la componente normal de rot F es el valor lmite del cociente entre
la circulacion de F y el area de la superficie circular. Para un valor pequeno pero fijo de r,
tenemos
1
(rot F(P0 )) n(P0 ) F dr (4.25)
Ar Cr
dy dx
Similarmente, si n = ( )i ( )j es una normal unitaria a C (verifique T n = 0),
ds ds
entonces (F n)ds puede expresarse en terminos de una integral doble que implique a
C
divF.
Del teorema de Green,
P Q i P Q i
h h
(F n)ds = P dy Qdx = dA = + dA
C C R x dy R x dy
Esto es,
(F n)ds = div FdA
C R
Teorema de Green en R3 : Es generalizar la ecuacion anterior al espacio tridimensional R3
y esta generalizacion es lo que llamaremos Teorema de la divergencia o de Gauss.
Este importante teorema da la relacion entre una integral
triple sobre una region solida D (en clase la llamamos pa-
pa) y una integral de superficie, donde la superficie S es la
frontera de D (frontera de D = S cascara de la papa). En
el enunciado del teorema, la superficie S es cerrada en el
sentido de que forma toda la frontera completa del solido
D. Ejemplos de superficies cerradas surgen de las regiones
limitadas o acotadas por esferas, elipsoides, cubos, tetrae-
dros, o combinaciones de estas superficies. Se supone que
D es una region solida sobre la cual se evalua una integral
triple, y que la superficie cerrada S esta orientada mediante
vectores normales unitarios dirigidos hacia el exterior, co-
mo se muestra en la figura. Con estas restricciones sobre S
y D, el teorema de la divergencia es como sigue.
190
Teorema 83 (Teorema de la Divergencia).
Suponga que D un solido acotado en el espacio tridimensional con una frontera suave por
partes S que esta orientada hacia arriba. Sea
un campo vectorial para el cual P , Q, y R son continuas y tienen primeras derivadas parciales
continuas en una region abierta del espacio tridimensional que contiene a D. Entonces
(F n)dS = (div F)dV. (4.26)
S D
Sobre S1 : Puesto que la normal hacia fuera apunta hacia abajo, describimos la superficie
como h(x, y, z) = g1 (x, y) z = 0. De tal modo
g1 g1
h x i + y j k 1
n= =r por lo que kn= r
khk
g1
2 2 2 2
1+ x + g
y
1
1+ g1
x + g1
y
191
Sobre S2 : Puesto que la normal hacia fuera apunta hacia abajo, describimos la superficie
como h(x, y, z) = z g2 (x, y) = 0. Por tanto,
g2 g2
i j+k
h x y 1
n= =r por lo que kn= r
khk g 2 g 2
2 2
g 2
1
g 2
1
1+ + 1+ +
x y x y
De la definicion de dS tenemos entonces
R(k n)dS = R(x, y, g2 (x, y))dA (4.29)
S2 R
Aunque demostramos (4.26) para una region especial D que tiene un lado vertical,
notamos que este tipo de region no se requiere en el teorema de l aDivergencia.
Una region acotada por una esfera o un elipsoide no tiene tampoco un lado vertical.
El teorema de la divergencia tambien se cumple para una region D acotada entre dos
superficies cerradas tal como las esferas concentricas Sa y Sb ; la superficie frontera S
de D es la union de Sa y Sb En este caso,
(F n)dS = div FdV
S D
se convierte en
(F n)dS + (F n)dS = div FdV
Sb Sa D
donde n apunta hacia fuera de D. En otras palabras, n apunta alejandose del origen
sobre Sb , pero n apunta hacia el origen sobre Sa .
Ejemplo 236. Sea D una region cerrada acotada por el hemisferio x2 + y 2 + (z 1)2 = 9,
1 z 4, y el plano z = 1 Verifique el teorema de la divergencia para el campo vectorial
F = xi + yj + (z 1)k.
La integral triple: Puesto que F = xi+yj+(z 1)k vemos que divF = 3. En consecuencia,
div FdV = 3dV = 3 dV = 54 (4.31)
Solucion. D D D
En el ultimo calculo, aprovechamos el hecho de que dV produce el volumen del hemis-
D
ferio.
192
La integral de superficie: Escribimos = + donde S1 es el hemisferio y S2
S S1 S2
es el plano z = 1. Si S1 es una superficie de nivel de h(x, y, z) = x2 + y 2 + (z 1)2 entonces
una normal unitaria que apunta hacia arriba es
h xi + yj + (z 1)k x y z1
n= =p = i+ j+ k
khk x2 + y 2 + (z 1)2 3 3 3
Ahora,
Fn=3
y por ello con la ayuda de coordenadas polares obtenemos
2 3
3
(F n)dS = 3 p dA = 9 (9 r2 )1/2 rdrd = 54
S1 R 9 x2 y 2 0 0
Sobre S2 , tomamos n = k de modo que Fn = z+1. Pero, como z = 1, (z+1)dS = 0.
S2
Por consiguiente, vemos que
(F n)dS = 54 + 0 = 54
S
Ejemplo 237. Evalue (F n)dS, donde S es el cubo unitario definido por 0 x 1,
S
0 y 1, 0 z 1 y F = xyi + y 2 zj + z 3 k.
Interpretacion fsica de la divergencia Vamos a suponer que P0 (x0 , y0 , z0 ) es cualquier
punto en el fluido y Sr es una pequena esfera de radio r centrada en P0 . Si Dr es la esfera
Sr y su interior, entonces el teorema de la divergencia produce
(F n)dS = divFdV (4.32)
Sr Dr
193