Funciones de varias variables

por Charlie

1 Funciones de dos variables

En Cálculo I se estudian funciones de una variable y = f (x). Ahora estudiaremos funciones de
varias variables. Estas funciones aparecen naturalmente en la naturaleza, en la economía y en los
más variados campos de la ciencia. Por ejemplo:

 La temperatura en un determinado punto de la tierra depende de la latitud, la longitud y

la altitud, es decir, depende de tres variables.

 El beneficio de un determinado producto depende del costo de la materia prima, del costo
de la mano de obra y en algunos casos de otros costos adicionales.

Veremos primero las funciones de 2 variables.

Definición 1. Dado un subconjunto D  R2 , una función f de dos variables es una regla que
asocia a cada par (x; y)2D un único valor real denotado por f (x; y). El conjunto D se llama
dominio de f y su imagen es el conjunto ff (x; y): (x; y)2DgR. Denotamos una tal función
como f : D ! R.

Cuando no mencionamos el dominio D, entendemos que el dominio de f es el subconjunto más

grande de R2 para el que se aplica la regla definida por f . También escribimos, a menudo, z = f (x;
y) para los valores asumidos por f en un punto cualquiera (x; y) de su dominio. Las variables x e
y son las variables independientes y z la variable dependiente.
p
Ejemplo 2. Determine el dominio e imágen de z = f (x; y) = 16 ¡ x2 ¡ y 2
El dominio de f es el conjunto de todos los puntos (x; y) para los que 16¡x2 ¡ y 2 0, luego

D = f(x; y): x2 + y 2 6 16g

Que es un disco de centro (0; 0) y radio 4. La imagen de f es el conjunto

p
fz: z = 16 ¡ x2 ¡y 2 y (x; y) 2 Dg = [0; 4]:

Ejemplo 3. Determine el dominio de

p
x¡1
z = f (x; y) =
x¡y

El dominio de f es el econjunto D = f(x; y): x > 1 y x =

/ yg.

Ejemplo 4. Determine el dominio e imágen de z = f (x; y) = ln(y ¡ x2)

La función ln(y ¡x2) se define solo para los puntos (x; y), tal que, y ¡ x2 >0 o y >x2. La parábola
y =x2 divide el plano R2 en dos regiones, una donde y <x2 y otra donde y >x2. Para comprobar
cuál es la región en la que y >x2, se toma un punto cualquiera para la prueba, por ejemplo el punto
T =(0; 2) y se comprueba fácilmente que este punto pertenece a la región que queremos y >x2. Como
y ¡ x2 >0, la imagen de la función ln(y ¡x2) es todo R, que es la imagen de la función logarítmica.

Veamos las regiones en Sagemath.

1
>>> var("x y z t a b c h k")
(x, y, z, t, a, b, c, h, k)
>>> region_plot(16-x^2-y^2>=0, (x, -5, 5), (y, -5, 5),aspect_ratio=1,
incol="lightblue",bordercol="black")

−4 −2 2 4

−2

−4

>>> region_plot([x-1 >=0, lambda x,y: x<y or x>y], (x, -5, 5), (y, -5,
5),aspect_ratio=1,incol="lightblue",bordercol="black")

−4 −2 2 4

−2

−4

2
>>> region_plot(y-x^2 > 0, (x, -5, 5), (y, 0, 5),incol="lightblue",
bordercol="black",borderstyle="dashed")

−4 −2 2 4

>>>

2 Gráficos y curvas de nivel

La gráfica de una función de dos variables z = f (x; y) con dominio D  R2 es el conjunto:

f(x; y; z) 2 R3: z = f (x; y); 8(x; y) 2 Dg

Este conjunto es una superficie en el espacio, normalmente llamada superficie z = f (x; y).

Para visualizar dicha superficie, utilizamos la misma técnica que se utilizó para trazar las superficies
cuadráticas. Primero, determinamos sus rasgos, es decir, las curvas de intersección de esta superficie
con los planos coordenados. Si no logramos una visualización razonable, determinamos las curvas
de intersección con planos paralelos a los planos coordenados, por ejemplo, z =k, o x=k o y =k.
Variamos k tantas veces como sea necesario para una mejor visualización de la superficie.

Si aún así no podemos visualizar la función, lo mejor es utilizar un ordenador con un buen programa
gráfico para trazar superficies (por ejemplo, SageMath, Maple, Matemathica, Winplot, Maxima, etc...).

Ejemplo 5. Entre las superficies cuadráticas, los paraboloides son gráficos de función.

Observación 6. Para comprobar si una superficie en el espacio es gráfica de una función, basta con
trazar una recta perpendicular al plano x ¡ y y comprobar cuántos son los puntos de intersección
de esta recta con la superficie. Si tiene más de un punto de intersección es porque a cada par (x; y)
corresponde a dos valores de z y, por lo tanto, la superficie no es gráfica de una función.

Como podemos ver, trazar el gráfico de una función, a menos que tengamos una computadora con
un buen programa gráfico, no es una de las tareas más fáciles. Además, otra forma de obtener
información sobre el comportamiento de una función es a través de sus curvas de nivel, una técnica
muy utilizada en cartografía para ver el mapa topográfico de una región y ver la variación de altitud
de un terreno.

