UNIVERSIDAD APEC

MÉTODOS MATEMÁTICOS - DERIVADAS PARCIALES

”¡Cómo es posible que la matemática, un producto del pensamiento humano independiente de

la experiencia, se adapte tan admirablemente a los objetos de la realidad!

Albert Einstein (1879-1955)

Funciones de varias variables

La mayor parte del cálculo elemental trata funciones escalares de una variable, f (x) : R → R.
En cambio, muchas de las magnitudes que aparecen más frecuentemente en ciencia son
vectoriales, por lo que son necesarias funciones vectoriales para representarlas.

Por ejemplo en el movimiento de una particula en 3D, el vector posición del objeto tendrı́a
tres coordenadas: x, y, z, y la variable serı́a el tiempo. La función apropiada para representar
ese movimiento seria de la forma F (t) : R → R3 .
Las funciones escalares también pueden tener varias variables, por ejemplo, una función
escalar de dos variables podrı́a representarse como la elevación sobre un plano xy, es decir, a
cada valor del plano xy le correponde un valor de z = f (x, y), donde f (x, y) : R2 → R.
Las curvas de nivel son una representación visual de todos los puntos para los que la función
toma un mismo valor. Son muy útiles por ejemplo en meteorologı́a, para representar las curvas
con la misma presión, o lı́neas isobaras.

Javier Garcı́a Maimó, PhD

2 Funciones multivariables y derivadas parciales

Los campos de vectores son útiles en problemas de mecánica de fluidos o de electricidad

y magnetismo por ejemplo. Un campo de vectores es una región en la que en cada posición se
tiene un vector, es decir, si tuvieramos un campo sobre un plano xy, en cada posición del plano
habrı́a un vector.

Derivadas parciales
Las derivadas parciales de una función multivariable son las derivadas respecto cada una de las
variables tomando las demás como constantes. Su interpretación geométrica es la pendiente de
la función al deplazarse desde un punto en la dirección de la variable respecto la que se deriva.
Por ejemplo, para la función f (x, y) = x2 y

∂f
= 2xy
∂x
∂f
= x2
∂y

Javier Garcı́a Maimó, PhD

3 Funciones multivariables y derivadas parciales

En la figura el eje en rojo representa el eje x y el verde el eje y. Evaluando las derivadas parciales
en el punto (1,-2) se obtiene
∂f
(1, −2) = 2xy|(1,−2) = −2
∂x
∂f
(1, −2) = x2 |(1,−2) = 1
∂y

Gradiente y Derivada direccional

El gradiente es la generalización de la derivada para funciones que dependen de más de una
variable, el gradiente ∇f es un vector cuyo módulo representa el ritmo de variación de la función,
o derivada direccional y que apunta en la dirección de máximo cambio de la función.
El gradiente para funciones de dos variables se define
 
∂f ∂f ∂f ∂f
∇f = i+ j= ,
∂x ∂y ∂x ∂y

Para la función f (x, y) = x2 y el vector gradiente es ∇f (x, y) = (2xy, x2 ). Una particularidad

del gradiente es que siempre es perpendicular a las curvas de nivel.

El gradiente indica la dirección de máximo cambio, si se quiere saber el cambio en cualquier

dirección, no necesariamente en la dirección de los ejes ni del gradiente, sino según la dirección
de un vector cualquiera ~v , es necesario definir la derivada direccional.

∇~v f = ~v · ∇f

Divergencia y rotacional
Dos operadores que se pueden aplicar a campos vectoriales son la divergencia ∇·F y el rotacional
∇×F . Los resultados de aplicar los operadores son magnitudes escalare que aportan información
valiosa sobre las caracterı́sticas del campo vectorial. Se puede ver una representación visual de
la divergencia y el rotacional en el siguiente video.

La regla de la cadena. Derivada implı́cita

La regla de la cadena es útil para derivar funciones compuestas:
df dg
si h(x) = (f ◦g)(x) = f (g(x)) entonces h0 (x) = . Por ejemplo, si h(x) = sin2 x entonces
dg dx
h0 (x) = 2 sin x cos x.
La derivada implı́cita es útil en expresiones en las que no es sencillo despejar la variable y.
Aunque este no es el caso, vamos a usar por simplicidad el ejemplo x2 + y 2 = 1

Javier Garcı́a Maimó, PhD

4 Funciones multivariables y derivadas parciales

d
x2 + y 2 = 0


dx
d(y 2 )
2x + =0
dx
para hacer la derivada tenemos que tener en cuenta que y = y(x) es una función de x, por lo
que debemos aplicar la regla de la cadena.

2x + 2yy 0 = 0
x
y0 = −
y
Cuando una función depende de x y de y hay que considerar que es necesario derivar también
respecto x la función y usando la regla de la cadena, es decir

df ∂f ∂f ∂y
= +
dx ∂x ∂y ∂x

2 cos x + y 2
 
Ejemplo: Calcular y0 para la expresión 2 sin x + xy 2 = 2, sol: −
2xy

Actividades

1. Definir las caracterı́sticas (dominio y rango) de las funciones que permitan modelar lo
siguiente:

a) la presión en cada punto de la atmósfera terrestre.

