Cap 2 Calculo IVRGG

Captulo 2
Diferenciaci
on
de Funciones Vectoriales
de Varias Variables
En este captulo se establece de forma general y natural el concepto de diferenciaci

on
para funciones vectoriales de varias variables, aplicando la definicion correspondiente
a las distintas clases de funciones de este tipo, descritas en la Seccion 1.3, llevando
en cada contexto a un tipo de derivada especfico. Se introduce en cada caso
la notacion convencional utilizada, se da la interpretacion fsica y geometrica para
algunos de los casos y se considera la aplicacion potencial respectiva. Para ilustrar
la teora se ejemplifica con algunos casos en dimensiones menores. En el contexto de
los campos escalares se trabaja con la Regla de la Cadena y se da un bosquejo de la
demostracion del Teorema de Taylor. Ademas, se analizan de forma conjunta todos los
conceptos diferenciales para diferentes sistemas alternos de coordenadas, y finalmente,
se establece como un tema complementario uno de los resultados fundamentales del
Analisis: el Teorema de la Funcion Implcita y uno de sus principales corolarios,
el Teorema de la Funcion Inversa.
81
82 Captulo 2. Diferenciacio
n
2.1 El Concepto General de Diferenciaci

on
En esta parte se estudia el concepto de derivada, iniciando con el caso mas sencillo:
las funciones reales de variable real, y de forma natural se aborda teoricamente el caso
general de las funciones vectoriales de varias variables.
Una funcion f : [a, b] R R se dice ser diferenciable en x0 (a, b) si existe

0
un n
umero real, denotado por f (x0 ), tal que
f (x0 + x) f (x0 )
lim = f 0 (x0 ). (2.1.1)
x0 x
Una forma alternativa de escribir la existencia del proceso de lmite dado en (2.1.1) es
f (x) f (x0 )
lim = f 0 (x0 ), con x = x0 + x, lo cual es equivalente a
xx0 x x0
f (x) RT (x)
lim = 0, (2.1.2)
xx0 x x0
donde RT (x) = f (x0 ) + f 0 (x0 )(x x0 ) es la expresion de la recta tangente a la grafica

de f (x) en el punto (x0 , f (x0 )), con pendiente f 0 (x0 ) (ver Figura 2.1).
y
6 RT (x)

f (x0 ) y = f (x)

-
x0 x
Figura 2.1 Aproximaci
on Lineal de una Funcion Diferenciable (Suave)
Geometricamente, la igualdad (2.1.2) significa que la recta tangente RT (x) es cercana

a f (x) (es una aproximacion lineal, de hecho la mejor) en el sentido que la diferencia
entre f (x) y f (x0 ) + f 0 (x0 )(x x0 ) se aproxima mas rapidamente a cero que (x x0 ),
conforme x tiende a x0 . Precisamente esta es la nocion de buena aproximacion
que se adapta a las funciones vectoriales de varias variables. En otras palabras,
si L : R R es la funcion lineal L(h) = f 0 (x0 )h, entonces la relacion (2.1.2) establece
que
f (x0 + h) f (x0 ) L(h)
lim = 0, (2.1.3)
h0 h
obteniendo as la definicion formal de diferenciabilidad.

2.1. El Concepto General de Diferenciacion 83
on 2.1.1 Una funcion f : [a, b]

Definici R R es diferenciable en x0 (a, b)
si existe una funcion lineal L : R R que satisface la ecuacion (2.1.3). En este caso,
tal funcion lineal se define en terminos de la derivada de la funcion f en x0 , la cual,
ademas de estar relacionada con la pendiente de la recta tangente, es precisamente
el coeficiente de expansion-contraccion de la transformacion lineal existente.
Nota: i) La diferenciabilidad de una funcion es una operacion local. Se dice que una funci
on
es diferenciable en un conjunto abierto s es diferenciable en todo punto de tal conjunto.
ii) La relaci
on de aproximaci
on (2.1.3) ciertamente no tiene sentido en el caso general de una
funcion vectorial f : Rn Rm , pero se puede formular de tal forma que s lo tenga. Esto es,
la Definici
on 2.1.1 se generaliza de forma inmediata para funciones vectoriales definidas en
dimensiones superiores utilizando el concepto de distancia.
on 2.1.2 Una funcion f : Rn Rm es diferenciable en x0 , donde

o
Definici
o
{ } denota el interior del dominio considerado, si existe una transformacion lineal
L : Rn Rm tal que
kf (x0 + h) f (x0 ) L(h)k
lim = 0. (2.1.4)
h0 khk
Nota: i) Observese que h es un vector en Rn, mientras que [f (x0 + h) f (x0)] L(h)
es un vector en Rm, razon por la cual es necesario considerar el cociente de las normas
involucradas, representando de esta forma la cercana o medida de las variaciones en cada
uno de los espacios respectivos.
ii) Aunque parezca una redundancia, cuando n = m = 1 la Definicion 2.1.2 ciertamente es

la Definici
on 2.1.1. Sin embargo, como ejercicio, argumentar porque no hay necesidad de
utilizar el valor absoluto en el proceso de lmite en la definicion para el caso real.
Un resultado conocido de la teora de operadores establece que todo operador lineal

entre espacios vectoriales de dimension finita tiene una representacion matricial que
depende de las bases de los dos espacios involucrados (Teorema 1.4.2). Por lo tanto,
el operador lineal L citado en la Definicion 2.1.2 es una matriz de orden mn, la cual
es u
nica, afirmacion que se demuestra en el siguiente resultado.
Teorema 2.1.1 S f : Rn Rm es diferenciable en x0 , entonces la trans-

o
formacion lineal L : Rn Rm que aproxima al cambio instantaneo de sus valores es

u
nica.
n
Demostraci
on. Supongamos que P : Rn Rm es otra transformacion lineal que satis-
face (2.1.4). Entonces, utilizando la funcion auxiliar definida por G(h) = f (x0 +h)f (x0 ),
se tiene que
kL(h) P (h)k kL(h) G(h) + G(h) P (h)k
lim = lim
h0 khk h0 khk
kL(h) G(h)k kG(h) P (h)k
lim + lim = 0.
h0 khk h0 khk
Ahora, si x Rn , entonces tx 0, cuando t 0. As, para x 6= 0, se tiene que

kL(tx) P (tx)k kL(x) P (x)k
0 = lim = , lo cual implica que kL(x) P (x)k = 0 y
t0 ktxk kxk
por lo tanto, L(x) = P (x), para toda x Rn . Es decir, L P , demostrando as la
unicidad de la transformaci
on lineal L.
2.1.1 C
alculo Diferencial sobre Trayectorias
En el caso particular de las funciones vectoriales que dependen u

nicamente de una
variable real, : [a, b] R Rm, conocidas como trayectorias, el concepto de
diferenciacion se reduce al del calculo diferencial de funciones reales de variable real,
ya que, en este caso, la operacion se aplica componente a componente. El objetivo de
esta seccion es establecer el concepto de derivada de una trayectoria, su interpretacion
fsica y geometrica, y su principal aplicacion en la determinacion de la longitud de
arco, curvatura y torsion de una curva1 .
on 2.1.3 La derivada de una trayectoria, (t) = (f1 (t), f2 (t), , fm (t))T ,

Definici
con fi : [a, b] R, i = 1, 2, , m, respecto al parametro t, en t = t0 se define,
a partir de la Definicion 2.1.2, entrada por entrada y se denota por

d def 1
0 ) (t)
(t = lim [(t0 + t) (t0 )]
dt t=t0
t0 t
T (2.1.5)
d d d
= f1 (t0 ), f2 (t0 ), , fm (t0 ) ,
dt dt dt
siempre que la operacion en cada una de las componentes este bien definida o exista,
d
donde () = (notacion que se debe a Sir Isaac Newton). Por lo tanto, es
dt
diferenciable en t0 (a, b) si y solo si lo son cada una de las funciones componentes
fi , i = 1, 2, . . . , m.
1
El an
alisis detallado de este tema corresponde a la Geometra Diferencial.
En particular, en el caso cuando la trayectoria esta en el espacio, (t) = (x(t), y(t), z(t))T
es diferenciable en (a, b) s el lmite
1
0 ) = lim
(t [(t0 + h) (t0 )]
h0 h
1
= lim [(x(t0 + h), y(t0 + h), z(t0 + h))T (x(t0 ), y(t0 ), z(t0 ))T ]
h0 h
z(t0 + h) z(t0 ) T

x(t0 + h) x(t0 ) y(t0 + h) y(t0 )
= lim , lim , lim
h0 h h0 h h0 h
existe para todo t0 (a, b). En tal caso, (t

0 ) = (x(t
0 ), y(t 0 ))T .
0 ), z(t
Ejemplo 24: a) (t) = (4t2 , tet , 12 sen t)T

(t) = (8t, (1 t)et , 2 cos t)T .
b) Si (t) = (ln(2 t2 ), |t 1|, 7t2/3 1)T , entonces (1)

= (2, @, 14/3)T , lo cual

implica que (1) no existe (basta con que una de las funciones componentes no sea
diferenciable para que la funcion vectorial tampoco lo sea).
2
Interpretaci
on Geom
etrica de la Derivada de una Trayectoria: La variacion
en los valores de la funcion, (t0 +h)(t0 ), es un vector anclado en el origen equipo-
lente al segmento de recta dirigido que va de p a q,
z
q = (t0 + h)
6
(t0 ) = p

1
: [(t0 + h) (t0 )]

h
-
y

x

1
de tal forma que el cambio instantaneo [(t0 + h) (t0 )] tambien es paralelo

h
a pq. Si h tiende a cero, el vector (t0 + h) se acerca cada vez mas al vector (t0 )
y la lnea recta que une a p con q (recta secante) se parece cada vez mas a la recta
tangente a la curva descrita por en el punto p = (t0 ). En el lmite, el vector
1
0 ) = lim [(t0 + h) (t0 )]
(t
h0 h
es paralelo al vector tangente a en p.

n
0 ) es paralelo a la recta tangente RT ( ) de la curva descrita

Por lo tanto, el vector (t
por en el punto (t0 ), cuya expresion esta dada por
0 ),
RT ( ) = (t0 ) + (t R, (2.1.6)
la cual representa la mejor aproximacion lineal a la trayectoria (t) cerca del punto
(t0 ). Geometricamente, se tiene que
RT ( )
z
6
p (t)

0)
(t

-
y

x
Nota: Al igual que en el caso de las funciones reales de variable real, tambien para las
trayectorias la diferenciabilidad implica continuidad, pero el recproco es falso.
Ahora, s se supone que la imagen de (t) = (x(t), y(t), z(t))T es la posicion de una
partcula en movimiento cuando t = tiempo transcurre, entonces

v(t) = (t) representa la velocidad,
r(t) = k(t)k
la magnitud del vector velocidad, representa
la rapidez de la partcula a lo largo de la curva en el tiempo t.

a(t) = (t) representa la aceleracion de la partcula en movimiento, donde

(t) x(t), y(t), z(t))T .
= (
4
Ejercicio 45. Sea (t) = (a cos t, a sen t, kt)T , con a, k constantes y t [0, ].

i) Trazar la curva descrita por .
ii) Calcular la velocidad para cualquier tiempo t y en particular
para t0 = 2 segundos.
iii) Calcular la rapidez.
iv) Calcular la aceleracion.
v) Determinar la expresion de la recta tangente en el punto (t0 ).

Ejercicio 46. Determinar (t)
y (0) para
T
a) (t) = ( sen 2t, cos 2t, 2t t2 )
T
b) (t) = (t2 cos(t2 + ), ln (1 + t), tet )
Ejercicio 47. Calcular los vectores velocidad y aceleracion as como la ecuacion de

la recta tangente en el punto especificado en cada caso
T
a) (t) = 2t, et , et , t0 = 0
3
T
b) (t) = sen (t 1), cos(t 1), 2t 2 , t0 = 1
Propiedades B
asicas de la Derivada de Trayectorias:
1). La derivada del producto de una funcion escalar s(t) y una trayectoria (t)
satisface la familiar regla del producto (Regla de Leibniz):
d

[s(t)(t)] = s(t)(t) + s(t)(t).
(2.1.7)
dt
2). La Regla de Leibniz tambien es valida para los otros productos definidos.
Por ejemplo, s 1 , 2 y 3 son trayectorias en R3, entonces
d
h1 (t), 2 (t)i = h1 (t),
2 (t)i + h
1 (t), 2 (t)i. (2.1.8)
dt
d
[1 (t) 2 (t)] = 1 (t) 1 (t) 2 (t).
2 (t) + (2.1.9)
dt
d
h1 (t), 2 (t) 3 (t)i = h
1 (t), 2 (t) 3 (t)i
dt
2 (t) 3 (t)i
+h1 (t), (2.1.10)
+h1 (t), 2 (t) 3 (t)i.
Nota: En la relaci
on (2.1.8) el orden de los factores en los terminos del lado derecho
no tiene importancia y de hecho se puede extender para cualquier par de trayectorias
1,2 : [a, b] R Rm .
n
3). El vector velocidad de una trayectoria de rapidez constante, es perpendicular

al vector aceleracion. Es decir, (t) (t)
(la velocidad es perpendicular a la
aceleracion).
Demostraci on. Sea (t) de rapidez constante. Esto es, k(t)k = k = constante.
2 d d
2
Entonces, k = k(t)k = h(t),

(t)i. As, por un lado h(t),
(t)i = k 2 = 0,
dt dt
d
y por el otro h(t),
(t)i
= 2h(t),
(t)i,
lo cual es compatible solo si (t) es el
dt

vector cero o si (t)
es perpendicular a (t).
4). An
alogo de la Regla de la Cadena: Si t = (s) es una funcion real
diferenciable y (t) es una trayectoria diferenciable, entonces
d d d
= . (2.1.11)
ds dt ds
Ejercicio 48. Demostrar las relaciones (2.1.7) (2.1.11), Propiedades 1), 2) y 4).
Longitud de Arco: La longitud de arco de una curva es un parametro natural,

com
unmente denotado por s, que establece sobre la curva una propiedad intere-
sante, las curvas reparametrizadas con su longitud de arco son de rapidez unitaria.
La formula para calcular la longitud de una curva dada se aproxima con la longitud
de una poligonal que cada vez se hace mas fina, lo cual lleva a un proceso de lmite
de sumas de Riemann y por consiguiente a la evaluacion de una integral real definida,
proceso que se resume en la siguiente definicion.
Definici on 2.1.4 Sea : [a, b] R R3 una trayectoria suave2 en el espacio,

con rapidez k(t)k.
Entonces, se define la longitud de la curva , comprendida entre
los puntos (a) y (b), por
Z b Z b q
2
2 + [z(t)]
2 dt.
def
` = `[] = k(t)kdt
= [x(t)]
+ [y(t)] (2.1.12)
a a
z
6
(b)
(a)

`
-
y

x
2 6= 0.
Que la curva sea suave significa que es diferenciable y que
b
S la curva esta en R , entonces (2.1.12) tiene la expresion ` =
Z q
2
2
[x(t)]
2 dt,
+ [y(t)]
Rm, se tiene que
a
mientras que en el caso general, cuando la curva esta en
Z b Z b
rh i2 h i2 h i2
`= k(t)kdt
= f1 (t) + f2 (t) + + fm (t) dt. (2.1.13)
a a
Ejemplo 25: i) La longitud de una circunferencia de radio r es igual a 2r.

El resultado se obtiene utilizando cualquier parametrizacion de la circunferencia dada. As,
sin perdida de generalidad
y
sea (t) = (r cos t, r sen t)T , t [0, 2]. Entonces,
6
'$
r Z 2 Z 2 p
-
x `= k(t)kdt
= r2 cos2 t + r2 sen2 t dt = 2r.
0 0
&%
ii) Determinar la longitud de un arco de la cicloide.
Como todos los arcos de la cicloide son iguales, es suficiente considerar el arco principal
descrito por (t) = (at a sen t, a a cos t)T , a > 0, 0 t 2.
y
6
-
0 2a x

Entonces, el vector velocidad de es (t) = (a a cos t, a sen t)T y la rapidez
q
k(t)k
= (a a cos t)2 + a2 sen2 t = a2 2a2 cos t + a2 cos2 t + a2 sen2 t
s
p t
= 2a2 (1 cos t) = 2a2 2 sen2 .
2
t
Como t [0, 2], se tiene que k(t)k
= 2a sen y por lo tanto,
2
Z 2 Z