Una curva de contorno es la curva de intersección de una superficie z = f (x; y) con el plano z = k.
Su proyección vertical en el plano x ¡ y es la llamada curva de nivel f (x; y)=k de la función f .
Al variar la constante k obtenemos un conjunto de curvas en el plano x ¡ y, llamado mapa de

3
contorno. El mapa de contorno es la representación bidimensional de una superficie tridimensional
z = f (x; y).

Ejemplo 7. Las curvas de nivel de la función f (x; y) = x2 + y 2 son las circunferencias concéntricas
x2 + y 2 = k.

Ejemplo 8. Las curvas de nivel de la función f (x; y) = x2 + y son las parábolas trasladadas
x2 + y = k, o y = k ¡ x2.

Ejemplo 9. Las curvas de nivel de la función f (x; y) = x2 ¡ y 2 son las hipérbolas x2 ¡ y 2 = k.

Ejemplo 10. Las curvas de nivel de la función f (x; y) = cos(x) ¡ sin(y) son las curvas cos(x) ¡
sin(y) = k.

>>> contour_plot(x^2+y^2,(x,-3,3),(y,-3,3),colorbar=True,cmap="jet",labels=True,
label_inline=True)
>>> contour_plot(x^2+y,(x,-3,3),(y,-3,3),colorbar=True,cmap="jet",labels=True,
label_inline=True)
>>> contour_plot(x^2-y^2,(x,-3,3),(y,-3,3),colorbar=True,cmap="jet",labels=True,
label_inline=True)
>>> contour_plot(cos(x)-sin(y),(x,-5,5),(y,-5,5),colorbar=True,cmap="jet",
labels=True,label_inline=True)

2.0
0

1.
1.0

0
0

00

1.
.5
.5

4 50
1.5
0
0

0 0
.5 .0
0 0 0
50 1.0
. 0
-0 .0
-1
1.0

0.5
-1
.0
0 -1.50 -1.50
0 0. 0
0 .0 50 .5 0.0
0 -0
1.50
0
.5
0

−0.5
1.00

1.00

−2
0
.5
0

−1.0
00
0 0 1.
.5 -1 .0
-0 .0 0 -1
0 .0
0
0

−1.5
.5

−4
-0
-1.50

-1.50

−2.0
−4 −2 0 2 4

>>>

3 Límites
Intuitivamente una función f : U  R2 ! R, z = f (x; y) tiene límite L en (a; b) 2 U , si podemos
hacer que f (x; y) esté tan próximo de L cuanto queramos haciendo (x; y) 2 U próximo suficiente
de (a; b), pero distinto de (a; b).

4
Esto significa, matemáticamente, que podemos hacer jf (x; y)¡Lj tan pequeño cuanto queramos,
haciendo k(x; y)¡(a; b)k lo suficientemente pequeño, pero no nulo. A su vez, hacer jf (x; y)¡Lj tan
pequeño cuanto queramos significa que dado cualquier  > 0, podemos hacer que jf (x; y)¡Lj < ,
cuando k(x; y)¡(a; b)k es lo suficientemente pequeño y no nulo, lo que se puede expresar como
0 < k(x; y) ¡ (a; b)k < , para algún  > 0. Geométricamente, esta última condición en el plano
x ¡ y significa que el punto (x; y) está en el llamado disco perforado D(a; b) = f(x; y): 0 <
(x ¡ a)2 +(y ¡ b)2 < 2g.

Definición 11. Sea f : U  R2 ! R, z = f (x; y) definida sobre un conjunto U  R2. Decimos

que el límite de f en (a; b) es L y excribimos

lim f (x; y)=L

(x;y)!(a;b)

si para todo número  > 0, hay un número  > 0, tal que,

jf (x; y)¡Lj < ; 8(x; y) 2 D(a; b) \ D:

Geométricamente, considerando el gráfico de la función f , dado  > 0 la condición jf (x; y)¡Lj < 
es equivalente a L ¡ < f (x; y)<L+, lo que significa que la gráfica de f para esos puntos está
entre los planos z = L¡ y z = L + .

Por otro lado, en el plano x ¡ y, los puntos (x; y) tales que 0 < k(x; y)¡(a; b)k < , son los puntos
del disco perforado D(a; b) de centro (a; b) y radio  (se excluye el centro (a; b) del disco habitual
D(a; b))

La definición establece entonces que si el límite de f en (a; b) es L, para cualquier  > 0, existe
un  > 0, tal que la función f lleva el disco perforado D(a; b) en un subconjunto del intervalo
(L¡; L+), es decir, a medida que disminuimos el valor de , el punto (x; y) se acerca a (a; b). Es
decir, el punto (x; y) se acerca a (a; b) por cualquier camino (curva) y no solo por los lados derecho
e izquierdo como sucede en una dimensión. Así, para demostrar que una función no tiene límite
en un punto (a; b), basta encontrar dos caminos que lleven (x; y) a (a; b) con diferentes límites.