b) la trayectoria de la partı́culas de un contaminante en un fluido.
c) el campo eléctrico de una distribución de cargas situada sobre un plano.
d) el potencial eléctrico de la distribución de cargas anterior.

2. ¿Qué forma tienen las curvas de nivel de la función f (x, y) = x2 y?

3. Calcular el gradiente de la función f (x, y) = x2 y + sin y, y evaluarlo en el punto (0,π)

4. Calcular la derivada direccional de la función f (x, y) = x2 y + sin y, en la dirección (1,2)

5. Calcular la divergencia y el rotacional del campo vectorial F (x, y, z) = (2xy, xey , 2z)
1
6. Calcular la derivada de + ln y + (xy)3 = 4x
x

Javier Garcı́a Maimó, PhD

5 Funciones multivariables y derivadas parciales

Aplicaciones. Cambios de variables en integrales

Integrales simples
Z 2
Este caso ya es conocido y vamos a ilustrarlo con un ejemplo: e2x+1 dx.
0
Parece natural intentar el cambio t = 2x + 1, de forma que dt = 2dx, y entonces
Z 2
1 5 t
Z
2x+1 1  t 5 1 5
e dx = e dt = e 1 = (e − e)
0 2 1 2 2
En general este tipo de cambios son de la forma
Z xb Z tb
dx
f (x)dx = fˆ(t) dt
xa ta dt

Caso multivariable. Integral doble

ZZ
En el caso multivariable se tiene f (x, y)dA. Donde la expresión de dA depende de las
coordenadas usadas, para coordenadas cartesianas dA = dxdy. En general si u = u(x, y) y
v = v(x, y) nos gustarı́a saber qué expresión se usa para transformar la integral en este caso
ZZ ZZ
f (x, y)dxdy = fˆ(u, v)???dudv

Por simplicidad estudiaremos el caso de la integral doble, y un ejemplo de un cambio de

coordenadas de cartesianas a polares1 . En este caso tenemos (x, y) → (r, θ)

x = r cos θ
y = r sin θ

y por tanto
∂x ∂x
dx = dr + dθ
∂r ∂θ
∂y ∂y
dy = dr + dθ
∂r ∂θ

   
dx cos θ −r sin θ dr
De forma matricial = . En general si u = u(x, y) y
dy sin θ r cos θ dθ
 
  ∂x ∂x  
dx  ∂u ∂v  du
v = v(x, y) tendremos =  ∂y ∂y  . La matriz del cambio de variables
dy dv
  ∂u ∂v
∂x ∂x
∂(x, y)  ∂u ∂v 
=  ∂y ∂y  se llama matriz Jacobiana.
∂(r, θ)
∂u ∂v

Cambios de variable en general

En general la expresión que permite hacer el cambio de variables es
ZZ ZZ  
 ∂(x, y) 
f (x, y)dxdy = ˆ
f (u, v) 
  dudv
∂(u, v) 
1
en coordenadas polares dA = rdrdθ

Javier Garcı́a Maimó, PhD

6 Funciones multivariables y derivadas parciales

ZZ
Ejemplo: Calcular xydxdy, donde R es la región limitada por 0 < x < 2 y x < y < x + 4,
R

a) usando las variables x, y

b) haciendo el cambio x = 2u, y = 2u + 4v, sol:80/3

c) Dibujar la región R con cada pareja de variables

Actividades
Z
7. Utilizar un cambio de variable para resolver la integral indefinida (x + 3)(x − 1)5 dx

x2 y
 
8. Dada la función F : → R2 F =R2 , Calcular su jacobiano y evaluarlo en
5x + sin y
el punto (1, π). Calcular su determinante.
ZZ
x+y
9. Calcular e x−y dxdy, donde R es la región limitada por x > 0, y > 0 y x + y = 2
R
ZZ
10. Calcular (x + y)2 ex−y dxdy, donde R es la región limitada por x + y = 1, x + y = 4,
R
x − y = −1 y x − y = 1

Javier Garcı́a Maimó, PhD

7 Funciones multivariables y derivadas parciales

Aplicaciones: Resolución de ecuaciones no lineales

El problema de la resolución de ecuaciones no lineales consiste en encontrar la solución de
f (x) = 0, donde f (x) es una función real de variable real 2 . Incluso en ecuaciones de una
sola incógnita sencillas, no siempre es posible despejar la variable, por ejemplo en el caso de la
ecuación ex − 5x = 0.

El método gráfico es la forma más sencilla de

obtener una solución aproximada de una ecuación,
consiste en graficar la función f (x) y determinar la
intersección de la función con el eje x. Este punto
representa el valor de x que anula la función, es
decir, si denotamos con ξ este valor de x, entonces
sustituyéndolo en la función se obtiene f (ξ) = 0.

Los métodos iterativos producen una sucesión

convergente a la solución del problema. El método gráfico suele usarse como primer paso
de un método iterativo, como el método de Newton-Raphson

Método de Newton-Raphson
El método de Newton es una aproximación de primer
orden a la solución de la ecuación no lineal, su
interpretación geométrica es trazar la tangente a
y = f (x) en el punto (xn , f (xn )) para obtener la
recta
y = f (xn ) + f 0 (xn )(x − xn )
que corta al eje x en