t
l[] = 2a sen dt = 4a sen u du = 4a( cos u) = 4a(cos + cos 0).
0 2 0 0
As, l[arco de la cicloide] = 8a.

2
n
x2 y 2
Ejercicio 49. Determinar la longitud de la elipse + 2 = 1.
a2 b
Ejercicio 50. Determinar la longitud de la astroide (hipocicloide) del Ejercicio 23.
y
6 2 2 2
a x3 + y 3 = a3

a
-
a x
A
a AU
(t) = (a cos3 t, a sen3 t)T , 0 t 2
" Z
#
2
Sugerencia: ` = cuatro veces la longitud de cada arco = 4 k(t)k
dt.
0
Ejercicio 51. a) Determinar la trayectoria tal que (1) = (1, 0, 1)T y

(t) = (1, 2t 2, cos t t sen t)T .
b) S una partcula se mueve sobre la curva (t) y en el tiempo t = 1 se sale por

la tangente donde estara en el tiempo t = 2?
Ejercicio 52. Calcular la longitud de la curva dada en el intervalo especificado:

3 T

i) (t) = sen 3t, cos 3t, 2t 2 , [0, 1].
ii) (t) = (t, t sen t, t cos t)T , [0, ].
Ejercicio 53. Una parametrizacion de la figura del ocho acostado o smbolo del
infinito, llamada Lemniscata de Bernoulli esta dada por
y
6 T

R
1 (t) = 2 cos t 2 , sen 2 t 2 , 0 t 2
-
2 R 2 x
1
a) Determinar la ecuacion de la recta tangente a la curva que forme con el eje x

un angulo de 45 .
b) Determinar la longitud de la curva dada.

on 2.1.5 Sea h : I1 I una funcion C1 y sea : I

Definici Rm una trayec-
toria C1 por piezas. Entonces, a la composicion = h : I1 Rm se le llama
reparametrizacion3 de . Esto es, ( ) = (h( )), donde t = h( ).
La funcion longitud de arco o parametro natural,

Z t
s = s(t) = k()k
d, (2.1.14)
0
ds
es una funcion continua creciente del parametro t, ya que = k(t)k
y lo mismo
dt
es cierto para la funcion inversa t = (s), por lo que es completamente valido describir
una curva mediante = (s) para 0 s `, donde ` es la longitud total de la curva
considerada.
Teorema 2.1.2 S : I Rm es de clase C1 y s(t) es el parametro natural, entonces

la reparametrizacion de con s es de rapidez unitaria.
Demostraci
on. Se deja como ejercicio.
Ejercicio 54. Demostrar el Teorema 2.1.2
Ejercicio 55. Calcular la longitud de la helice circular (t) = (a cos t, a sen t, kt)T ,
t 0, 4

. Reparametrizar con su longitud de arco y verificar por calculo directo
el Teorema 2.1.2.
De acuerdo al Teorema 2.1.2, el vector tangente unitario a la curva esta dado por
d

T= = . (2.1.15)
ds kk

Entonces, por la Propiedad Basica 3) de la Derivada de una Trayectoria, el vector

normal unitario esta dado por
d2
2
N = ds
2 , (2.1.16)
d
k 2k
ds
3
Es importante observar que existen reparametrizaciones que preservan la orientacion de la
curva y existen tambien reparametrizaciones que invierten la orientacion. Por ejemplo, una
on directa que invierte la orientacion es h( ) = a + b , [a, b].
reparametrizaci
n
mientras que el vector binormal unitario esta dado por

d d2
2
B=TN= ds ds = , (2.1.17)
2
d k
k

k 2k
ds
2
3
d d dt
ya que 2 = (
)
. Por consiguiente,
ds ds ds
)
( 2
kk (
)

N=BT= = . (2.1.18)
k
k
kk
kk
k k

La terna de vectores ortonormales LI (T, N, B), llamada triedro intrnseco o aparato

de Frenet, constituye un sistema movil de referencia a lo largo de la curva .
El plano generado por T y N se llama plano osculante, mientras que el plano generado
por T y B se llama plano rectificante y el plano generado por N y B se llama plano
normal.
Ahora, para curvas de clase C3 , la rapidez instantanea de cambio del vector tangente
dT d2 d2
unitario respecto a la distancia recorrida a lo largo de la curva, = 2 =k 2kN
ds ds ds
proporciona una medida de la curvatura4 de una curva en cada uno de sus puntos.
d2
Introduciendo la notacion = k 2 k, esta derivada se escribe como
ds
dT
= N. (2.1.19)
ds
En terminos de los tres vectores LI T, N, B se tiene que

dN dB
= T + N + B, = T + N + B,
ds ds
de tal forma que al efectuar el producto interior con T, N, B lleva a
dN dN dN
=T , =N , =B ,
ds ds ds
dB dB dB
=T , =N , =B .
ds ds ds
4
La curvatura de una curva mide que tanto se dobla la curva en cada punto o que tanto la
curva se desva de ser recta.
Las propiedades del triedro (T, N, B) implican que
dT dN dB
T = 0, N = 0, B = 0,
ds ds ds
dN dT dB dT db dN
T = N , T = B , N = B ,
ds ds ds ds ds ds
de donde = = 0 y = , cuyo valor com un se denota por y representa la

5 dB
torsion de la curva en cada punto. Esto es, = = = N. Ademas,
ds
dT dT
se tiene que N = , B = 0 y por consiguiente, = , = 0.
ds ds
Ejercicio 56. Mostrar que la curvatura de una lnea recta es cero, mientras que la
curvatura de una circunferencia de radio a es constante igual a 1/a.
As, la evolucion del triedro de Frenet conforme se mueve sobre la curva satisface
la dinamica gobernada por el sistema establecido en el siguiente resultado.
Teorema 2.1.3 Las derivadas con respecto al parametro longitud de arco de los
vectores del triedro de Frenet satisfacen el sistema
dT

= N,

ds

dN

= T + B, (2.1.20)

ds

dB

= N,

ds
el cual se conoce como las formulas de Frenet-Serret6 . Los coeficientes y son
llamados la curvatura y la torsion de la curva, respectivamente.
5
La torsi
on es una medida del giro o enroscamiento de una curva en cada punto.
6
En el Siglo XIX, los matem aticos franceses Jean Frederic Frenet (1816-1900) y Joseph Alfred
Serret (1819-1885) fueron los primeros en formalizar el concepto de curvatura de una curva en
cada punto y su inverso = 1/, el radio de curvatura, as como la torsi
on de tal curva.
n
Nota: Las ecuaciones de Frenet-Serret (2.1.20) forman un sistema de tres ecuaciones

diferenciales vectoriales ordinarias de primer orden en T, N y B, por lo que es natural
preguntarse si dadas dos funciones escalares reales y existen las soluciones T, N y B del
sistema de Frenet-Serret (2.1.20) y por tanto, si existe una curva que tenga la curvatura y
la torsion especificadas. La respuesta es afirmativa y se demuestra en Geometra Diferencial
(Teorema de Existencia y Unicidad de Curvas en el Espacio).
Ejercicio 57. Demostrar que la curvatura y la torsion de una curva parametrizada

...
k
k
h,

i
por (t) estan dadas por = y = .
3
kk k 2
k
Ejercicio 58. Determinar la curvatura y la torsion de la helice
(t) = (a cos t, a sen t, bt)T , a, b > 0.
Ejercicio 59. Si a es una constante diferente de cero, calcular la curvatura y la

torsion de la curva descrita por (t) = (a(3t t3 ), 3at2 , a(3t + t3 ))T .
Ejercicio 60. Mostrar que la curvatura de una curva explcita en el plano, y = f (x),
|f 00 (x)|
esta dada por (x) = .
[1 + (f 0 (x))2 ]3/2
2.1.2 C
alculo Diferencial sobre Campos Escalares
En la Subseccion 1.3.2 se define y se proporcionan algunos ejemplos de funciones

reales de varias variables, f : Rn R, llamadas campos escalares reales, cuyas
graficas son hipersuperficies en Rn+1. La Definicion 2.1.2 establece que significa que
o
xn+1 = f (x1 , x2 , , xn ) = f (x) sea diferenciable en x0 = (x01 , x02 , , x0n )T .
Sin embargo, es mas conveniente e ilustrativo llegar a este concepto utilizando solo la
nocion de derivada parcial, ya que precisamente en terminos de las derivadas parciales
se definen los dos conceptos diferenciales mas importantes de los campos escalares:
el de vector gradiente y el de derivada direccional.
on 2.1.6 Sea f :
Definici Rn R un campo escalar real de n variables.
Las derivadas parciales de f con respecto a la primera, segunda,..., nesima variable,
f f f
denotadas por , , , , tambien son funciones reales de n variables
x1 x2 xn
definidas en el punto x = (x1 , x2 , , xn )T
e i por
f f (x1 , x2 , , xi + h, , xn ) f (x1 , x2 , , xn )
(x) = lim
xi h0 h (2.1.21)
f (x + hni ) f (x)
= lim ,
h0 h
si el lmite existe, donde i = 1, 2, , n y ni es el iesimo vector de la base natural

de Rn (ver relacion (??)).
f
Nota: Es importante se
nalar que es la derivada ordinaria de la funcion f con respecto
xi
f
a la variable xi , manteniendo las otras (n 1) variables fijas. Esto es, para calcular (x0 )
xi
se deriva f (x01 , x02 , , xi , , x0n ) = g(xi ) en la forma ordinaria respecto a la variable
xi ,
f dg
considerando como constantes al resto de las variables. Es decir, (x0 ) = .
xi dxi xi =x0i
De aqu el nombre derivada parcial. Por ejemplo,
f f (x + x, y) f (x, y)
fx (x, y) = lim ,
x x0 x
para (x, y) en el interior del dominio de f
donde el lmite existe
f f (x, y + y) f (x, y)
fy (x, y) = lim .
y y0 y
n
De forma similar a la interpretacion que tiene la derivada de una funcion de una

f
variable, (x0 ) indica la rapidez de cambio de la funcion f en el punto x0 en la
xi
direccion del vector ni . Esto se debe a que el incremento de f en el intervalo de tiempo
[0, t] en esa direccion es f (x0 + tni ) f (x0 ) y que la velocidad media del incremento
f (x0 + tni ) f (x0 )
en dicho intervalo sea . Por lo tanto, las derivadas parciales de
t
un campo escalar f miden la razon de cambio de la funcion en la direccion de cada
uno de los ejes coordenados. En el caso de una funcion real de dos variables, la cual
representa a la superficie explcita z = f (x, y), las derivadas parciales proporcionan
las pendientes de las rectas tangentes a las curvas secciones de f (x, y) paralelas a los
planos coordenados xz y yz, respectivamente.
f f
Ejemplo 26: a) Las derivadas parciales y de la funcion real de dos variables
x y
y f y f x
f (x, y) = tan1 son = 2 y = 2 .
x x x + y2 y x + y2
f
Observaci
on: i) Al derivar las funciones , i = 1, 2, , n (cuando son derivables en
xi
el sentido de la Definici
on 2.1.6), se obtienen las derivadas parciales de segundo orden:
2f 2f 2f
, , , ,
x21 x2 x1 xn x1
2f 2f 2f
, , , ,
x1 x2 x22 xn x2 (2.1.22)
.. .. ... ..
. . .
2f 2f 2f
, , , .
x1 xn x2 xn x2n
ii) S las derivadas parciales de primer orden son continuas, entonces el arreglo (2.1.22)
es simetrico. Esto es,
2f 2f
= i, j. (2.1.23)
xi xj xj xi
Es decir, se tiene igualdad en las parciales mixtas7 .
7
El problema es cuando no hay continuidad de las primeras derivadas parciales, porque entonces
cabe la posibilidad de que las parciales mixtas no existan o tomen valores diferentes.
y
1
b) Al calcular las segundas derivadas parciales de la funcion f (x, y) = tan ,
x
se obtiene
2f

y 2xy
= = ,
x2 x x + y2
2
(x2 + y 2 )2
2f y 2 x2 2f

x y
= = 2 = y = ,
xy x x2 + y 2 2 2
(x + y ) x2 + y 2 yx
2f

x 2xy
= = .
y 2 y x + y2
2
(x2 + y 2 )2
Definici
on 2.1.7 Ademas de las derivadas parciales de primer y segundo orden,
se definen las derivadas parciales de tercer orden, . . . , derivadas parciales de orden
superior (cuando es posible en terminos de la Definicion 2.1.6) como

( derivada parcial previa ) . (2.1.24)
xi
c) Las primeras y segundas derivadas parciales de la funcion norma euclideana en Rn
con respecto a cada una de las variables xi , respectivamente, son
n
xi 2 1 x2i kxk2 x2i 1 X 2
kxk = y kxk = = = xk ,
xi kxk x2i kxk kxk3 kxk3 kxk3 k6=i
2 /4t
d) Sea v(r, t) = tn er . Determinar el
valor de
la constante n tal que v satisfaga la
v 1 v
ecuacion diferencial parcial = 2 r2 .
t r r r
Soluci
on.
r2 r2 r2

v 2 2 n1 r2 /4t n2 r2 /4t
= n tn1 er /4t + tn er /4t =t e n+t 2 =t e nt + .
t 4t2 4t 4

v n r2 /4t 2r v 1 2
=t e r2 = r3 tn1 er /4t
r 4t r 2

v 3 2 1 2 2r
r2 = r2 tn1 er /4t + r3 tn1 er /4t
r r 2 2 4t
r2

3 2 1 2 2 3
= r2 tn1 er /4t + r4 tn2 er /4t = r2 tn2 er /4t t + .
2 4 2 4

v 1 2 v 3
Por lo tanto, = 2 r n= .
t r r r 2
n
e) Las primeras y segundas derivadas parciales, respecto a cada una de las variables,
n X n
de la siguiente forma cuadratica en R , f (x) = x Ax =
X
n T
aij xi xj , donde A
i=1 j=1
es una matriz simetrica de orden n, respectivamente, son
n
f X 2f
=2 aij xj y = 2aij .
xi xj xi
j=1
Definici
on 2.1.8 El Operador de Laplace (o simplemente Laplaciano) ndimensional
se define por
n
2 2 2
def 2 X 2
= + + + = . (2.1.25)
x21 x22 x2n i=1
x2i
f ) S u(x) = eha,xi , x, a Rn , con kak = k. Demostrar que la funcion u satisface la

Ecuacion de Helmholtz ndimensional u k 2 u = 0.
u 2u
Soluci
on. = eha,xi ai y = eha,xi a2i , as que
xi x2i
n n n
X 2 u X ha,xi 2 ha,xi
X
2 = e a i = e a2i = ukak2 = k 2 u.
i=1
x i i=1 i=1
g) Si f (x) = kxk(2n) , 0 6= x Rn y R. Demostrar que la funcion f satisface

la Ecuacion de Laplace ndimensional u = 0.
2n
on. Como f (x) = x21 + x22 + + x2n
Soluci 2
, entonces

f 2n n
= (2xi ) x21 + x22 + + x2n 2
= (2 n) xi kxkn
xi 2
2f n n n2
= (2 n) x21 + x22 + + x2n (2 n) 2x2i x21 + x22 + + x2n

2
2
x2i 2
n 2 n2
= (2 n)kxk n (2 n)xi kxk .
As,
n
" n #
X 2f n2
X 2 2

n2

2 2

2 = (2 n)kxk kxk n x i = (2 n)kxk nkxk nkxk ,
i=1
xi i=1
de donde se tiene que f = 0.

2
Ejercicio 61. En cada uno de los siguientes ejercicios, calcular todas las derivadas
de primer orden:
1
a) f (x, y) = cos x2 , y 6= 0,
y

1 x+y
b) f (x, y) = tan ,
1 xy
loga (x2 + y 2 + z 2 )
c) f (x, y) = ,
x2 + y 2 + z 2
x2 z 2 2 3 2
+ 4x1/2 y 1/3 z 1/4 ex y z cos x2 4y 2 8z 2 ,

d) f (x, y, z) = 2 2
x y
w x+y+z 2 2
z
e) f (x, y, z, w, t) = sen e w2 t2 t exy +w + 7cos1 ,
t x
f ) f (x) = ln kxk4 , x Rn .

Ejercicio 62. Un automovil N viaja por una carretera hacia el Norte a 90 km/h,
mientras que otro automovil T viaja por una carretera hacia el Oeste a 80 km/h.
Cada auto se acerca a la interseccion de estas carreteras. Que tan rapido cambia
la distancia que hay entre los autos cuando el auto N esta a 0.3 km del punto de
interseccion, mientras que el auto T esta a 0.4 km de tal punto de cruce?
Ejercicio 63. Verificar que zxy = zyx , donde
a) z = 2x3 + 5x2 y 6y 2 + xy 4 ,
b) z = sen xy + tan1 xy,
c) tambien calcular zxx y zyy en ambos incisos.
n
!
2u X
Ejercicio 64. Calcular de la funcion u = sen kxk .
x2i k=1
yz
p
Ejercicio 65. S f (x, y, z) = xx + ln (x + z) tan1 [cos xy ln (x + y + z)].
f
Determinar (1, 2, 0).
y
n

xy (x2 y 2 )
si (x, y) 6= (0, 0),

x2 + y 2

Ejercicio 66. Sea f (x, y) =

0 si (x, y) = (0, 0).
f f
a) S (x, y) 6= (0, 0), calcular y .
x y
f f
b) Mostrar que (0, 0) = 0 = (0, 0).
x y
2f 2f
c) Mostrar que (0, 0) = 1 y (0, 0) = 1.
xy yx
Cual es el error? Porque las parciales mixtas no son iguales?
z z
Ejercicio 67. Calcular y de las siguientes funciones definidas implcitamente:
x y

x z p
a) = z xy xy 2 z 3 + ln x2 + y 2 + z 2 .
z y
x 4x 8y 2 2

b) e z cos(x + y + z) + xyz sec x yz = 0.
x2 y 2 + z 2
Ejercicio 68.
2 k2 t
i) Verificar que la funcion u = e sen kx es una solucion de la ecuacion de
2
calor 1dimensional, ut = uxx .
ii) Mostrar que u(x, t) = sen(x at) + ln (x + at) satisface la ecuacion de onda
1dimensional, utt = a2 uxx .
1 x2 +y 2
iii) Mostrar que u(t, x, y) = e 4kt satisface la ecuacion de calor 2dimensional,
2 4kt
u u 2u
=k + y tambien u(t, x, y) = exp [(m2 + n2 )kt] sen mx cos ny
t x2 y 2
para cualquier eleccion de m y n.
iv) Demostrar que s f : Rn R satisface la ecuaci

on de
Laplace ndimensional,
x
f = 0, entonces tambien la funcion (x) = f es una solucion.
kxk
2u
v) Dada z = u(x, y)eax+by , con = 0, encontrar valores de las constantes
xy
2z z z
a y b tales que + z = 0.
xy x y
El Vector Gradiente: S f es un campo escalar de n variables, entonces el valor

f f f
de las derivadas parciales , , , en un determinado punto represen-
x1 x2 xn
tan, respectivamente, las razones de cambio de f con respecto a la distancia en las
direcciones de los vectores de la base natural, n1 , n2 , , nn . Considerar el vector
f 1 f 2 f n
v= n + n + + n , (2.1.26)
x1 x2 xn
cuyas componentes escalares en las direcciones n1 , n2 , , nn en cualquier punto p
sean precisamente las respectivas razones en las cuales f cambia con respecto a la
distancia en esas direcciones en p. Para determinar la componente de v en cualquier
otra direccion en p, se considera al vector r como la posicion del punto p y se indica
un desplazamiento diferencial desde p en cualquier direccion por
dr = dx1 n1 + dx2 n2 + + dxn nn . (2.1.27)
Geometricamente,
6 p(x1 , x2 , , xn )
R dr
@
@
r
-

La componente escalar de v en esa direccion especfica se obtiene mediante el producto

dr
interior de v y el vector unitario u = , donde ds = kdrk, dando como resultado
ds
f dx1 f dx2 f dxn df
hv, ui = + + + = . (2.1.28)
x1 ds x2 ds xn ds ds
As, la componente de v en cualquier direccion es la razon de cambio de f con respecto

a la distancia en esa direccion y es importante se
nalar que v es una funcion vectorial
que esta asociada con la funcion escalar f en una forma independiente del sistema
coordenado. El vector v es llamado el vector gradiente de f .
on 2.1.9 S f : Rn R es diferenciable en x0 , entonces el vector

o
Definici
gradiente de f en x0 = (x01 , x02 , , x0n )T es el vector en el espacio dominio Rn
definido por
f f f
grad(f ) = , , , , (2.1.29)
x1 x2 xn x=x0
en coordenadas cartesianas.
n
Tambien se denota al vector gradiente por f (x)|x=x0 o simplemente por f (x0 ),

el cual por definicion (convencion) es un vector renglon, donde precisamente f (x0 )
define a la transformacion lineal a la que hace referencia la Definicion 2.1.2. Esto es,
f
L(h) = hf (x0 ) , hi. Otra notacion frecuentemente utilizada es (x0 ).
x
El operador , llamado operador del, representa a un operador diferencial vectorial
de la forma

1 2 n
= , , , = n + n + + n , (2.1.30)
x1 x2 xn x1 x2 xn
en coordenadas rectangulares, ya que tiene diferente representacion en otros sistemas

de coordenadas (lo cual se establece en la Seccion 2.3), de tal forma que al aplicar
a un campo escalar diferenciable f , se obtiene el gradiente del campo escalar dado,

1 2 n f 1 f 2 f n
f = n + n + + n f= n + n + + n .
x1 x2 xn x1 x2 xn
En general, s r representa al vector posicion de un punto p y s dr representa un

desplazamiento diferencial de longitud ds desde p, entonces la ecuacion (2.1.30) provee
la siguiente caracterstica fundamental
dr df
h , gradf i = (2.1.31)
ds ds
o en terminos diferenciales
df = hf, dri. (2.1.32)
Ejemplo 27: Calcular el gradiente de la funcion norma euclideana en Rn.

Soluci
on. En el Ejemplo 26 c) se determinaronlas primeras derivadas parciales para
f f f 1
esta funcion, as que f (x) = , , , = xT . En este caso, el gradiente
x1 x2 xn kxk
on del vector x. Esto es, kf (x)k = 1.
es un vector unitario en la direcci 2
Nota: S f es un campo escalar en el espacio 3dimensional, f : R3 R , es decir, si f

es una funci
on escalar de las tres variables rectangulares x, y, z, entonces es com
un tambien
f f f
que el gradiente se denote por f = i+ j+ k, el cual resulta ser el vector que
x y z
apunta en la direcci
on en la cual la razon de cambio de los valores de f es n umericamente
mas grande. Es decir, f proporciona la direccion de maximo crecimiento de la funcion f .
La demostraci
on de esta afirmaci
on es casi trivial, pero importante, la cual se incluye mas
adelante.
xz yz
Ejemplo 28: i) Sea f (x, y, z) = ex+y + sen xz + . Entonces,
x2 + y 2
f (1, 0, 0) = (e, e, 2).
ii) S la temperatura en cualquier punto p(x, y, z) es directamente proporcional a

la distancia al cuadrado desde el origen, entonces T = c (x2 + y 2 + z 2 ) y el gradiente
tiene la forma T = 2c(xi + yj + zk), un vector radial que sale del origen para c > 0
(calentamiento), y por lo tanto, perpendicular a las superficies esfericas equitermales,
x2 + y 2 + z 2 = constante.
p
La magnitud del gradiente esta dada por kT k = 2cr, donde r = x2 + y 2 + z 2 .
dT
Como T = cr2 , entonces tenemos que kT k = , la cual es la derivada de la
dr
temperatura en la direccion radial, que en este caso representa la razon maxima de
crecimiento de la funcion temperatura.
2
Recordemos que las curvas de nivel de una funcion real de dos variables z = f (x, y)
son curvas sobre el plano xy de la forma f (x, y) = c, con c una constante. As,
por cada punto (x0 , y0 ) D = Dom(f ) pasa una curva de nivel de f , a saber,
la curva f (x, y) = f (x0 , y0 ). El vector gradiente, f |x=x0 , es normal a la curva
de nivel correspondiente en el punto x0 = (x0 , y0 ).
* f

(x0 , y0 )
f = c1
f = c2
..
.
f =c
La demostracion de esta afirmacion queda

incluida en el caso mas general de funciones de
tres variables, que se demostrara mas adelante
y este a su vez facilmente se generaliza cuando
la dimension n > 3.
n
x2
Ejemplo 29: y g(x, y) = x2 + 2y 2 . Entonces, las curvas de
S f (x, y) =
y
nivel de f cortan perpendicularmente a las curvas de nivel de g. En efecto, para todo
punto (x0 , y0 ) Dom(f ) Dom(g), se tiene que
2x0 x20

f (x0 , y0 ) = , 2 , g(x0 , y0 ) = (2x0 , 4y0 )
y0 y0
x20

2x0
hf, gi = (2x0 ) (4y0 ) = 0. Por lo tanto, f g.
y0 y02 2
Nota: Las superficies de nivel de una funcion real de tres variables f : R3 R ,

son superficies de la forma f (x, y, z) = k, con k = constante. En este caso, por cada punto
(x0 , y0 , z0 ) pasa una y solo una superficie de nivel de f .
Teorema 2.1.4 El vector gradiente f |x=x0 es un vector normal a la superficie

de nivel correspondiente, f (x) = f (x0 ).
on. Sea S una superficie de nivel de una funcion f con derivadas parciales de
Demostraci
primer orden continuas, la cual se describe por f (x, y, z) = k. Sea p(x0 , y0 , z0 ) un punto
sobre S y sea C una curva arbitraria sobre la superficie S pasando a traves del punto p, la cual
se puede parametrizar mediante una funcion vectorial continua (t) = (x(t), y(t), z(t))T .
z
6
f (x0 , y0 , z0 )
f
@
C ..... @
S .....
p@ @
p @
@ 0)
R (t @
- @
y

f =k
x

Sea t0 el valor del par

ametro correspondiente a p. Es decir, (t0 ) = (x0 , y0 , z0 ) = p. Como
la curva C est
a sobre la superficie, entonces cualquier punto (t) = (x(t), y(t), z(t)) debe
on de S. Esto es, f ((t)) = f (x(t), y(t), z(t)) = k. Si x, y, z son funciones
satisfacer la ecuaci
diferenciables de t y f tambien es diferenciable respecto a sus argumentos, entonces la
d f dx f dy f dz
Regla de la Cadena (Seccion 2.2) implica que f ((t)) = + + = 0.
dt x dt y dt z dt

dx dy dz

Como la derivada de una trayectoria es (t) = , , , entonces la relacion anterior
dt dt dt
puede escribirse en terminos del producto interior como hf ((t)), (t)i
= 0. En particular,
cuando t = t0 , (t0 ) = (x0 , y0 , z0 ) = p y as, hf (p), (t
0 )i = 0. Por lo tanto, el vector
gradiente de f en p, f (x0 , y0 , z0 ), es perpendicular al vector tangente (t
0 ) a la curva C
que pasa a traves de p sobre la superficie S. Pero, como la curva elegida es arbitraria,
entonces f es perpendicular a las superficies de nivel de la funcion f en cada punto del
dominio.
De este resultado se deduce que el vector gradiente proporciona la direccion en la

cual ocurre el incremento mas grande de la funcion f (esto es, la direccion de maximo
crecimiento). Estas dos propiedades del gradiente son casi equivalentes, ya que si
se mueve a partir del punto p sobre una superficie de nivel, el valor de la funcion f
no cambia, por lo que resulta razonable que s se mueve en la direccion perpendicular a
tal superficie, entonces se obtiene el maximo crecimiento y por consiguiente, el mnimo
crecimiento se obtiene en la direccion opuesta, f .
Ejemplo 30: I) El gradiente para cada uno de los campos escalares siguientes,
respectivamente esta dado por
a) f (x, y, z) = axy + byz + cxz, f (x, y, z) = (ay + cz, ax + bz, cx + by)T .

1 2 T
b) g(x, y, z) = , g(x, y, z) = 2 (x, y, z) .
x2 + y 2 + z 2 2 2 2
(x + y + z )
1 2
c) En general, s g(x) = 2 , entonces g(x) = xT .
kxk kxk4

1 r
II) Probar que = 3, donde r = krk.
r r
v
u n 2
uX
on. Sea r = xT , x IRn . Entonces, krk = r = t
Soluci xi . Por consiguiente,
i=1

r xi 1 1 1 xi xi
= , = 2 (r) = 2 krk = ,
xi krk xi r r xi krk krk3
Xn
1 1 r r
de donde = ni = 3 = r3 .
r xi r krk
i=1 2
n
Ejercicio 69. S r y son las coordenadas polares en el plano xy, determinar grad r
y grad .
Ejercicio 70. Determinar un vector normal unitario a la superficie x3 xyz +z 3 = 1

en el punto (1, 1, 1)T .
La Derivada Direccional: Sea f : Rn R un campo escalar. Sean v y x

dos vectores fijos en Rn , con kvk = 1. Consideremos la funcion de R a R definida por
t 7 f (x + tv). El conjunto de puntos de la forma x + tv, t R, es precisamente la
recta `(t) que pasa por el punto x en la direccion del vector v (tal como se muestra
en la figura siguiente).
6 `(t) = x + tv

x 3

v:

tv

: -
IRn

La funcion t 7 f (x + tv) es la funcion f restringida a la recta `. Ahora, con que