Observación 12. La función f no necesita estar definida en el punto (a; b), pero para garantizar
la consistencia de la definición la función debe estar bien definida en puntos (x; y) muy próximos de
(a; b). Cuando este es el caso decimos que el punto (a; b) es un punto de acumulación del dominio
de la función. Formalmente esto significa que dado cualquier  > 0, el disco perforado D(a; b)
siempre tiene puntos del dominio de f .

p 
sin x2 + y 2
Ejemplo 13. Sea z = f (x; y) = p . La función f se define en todo R2 con la excepción
x2 + y 2
del origen. Observe su gráfico y vea que visiblemente nos parece una superficie suave definida en
todos los puntos del plano, incluso en el punto (0; 0), que es el punto donde la función no está
definida. Veamos qué sucede si tomamos un punto cualquiera cerca del origen y hacemos que este
punto se acerque al origen.

Tenga en cuenta que las curvas de nivel son circunferencias concéntricas con centro en el origen
con ecuación
p 
sin x2 + y 2
p =k
x2 + y 2

5
 p
Así, si tomamos circunferencias con radio r = x2 + y 2 , cuando r tiende a cero

p 
sin x2 + y 2 sin(r)
p = !1
x +y
2 2 r

Esto significa que cualquier punto (x; y) que tienda a (0; 0) por cualquier camino, el límite siempre
será 1.

>>> contour_plot(sin(sqrt(x^2+y^2))/(sqrt(x^2+y^2)),(x,-3,3),(y,-3,3),
colorbar=True,cmap="jet",labels=True,label_inline=True)

1.0
3
00
0.
0
.0
0
0.8
20
2 0. 0.
40

0.6

0.4

0.2

−1
0.80

0.0
0.6
0

−2
−0.2
0
.0
0

00
0.
−3
−0.4
−3 −2 −1 0 1 2 3

>>>

x2 ¡ y 2
Ejemplo 14. Calcula lim (x;y)!(0;0) x2 + y2

La función no está definida en (0; 0). Observando su gráfico vemos un agujero cercano al origen,
pero también un comportamiento irregular cerca del origen. Calculamos el límite por caminos

 Camino por el eje x (y = 0):

x2 ¡ 02 x2
lim = lim 2 = 1
(x;0)!(0;0) x + 0
2 2 x!0 x

 Camino por el eje y (x = 0):

02 ¡ y 2 ¡y 2
lim = lim 2 = ¡1
(0;y)!(0;0) 0 + y
2 2 y!0 y

Vemos, por lo tanto, que el límite tiene diferentes valores para diferentes caminos y, por lo tanto,
el límite no existe.

6
>>> contour_plot((x^2-y^2)/(x^2+y^2),(x,-3,3),(y,-3,3),colorbar=True,cmap="jet",
labels=True,label_inline=True)

1.00
3

0
.2

.5
-0

.7
0.75

-0
-0
2 0.
25 0
.0
0
0.50
50
0.

1
0.25
0.7
5

0 0.00

0.7
5

−0.25
−1

50
0.
−0.50
0.
25
−2
5
0

.7

5 −0.75
.5

.2
-0

0
-0

.0 -0
0

−3
−1.00
−3 −2 −1 0 1 2 3

>>> f(x,y)=(x^2-y^2)/(x^2+y^2)
>>> f
x2 ¡ y 2
(x; y) 7!
x2 + y 2
>>> f(0,x)
¡1
>>> f(0,y)
¡1
>>> f(x,x)
0
>>>
Vemos, por lo tanto, que el límite tiene diferentes valores para diferentes caminos y, por lo tanto,
el límite no existe.

x2y
Ejemplo 15. Calcula lim (x;y)!(0;0) x4 + y2

La función no está definida en (0; 0). Observando su gráfico vemos un agujero cercano al origen,
pero un comportamiento irregular cerca del origen. Calculamos el límite por caminos

 Camino por el eje x (y = 0):

x2
0
lim = lim 0 = 0
(x;0)!(0;0) x4 + 02 x!0

 Camino por el eje y (x = 0):

0
y
lim = lim 0 = 0
(0;y)!(0;0) 04 + y 2 y!0
 Camino lineal y = kx:
x2(kx) kx
lim = lim = 0;
(x;kx)!(0;0) x4 + (kx)2 x!0 x2 + k

7
 para cualquier valor de k. Por lo tanto, el límite es cero por cualquier recta que pase por el
origen.

 Camino cuadrático y = x2:

x2 x2 x4 1
lim = lim =
(x;0)!(0;0) x4 + (x2)2 x!0 2x4 2

Por lo tanto, cuando un punto (x; y) va al origen sobre una parábola, el límite es igual a 1/2.

Tenemos diferentes valores para el límite cuando tomamos caminos diferentes. Por lo tanto, el
límite no existe.
>>> contour_plot((x^2*y)/(x^4+y^2),(x,-1,1),(y,-1,1),colorbar=True,cmap="jet",
labels=True,label_inline=True)

0.60
1.0 30

45
5

0.45
0.1

0.
0.
0.
54

0.5 0.30
0.
03

0.15

0.15

0.0 0.00 0.00

-0.15