f (xn )
x = xn −
f 0 (xn )

que es precisamente el valor de xn+1 de la fórmula

de Newton-Raphson, también conocido como método
de la tangente.

f (xn )
xn+1 = xn −
f 0 (xn )

Ejemplo: Encontrar la raı́z positiva real de la ecuación no lineal 3x4 − 4x3 + 9x2 + 3x − 1 = 0
tomando como punto inicial x0 = 1

n x f (x) f 0 (x) εr
0 1.0000 10.0000 21.0000 -
1 0.5238 2.6918 10.8607 90.91
2 0.2760 0.4466 7.30567 89.81
3 0.2148 0.0266 6.43208 28.46
4 0.2107 0.0001 6.37206 1.96
5 0.2107 0.0000 6.37178 0.00

Se puede ver un video con la interpretación geométrica del método de Newton en el siguiente
enlace. Una hoja dinámica para practicar se puede encontrar aqui.
2
la solución de una ecuación también se llama raı́z o cero de una función

Javier Garcı́a Maimó, PhD

8 Funciones multivariables y derivadas parciales

Actividades
11. Utilizar tres iteraciones del método de de Newton para encontrar la raı́z que se encuentra
en el intervalo [2,3] de x2 − ln x − 4 = 0
12. El balance de masa de un contaminante en un lago se expresa:
dc √
V = W –Qc–kV c
dt
tomando V = 1 × 106 , Q = 1 × 105 y W = 1 × 106 , y k = 0.25, use tres iteraciones del
método de Newton para calcular la concentración en el estado estacionario tomando como
punto inicial c = 4.

El método de Newton-Raphson multidimensional

El método de Newton-Raphson para sistemas consiste en una extensión del método para
ecuaciones al caso multivariable. Por ejemplo, sean las funciones
x2 − 2x − y = −0.5
x2 + 4y 2 = 4
Para encontrar la solución es necesario resolver el sistema de dos ecuaciones no lineales con dos
incógnitas que definen dos curvas en el plano XOY : una parábola y una elipse. Los puntos de
intersección de las dos curvas son las soluciones del sistema.

Al tener una derivada de una función de varias variables la derivada de la expresión iterativa
unidimensional se transforma en la matriz JacobianaLa, de manera que fórmula iterativa del
método de Newton Raphson para sistemas es:
x(k+1) = x(k) − J−1 (x(k) )F(x(k) ),
donde el vector x es el vector cuyas coordenadas son las variables buscadas (x1 = x, x2 = y
en este caso), J representa la matriz Jacobiana asociada al sistema y F es el vector cuyas
coordenadas son las funciones fi del sistema.
 
f1 (x)
 f2 (x) 
F= 
 .... 
fn (x)
La matriz Jacobiana asociada al sistema es
 ∂f1 (x) ∂f1 (x) ∂f1 (x) 
∂x1 ∂x2 ... ∂xn
∂f2 (x) ∂f2 (x) ∂f2 (x)
...
 
J=
 ∂x1 ∂x2 ∂xn 

 ... ... ... ... 
∂fn (x) ∂fn (x) ∂fn (x)
∂x1 ∂x2 ... ∂xn

Javier Garcı́a Maimó, PhD

9 Funciones multivariables y derivadas parciales

Formamos el vector F y la matriz Jacobiana J para nuestro caso:

x2 − 2x − y + 0.5
   
2x − 2 −1
F= , J=
x2 + 4y 2 − 4 2x 8y

De la gráfica vemos que podemos tomar como estimación inicial P0 el punto (2.0000, 0.2500)
Calculamos F0 y J0 :
     
0.25 2 −1 0.25 0.125
F0 = , J0 = , J−1
0 =
0.25 4 2 −0.5 0.25

Esto permite calcular P1 :

        
x1 2.0000 0.25 0.125 0.25 1.9063
P1 = = − =
y1 0.2500 −0.5 0.25 0.25 0.3125

Hacemos una iteración más. Calculamos F1 , J1 y J−1

1 :
     
0.0088 1.8125 −1 −1 0.2996 0.1199
F1 = , J0 = , J0 =
0.0244 3.8125 2.5 −0.4569 0.2172

De forma que:
      
1.9063 0.2996 0.1199 0.0088 1.9007
P2 = − =
0.3125 −0.4569 0.2172 0.0244 0.3112

Del mismo modo se puede continuar iterando para obtener aproximaciones más precisas. El
número de iteraciones dependerá de la tolerancia requerida en el problema. En este caso después
de la segunda iteración se tiene εx = 0.29% y εy = 0.41% de error. Con una iteración más se
obtiene P3 = (1.9007, 0.3113) y errores εx = 0.00% y εy = 0.002%.

Actividades

13. Determine las raı́ces de las siguientes ecuaciones no lineales simultáneas, por medio de dos
iteracione del método de Newton-Raphson:

y = −x2 + x + 0.75
y + 5xy = x2

utilizando valores iniciales x = 1.2, y = 1.2

Javier Garcı́a Maimó, PhD