rapidez cambia la funcion f a lo largo de la recta ` en el punto x? Esta pregunta
es facil de contestar introduciendo el concepto de derivada direccional.
on 2.1.10 Si f :
Definici Rn R, entonces la derivada direccional de f en x
en la direccion de un vector unitario v se define y se denota por
def d

Dv f (x) = f (x + tv) . (2.1.33)
dt t=0
Nota: De la definici
on, la derivada direccional tambien se puede establecer directamente
f (x + tv) f (x)
mediante la f
ormula Dv f (x) = lim , si el lmite existe.
t0 t
Teorema 2.1.5 S f : Rn R es diferenciable en x , entonces la derivada

o
direccional de f en x en la direccion de un vector unitario v esta dada por

n
f f f X f
Dv f (x) = hf (x), vi = v1 + v2 + + vn = vi , (2.1.34)
x1 x2 xn i=1
xi
donde v = (v1 , v2 , , vn )T Rn con kvk = 1.

Demostraci
on. Sea `(t) = x + tv, de modo que f (x + tv) = f (`(t)). Entonces, aplicando
la regla de la cadena,
d f dx1 f dx2 f dxn
f (`(t)) = + + + = hf (`(t)), `(t)i.
dt x1 dt x2 dt xn dt
Pero, como `(0) = x y
`(0) = v, entonces

d
Dv f (x) = f (x + tv) = hf (x), vi.
dt t=0
Nota: hf (x), vi es la raz

on de cambio de la funcion f en el punto x, en la direccion del
vector v y se pide que kvk = 1 para que el resultado sea uniforme ademas de establecer una
escala apropiada sobre la recta `(t).
Ejemplo 31: Deteminar la derivada direccional de la funcion

g(x, y, z) = xeyz + xyez xyz
en el punto x0 = (1, 1, 1)T en la direccion del vector v = (3, 6, 2)T .
v 1
Soluci
on. El vector unitario en la direccion de v es u = = (3, 6, 2)T , mientras que
kvk 7
el gradiente de la funci
on g es

yz z yz yz z xz yz z xy
g(x) = e + ye , xze + xe , xye + xye ,
2 xyz 2 xyz 2 xyz

1 1 1
de donde g(x0 ) = 2e , 2e , 2e . Por lo tanto,
2 2 2
1
Dv g(x0 ) = hg(x0 ), ui = 2e .
2 2
Teorema 2.1.6 Sea f (x) 6= 0. Entonces, el vector gradiente f apunta en la

direccion a lo largo de la cual la funcion f crece mas rapidamente.
Demostraci
on. Sea v un vector unitario cuyo punto de aplicacion esta sobre la superfi-
cie de nivel f (x) = f (x0 ). Entonces, la razon de cambio de f en x0 en la direccion v es
hf (x0 ), vi = kf (x0 )kkvk cos , donde 0 es el angulo entre v y f (x0 ).
axima cuando = 0. Esto es, cuando v y f (x0 ) son paralelos.
Esta cantidad es m
Por lo tanto, f (x0 ) representa la direccion de maximo crecimiento de la funcion f , a partir
del punto x0 , el cual apunta en la direccion normal a las superficies de nivel de f , as que
kf (x0 )k representa la raz
on maxima de crecimiento.
n
on 2.1.11 S f (x0 ) 6= 0, entonces el hiperplano tangente a la superficie de

Definici
nivel f (x) = k, en el punto x0 , se define como
hf (x0 ), x x0 i = 0. (2.1.35)
Por ejemplo, s f es una funcion real de tres variables, entonces el plano tangente a
la superficie f (x, y, z) = k en el punto x0 = (x0 , y0 , z0 )T esta dado por
f f f
(x0 , y0 , z0 )(x x0 ) + (x0 , y0 , z0 )(y y0 ) + (x0 , y0 , z0 )(z z0 ) = 0. (2.1.36)
x y z
x2 z2
Ejemplo 32: Determinar la ecuacion del plano tangente al elipsoide +y 2 + = 3
z 4 9
en el punto x0 = (2, 1, 3)T .
6
x2 z2
+ y2 + =3
4 9
.
...
-
. . . . y
. . . . . . . . . .. . .

..
..
..

x ..
..

3x + 6y 2z 18 = 0
Soluci on. El elipsoide considerado es una superficie de nivel (con k = 3) de la funcion de

x2 z2 f x f f 2z
tres variables f (x, y, z) = + y 2 + . Como = , = 2y, = , entonces
4 9 x 2 y z 9
f f f 2
(x0 ) = 1, (x0 ) = 2, (x0 ) = y por consiguiente, la formula (2.1.36) implica
x y z 3
2
1 (x 2) + 2 (y 1) (z + 3) = 0, de donde finalmente se tiene la ecuacion del plano,
3
3x + 6y 2z 18 = 0. 2
Ejercicio 71. Sea f (x, y, z) = z 2 tan1 (x + y) zexy y sean p = (0, 0, 4)T y

q = (3, 6, 2)T dos puntos en el espacio 3dimensional.
a) Calcular la derivada direccional de f en p en la direccion de p a q.
b) Encontrar un vector unitario en la direccion del maximo crecimiento de f en
p y calcular la tasa de crecimiento de f en esa direccion.
c) Encontrar un vector unitario en la direccion en la que f disminuye mas rapida-
mente en p y calcular la razon de cambio de f en esa direccion.
Ejercicio 72. Se dice que dos superficies son ortogonales en alguno de sus puntos
de interseccion p(x0 , y0 , z0 ) si sus vectores normales en p son perpendiculares.
a) Demostrar que las graficas de F (x, y, z) = 0 y de G(x, y, z) = 0 son ortogonales
en p si y solo si Fx Gx + Fy Gy + Fz Gz = 0, en el dominio com
un donde F y G
tienen derivadas parciales.
b) Mostrar que la esfera x2 + y 2 + z 2 = a2 y el cono x2 + y 2 z 2 = 0 son
ortogonales en todos sus puntos de interseccion.
c) Mostrar que todo plano tangente al cono z 2 = x2 + y 2 pasa por el origen.
Ejercicio 73. Mostrar que la esfera x2 + y 2 + z 2 8x 6y 8z + 24 = 0 y el elipsoide

3x2 + 2y 2 + z 2 = 9 son tangentes uno a otro en el punto (1, 1, 2)T .
Ejercicio 74. Determinar la ecuacion del plano tangente a la superficie dada, en el

punto indicado:
a) xy 2 z 3 = 12, (3, 2, 1)T ,
b) z + 1 = xey cos z, (1, 0, 0)T .

xyz
Ejercicio 75. S f (x, y, z) = ,
x2 + y2 + z2
i) Determinar la expresion del plano tangente a la superficie de nivel correspon-
diente al punto x0 = (1, 1, 1)T ,
ii) Hacia donde f crece mas rapido?, a partir de x0 ,
iii) Calcular la razon de cambio mnima de f en x0 .
Ejercicio 76. Supongamos que la temperatura en un punto (x, y, z)T en el espacio

80
es dada por T (x, y, z) = , donde T esta medida en o C y x, y, z
1 + x + 2y 2 + 3z 2
2
en metros. En cual direccion, a partir del punto (1, 1, 2)T , se incrementa mas
rapidamente la temperatura? Cual es la razon maxima de crecimiento?
Ejercicio 77. La superficie de un lago puede ser referida sobre el plano xy de tal
modo que la profundidad en cada punto (x, y) esta dada por P (x, y) = 3002x2 3y 2
(en metros). S un individuo se encuentra nadando en ese lago, precisamente en el
punto (4, 9), cual es la profundidad en dicho punto?, en cual direccion debe nadar
para que la profundidad disminuya lo mas rapido posible?
n
Ejercicio 78. Supongase que z = f (x, y) es una funcion diferenciable para todo
punto (x, y) y sea p(2, 1). S la derivada direccional de z en p en la direccion haca
el punto a(0, 1) vale 7, mientras que en la direccion haca el punto b(1, 1) vale 3,
cual es el valor de la derivada direccional de z en el mismo punto p, pero en la
direccion haca c(6, 2)?
Ejercicio 79. Supongamos que Usted va a escalar una monta

na cuya forma es dada
por la ecuacion z = 1000 0.001x2 0.02y 2 , donde el eje y apunta haca el Norte y el
eje x haca el Este (todas las distancias en metros). S Usted esta parado en el punto
con coordenadas (60, 100, 764)T .
a) Asciende o desciende cuando se mueve en la direccion Noroeste?

b) En que direccion debera Usted dirigirse inicialmente para alcanzar la cima
de la monta
na de la forma mas rapida?
c) S Usted esta subiendo en tal direccion, en que angulo sobre la horizontal
debera Usted subir?
d) En que direccion se debe Usted mover para no variar su altura?
Ejercicio 80. Los Biologos marinos han determinado que cuando un tiburon detecta
la presencia de sangre en el agua, nada en la direccion en la que concnetracion de la
sangre se incrementa mas rapidamente. Despues de un derrame de sangre y basandose
sobre ciertas pruebas del agua del mar, han determinado que la concentracion de
sangre en partes por millon en un punto p(x, y) sobre la superficie es aproximada
2 +2y 2 )/104
por C(x, y) = e(x , donde x y y son medidas en metros en un sistema de
coordenadas rectangulares con la fuente de sangre en el origen.
a) Identificar las curvas de nivel de la funcion concentracion C, bosquejar algunos

miembros de esta familia y graficar una porcion de la superficie, correspondiente
a un dominio alrededor del origen.
b) S un tiburon se encuentra en el punto p0 (x0 , y0 ) cuando por primera vez detecta
la presencia de sangre en el agua, determinar la trayectoria (direccion en cada
punto) que seguira el tiburon hacia la fuente de la sangre (el origen).
Ejercicio 81. Verificar que Dv hv, xi = 1 en cualquier direccion.

2.1.3 La Matriz Jacobiana
En esta parte se define la derivada general de una funcion vectorial de varias variables
de una forma natural teniendo como base solida el concepto de derivada parcial.
S f : Rn Rm, entonces se escribe f (x) = (f1(x), f2(x), , fm(x)) T
, con
fi : R . La transformacion lineal L a la que hace referencia la Definicion 2.1.2,
m
nica, es una matriz de orden m n cuyas entradas

la cual por el Teorema 2.1.1 es u
son las derivadas parciales de las funciones componentes respecto a cada uno de sus
argumentos tal como se establece en el siguiente resultado.
Teorema 2.1.7 La transformacion lineal L : Rn Rm, denotada por Df (x0),

es una matriz m n cuyos elementos
`ij son precisamente las derivadas parciales de
fj
las funciones componentes, , i = 1, 2, , n, j = 1, 2, , m. Esto es,
xi x=x0

f1 f1 f1
(x0 ) (x0 ) (x0 )
x1 x2 xn f1 (x0 )
f2 f2 f2

(x0 ) (x0 ) (x0 ) f2 (x0 )

Df (x0 ) = x1 x2 xn .. , (2.1.37)
.. .. .. ..
.

. . . .

fm fm fm fm (x0 )
(x0 ) (x0 ) (x0 )
x1 x2 xn
la cual es conocida como la matriz jacobiana de f en x0 .
on. Utilizando la propiedad |xk | kxk, k = 1, 2, , n, para x

Demostraci Rn,
de la definici
on de diferenciabilidad, relacion (2.1.4), se tiene que
|fj (x0 + h) fj (x0 ) (L(h))j |

lim = 0, j = 1, 2, , m,
h0 khk
donde (L(h))j es la jesima componente del vector columna L(h). Ahora, tomando
h = tni = (0, , t, , 0) , t R, se tiene que el lmite anterior es equivalente a
T
fj (x0 + tni ) fj (x0 )

i

lim L(n ) j = 0,
t0 t
n
el cual a su vez implica

fj (x0 + tni ) fj (x0 )
= L(ni ) j ,

lim
t0 t
pero el lmite de la izquierda es la derivada parcial de fj con respecto a la iesima variable,

evaluada en el punto x0 (Definici
on 2.1.6, relacion 2.1.21). Por lo tanto,

fj
`ij = L(ni ) j =

.
xi x=x0
Nota: La matriz jacobiana Df (x0 ) esta definida en cada punto en el que las mn derivadas
fj
parciales existen, la cual contiene la informacion de los cambios de la funcion f a lo
xi
largo de los ejes principales (direcci
on de cada vector de la base). Esta es la razon por la
que se le nombra derivada general o total8 de la funcion vectorial de varias variables f .
Ejemplo 33: Sea f : R3 R2 dada por

T
f (x, y, z) = exy + yz, xy 2 z ln (x2 y 2 z 2 ) .
Calcular la matriz jacobiana de f para cualquier x en el interior de su dominio de

definicion, = R3 \ {0}, y en particular, para x0 = (1, 1, 1)T .
on. Como f1 (x, y, z) = exy + yz

Soluci y f2 (x, y, z) = y 2 xz ln (x2 y 2 z 2 ). Entonces,

exy exy + z y

Df (x) =

,

y 2 z ln (x2 y 2 z 2 ) + 2y 2 z 2yxz ln (x2 y 2 z 2 ) + 2xyz y 2 x ln (x2 y 2 z 2 ) + 2y 2 x
!
1 0 1
de donde se obtiene mediante evaluacion directa Df (x0 ) = .
2 2 2 2
Ejercicio 82. Sea T : Rn Rm una transformacion lineal dada por la matriz mn,
n
X
A = (aij ) , de tal forma que Ax tiene componentes yi = aij xj . Usar la desigualdad
j=1
qX X
de Schwartz para probar que kAxk Mkxk, donde M = a2ij .
8
Observaci
on: El termino derivada total tambien es utilizado en Fsica en otro contexto,
lo cual debe ser claro a partir del mismo contexto.
La diferenciabilidad de una funcion implica su continuidad. Sin embargo, el recproco

no necesariamente se cumple (dar ejemplos).
Teorema 2.1.8 Sea f : Rn Rm una funcion diferenciable en x0 , entonces

o
f satisface la condicion de Lipschitz
kf (x) f (x0 )k < Mkx x0 k, (2.1.38)
para alguna constante M y x B (x0 ), pero x 6= x0 . Por consiguiente, f es continua

en x0 .
Demostraci
on. La Definici
on 2.1.2 de diferenciabilidad, establecida en la relacion (2.1.4)
y el hecho que kDf (x0 )hk Kkhk, donde (por el Ejercicio 82) K es la raz cuadrada de
la suma de los cuadrados de los elementos de la matriz Df (x0 ), es decir, K = kDf (x0 )k,
eligiendo = 1, se tiene que
1 > 0 3 0 < khk < 1 kf (x0 + h) f (x0 ) Df (x0 )hk < khk = khk.
Ahora, si khk < 1 , entonces
kf (x0 + h) f (x0 )k = kf (x0 + h) f (x0 ) Df (x0 )h + Df (x0 )hk

kf (x0 + h) f (x0 ) Df (x0 )hk + kDf (x0 )hk
< khk + Kkhk = (1 + K)khk.
Tomando M = 1 + K se obtiene la condicion de Lipschitz. Por otro lado, si se elige 0 > 0

y = min{1 , 0 /(1 + K)}, entonces khk < , lo cual implica que
kf (x0 + h) f (x0 )k < (1 + K)0 /(1 + K) = 0 .
Por lo tanto, lim f (x0 + h) = f (x0 ), de donde se concluye que f es continua en x0 .

h0
Rn Rm no es continua en x0 , entonces no es
o
Corolario 2.1.1 S f :
diferenciable9 en x0 .
Por lo tanto, una funcion que no es continua en alg

un punto de su dominio no puede
ser aproximada por ninguna transformacion lineal en tal punto.
9
En el sentido de la Definici
on 2.1.2, ya que existe un concepto mas general de diferenciacion en
el contexto de la Teora de Distribuciones o Funciones Generalizadas.
n
Un criterio importante de diferenciabilidad es que todas las derivadas parciales de

o
primer orden existan y sean continuas. Es decir, que la funcion f sea de clase C1 (),
lo cual significa que cada una de sus funciones componentes lo sea.
Teorema 2.1.9 Sea f : Rn Rm . Supongamos que las derivadas parciales

fj o
existen y son continuas en una vecindad de un punto x . Entonces, f es
xi
diferenciable en x.
Demostraci
on. Es suficiente considerar el caso m = 1, ya que en general el criterio se
cumple componente a componente y las componentes de una funcion vectorial son campos
escalares. As, sea f : Rn R. En este caso, debemos verificar la relacion (2.1.4).
Es decir,
n
X f
|f (x + h) f (x) (x)hi |
xi |f (x + h) f (x) f (x)h|
i=1
lim = lim = 0.
h0 khk h0 khk
Para empezar,
f (x + h) f (x) = f (x1 + h1 , x2 + h2 , . . . , xn + hn ) f (x1 , x2 , . . . , xn )

= [f (x1 + h1 , x2 + h2 , . . . , xn + hn ) f (x1 , x2 + h2 , . . . , xn + hn )]
| {z }
f
(w1 )h1
x1
+ [f (x1 , x2 + h2 , . . . , xn + hn ) f (x1 , x2 , . . . , xn + hn )]
| {z }
f
(w2 )h2
x2
..
.
+ [f (x1 , x2 , . . . , xn1 + hn1 , xn + hn ) f (x1 , x2 , . . . , xn1 , xn + hn )]
| {z }
f
(wn1 )hn1
xn1
+ [f (x1 , x2 , . . . , xn1 , xn + hn ) f (x1 , x2 , . . . , xn1 , xn )],
| {z }
f
(wn )hn
xn
donde por el Teorema del Valor Medio
w1 = (c1 , x2 + h2 , x3 + h3 , . . . , xn + hn )T , c1 (x1 , x1 + h1 ),
w2 = (x1 , c2 , x3 + h3 , . . . , xn + hn )T , c2 (x2 , x2 + h2 ),
..
.
wn = (x1 , x2 , . . . , xn1 , cn )T , cn (xn , xn + hn ).
As,

n
X f f f f f
f (x + h) f (x) (x)hi = (w1 ) (x) h1 + (w2 ) (x) h2

xi x1 x1 x2 x2
i=1
f f
+ + (wn ) (x) hn
xn xn

f f f f

(w1 ) (x) |h1 | +
(w2 ) (x) |h2 |
x1 x1 x2 x2
f f
+ + (wn ) (x) |hn |
xn xn

f f f f

(w1 ) (x) +
(w2 ) (x)
x1 x1 x2 x2
f f
+ + (wn ) (x) khk.
xn xn
ltimo paso se debe a que |hi | khk, i = 1, 2, . . . , n. Por lo tanto,
El u
n

X f
f (x + h) f (x) (x)hi

xi
f f

i=1
lim lim (w1 ) (x) +
h0 khk h0 x 1 x1
f f
+ (wn ) (x) .
xn xn
Cuando h 0, wk x, k = 1, 2, . . . , n. Como las parciales son continuas, el lado

derecho tiende a cero cuando h 0, as que el lado izquierdo tambien tiende a cero, lo cual
se quera probar para asegurar la diferenciabilidad del campo escalar f .
Ejercicio 83. Sea f (x, y, z) = (ex , cos y, sen z) .
i) Calcular la matriz jacobiana Df .

ii) En general, bajo que condiciones Df es una matriz diagonal?
Ejercicio 84. Sea f : R3 R4 dada por

!
yz (x2 +y 2 +z 2 ) 1 x2 4xz 4yz + y 2
f (x, y, z) = x e ,p , x cos(xz) + y sen (yz), .
x2 + y 2 + z 2 5x 6y 13z + xyz
Determinar la matriz jacobiana de f en el punto x0 = (0, 1, 0).

n
2.1.4 C
alculo Diferencial de Campos Vectoriales
El concepto general de diferenciacion adquiere un matiz especial cuando las funciones

sobre las que se define son campos vectoriales. El objetivo de esta parte es trabajar de
la forma mas natural posible dos de las operaciones diferenciales mas importantes del
calculo, la divergencia y el rotacional de un campo vectorial, describiendo rapidamente
el significado fsico y geometrico respectivo.
En general, un campo vectorial en Rn es una funcion que en cada punto del dominio de
definicion asigna un vector (magnitud y direccion especficas). Un ejemplo concreto
de campo vectorial es el campo vectorial gradiente de un campo escalar diferenciable,
f : Rn R, ya que f (x) asocia a cada punto del interior del dominio de
f e Rn Rn , donde
un vector anclado en ese punto. Es decir, f : e es el
f
dominio reducido para el que existen todas y cada una de las derivadas parciales .
xi
Geometricamente,
f (p3 ) f (p2 )
f (p4 ) 6
I
@
@ p3 p2
@p4
p1 f (p1 )
:

...
pn -

f (pn )
) n

IR
Definici
on 2.1.12 Un campo vectorial V en Rn se denomina campo vectorial
conservativo si es el gradiente de alg
un campo escalar. Esto es, si existe una funcion
f : Rn R tal que f = V. En este caso, se dice que f es una funcion de
potencial para V.
Nota: i) M
as adelante, se establece un criterio que ayuda a determinar cuando un campo
vectorial es conservativo.
ii) En particular, el caso de mayor interes, por su importancia y significado fsico-geometrico,
es el de los campos vectoriales tridimensionales. En lo que sigue, se consideraran campos
vectoriales en R3 o R2, a menos que se mencione lo contrario (en la Subseccion 1.3.5
se ilustran algunos ejemplos cl
asicos en R3).
Definici
on 2.1.13 Un campo vectorial en el espacio euclideano tridimensional es una
funcion, V : R3 R3, que a cada punto de cierta region le asigna un vector
3-dimensional. Esto es,
V(x) = V(x, y, z) = (P (x, y, z), Q(x, y, z), R(x, y, z))T = P (x)i+Q(x)j +R(x)k. (2.1.39)
Ejemplo 34: De acuerdo a la Ley de Gravitacion de Newton, la fuerza gravitacional

sobre una unidad de masa localizada en un punto r = (x, y, z)T producida por una
masa M localizada en el origen de R3 es F = 2 v, donde G es la constante
GM
p r
universal de gravitacion, r = krk = x2 + y 2 + z 2 es la distancia desde el punto
(x, y, z)T al origen
y v = r/r es el vector unitario en la direccion radial. Notese que
GM
F = = . Es decir, F es menos el gradiente del potencial gravita-
r
GM
cional = . Por lo tanto, el campo gravitacional F es un campo vectorial
r
conservativo.
z
6
..
...
m = 1
..
F
..
..
..
M
.}. .
. ...............
-
...
.. . y
.

Superficie de Nivel de

x

+
Ejercicio 85. Las lneas de flujo (o lneas de corriente) de un campo vectorial son las
trayectorias que sigue una partcula cuyo campo de velocidades es el campo vectorial
dado. Esto es, los vectores definidos por un campo vectorial son tangentes a las lneas
de flujo en cada punto del dominio de definicion.
a) Bosquejar algunas de las lneas de flujo para el campo vectorial 2-dimensional

F (x, y) = yi xj. De su bosquejo, Podra usted determinar las ecuaciones de
las lneas de flujo?
b) Determinar la ecuacion para la lnea de flujo que pasa a traves del punto
x0 = (1, 1)T , resolviendo el sistema de ecuaciones diferenciales lineales
dx

= y,
dt

dy = x.

dt
n
Divergencia y Rotacional: La divergencia y el rotacional son dos conceptos

diferenciales muy importantes que caracterizan de una forma especial a los campos
vectoriales.
Definici
on 2.1.14 Sea V un campo vectorial en R3. Se define la divergencia y el
rotacional de V mediante la operacion producto interior y producto cruz del operador
y V, respectivamante, de la forma siguiente:
V V V def
h, Vi V = i +j +k = div V, (2.1.40)
x y z
V V V def
V =i +j +k = rot V. (2.1.41)
x y z
Si V = P i + Qj + Rk, entonces la divergencia tiene la expresion

P Q R
div V = + + , (2.1.42)
x y z
mientras que el rotacional,

R Q P R Q P
rot V = i+ j+ k
y z z x x y

i j k

(2.1.43)

=

x y z

P Q R
Nota: La notaci
on para el rotacional algunas veces vara en la literatura en ingles, donde
se utiliza curl en lugar de rot
Ejemplo 35: La divergencia y el rotacional del campo vectorial

T
xz x2 y 2 2 2
V(x, y, z) = , ze , x yz ,
y
respectivamente, es
z 2 2
div V = 2yzex y + 2x2 yz,
y T
2 2 x2 y 2 x 2 x2 y 2 xz
rot V = x z e , 2xyz , 2xze + 2 .
y y 2
Teorema 2.1.10 Si f : R3 R es una funcion de clase C2 (), entonces

rot (f ) = 0.
Es decir, el rotacional de cualquier campo gradiente es el vector cero.

f f f
on. Como f =
Demostraci , , , entonces
x y z

i j k

2 2
2 2
2 2

= f f i+ f f j+ f f k.

f = x y z yz zy zx xz xy yx

f f f

x y z
Cada componente del vector del lado derecho es cero si el orden de la diferenciacion no tiene
importancia. Es decir, cuando las derivadas parciales mixtas son identicas, lo cual es valido
cuando f C2 (). Por lo tanto, rot f = 0.
As, se concluye que si V es un campo vectorial conservativo, entonces rot (V) = 0.

Por lo tanto, este resultado proporciona un criterio para verificar cuando un campo
vectorial no es conservativo. Esto es, si rot (V) 6= 0, entonces V no es conservativo.
Definici
on 2.1.15 Si rot V = 0, entonces el campo vectorial V se dice ser
irrotacional.
T
y x
Ejemplo 36: Mostrar que el campo vectorial V(x, y, z) = , , z2
x2 + y 2 x2 + y 2
es irrotacional. Es conservativo?
Soluci
on. Como

i j k

x

= y

V = x y z k = 0,
y x

x x2 + y 2 y x2 + y 2

x2 + y 2 z2
x2 + y 2
entonces el campo vectorial dado es irrotacional. Ademas, este campo

vectorial es conser-
1 x 1
vativo, ya que es el gradiente del campo escalar (x, y, z) = tan + z3.
y 3 2
Afirmaci
on: Todo campo conservativo (expresable como el gradiente de un potencial
escalar) es irrotacional y viceversa! (Mostrar el recproco, i.e., V = 0 V = f ).
n
Teorema 2.1.11 S V : R3 R3 es un campo vectorial de clase C2 (), entonces

div rot V = ( V) = 0.
Es decir, la divergencia de cualquier campo rotacional es cero.
Demostraci
on. Para el campo vectorial considerado, se verifica que

x y z
( V) = ,

x y z

P Q R
de donde, de la teora de determinantes, se sigue que div rot V = 0.
Nota: El rotacional, V, de un campo vectorial V esta relacionado con las rotaciones o

giros que ocurren en el campo de direcciones, mientras que la divergencia, V, esta rela-
cionada con la variaci
on de la magnitud (expansiones o contracciones) del campo vectorial.
Definici
on 2.1.16 Si div V = 0, entonces se dice que el campo vectorial V es
incompresible o solenoidal.
xy T
Ejemplo 37: A) Calcular la divergencia de V = x2 y, zx, .
z
xy xy
Soluci
on. div V = 2xy +0 = 2 (2z 2 1). Este campo vectorial no puede ser el rota-
z2 z
cional de otro campo vectorial, ya que de serlo su divergencia tendra que ser identicamente
igual a cero. Tampoco es conservativo (verificarlo!).
B) Calcular el rotacional de V(x, y, z) = (x, xy, z)T .
Soluci
on.
i j k

V = = (0 0)i + (0 0)j + (y 0)k.

x y z

x xy z
Por lo tanto, V = (0, 0, y)T 6= 0 si y 6= 0, lo cual implica que V no puede ser el

gradiente de ning
un campo escalar, ya que tendra rotacional cero. Tampoco es el rotacional
de otro campo vectorial, ya que de ser as, entonces tendra divergencia cero y en este caso
div V = 2 + x.
2
Observaci
on: i) Si un campo vectorial V representa el flujo de un fluido, entonces
rot V = 0 significa, fsicamente, que el fluido esta libre de rotaciones (es irrotacional).
Es decir, est
a libre de remolinos o vortices, lo cual esta estrechamente relacionado con la
viscosidad y compresibilidad del fluido. En el Captulo 4, una vez que se establezca el
Teorema de Stokes, se dar
a detalladamente el significado fsico preciso del rotacional de un
campo vectorial.
ii) Ahora, si V es el campo de velocidades de un gas (o un fluido cualquiera), entonces

se puede adelantar que div V representa la razon de expansion del gas por unidad de
volumen. Cuando se discuta el Teorema de Gauss, se dara el significado fsico completo de
la divergencia de un campo vectorial.
Por u
ltimo, otro operador diferencial muy importante y utilizado que aparece cuando
se calcula la divergencia de un campo vectorial gradiente es el operador de Laplace,
2f 2f 2f
2 f = (f ) div (f ) = + + . Sin embargo, es mas com un que el
x2 y 2 z 2
operador 2 = se denote por , Definicion 2.1.8, el cual en el caso n-dimensional
tiene la forma 2.1.25. Cuando este operador se anula, f = 0, la relacion obtenida es
conocida como Ecuacion de Laplace ndimensional, la cual es una ecuacion diferencial
parcial lineal de segundo orden y a la funcion f se le llama funcion armonica.
Por otro lado, s V = P i + Qj + Rk es un campo vectorial de clase C2 , entonces 2 V

se define componente a componente. Esto es,
T
2 V = 2 P i + 2 Qj + 2 Rk = 2 P, 2 Q, 2 R . (2.1.44)
Nota: El operador 2 se conoce como el Laplaciano, tiene diferente representaci

on
para diferentes sistemas de coordenadas y esta presente en una gran variedad de problemas
fsicos (realizar una amplia descripcion en clase).
Ahora, si V = P i + Qj + Rk es un campo vectorial diferenciable en R3 ,

o
entonces la divergencia y el rotacional de V tienen una relacion directa con la Matriz

Jacobiana de V,
P P P

x y z
Q Q Q
DV = ,

x y z
R R R
x y z
n
en el siguiente sentido
div V = tr DV, (2.1.45)
donde tr denota la traza de la matriz10 . Ademas, de la Teora de Matrices, se tiene
que toda matriz real cuadrada A puede expresarse como
1 1
A + AT + A AT = S + SS,

A=
2 2
T
1 1
A AT

donde S = 2
A + A es una matriz simetrica (sij = sji , i 6= j) , SS = 2
es una matriz semisimetrica o antisimetrica (elementos por arriba y por abajo de
la diagonal son iguales, pero opuestos en signo, es decir ssij = ssji , i 6= j) y AT
denota a la transpuesta de la matriz A. Cuando A es la Matriz Jacobiana, la parte
semisimetrica es
rot V $
P Q P R
0 y

x z

x
6Q P Q R
1h T
i 1
2
DV (DV) =
2 x

y 0 z

y
R P R Q
x

z y

z 0
&
rot V
Los elementos de rot V son los elementos de esta matriz, en el sentido que indica
la flecha, con los respectivos signos heredados de la posicion correspondiente en la
matriz.
Nota: Si la Matriz Jacobiana es simetrica, entonces cada elemento de la matriz SS es cero

y en este caso rot V = 0.
Ejemplo 38: a) Sea F (x, y, z) = (f (x), g(y), h(z))T , entonces

0
f (x) 0 0
DF = 0 0
g (y) 0 .

0 0 h0 (z)
Esto implica que div F = f 0 (x) + g 0 (y) + h0 (z) y rot F = 0.
10
Recuerde que la traza de una matriz se define como la suma de los elementos de la diagonal
principal de la matriz.
T
b) Sea V(x, y, z) = (xy 2 z 2 , z 2 sen y, x2 ey ) , entonces

y2z2 2xyz 2 2xy 2 z

DV =
0 z 2 cos y 2z sen y
,
2xey 2
x e y 0
de donde div V = z 2 (y 2 + cos y). Ahora,

0 2xyz 2 2xy 2 z 2xey
1h i 1
DV (DV)T = 0

2xyz 2 2z sen y x2 ey
,
2 2
2xey 2xy 2 z x2 ey 2z sen y 0
T
de donde rot V = (x2 ey 2z sen y, 2xy 2 z 2xey , 2xyz 2 ) .
c) Establecer un campo vectorial con divergencia y rotacional iguales a cero.

T
y x
Soluci
on. Sea V(x, y, z) = , , 0 , entonces
x2 + y 2 x2 + y 2

2xy y 2 x2 0
1
y 2 x2 2xy 0 ,

DV = 2
(x + y 2 )2
0 0 0
de donde div V = 0 y rot V = 0, ya que DV es simetrica.
d) Divergencia y rotacional de un campo gradiente.

on. Sea V(x, y, z) = =
Soluci , , , entonces
x y z
2
2 2

x2 yx zx
2
2 2
DV = .
xy
2 y 2 zy
2 2

xz yz z 2
2 2 2
De aqu, div V = + + 2 = , la divergencia del gradiente de un campo
x2 y 2 z
escalar es el Laplaciano del campo escalar. Por lo tanto, el gradiente de una funci
on
armonica tiene divergencia nula. Por otro lado, rot(grad ) = 0 para todo campo escalar
de clase C2 , ya que en este caso, las parciales mixtas son iguales lo cual implica que la
Matriz Jacobiana DV es simetrica.
n
e) Divergencia y rotacional de un campo rotacional.
on. Si F = (f1 , f2 , f3 )T , entonces rot F es de nuevo un campo vectorial y por tanto

Soluci
se puede calcular divergencia y rotacional. De la formula (2.1.43), para el calculo del
rotacional, se tiene que
2 f3 2 f2 2 f3 2 f2 2 f3 2 f2

xy xz y 2 yz zy z 2

2 f1 2 f3 2 f1 2 f3 2
f1 2 f3
D(rot F ) = .

xz x2 yz yx z 2 zx
2 f2 2 f1 2 f2 2 f1 2 f2 2 f1

x2 xy yx y 2 zx zy
Si el campo es de clase C2 , entonces las parciales mixtas son iguales y por lo tanto,
div (rot F ) = 0 y rot (rot F ) = grad (div F ) F , donde F = (f1 , f2 , f3 )T .
Nota: Esta identidad relaciona los cuatro operadores diferenciales principales del Calculo
Vectorial: el gradiente, la divergencia, el rotacional y el Laplaciano de un campo vectorial
(Formula 10, de la tabla del Ejercicio 90). 2
Ejercicio 86. Sea f un campo escalar en R3 y sea F un campo vectorial en R3.

Establecer cada una de las expresiones siguientes como Campo Vectorial (CV),
Campo Escalar (CE) o No Tiene Sentido la operacion (NTS).
a) rot f, b) grad f, c) div F , d) rot (grad f ), e) grad F ,

f ) grad (div F ), g) div (grad f ), h) grad (div f ), i) rot (rot F ),
j) div (div F ), k) (grad f ) (div F ), l) div (rot (grad f )).
Ejercicio 87. En cada caso, verificar la identidad dada, donde r = xi + yj + zk y

r = krk.

r 1 r
i) r = , ii) r = 3, iii) r = 0, iv) = 3 , r 6= 0,
r r r
r
v) (rr) = 4r, vi) ln r = 2
, r 6= 0, vii) r = 12r, viii) (rn r) = 0,
3
r
1
ix) (rn ) = nrn2 r, x) = 0, r 6= 0, xi) rn = n(n + 1)rn2 .
r
Ejercicio 88. Determinar DF , div F y rot F para el campo vectorial

T
xy 2 2
F (x, y, z) = xyz, 2 , 3x y 7xz .
x + y2
Ejercicio 89. Considerar el campo vectorial
V(x, y, z) = 2xyi + (x2 + 2yz)j + (y 2 + 1)k.

Es V un campo vectorial gradiente? Si as es, encontrar una funcion potencial
para V.
Ejercicio 90. Si f, g son campos escalares y F , G son campos vectoriales en R3,

entonces f F y F G tambien son campos vectoriales definidos por
(f F )(x, y, z) = f (x, y, z)F (x, y, z), (F G)(x, y, z) = F (x, y, z) G(x, y, z).

Ademas, hF , Gi(x, y, z) = hF (x, y, z), G(x, y, z)i define un campo escalar. Bajo
estas consideraciones y asumiendo que las derivadas parciales apropiadas existen y
son continuas, probar las identidades de la siguiente tabla:
rmulas de Diferenciacio
Fo n Vectorial
1. grad (f + g) = grad f + grad g, R, constante
2. grad (f g) = (grad f ) g + f (grad g)
f 1
3. grad ( ) = 2 [(grad f ) g f (grad g)]
g g
4. div (F + G) = div F + div G
5. rot (F + G) = rot F + rot G
6. div (f F ) = f div F + f F
7. rot (f F ) = f rot F + (f ) F
8. div (F G) = G rot F F rot G
9. div (f g) = 0
10. rot rot F = grad (div F ) F
11. grad (F G) = (F )G + (G )F + F rot F
12. (f g) = f g + g f + 2 (f g)
Ejercicio 91. Sea v un campo vectorial en R3.

a) Utilizando la formula 11 de la Tabla del Ejercicio 90, mostrar que
(v v) = 2v v + 2v ( v)
y deducir que la relacion v v = 21 v 2 es valida cuando el campo vectorial es

irrotacional, donde v = kvk.
b) Verificar el resultado general del inciso a) cuando v = xyi + (y 2 + 1)j.

n
2.2 La Regla de la Cadena
En esta seccion se considera el concepto general de diferenciacion para la composicion

de funciones, lo cual lleva a uno de los principales resultados del calculo diferencial,
la regla de la cadena. En general, el procedimiento para diferenciar composicion de
funciones vectoriales de varias variables es similar, aunque relativamente mas compli-
cado, al caso de funciones reales de una variable. Esta analoga se establece utilizando
notacion matricial. Una aplicacion inmediata de la regla de la cadena se encuentra en
el establecimiento del teorema de Taylor generalizado. Ademas, la regla de la cadena
tambien se utiliza para obtener la representacion de cualquier operador diferencial
en otros sistemas coordenados de referencia.
2.2.1 F
ormula General y Casos Especiales
La formula general de la regla de la cadena en termino de las matrices jacobianas de

las funciones que intervienen en la funcion composicion es analoga al caso simple de la
composicion de funciones reales de una variable tal como se establece en el siguiente
resultado.
Teorema 2.2.1 Sean f : Rn Rm y g:V Rm Rr funciones dadas

tales que f manda a en V (o en una parte de V), de tal forma que la composicion
h = g f : Rr esta bien definida. Si f es diferenciable en x0 y g es diferenciable
def
en y 0 = f (x0 ), entonces h es diferenciable en x0 y
Dh(x0 ) D (g f ) (x0 ) = Dg(f (x0 ))Df (x0 ). (2.2.1)
Demostraci
on. Por verificar que
kg(f (x)) g(f (x0 )) Dg(f (x0 )) Df (x0 )(x x0 )k
lim = 0.
xx0 kx x0 k
Es decir, g f satisface la definici

on general de diferenciabilidad (2.1.4).
Reescribiendo el numerador de la relacion anterior y aplicando la desigualdad del triangulo,

se tiene que
kg(f (x)) g(f (x0 )) Dg(f (x0 )) (f (x) f (x0 ))
+Dg(f (x0 )) [(f (x) f (x0 )) Df (x0 )(x x0 )] k
kg(f (x)) g(f (x0 )) Dg(f (x0 )) (f (x) f (x0 )) k
+kDg(f (x0 )) [f (x) f (x0 ) Df (x0 )(x x0 )] k.
2.2. La Regla de la Cadena 127
Como kDg(f (x0 ))hk Mkhk para alguna constante M (Ejercicio 82), entonces la ex-
presion anterior es menor o igual que
kg(f (x)) g(f (x0 )) Dg(f (x0 )) (f (x) f (x0 )) k + Mkf (x) f (x0 ) Df (x0 )(x x0 )k.
(2.2.2)
Ahora, como f es diferenciable en x0 , entonces dado > 0 existe 1 > 0 tal que
0 < kx x0 k < 1 implica
kf (x) f (x0 ) Df (x0 )(x x0 )k
< .
kx x0 k 2M
kx x0 k
Esto hace que el segundo termino de la desigualdad (2.2.2) sea menor que .
2
Para el primer termino, como f es diferenciable, entonces el Teorema 2.1.8 implica que
kf (x) f (x0 )k < M1 kx x0 k para alguna constante M1 si x esta cerca de x0 , digamos
0 < kx x0 k < 2 . Ahora, sea 3 > 0 tal que 0 < kf (x) f (x0 )k < 3 implique que
kf (x) f (x0 )k
kg(f (x)) g(f (x0 )) Dg(f (x0 ))(f (x) f (x0 ))k < .
2M1
3
As, si kx x0 k < y kx x0 k < 2 , entonces
M1
kf (x) f (x0 )k kx x0 k
kg(f (x)) g(f (x0 )) Dg(f (x0 ))(f (x) f (x0 ))k < < .
2M1 2
Por consiguiente, si = min {1 , 2 , 3 /M1 }, entonces
kg(f (x)) g(f (x0 )) Dg(f (x0 )) (f (x) f (x0 )) k + Mkf (x) f (x0 ) Df (x0 )(x x0 )k
kx x0 k kx x0 k
es menor que + , lo cual implica que
2 2
kg(f (x)) g(f (x0 )) Dg(f (x0 ))Df (x0 )(x x0 )k
< + = ,
kx x0 k 2 2
para 0 < kx x0 k < . Por lo tanto, g f es diferenciable en x0 y su derivada esta dada
por la ecuaci
on (2.2.1).
Algunos Casos Especiales de la Regla de la Cadena:

(
Caso I: y = f (x) dy dy dx df dg
= = (g(t)) (t).
x = g(t) dt dx dt dx dt
(2.2.3)
Caso II:

z = f (x, y)

dz f dx f dy
x = g(t) = + . (2.2.4)
dt x dt y dt
y = h(t)

n
Caso III:
z f x f y
= + .
z = f (x, y)
s x s y s
x = g(s, t) (2.2.5)
z f x f y

y = h(s, t)

= + .
t x t y t
Caso IV:
( n
u = u(x) = u(x1 , x2 , . . . , xn ) u X u xj
= , i = 1, 2, . . . , l.
xj = xj (t1 , t2 , . . . , tl )) ti j=1
x j t i
(2.2.6)
Nota: Todos estos casos simples de la regla de la cadena se estudian en un primer curso de
calculo vectorial y se resumen de manera adecuada en la definicion mas general dada por
la ecuacion (2.2.1).
Ejemplo 39: i) La temperatura de una placa circular delgada se representa por

un campo escalar , siendo (x, y) la temperatura en el punto (x, y). Introduciendo
las coordenadas polares x = r cos , y = r sen , la temperatura se convierte en una
funcion de r y , determinada por la ecuacion (r, ) = (r cos , r sen ). Expresar

las derivadas parciales y en funcion de las derivadas parciales y .
r x y
Soluci
on. Utilizando la Regla de la Cadena, relacion (2.2.5), se tiene que
x y
= + = cos + sen
r x r y r x y
x y
= + = r sen + r cos
x y x y
Equivalentemente, se puede manipular la Regla de la Cadena en forma matricial de la

siguiente forma: sea (r, ) = (P (r, )), donde P es la transformacion polar. Es decir,
P (r, ) = (x(r, ), y(r, ))T = (r cos , r sen )T .
Por consiguiente, como : R2 R, entonces

D(P (r, )) = (P (r, )) = (P (r, )), (P (r, )) ,
x y
mientras que para P : R2 R2 se tiene que

x x

(r, ) (r, ) " #
r cos r sen
DP (r, ) = = .

y y
sen r cos
(r, ) (r, )
r
Por lo tanto, para : R2 R
" #
cos r sen

D(r, ) = DDP = ,
x y sen r cos

= cos + sen , r sen + r cos
x y x y

= (r, ), (r, ) (r, ).
r
dz x
ii) Determinar , si z = , con x = et , y = ln t.
dt y
Soluci
on. Aplicando la relaci
on (2.2.4), se tiene que
et et

dz z dx z dy 1 t x 1 1 t 1 1
= + = e 2 = e = 1 .
dt x dt y dt y y t ln t (ln t)2 t ln t t ln t
iii) Si u(x, y, z) = k (x2 + y 2 + z 2 ), con x = R cos cos , y = R cos sen ,

u u
z = R sen , Rfija, entonces =0 y = 0.

Soluci
on. Aplicando (2.2.6), se tiene que
u u x u y u z
= + +
x y z
= 2kx(R sen cos ) + 2ky(R sen sen ) + 2kz(R cos )
= 2k [(R cos cos )(R sen cos ) + (R cos sen )(R sen sen )
+(R sen )(R cos )]
= 2k R2 cos sen cos2 R2 cos sen sen2 + R2 sen cos

= 2kR2 cos sen cos2 + sen2 + cos sen = 0

u u x u y u z
= + +
x y z
= 2kx(R cos sen ) + 2ky(R cos cos ) + 2kz( 0 )
= 2k [(R cos cos )(R cos sen ) + (R cos sen )(R cos cos )]
= 2k R2 cos2 cos sen + R2 cos2 cos sen = 0.

n
En forma matricial, como u : R3 R, u(x, y, z) = k(x2 + y 2 + z 2 ), entonces

u u u
Du = u = , , = (2kx, 2ky, 2kz),
x y z
mientras que T : R2 R3 , T (, ) = (R cos cos , R cos sen , R sen ), implica

R sen cos R cos sen

DT =
R sen sen R cos cos .

R cos 0
on u T : R2 R, definida por u T = u(T (, )), tiene derivada

La composici
u(, ) = u(T (, ))DT (, )

R sen cos R cos sen

= (2kR cos cos , 2kR cos sen , 2kR sen )
R sen sen R cos cos

R cos 0

u u
= (0, 0) = , .

2
w w dz
Ejercicio 92. Utilizar la regla de la cadena para determinar , y , donde
t s dt
a) w = x tan1 (xy), x = t2 , y = set ,
b) w = xey + yex , x = et , y = st2 .

p x
c) x2 + y 2 + z 2 e z = xyz ln(x2 y 2 z 2 ), x = et cos t, y = et sen t.
Ejercicio 93. Sea f : Rn Rm una funcion vectorial diferenciable y sea

u(x) = lnhf (x), f (x)i un campo escalar de n variables. Calcular Du(x).
Ejercicio 94. Sea u = f (x, y), con x = es cos t, y = es sen t. Mostrar que
2 2 " 2 2 #
u u 2s u u
a) + =e + ,
x y t s
2u 2u
2
2u

2s u
b) + =e + 2 .
x2 y 2 t2 s
Ejercicio 95. Calcular ws y wt , donde w = ln(x2 + y 2 + z 2 ) y x = (s t)cos2 t,

y = (s + t)sen2 t y z = 2 st es/t .
p
Ejercicio 96. Supongase que u = f (r), donde r = x2 + y 2 + z 2 . Mostrar que
2u 2u 2u d2 u 2 du
u + + = + .
x2 y 2 z 2 r2 r dr
Ejercicio 97. Sean f : R3 R2 y g : R3 R3 dos funciones vectoriales de

variable vectorial, definidas como sigue:
f (x, y, z) = (x2 + y + z)i + (2x + y + z 2 )j,

g(u, v, w) = uvw2 i + w2 sen vj + u2 ev k.
a) Calcular cada una de las Matrices Jacobianas, Df (x, y, z) y Dg(u, v, w),
b) Determinar la funcion compuesta h(u, v, w) de f con g. Es decir,
h(u, v, w) = f (g(u, v, w) (f g)(u, v, w).
c) Determinar la Matriz Jacobiana de la funcion composicion, Dh(u, 0, w). Que

es mas facil, aplicar la Regla de la Cadena o calcular esta derivada directamente
de b)?
Z x
Ejercicio 98. Discutir porque si y(x) = h(t) sen t(x t) dt, entonces
a
Z x
dy
= y 0 (x) = h(t) cos(x t) dt.
dx a
Z B(x)
d
En general. Sabra Usted como calcular la para (x) = f (x, t) dt?
dx A(x)
Discutir su respuesta, ya que aqu se aplica la Regla de la Cadena.
Ejercicio 99. Considere un recipiente conico circular recto el cual se va a llenar con
agua. Si mientras se llena, en cierto momento la altura del lquido es de 30 cm y esta
aumentando 2 cm por segundo, y en el mismo intante, el radio de la superficie libre
es de 20 cm y esta aumentando 1 cm por segundo, con que rapidez esta aumentando
el volumen del agua contenida en el cono en ese preciso instante?
n
2.2.2 El Teorema de Taylor
Bajo condiciones especficas de diferenciabilidad, una funcion de varias variables

tambien se puede representar mediante una serie de potencias de orden determinado,
para efectos de aproximacion similar al caso de las expansiones familiares correspon-
dientes a las funciones de una variable. La generalizacion del Teorema de Taylor
practicamente proporciona una aproximacion de orden superior para una funcion
f : Rn R, sin restringirse solo a la primera derivada. Sin embargo, cabe
aclarar que el resultado se establece u
nicamente para aproximaciones hasta de tercer
orden debido a que para ordenes mayores la notacion y manipulacion del resultado
se vuelve una tanto tediosa y en terminos vectoriales es limitada.
Teorema 2.2.2 Sea f : Rn R un campo escalar de clase C3 (). Entonces,

n n n
X f 1 X X 2f
f (x0 +h) = f (x0 )+ (x0 ) hi + (x0 ) hi hj +R2 (h, x0 ), (2.2.7)
i=1
xi 2 i=1 j=1 xi xj
R2 (h, x0 )
donde 0 cuando h 0.
khk2
Demostraci
on. La regla de la cadena implica que
n
d X f
f (x0 + th) = hf (x0 + th), hi = (x0 + th) hi .
dt xi
i=1
Integrando ambos lados de esta expresion y aplicando el Teorema Fundamental del Calculo,
sobre el intervalo [0, 1], se obtiene
Z 1 n
X f
f (x0 + h) f (x0 ) = (x0 + th) hi dt.
0 xi
i=1
f
Integrando por partes el lado derecho, tomando u = y v = (t 1), se tiene que
xi
n Z 1 n n X n Z 1
X f X f X 2f
(x0 + th) hi dt = (x0 ) hi + (1 t) hi hj dt,
0 xi xi 0 xi xj
i=1 i=1 i=1 j=1
llegando as a la igualdad
n

X f
f (x0 + h) = f (x0 ) + (x0 )hi + R1 (h, x0 ),

xi

i=1
n n 1 (2.2.8)
2f
XX Z
(1 t)

R1 (h, x0 ) = (x0 + th) hi hj dt,

0 xi xj

i=1 j=1
la cual corresponde a la F
ormula de Taylor de primer orden.
Integrando de nuevo por partes, la expresion R1 (h, x0 ), con
2f (t 1)2
u= (x0 + th) hi hj y v= ,
xi xj 2
n X n
X 1 2f
se obtiene R1 (h, x0 ) = (x0 ) hi hj
2 xi xj
i=1 j=1
n X n X n Z 1
X (t 1)2 3f
+ (x0 + th) hi hj hk dt.
0 2 xi xj xk
i=1 j=1 k=1
As,
n n n

X f 1 X X 2f
f (x0 + h) = f (x0 ) + (x0 ) hi + (x0 ) hi hj + R2 (h, x0 ),

xi 2 xi xj

i=1 i=1 j=1

n X
n X
n Z
(t 1)2
1

3f
X

R2 (h, x0 ) = (x0 + th) hi hj hk dt.

0 2 xi xj xk

i=1 j=1 k=1
(2.2.9)
El integrando es una funci

on continua de la variable t y por lo tanto esta acotada en una
vecindad de x0 . Por lo tanto, para una constante M 0, para khk peque
no, se tiene que
3 |R2 (h, x0 )|
|R2 (h, x0 )| khk M. En particular, khkM 0 cuando h 0, a lo cual
khk2
se deseaba llegar.
Ahora, el Segundo Teorema del Valor Medio para Integrales:

Z b Z b
f (x)g(x)dx = f (c) g(x)dx,
a a
para f y g continuas y g 0 en [a, b], con c (a, b)), implica una forma explcita
para el residuo R2 (h, x0 ), a saber
n X
n X
n 1
(t 1)2 3f
X Z
R2 (h, x0 ) = (x0 + th) hi hj hk dt
i=1 j=1 k=1 0 2 xi xj xk
n X n X n (2.2.10)
X 1 3f
= (cijk ) hi hj hk ,
i=1 j=1 k=1
3! xi xj xk
donde cijk esta sobre el segmento de recta que une x0 con x0 +h. La expresion (2.2.10)
se denomina forma de Lagrange para el residuo.
n
No es dficil inferir la forma del Teorema de Taylor para ordenes mayores que dos.
Por ejemplo, la Formula de Taylor de tercer orden es
n n n
X f 1 XX 2
f (x0 + h) = f (x0 ) + (x0 ) hi + (x) hi hj
i=1
xi 2 i=1 j=1 xi xj
n n n
1 XXX 3f
+ (x0 ) hi hj hk + R3 (h, x0 ),
3! i=1 j=1 k=1 xi xj xk
(2.2.11)
R3 (h, xo )
donde 0 cuando h 0 y as sucesivamente. La formula general puede
khk3
probarse por induccion usando el metodo de demostracion anterior (intentarlo!).
Un buen ejercicio es la manipulacion del Teorema de Taylor general para funciones
vectoriales (trate de escribir la versi
on general!).
Ejemplo 40: Determinar la formula de Taylor de segundo orden para la funcion

f (x, y) = cos(xy) + sen(xy) alrededor del origen, x0 = (0, 0).
Soluci
on. Para la funci
on dada, se tiene
f (x, y) = cos(xy) + sen(xy), f (0, 0) = 1,

f f
(x, y) = y sen(xy) + y cos(xy), (0, 0) = 0,
x x
f f
(x, y) = x sen(xy) + x cos(xy), (0, 0) = 0,
y y
2
f 2
f
2
(x, y) = y 2 cos(xy) y 2 sen(xy), 2
(0, 0) = 0,
x x
2f 2
f
(x, y) = x2 cos(xy) x2 sen(xy), (0, 0) = 0,
y 2 y 2
2f 2f
(x, y) = sen(xy) xy cos(xy) + cos(xy) xy sen(xy), (0, 0) = 1.
xy xy
As, aplicando la f
ormula 2dimensional del Teorema de Taylor,
f f 1 2f 2 2f
f (x0 + h) = f (x0 ) + (x0 ) h1 + (x0 ) h2 + (x0 ) h 1 + (x0 ) h1 h2
x y 2 x2 xy
1 2f
+ (x0 ) h22 + R2 (h, x0 ),
2 y 2
se tiene que f (h) = f (h1 , h2 ) = 1 + h1 h2 + R2 (h, 0).

2
Ejercicio 100. Determinar la formula de Taylor de segundo orden para las funciones
dadas, respecto del punto indicado:
2 y 2
a) f (x, y) = ex cos(xy), x0 = 0, y0 = 0.
ln (3x2 + 2y 2 + z 2 )
b) f (x, y, z) = , x0 = (0, 0, 1).
x2 + y 2 + 2z 2 1
Z x2 y
2
c) f (x, y) = es ds, x0 = (0, 0).
0
[Sugerencia: Ver la parte general del Ejercicio 98].
Ejercicio 101. a) Mostrar que la formula (2.2.7) del Teorema de Taylor puede
expresarse mediante la formula mas compacta
f (x0 + h) = f (x0 ) + f (x0 )h + H f (x0 )h + R2 (h, x0 ), (2.2.12)
donde H f (x0 ) : Rn IR, es el Hessiano de f en x0 , definido por

n n
1 X X 2f
(h1 , h2 , . . . , hn ) 7 (x0 ) hi hj . (2.2.13)
2 i=1 j=1 xi xj
Por ejemplo, para funciones reales de dos variables f (x, y), el Hessiano se escribe
como
2f 2f

x2 xy
!
1 h1
H f (x, y)h = (h1 , h2 ) .
2
2f
h2
2f
yx y 2
2 y 2
b) Calcular el Hessiano, en el origen, de la funcion f (x, y) = x2 2xy +2y 2 +e1+x .
c) Aproximar el valor de la funcion f (x, y) en el punto (0.5, 0.5).
Ejercicio 102. Aproximar cerca del origen con una funcion cuadratica a la funcion
sen(x2 + y 2 )
sinc, f (x, y) = (el sombrero de dama antinguo). Graficar (en Maple,
x2 + y 2
Mathematica, etc.) ambas funciones y comparar geometricamente.
n
2.3 Transformaci
on de Coordenadas
Una alternativa viable en el estudio de alg

un problema especfico es cambiar el sistema
de referencia original, con el fin de transformarlo a un formato en el que sea posible
un mejor analisis del mismo. En particular, en muchos contextos del Calculo Vecto-
rial muchas veces es imprescindible introducir un cambio de variable mediante una
transformacion de coordenadas apropiada, que permita reducir un determinado prob-
lema a una version mas simple en la que sea posible obtener informacion de interes
de una forma mas facil. Este es uno de los temas mas importantes del Calculo Vec-
torial, fundamental para el teorema de cambio de variable en la integracion m
ultiple
(Captulo 3). El objetivo, ademas de obtener los elementos geometricos de distancia,
area y volumen, es determinar la representacion de los principales operadores difer-
enciales cuando se introduce un nuevo sistema de coordenadas, lo cual por supuesto
se puede realizar directamente aplicando la regla de la cadena. Sin embargo, se opta
por adoptar el enfoque geometrico en el que solo se utilizan las propiedades elemen-
tales de las curvas naturales del nuevo sistema. Se trabaja inicialmente en el es-
pacio y rapidamente se generaliza a dimensiones mayores cuando aplica, quedando
por supuesto includos los casos mas simples en dimensiones menores.
Supongamos que las coordenadas rectangulares x, y, z del espacio R3 son expresadas

en terminos de nuevas coordenadas u1 , u2 , u3 mediante las ecuaciones
x = x(u1 , u2 , u3 ), y = y(u1 , u2 , u3 ), z = z(u1 , u2 , u3 ), (2.3.1)
y que, recprocamente, estas relaciones pueden invertirse para expresar u1 , u2 , u3

en terminos de x, y, z, con x, y, z y/o u1 , u2 , u3 apropiadamente restringidas. Esto es,
introducir una transformacion suave T : R3 R3 unoauno (invertible).
Entonces, al menos en alguna region, cualquier punto con coordenadas (x, y, z) tiene
sus correspondientes coordenadas (u1 , u2 , u3 ).
S se considera que una partcula se mueve a partir de un punto p de forma que

u2 y u3 se mantienen constantes (fijas) y u
nicamente u1 vara, entonces se genera una
curva en el espacio a la que se le denomina u1 curva. Similarmente, se obtienen las
otras dos curvas 1parametricas: la u2 curva y la u3 curva, en cada punto. Ademas,
si u
nicamente una de las coordenadas se mantiene constante, entonces se determinan
sucesivamente tres superficies parametrizadas que pasan a traves de un punto del
espacio, las cuales se intersectan en las curvas 1parametricas.
2.3. Transformaci
on de Coordenadas 137
En la Subseccion 1.3.6 se ilustro el caso de las coordenadas cilndricas y esfericas.

Por ejemplo, para las coordenadas cilndricas (r, , z), se tiene
x = x(r, , z) = r cos , y = y(r, , z) = r sen , z = z(r, , z) z,

con r > 0 y 0 < 2.
S u1 = r, u2 = , u3 = z, entonces en cualquier punto p (localmente) la u2 curva es

una circunferencia de radio r mientras que las u1 y u3 curvas son, respectivamente,
segmentos de lnea recta. z
6 uz
HH 6
H : u

H
HH
j
uH
rH
z
...
-
H
. . ... y
HH r ...
.

H ...
. .
......

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .H
.H
j
x
Figura 2.2 Sistema Coordenado Cilndrico

Un caso especial de coordenadas es cuando las curvas 1parametricas son mutua-

mente perpendiculares en cada punto del espacio, las cuales son llamadas coordenadas
curvilneas ortogonales. Esto es,
z u3
6 @
I u2
@
p@

u1
-
y

x

S r es el vector que representa al vector posicion de un punto p en el espacio,

r = xi + yj + zk, (2.3.2)
entonces el vector tangente a la uk curva en el punto p esta dado por
r r dsk
Uk = = , (2.3.3)
uk sk duk
r
donde sk es la longitud de arco de la uk curva. Como es un vector unitario,
sk
entonces Uk = hk uk , con uk el vector tangente unitario a la uk curva en la direccion
dsk
del incremento de la longitud de arco y hk = la magnitud de Uk . Esto es,
duk
n
U1 = h1 u1 , U2 = h2 u2 , U3 = h3 u3 , (2.3.4)
donde uk , k = 1, 2, 3, es el vector tangente unitario a la uk curva y

ds1 r ds 2 r ds3 r
h1 = = , h2 = = , h3 = = . (2.3.5)
du1 u1 du2 u2 du3 u3
Las ecuaciones (2.3.5) pueden escribirse en la forma diferencial

ds1 = h1 du1 , ds2 = h2 du2 , ds3 = h3 du3 , (2.3.6)
donde los coeficientes h1 , h2 , h3 son por naturaleza los factores de escala que establecen
la relacion que existe entre los diferenciales de distancia y los diferenciales de los
parametros coordenados. Cabe se
nalar que el calculo de Uk , de la ecuacion (2.3.3),
determina tanto el factor de escala, hk = kUk k, como el vector unitario uk = Uk /hk .
Por otro lado, es recomendable adoptar siempre la convencion de formar un sistema
coordenado derecho u1 , u2 y u3 mediante un orden apropiado al enumerar las nuevas
coordenadas (ver la nota al final de esta pagina).
Ahora, s s es el parametro natural o longitud de arco a lo largo de una determinada

curva C, entonces el vector
dr r du1 r du2 r du3 du1 du2 du3

= + + = U1 + U2 + U3 . (2.3.7)
ds u1 ds u2 ds u3 ds ds ds ds
tiene longitud unitaria. Ademas, si las curvas coordenadas son ortogonales, entonces
hU1 , U2 i = hU2 , U3 i = hU3 , U1 i = 0, (2.3.8)
de donde se sigue que hU1 , U1 idu1 2 + hU2 , U2 idu2 2 + hU3 , U3 idu3 2 = ds2 . Utilizando
las relaciones (2.3.5), esto es hUk , Uk i = kUk k2 = h2k , k = 1, 2, 3, se obtiene el
elemento de longitud de arco de la curva C,
ds2 = h21 du1 2 + h22 du2 2 + h23 du3 2 . (2.3.9)
Por otra parte, para encontrar el elemento de volumen, observese que los vectores
U1 du1 , U2 du2 y U3 du3 son mutuamente perpendiculares y que tienen por longitudes
los diferenciales de longitud de arco ds1 , ds2 y ds3 , respectivamente. As, el elemento
de volumen dV esta dado por el volumen del paraleleppedo rectangular determinado
2.3. Transformaci
por estos vectores, dV = hU1 du1 U2 du2 , U3 du3 i = hu1 u2 , u3 ih1 h2 h3 du1 du2 du3 ,
y por lo tanto,
dV = h1 h2 h3 du1 du2 du3 . (2.3.10)
Nota: El hecho que hU1 U2 , U3 i es positivo y que hu1 u2 , u3 i = +1 es una consecuencia

de la suposici
on de que se tiene un sistema derecho.
n
Sobre la superficie coordenada u1 = constante, se define el vector elemento de area,

d 1 , mediante el producto vectorial U2 du2 U3 du3 , de tal forma que
d 1 = u1 h2 h3 du2 du3 . (2.3.11)
Similarmente, se obtienen las expresiones correspondientes para las otras superficies

coordenadas, de donde, en especfico, se obtienen los elementos de area sobre las
superficies coordenadas
d1 = h2 h3 du2 du3 , d2 = h3 h1 d u3 du1 , d3 = h1 h2 du1 du2 . (2.3.12)
A continuacion se determinan las expresiones para los principales operadores diferen-

ciales del Caculo Vectorial: el gradiente, la divergencia, el rotacional y el Laplaciano
cuando se introduce un nuevo sistema de coordenadas. Primero, para determinar el
gradiente de un campo escalar f , es necesario considerar la relacion
f f f
df = hf, dri = dx + dy + dz, (2.3.13)
x y z
la cual es una cantidad intrnseca (invariante bajo transformacion de coordenadas).
En terminos de las coordenadas u1 , u2 , u3 , se tiene entonces que
f f f
df = du1 + du2 + du3 .
u1 u2 u3
Como dr = U1 du1 + U2 du2 + U3 du3 = h1 u1 du1 + h2 u2 du2 + h3 u3 du3 , escribiendo
f = 1 u1 + 2 u2 + 3 u3 , con k s escalares por determinar, sustituyendo estas
expresiones en (2.3.13) se tiene que
f f f
du1 + du2 + du3 = h1 1 du1 + h2 2 du2 + h1 3 du3 ,
u1 u2 u3
1 f
y como las variables u1 , u2 , u3 son independientes, entonces k = , k = 1, 2, 3.
hk uk
Por lo tanto, la expresion para el gradiente en el nuevo sistema de coordenadas toma
la forma
1 f 1 f 1 f
f = u1 + u2 + u3 . (2.3.14)
h1 u1 h2 u2 h3 u3
Esto es, la nueva representacion del operador tiene la expresion
1 1 1
= u1 + u2 + u3 . (2.3.15)
h1 u1 h2 u2 h3 u3
2.3. Transformaci
En particular, el gradiente de cada una de las variables independientes u1 , u2 , u3 ,

respectivamente, es
u1 u2 u3
u1 = , u2 = , u3 = . (2.3.16)
h1 h2 h3
Ahora, de estos resultados especiales, se obtienen las representaciones del resto de los
operadores diferenciales. Primero, como uk = 0 (Teorema 2.1.10), se concluye
que u1 /h1 , u2 /h2 y u3 /h3 tienen rotacional cero,
u1 u2 u3
= = = 0. (2.3.17)
h1 h2 h3
Segundo, como u1 = u2 u3 , entonces se sigue que
u1 u2 u3
= = (u2 ) (u3 ),
h2 h3 h2 h3
y analogamente se obtienen las otras dos expresiones. Ahora, haciendo referencia a la
Tabla del Ejercicio 90 (Formula 9), se muestra que u1 /(h2 h3 ), u2 /(h3 h1 ) y u3 /(h1 h2 )
tienen divergencia cero,
u1 u2 u3
= = = 0.
h2 h3 h3 h1 h1 h2
Para encontrar una expresion para la divergencia de un campo vectorial F en las

coordenadas u1 , u2 , u3 ,
F = f1 u1 + f2 u2 + f3 u3 , (2.3.18)
primero se escribe (Formula 4, Tabla del Ejercicio 90)
F = (f1 u1 ) + (f2 u2 ) + (f3 u3 )
y despues se escribe el primer termino en la forma (Formula 6, Tabla del Ejercicio 90)
u1 1
(f1 u1 ) = (h2 h3 f1 ) = (h2 h3 f1 ).
h2 h3 h1 h2 h3 u1
Finalmente, tratando los otros terminos de forma similar, se obtiene

1
div F = F = (h2 h3 f1 ) + (h3 h1 f2 ) + (h1 h2 f3 ) . (2.3.19)
h1 h2 h3 u1 u2 u3
n
En particular, al combinar (2.3.14) y (2.3.19), se determina la expresion para el

operador Laplaciano,

1 h2 h3 h3 h1 h1 h2
= + + . (2.3.20)
h1 h2 h3 u1 h1 u1 u2 h2 u2 u3 h3 u3
Similarmente, para determinar la expresion del rotacional de F primero se expresa

como F = (f1 u1 ) + (f2 u2 ) + (f3 u3 ) y despues se calcula el
rotacional de cada uno
de los
sumandos.
El primer termino se puede expresar como
u1
(f1 u1 ) = (h1 f1 ) , pero de acuerdo a (2.3.17), el segundo factor
h1
en los corchetes tiene rotacional cero. Entonces, utilizando de nuevo la Tabla del
Ejercicio 90 (Formula 7), se obtiene

u1 u3 u2
(f1 u1 ) = (h1 f1 ) = (h1 f1 ) + (h1 f1 )
h1 h1 h2 u
2 h3 h1 u3
1
= h2 u2 h3 u3 (h1 f1 ).
h1 h2 h3 u3 u2
Analogamente, considerando los otros dos terminos y combinando las componentes

correspondientes, el resultado se determina utilizando el determinante

h1 u1 h2 u2 h3 u3

1
rot F = F = . (2.3.21)
h1 h2 h3 u1 u2 u3

h1 f1 h2 f2 h3 f3
Nota: En las coordenadas rectangulares, (u1 , u2 , u3 ) (x, y, z), se tiene que los factores de
escala h1 = h2 = h3 = 1, de manera que todos los resultados obtenidos en esta seccion se
reducen a las f
ormulas conocidas originalmente con las que empieza el estudio del Calculo
Diferencial Vectorial.
Ejemplo 41: Derteminar todas las cantidades diferenciales involucradas cuando se

trabaja con el sistema de coordenadas cilndricas circulares (r, , z),
x = r cos , y = r sen , z = z, (2.3.22)
con r 0 y 0 2.
2.3. Transformaci
Soluci
on. El orden en el que se consideran las nuevas coordenadas es r, y z, ya que
as forman un sistema derecho. Entonces, el vector posicion tiene la forma r = r cos i +
r sen j + zk y los vectores tangentes del sistema son
r r
Ur = = cos i + sen j, U = = r sen i + r cos j, Uz = k.
r
De la ecuaci
on (2.3.5), se sigue que
hr = 1, h = r, hz = 1, (2.3.23)
de donde
ur = cos i + sen j, u = sen i + cos j, uz = k. (2.3.24)
Las ecuaciones (2.3.9), (2.3.10) y (2.3.12), respectivamente, implican
q
ds = [dr]2 + r2 [d]2 + [dz]2 ,
dV = r ddrdz, (2.3.25)
dr = r ddz, d = drdz, dz = z ddr.
Finalmente, las ecuaciones (2.3.14), (2.3.19) (2.3.20) y (2.3.21) implican los resultados
f 1 f f
f = ur + u + uz ,
r r z
1 1 f fz
F = (rfr ) + + ,
r r r z
1 2f 2f

1 f
f = r + 2 2 + 2,
r r r r z
(2.3.26)

ur ru uz

1
F = .
r r z

fr f fz
2
n
Ejercicio 103. En las coordenadas esfericas, ordenadas positivamente (, , ), con

la distancia desde el origen, el angulo conico medido desde el eje z positivo y
el angulo polar medido sobre el plano xy, la transformacion es de la forma
x = sen cos , y = sen sen , z = cos , (2.3.27)
donde 0, 0 y 0 2.
a) Graficar las curvas y superficies coordenadas del sistema esferico.
b) Calcular, en este sistema de coordenadas, el diferencial de longitud de arco, ds,

el diferencial de volumen, dV y los diferenciales de area sobre las superficies
coordenadas, d , d y d .
c) Verificar que
f 1 f 1 f
f = u + u + u ,
sen
1 2 1 1 f
F = 2
( f ) + (f sen ) + ,
sen sen

u u sen u
(2.3.28)

1
F = 2 .
sen

f f sen f
d) En particular, mostrar que cuando la ecuacion de Laplace se escribe en el sistema

de coordenadas esfericas tiene la expresion
2f

1 2 f 1 f 1
f = 2 + 2 sen + 2 =0
sen sen2 2
o
2f 2 f cot f 1 2f 1 2f
+ + + + = 0.
2 2 2 2 2 sen2 2
(2.3.29)
2.3. Transformaci
Ejercicio 104. Las coordenadas cilndricas elpticas se definen por
x = a cosh u cos v, y = a sen hu sen v, z = z, (2.3.30)
con u 0 y 0 v 2.
a) Demostrar que este sistema de coordenadas (u, v, z) es ortogonal, verificando

que se satisface (2.3.8).
b) Mostrar que en el plano xy las curvas uconstante son elipses con semiejes
a cosh u, a senh u, en la direccion x, y, respectivamente; tambien que la curva
vconstante es la mitad de una rama de la hiperbola con semiejes a cos v y
a sen v. En particular, mostrar que el lugar geometrico correspondiente a u = 0
degenera en el segmento (a, a) del eje x, mientras que el lugar geometrico
correspondiente a v = 0 y v = es, respectivamente, el exterior positivo y
negativo de este segmento; tambien que el lugar geometrico correspondiente a
v = /2 y v = 3/2 es, respectivamente, la porcion positiva y negativa
del eje y. Bosquejar y se
nalar en un solo diagrama las curvas u = 0, 1 y
v = 0, /4, /2, 3/4, 3/2, 5/4 y 7/4.
c) Verificar que en este sistema coordenado, se tiene

p
hu = hv = a cosh2 u cos2 v, hz = 1, (2.3.31)
i senh u cos v + j cosh u sen v
u1 = p , (2.3.32)
cosh2 u cos2 v
i cosh u sen v + j senh u cos v
u2 = p , (2.3.33)
cosh2 u cos2 v

1 f f f
f = p u1 + u2 +k , (2.3.34)
2
a cosh u cos2 v u v z
2
2f 2f

1 f
f = 2 + + . (2.3.35)
a cosh2 u cos2 v u2 v 2 z 2
d) Mostrar que para valores grandes de u se sigue que
u1 i cos v + j sen v, u2 i sen v + j cos v

y
u1 cos v u2 sen v i, u1 sen v + u2 cos v j.
n
Ejercicio 105. Las coordenadas cilndricas parabolicas se definen por las ecuaciones
1 2
u v2 ,

x= y = uv, z = z, (2.3.36)
2
donde u y v 0.
a) Demostrar que este sistema es ortogonal.
b) Mostrar que en el plano xy una curva vconstante es una parabola simetrica

con respecto al eje x y abierta hacia la derecha, mientras que una uconstante
es una mitad de una parabola similar abierta hacia la izquierda. En particular,
mostrar que el lugar geometrico correspondiente a v = 0 es el eje x positivo,
mientras que el lugar geometrico correspondiente a u = 0 es el eje x negativo y
que el eje y positivo es dado por u = v y el eje y negativo por u = v. Bosquejar
nalar en un solo dibujo las curvas u = 0, 1 y v = 0, 1.
y se
c) Realizar los calculos necesarios para mostrar que el Laplaciano en este sistema
coordenado es de la forma
2f 2f 2f

1
f = 2 + + . (2.3.37)
u + v2 u2 v 2 z 2
Ejercicio 106. En una traslacion y rotacion de ejes, en la cual el origen de las nuevas
variables del espacio x0 , y 0 , z 0 se toma en el punto (a, b, c) en el espacio xyz y en el
cual las direcciones de los ejes x0 , y 0 y z 0 son especificados por las direcciones de los
cosenos directores (l1 , m1 , n1 ), (l2 , m2 , n2 ) y (l3 , m3 , n3 ) relativos a los ejes originales,
la trasformacion de coordenadas es de la forma
x = a + l1 x0 + l2 y 0 + l3 z 0 ,
y = b + m1 x0 + m2 y 0 + m3 z 0 ,
x = c + n1 x0 + n2 y 0 + n3 z 0 .
Demostrar que la ecuacion de Laplace es invariante bajo este tipo de transformacion.

Es decir, la ecuacion de Laplace es de la forma
2f 2f 2f
+ + = 0,
x0 2 y 0 2 z 0 2
en terminos de las nuevas variables.
2.3. Transformaci
Ejercicio 107. Si r, y z son las coordenadas cilndricas circulares, entonces evaluar

las cantidades siguientes:
a) , b) rn , c) (rn cos n) .
Ejercicio 108. Si u , u y u son los vectores tangentes unitarios a las curvas

coordenadas del sistema esferico, entonces mostrar que
i = (u sen + u cos ) cos u sen ,
j = (u sen + u cos ) sen u cos ,
k = u cos u sen .
Ejercicio 109. Si , y son las coordenadas esfericas, entonces evaluar las

cantidades siguientes:
1
a) , b) , c) r + 2 cos .
r
Ejercicio 110. Mostrar que los vectores tangentes unitarios en coordenadas esfericas
satisfacen las relaciones siguientes:
u u
= u , = u sen ,

uv a u
= u , = u cos ,

u u
= 0, = u sen u cos .

Ejercicio 111. Usando los resultados del Ejercicio 110 y escribiendo el vector posicion
en la forma r = rur , obtener expresiones para las componentes de la aceleracion a lo
largo de las curvas coordenadas de las coordenadas esfericas como sigue:
a = 2 2 sen2 ,
a = 2 + 2 sen cos ,
a = 2 sen + sen + 2 cos .

n
Ejercicio 112. Las coordenadas elipsoidales (u, , ) son definidas por

x = au sen cos , y = bu sen sen , z = cu cos ,
donde u > 0, 0 , 0 < 2 y a, b, c son constantes positivas.
a) Verificar que el sistema es curvilneo ortogonal.
b) Mostrar que el elemento de volumen en este sistema coordenado es de la forma

dV = abc u2 sen dudd.
c) Determinar las expresiones para la divergencia y el Laplaciano en este sistema

de coordenadas.
d) Mostrar que
ux uy uz
(u, , )
hu, i = det x y z

(x, y, z)
x y z
def 1
= Jacobiano de la Transformacion Inversa = .
hu h h
Ejercicio 113. Las coordenadas toroidales se definen por
x = (R0 + u cos ) cos , y = (R0 + u cos ) sen , z = u sen ,

donde R0 = cte., 0 u a R0 y 0 2, 0 2. Trabajar los
incisos a) d) del Ejercicio 112 para este sistema de coordenadas especiales.
Ejercicio 114. (Primer Proyecto de Investigacio

n Simple). Insvestigar so-
bre alg
un sistema de coordenadas de interes particular (distinto de los previamente
abordados en esta seccion), trabajar todos los conceptos diferenciales sobre el nuevo
sistema de coordenadas y describir su aplicacion directa en alg
un problema fsico y/o
geometrico interesante.
Ejercicio 115. (Segundo Proyecto de Investigacio

n Simple). Investigar
sobre el Teorema de la Funcion Implcita (enunciado y demostracion), en particular,
sobre uno de sus principales corolarios: el Teorema de la Funcion Inversa. Es muy
importante que el estudiante comprenda que es lo que establece (o en otras palabras,
que aprenda a utilizarlo) y de ser posible que tenga una nocion general sobre el
bosquejo de su demostracion (de preferencia con su propio lenguaje). En la siguiente
seccion se analiza este resultado desde un enfoque muy particular del autor.

Cap 2 Calculo IVRGG

Cargado por

Copyright:

Formatos disponibles

Cap 2 Calculo IVRGG

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Cap 2 Calculo IVRGG

Cargado por

Copyright:

Formatos disponibles

Captulo 2

En este captulo se establece de forma general y natural el concepto de diferenciaci

2.1 El Concepto General de Diferenciaci

Una funcion f : [a, b] R R se dice ser diferenciable en x0 (a, b) si existe

donde RT (x) = f (x0 ) + f 0 (x0 )(x x0 ) es la expresion de la recta tangente a la grafica

Geometricamente, la igualdad (2.1.2) significa que la recta tangente RT (x) es cercana

obteniendo as la definicion formal de diferenciabilidad.

on 2.1.1 Una funcion f : [a, b]

on 2.1.2 Una funcion f : Rn Rm es diferenciable en x0 , donde

ii) Aunque parezca una redundancia, cuando n = m = 1 la Definicion 2.1.2 ciertamente es

Un resultado conocido de la teora de operadores establece que todo operador lineal

Teorema 2.1.1 S f : Rn Rm es diferenciable en x0 , entonces la trans-

formacion lineal L : Rn Rm que aproxima al cambio instantaneo de sus valores es

Ahora, si x Rn , entonces tx 0, cuando t 0. As, para x 6= 0, se tiene que

En el caso particular de las funciones vectoriales que dependen u

on 2.1.3 La derivada de una trayectoria, (t) = (f1 (t), f2 (t), , fm (t))T ,

existe para todo t0 (a, b). En tal caso, (t

Ejemplo 24: a) (t) = (4t2 , tet , 12 sen t)T

b) Si (t) = (ln(2 t2 ), |t 1|, 7t2/3 1)T , entonces (1)

es paralelo al vector tangente a en p.

0 ) es paralelo a la recta tangente RT ( ) de la curva descrita

Ejercicio 47. Calcular los vectores velocidad y aceleracion as como la ecuacion de

3). El vector velocidad de una trayectoria de rapidez constante, es perpendicular

Longitud de Arco: La longitud de arco de una curva es un parametro natural,

Definici on 2.1.4 Sea : [a, b] R R3 una trayectoria suave2 en el espacio,

Ejemplo 25: i) La longitud de una circunferencia de radio r es igual a 2r.

ii) Determinar la longitud de un arco de la cicloide.

As, l[arco de la cicloide] = 8a.

Ejercicio 50. Determinar la longitud de la astroide (hipocicloide) del Ejercicio 23.

Ejercicio 51. a) Determinar la trayectoria tal que (1) = (1, 0, 1)T y

b) S una partcula se mueve sobre la curva (t) y en el tiempo t = 1 se sale por

Ejercicio 52. Calcular la longitud de la curva dada en el intervalo especificado:

ii) (t) = (t, t sen t, t cos t)T , [0, ].

a) Determinar la ecuacion de la recta tangente a la curva que forme con el eje x

b) Determinar la longitud de la curva dada.

on 2.1.5 Sea h : I1 I una funcion C1 y sea : I

La funcion longitud de arco o parametro natural,

Teorema 2.1.2 S : I Rm es de clase C1 y s(t) es el parametro natural, entonces

Ejercicio 54. Demostrar el Teorema 2.1.2

Entonces, por la Propiedad Basica 3) de la Derivada de una Trayectoria, el vector

mientras que el vector binormal unitario esta dado por

La terna de vectores ortonormales LI (T, N, B), llamada triedro intrnseco o aparato

En terminos de los tres vectores LI T, N, B se tiene que

Las propiedades del triedro (T, N, B) implican que

de donde = = 0 y = , cuyo valor com un se denota por y representa la

Nota: Las ecuaciones de Frenet-Serret (2.1.20) forman un sistema de tres ecuaciones

Ejercicio 57. Demostrar que la curvatura y la torsion de una curva parametrizada

Ejercicio 58. Determinar la curvatura y la torsion de la helice

(t) = (a cos t, a sen t, bt)T , a, b > 0.

Ejercicio 59. Si a es una constante diferente de cero, calcular la curvatura y la

En la Subseccion 1.3.2 se define y se proporcionan algunos ejemplos de funciones

si el lmite existe, donde i = 1, 2, , n y ni es el iesimo vector de la base natural

De forma similar a la interpretacion que tiene la derivada de una funcion de una

Es decir, se tiene igualdad en las parciales mixtas7 .

f ) S u(x) = eha,xi , x, a Rn , con kak = k. Demostrar que la funcion u satisface la

g) Si f (x) = kxk(2n) , 0 6= x Rn y R. Demostrar que la funcion f satisface

de donde se tiene que f = 0.

Ejercicio 63. Verificar que zxy = zyx , donde