−0.15

5
.4
−0.5 -0
−0.30
0
.3
-0

-0.

-0

-0
15

.3

−0.45
.4
0

−1.0
−0.60
−1.0 −0.5 0.0 0.5 1.0

>>> f(x,y)=(x^2*y)/(x^4+y^2)
>>> f
x2 y
(x; y) 7!
x4 + y 2
>>> f(0,y)
0
>>> f(x,0)
0
>>> f(x,x)
x3
x4 + x2
>>> limit(f(x,x),x=0)
0
>>> f(x,x^2)
1
2

8
>>>

2x2y
Ejemplo 16. Calcula lim (x;y)!(0;0) 3x2 + 3y2

La función no está definida en (0; 0). Observando su gráfico vemos un agujero cercano al origen,
pero también un comportamiento más regular cerca del origen que las anteriores funciones. En
este caso el gráfico de f parece una superficie completamente plana, lo que nos lleva a creer que
el límite por cualquier camino hacia el origen dará el valor cero.
Calculamos el límite por caminos

 Camino por el eje x:

2x2 0
lim = lim 0 = 0
(x;0)!(0;0) 3x2 + 3
02 x!0

 Camino por el eje y:

2
0y 2
lim = lim 0 = 0
(0;y)!(0;0) 3
02 + 3y 2 y!0

 Camino lineal y = kx:

2x2(kx) 2kx
lim = lim =0
(x;0)!(0;0) 3x2 + 3(kx)2 x!0 3 (1 + k 2)

para cualquier valor de k. Por lo tanto, el límite es cero en cualquier recta que pase por el
origen.

 Camino cuadrático y = x2:

2x2 x2 2x2
lim = lim =0
(x;0)!(0;0) 3x2 + 3(x2)2 x!0 3(1 + x2)

Evidentemente, esto no prueba nada, pues no podemos hacer esto para todos los caminos que
van para el origem.
Una técnica para el cálculo del límite que podemos utilizar en este caso, es la de la sustitución de
variables utilizando coordenadas polares. Sean, entonces, (r; ) las coordenadas polares del punto
(x; y), es decir, x=r cos  e y =r sen . Así,

2x2y 2r3cos2 sin  2

= = r cos2 sin ;
3x + 3y
2 2 3r2 3
p
donde r = x2 + y 2 y claro que r ! 0 cuando (x; y) ! (0; 0). Por lo tanto

2x2y 2
lim = lim r cos2 sin  = 0
(x;y)!(0;0) 3x2 + 3y 2 r!0 3

pues jcos2 sin j 6 1 para todos los valores de .

Otra técnica que se puede aplicar es la del Teorema del Sándwich. Tenga en cuenta que
p p
2x2y jyj y2 x2 + y 2
6 2 = 2 6 2
3x2 + 3y 2 3 3 3

Así, como la última expresión va a cero cuando (x; y)!0, entonces, f (x; y)!0.

9
>>> contour_plot((2*x^2*y)/(x^2+y^2),(x,-1,1),(y,-1,1),colorbar=True,cmap="jet",
labels=True,label_inline=True)
1.00
1.0

0
5

0.75
0.5
0.75

0.2
0.2
5
5
.7 0.50
0.5 0

0.50
0.25

0.0 0.00 0.00

-0.25

−0.25

−0.5 −0.50
-0
.7
5

-0.75
0

-0.5
-0.2
-0.5

−0.75
5

−1.0
−1.00
−1.0 −0.5 0.0 0.5 1.0

>>>

4 Continuidad
La definición de continuidad para funciones de varias variables es análoga a la de función de una
variable:

Definición 17. Sea f : U  R2 ! R, z = f (x; y) definida sobre un conjunto U  R2. Decimos

que f es continua en (a; b) 2 U, si para todo número  > 0, hay un número  > 0, tal que,

jf (x; y)¡ f (a; b)j <  8(x; y) 2 D(a; b) \ D:

Cuando el punto (a; b) es un punto de acumulación de D, esta condición es equivalente a que

lim f (x; y)=f (a; b)

(x;y)!(a;b)

Observación 18. Cuando el punto (a; b) es un punto de acumulación de D, se requiere el cum-

plimiento de tres condiciones

1. El punto (a; b) debe estar en U , el dominio de f , para que esté bien definido f (a; b)

2. Debe existir el límite lim (x;y)!(a;b) f (x; y)

3. Debe verificarse la igualdad

lim f (x; y)=f (a; b):
(x;y)!(0;0)

Cuando el punto (a; b) 2 U , no es un punto de acumulación de U , decimos que es un punto aislado

de U . En este caso la función siempre es continua en (a; b).
p 
sin x2 + y 2
La función z = f (x; y) = p no está definida en el origen, luego no puede ser continua en
x2 + y 2
ese punto, aunque el límite de f existe y es 0 y su gráfica parece no tener un agujero en el origen.

10
Ahora bien, si definimos una nueva función con el valor del límite, por ejemplo, la función

8 p 
>
< sin x2 + y 2
p ; (x; y) =
/ (0; 0)
g(x; y) = x2 + y 2
>
:
0 (x; y) = (0; 0)

esta función es continua en el origen, pues esta función está definida en ese punto, tiene el mismo
límite que f en ese punto y vale lim (x;y)!(0;0)g(x; y) = 0 = g(0; 0). En este caso decimos que el
origen es una discontinuidad removible.
x2 ¡ y 2
Por otro lado, la función z = f (x; y) = x2 + y2 tampoco está definida en el origen, luego no es continua
en ese punto. La diferencia con el ejemplo anterior es que para esta función el límite en el origen
no existe.

Si definimos la función g como

8
< x2 ¡ y 2
; (x; y) =
/ (0; 0)
g(x; y) = x2 + y 2
:0 (x; y) = (0; 0);

esta función no es continua en el origen, pues si bien esta función está definida en ese punto, tiene
el mismo límite que f en ese punto, luego en este caso lim (x;y)!(0;0)g(x; y) no existe. En este caso
decimos que el origen es una discontinuidad no removible.

Todas las definciones y propiedades presentadas, excepto aquellas referidas a las gráficas de fun-
ciones, pueden extenderse al caso de funciones de más de dos variables.

Definición 19. Dado un subconjunto D  Rn, una función f de n variables es una regla que
asocia a cada punto x 2 D un único valor real denotado por y = f (x). El conjunto D se llama
dominio de f y su imagen es el conjunto ff (x): x2DgR. Denotamos una tal función como
f : D ! R.

En coordenadas escribimos los puntos de x 2 Rn, como x = (x1; : : : ; xn) y a menudo y = f (x), como
y = f (x1; : : : ; xn) para los valores asumidos por f en un punto cualquiera x = (x1; : : : ; xn) de su
dominio. Las variables x1; : : : ; xn son las variables independientes y y la variable dependiente.

En el caso n = 3, los puntos de R3 se suelen escribir en coordenadas como (x; y; z) y una función
f como w = f (x; y; z).

La nominación de las variables si bien matemáticamente no es relevante, suele estar influenciada

por el contexto. Por ejemplo, si se quiere estudiar una función que expresa la temperatura en una
placa rectangular dada, lo mejor sería nominarla con la letra T y las coordendas que determinan
un punto arbitrario de la placa por (x; y), luego tenemos una función T = T (x; y). Si consideramos
también el tiempo, entonces podemos elegir la nominación para las variables como (x; y; t) y la
función sería T = T (x; y; t).

5 Derivadas parciales

Sea z = f (x; y) una función definida en D  R2.

11
Definición 20. La derivada parcial de f respecto a x en (a; b) se definen como

@f f (a + t; b) ¡ f (a; b)
fx(a; b) = (a; b) = lim ;
@x t!0 t
si este límite existe.
La derivada parcial de f respecto a y en (a; b) se definen como

@f f (a; b + t) ¡ f (a; b)
fy(a; b) = (a; b) = lim ;
@y t!0 t
si este límite existe.

Observación 21. Observe que en la definición, si consideramos la función parcial '(x) = f (x; b),
tenemos que

@f f (a + t; b) ¡ f (a; b) '(a + t) ¡ '(a) d'

(a; b) = lim = lim = ' 0(a) = (a):
@x t!0 t t!0 t dx

@f
Esto significa que para calcular la derivada parcial de f respecto a x, @x (a; b), se considera y
constante (y = b) y se deriva f normalmente como una función de una variable x, que se evalua
luego en el punto x = a.
@f
Análogamente, para calcular la derivada parcial de f respecto a y, @y (a; b), se considera x constante
(x = a) y se deriva f normalmente como una función de una variable y, que se evalua luego en el
punto y = b.

@z @z
Ejemplo 22. Calcule las derivadas parciais @x
e @y
de z = f (x; y) = x2y 3 sin(xy), en un punto
genérico (x; y).
Considerando y constante y derivando respecto a x tenemos que

@z
= (x2) 0y 3 sin(xy) +x2(y 3sin(xy)) 0 =2xy 3 sin(xy) +x2y 4cos(xy):
@x

Considerando ahora x constante y derivando respecto a y tenemos que

@z
= (y 3) 0x2 sin(xy) + y 3(x2 sin(xy)) 0 =3 y 2 x2 sin(xy) + y 3x3cos(xy):
@y

>>> var("x y z")

>>> f(x,y) = x^2*y^3*sin(x*y)
>>> f
>>> diff(f(x,y),x)
>>> diff(f(x,y),y)
>>> f(x,y).diff(x)
>>> f(x,y).diff(y)
>>> derivative(f(x,y),x)
>>> derivative(f(x,y),y)
>>> diff(f(x,y),x)(x=2,y=5)
>>> diff(f,y)
>>> diff(f,y)(2,3)

12
>>>

En general,

Definición 23. Sea f : U  Rn ¡! R y a = (a1; : : : ; an) 2 U. Para cada i = 1; : : : ; n, la derivada

parcial de f en a respecto a la i ¡ ésima coordenada se definen como

@f f (a + tei) ¡ f (a) f (a1; : : : ; ai + t; : : : an) ¡ f (a1; : : : ; ai; : : : ; an)

(a) = lim = lim
@xi t!0 t t!0 t

si este límite existe.

>>> f=function("f")(x,y)
>>> f
f (x; y)
>>> diff(f,x)
@
f (x; y)
@x
>>> diff(f,y)
@
f (x; y)
@y
>>>

Interpretación geométrica de la derivada parcial

Consideremos una función z = f (x; y), un punto (a; b) en el dominio de f y c = f (a; b). Geométri-
camente, en el espacio R3, el punto P = (a; b; c) se encuentra sobre la superficie z = f (x; y) que es
la gráfica de f .

 Cuando derivamos respecto de x, mantenemos la variable y fija (y = b) y consideramos la

función parcial '(x) = f (x; b).

 El gráfico de esta función de una variable x, coincide geométricamente con la proyección

sobre el plano x ¡ z, de la curva C1 que resulta de la intersección la superficie z = f (x; y)
con el plano y =b.

 La derivada ' 0(a) es la pendiente de la recta tangente a la gráfica de ' en el punto (a; c). En
el espacio, esta recta se encuentra en el plano x ¡ z y si la trasladamos paralelamente al eje
y hasta el punto P obtenemos una recta T1 tangente a la superficie z = f (x; y) en el punto P .

@f
 La derivada parcial @x (a; b) representa, entonces, la tangente del ángulo , que es el ángulo
que la recta tangente T1 forma con la recta z = 0; y = b paralela al eje x, es decir, es la
pendiente de la recta tangente al gráfico de la curva z = f (x; b).

 Análogamente, si C2 es la curva que resulta de la intersección de la superficie z = f (x; y) con

el plano x=a y T2 es la recta tangente a la curva C2 y a la gráfica de f en el punto P , la
@f
derivada parcial @y (a; b) representa, la tangente del ángulo , que es el ángulo que la recta
tangente T2 forma con la recta z = 0; x = a paralela al eje y, es decir, es la pendiente de la
recta tangente al gráfico de la curva z = f (a; y).

13
Derivadas direccionales

Sea f : U  R2 ¡! R; z = f (x; y), la derivada direccional de f en (a; b) en la dirección de v = (h; k)

se define como

@f f (a + th; b + tk) ¡ f (a; b)

(a; b) = lim
@v t!0 t

Observación 24. Observe que en la definición, si consideramos la función '(t) = f (x + th; b + tk),
tenemos que

@f f (a + th; b + tk) ¡ f (a; b) '(t) ¡ '(0) d'

(a; b) = lim = lim = ' 0(0) = (0):
@v t!0 t t!0 t dt

Esto significa que para calcular la derivada direccional de f en (a; b) en la dirección v = (h; k),
@f d
@v
(a; b), se calcula la derivada dt f (x + th; b + tk) normalmente como una función de una variable
t, que se evalua luego en el punto t = 0.

Definición 25. Sea f : U  Rn ¡! R, la derivada direccional de f en a = (a1; a2; : : : ; an) 2 U

en la dirección de v = (v1; v2; : : : ; vn) se define como

@f f (a + tv) ¡ f (a) f (a1 + tv1; : : : ; ai + tvi; : : : ; an + tvn) ¡ f (a1; : : : ; ai; : : : ; an)

(a) = lim = lim
@v t!0 t t!0 t

Derivadas de orden superior

@f @f
Como vimos en los ejemplos de la sección anterior, las derivadas parciales @x y @y de una función
z = f (x; y) son también funciones de dos variables. Así, podemos volver a considerar sus deri-
vadas parciales, llamadas derivadas parciales de segundo orden de z = f (x; y), definidas de la
siguiente manera:
       
@f @f @f @f
@ @x 2
@ f @ @y 2
@ f @ @y 2
@ f @ @x @ 2f
= ; = ; = ; =
@x @x2 @y @y 2 @x @x@y @y @y@x

Del mismo modo podemos tener derivadas parciales de 3er orden, 4to orden, . . .

Ejemplo 26. Determine las derivadas parciales de primer y segundo orden de

z = f (x; y)=x3 +2y 3 +3x2y 2

>>> var("x y z")

f(x,y)=x^3 +2*y^3 +3*x^2*y^2
>>> f(2,3)
170
>>> diff(f(x,y),x)
6 x y 2 + 3 x2
>>> diff(f(x,y),y)
6 x2 y + 6 y 2

14
>>> diff(f,x)
(x; y) 7! 6 x y 2 + 3 x2
>>> diff(f(x,y),x,x)
6 y2 + 6 x
>>> diff(f(x,y),x,y)
12 x y
>>> diff(f(x,y),y,x)
12 x y
>>> diff(f(x,y),y,y)
6 x2 + 12 y
>>>

Notemos que en este ejemplo

@ 2z @ 2z
= = 12xy
@x@y @y@x

Teorema 27. (Clairaut) Si z = f (x; y) tiene derivadas parciales mixtas continuas, entonces
son iguales, es decir,
@ 2z @ 2z
= :
@x@y @y@x

Plano tangente
En Cálculo I si una recta es tangente al gráfico de una función de una variable en un punto P ,
a medida que nos acercamos a este punto el gráfico de la función y la recta prácticamente se
confunden. Este concepto geométrico también se aplica en el caso de funciones de dos variables.
En este caso, el concepto de tangencia de un plano a una superficie que es gráfica de una función
de dos variables es el mismo: A medida que nos acercamos al punto de tangencia, el plano y la
superficie se confunden. También se puede pensar que, a medida que hacemos zoom en un punto
de la superficie, la superficie se vuelve cada vez más plana.
Un ejemplo de una superficie que no tiene un plano tangente en el origen es el cono z 2 =x2 + y 2.
Cualquier plano que pase por el origen es perfectamente distinguible del cono, por mucho que nos
acerquemos al origen.
Supongamos, entonces, que una superficie S, gráfica de una función z = f (x; y), tiene un plano
tangente en un punto Q = (a; b; c), donde c = f (a; b).
Ya vimos que las derivadas parciales de z = f (x; y) son las pendientes de las rectas tangentes rx y
ry a las curvas x y y que resultan de la intersección de la superficie S con planos verticales x=a
e y =b en las direcciones de los ejes x e y, respectivamente. Por lo tanto, estas dos rectas tangentes
definen un plano y, por supuesto, si S tiene un plano tangente, este plano contiene estas dos rectas.
Las curvas x y y tienen parametrizaciones naturales:

8 8
>
< x=x >
< x=a
x= y=b y= y=y
>
: z = f (x; b) >
: z = f (a; y)

15
luego
8 8
>
< x0 = 0
0
< x =1
0
> >
y =0
0 y =1
x(a) = y (b) =
0 0
>
: z 0 = @f (a; b) >
>
: z 0 = @f (a; b)
@y
@x

   
@f @f
Así, los vectores ux = 1; 0; @x (a; b) y uy = 0; 1; @y (a; b) dan las direcciones de las rectas rx y
ry, respectivamente.
El vector n, normal al plano tangente en el punto P , es

i j k  
@f @f @f @f @f
n = ux  uy = 1 0 @x
(a; b) = k ¡ (a; b) i ¡ (a; b) j = ¡ (a; b); ¡ (a; b); 1 ;
@f
@x @y @x @y
0 1 @y
(a; b)

donde i = (1; 0; 0); j = (0; 1; 0); k = (0; 0; 1) son los vectores básicos de R3.
Así, si P = (x; y; z) es un punto cualquiera del plano, la ecuación del plano tangente se obtiene por
el producto escalar del vector n y el vector del plano P ¡ Q = (x ¡ a; y ¡ b; z ¡ c),
 
@f @f
n
(P ¡Q) = ¡ (a; b); ¡ (a; b); 1
(x ¡ a; y ¡ b; z ¡ c) = 0
@x @x

esta ecuación puede escribirse como

@f @f
z = f (a; b) + (a; b) (x ¡ a) + (a; b) (y ¡ b)
@x @y

Este plano se denomina plano tangente a la gráfica de la función f en el punto Q = (a; b; c).
8
< xy
; (x; y) =
/ (0; 0)
Ejemplo 28. f (x; y) = : x2 + y2
0 (x; y) = (0; 0)

En los puntos (x; y) =

/ (0; 0) la función tiene derivadas parciales.
Veamos la función parcial '(x) = f (x; y), entonces

@f (xy) 0(x2 + y 2) ¡ xy(x2 + y 2) 0 y(x2 + y 2) ¡ 2x2y (y 2 ¡ x2)y

= ' 0(x) = = = 2
@x (x2 + y 2)2 (x2 + y 2)2 (x + y 2)2

Veamos la función parcial (x) = f (x; y), entonces

@f (xy) 0(x2 + y 2) ¡ xy(x2 + y 2) 0 x(x2 + y 2) ¡ 2xy 2 (x2 ¡ y 2)x

= 0
(y) = = = 2
@y (x2 + y 2)2 (x2 + y 2)2 (x + y 2)2

>>> diff(x*y/(x^2+y^2),x).factor()
>>> diff(x*y/(x^2+y^2),y).factor()
>>>
En el punto (0; 0), la función tiene derivadas parciales.
En este caso la función parcial es '(x) = f (x; 0) = 0, luego

@f '(0 + t) ¡ '(0)
(0; 0) = ' 0(0) = lim = 0:
@x t!0 t

16
De la misma forma, la función parcial es (y) = f (0; y) = 0, luego

@f (0 + t) ¡ (0)
(0; 0) = 0
(0) = lim = 0:
@y t!0 t

Veamos ahora lo que pasa en el origen (0; 0) y en la dirección v = (h; k), con h; k =
/ 0.

@f f (th; tk) ¡ f (0; 0) (th)(tk) hk

(0; 0) = lim = lim = 2
@v t!0 t t!0 (th)2 + (tk)2 h + k2

8
< x2 y
Ejemplo 29. g(x; y) = : x2 + y2 ; (x; y) =
/ (0; 0)
0 (x; y) = (0; 0)

En los puntos (x; y) = / (0; 0) la función tiene derivadas parciales (verificar como en el ejemplo
anterior). En el origen de nuevo las funciones parciales son '(x) = g(x; 0) = 0 y (y) = g(0; y) = 0,
@g @g
luego @x (0; 0) = @y (0; 0) = 0.

Veamos las derivadas direccionales. En este caso, en los puntos (x; y) = / (0; 0), en la dirección de
un vector v = (h; k), con h; k =
/ 0, consideramos la función parcial  (t) = g(x + th; y + tk), entonces

@f f (a + th; b + tk) ¡ f (a; b)  (t) ¡  (0)

(a; b) = lim = lim =  0(0):
@v t!0 t t!0 t

Calcularemos  0(0).
>>> g(x,y)=x^2*y/(x^2+y^2)
>>> g
>>> var("t h k")
>>> p(t)=g(x+t*h,y+t*k)
>>> p(t)
>>> diff(p(t),t)
>>> diff(p(t),t)(t=0).factor()
>>>
Veamos ahora lo que pasa en el origen (0; 0) y en la dirección v = (h; k), con h; k =
/ 0.

@f f (th; tk) ¡ f (0; 0) (th)2(tk)

(0; 0) = lim = lim =
@v t!0 t t!0 (th)2 + (tk)2

Diferenciabilidad
Una función f : U  R2 ¡! R, es diferenciable en (a; b) 2 U , si:

@f @f
1. Existen las derivadas parciales @x
(a; b), @y (a; b)

2. Para (h; k) 2 R2 tal que (a + h; b + k) 2 U podemos escribir

@f @f
f (a + h; b + k) = f (a; b) + (a; b) h + (a; b) k + r(h; k);
@x @y
donde
r(h; k)
lim = 0:
(h;k)!(0;0) k(h; k)k

17
Una función f : U  R3 ¡! R, es diferenciable en (a; b; c) 2 U , si:

@f @f @f
1. Existen las derivadas parciales @x
(a; b; c), @y (a; b; c), @z (a; b; c)

2. Para (h; k; l) 2 R2 tal que (a + h; b + k; c + l) 2 U podemos escribir

@f @f @f
f (a + h; b + k; c + l) = f (a; b; c) + (a; b; c) h + (a; b; c) k + (a; b; c) l + r(h; k; l);
@x @y @z

donde
r(h; k; l)
lim = 0:
(h;k;l)!(0;0;0) k(h; k; l)k
En general,

Definición 30. Una función f : U  Rn ¡! R, es diferenciable en a = (a1; : : : ; an) 2 U, si:

@f
1. Existen las derivadas parciales @xi
(a), para i = 1; : : : ; n

2. Para v = (v1; : : : ; vn) 2 Rn tal que a + v 2 U podemos escribir

@f @f
f (a + v) = f (a) + (a) v1 + : : : + (a) vn + r(v);
@x1 @xn
donde
r(v)
lim = 0:
v!0 kvk

Observación 31.

1. La segunda condición en la definición de diferenciabilidad, expresa el hecho que

@f @f
A(v1; : : : ; vn) = f (a) + (a) v1 + : : : + (a) vn
@x1 @xn

es la mejor aproximación lineal a f (a + v), en el sentido que el resto, que es la diferencia

r(v) = f (a + v) ¡ A(v);

r(v)
es un infinitesimo de orden superior a v, pues la condición lim v!0 kvk = 0, implica que
lim v!0r(v) = 0, pero también que r(v) tiende a cero, mucho más rápido que v.

2. La función A(v) es una funición afin 1. En el caso n = 2, la gráfica de esta función es un

plano que pasa por el punto (a; b; f (a; b)).

1. Una función A: Rn ! Rm es afin si tiene la forma A(v) = a + T (v), donde T : Rn ! Rm es una transformación
lineal. A su vez, T : Rn ! Rm es una transformación lineal si se cumple:

a. T (u + v) = T (u) + T (v), para u; v 2 Rn

b. T (u) = T (u), para u 2 Rn y  2 R.

En el caso m = 1, una transformación lineal T : Rn ! R tiene siempre la forma T (v1; : : : ; vn) = a1v1 +


+ anvn.

18
Teorema 32. Si f : U  Rn ! R es diferenciable en a 2 U, entonces f es continua en a.

Demostración. Recordemos que la continuidad se expresa por la condición

lim f (x) = f (a)

x!a

Esta condicionequivale a la siguiente, considerando el cambio de variable x $ a + v, entonces con

este cambio de variable, si x ! a, entonces v ! 0, luego

lim f (a + v) = f (a):
v!0

Veamos si esta condición se cumple. Como f es diferenciables sabemos que

  
@f @f @f
lim f (a + v) = lim f (a) + (a) v1 + : : : + (a) vn + r(v) = lim f (a) + lim (a) v1 + : : : +
v!0 v!0 @x1 @xn v!0 v!0 @x1
  
@f r(v) r(v)
(a) vn + lim r(v) = f (a) + 0 + lim kvk = f (a) + lim lim kvk = f (a) + 0
0 = f (a)
@xn v!0 v!0 kvk v!0 kvk v!0

(f  g) 0(x) = f 0(g(x))g 0(x)

19