Cálculo de Campos Vectoriales

Calculus. Early trascendentals. Seventh Edition.

James Stewart.
Presentación con diapositivas Beamer
elaboradas por
Luis Angel Zaldívar Cruz

Departamento de Ciencias Básicas

Instituto Tecnológico de Tehuacán

30 de mayo de 2020
Competencias del tema

En este tema el estudiante. . .

I Aprenderá a determinar las integrales de línea.
I Utilizará las integrales de línea para calcular el trabajo realizado por un
campo de fuerzas al mover un cuerpo a lo largo de una curva.
I Aprenderá a resolver integrales de superficie (que pueden usarse para
hallar la rapidez de un fluido por una superficie).
I Comprenderá que la conexión entre estos nuevos tipos de integrales
simples, dobles y triples que ya hemos visto están dadas por las versiones
de dimensiones más altas del teorema fundamental del cálculo: el teorema
de Green, el teorema de Stokes y el teorema de la divergencia.
I Utilizará las TIC’s para resolver integrales de línea.
I Utilizará las TIC´s para estudiar el Teorema de Green, el Teorema de
Stokes y el Teorema de la divergencia.
Campos vectoriales
Las flechas de la Figura 1 son vectores velocidad que indican la rapidez y dirección del
viento en los puntos que están 10 m por arriba de la superficie en el área de la bahía
de San Francisco. A primera vista, se observa que las flechas más largas en el inciso a)
indican que la mayor rapidez del viento en este tiempo ocurrió cuando todos los
vientos atravesaron la bahía por el Golden Gate Bridge. El inciso b) muestra los muy
diferentes patrones de viento 12 horas antes. Imagine un vector velocidad del viento
asociado con cada punto en el aire. Éste es un ejemplo de un campo vectorial de
velocidad.

Figura 1: Campos vectoriales de velocidad que muestran los patrones de viento en la

bahía de San Francisco.
Campos vectoriales
I Otros ejemplos de campos vectoriales de velocidad se ilustran
en la Figura 2: corrientes oceánicas y el flujo que se encuentra
en un automóvil.

Figura 2: Campos vectoriales de velocidad.

I Otro tipo de campo vectorial, llamado campo de fuerza, asocia

un vector fuerza con cada punto de una región. Un ejemplo es
el campo de fuerza gravitacional.
Campos vectoriales
En general, un campo vectorial es una función cuyo dominio es un
conjunto de puntos en R2 (o R3 ) y cuyo rango es un conjunto de
vectores en V2 o (V3 ).
Definición 1
Sea D un conjunto en R2 (una región plana). Un campo vectorial
sobre R2 es una función F que asigna a cada punto (x, y) en D un
vector bidimensional F(x, y).
I La mejor manera de representar un campo vectorial es dibujar
la flecha que representa al vector F(x, y) que inicie en el punto
(x, y).
I Naturalmente, es imposible hacerlo para todos los puntos
(x, y), pero podemos conseguir una representación razonable
de F trazando la flecha para algunos puntos representativos en
D como en la Figura 3.
Campos vectoriales

Figura 3: Campo vectorial en R2 .

Puesto que F(x, y) es un vector bidimensional, podemos expresarlo

en términos de sus funciones componentes P y Q como sigue:
F(x, y) = P (x, y)i + Q(x, y)j = (P (x, y), Q(x, y))
o bien, simplificando,
F = P i + Qj.
Observe que P y Q son funciones escalares de dos variables y,
algunas veces, se les llama campos escalares para distinguirlos de
los campos vectoriales.
Campos vectoriales
Definición 2
Sea E un subconjunto de R3 . Un campo vectorial sobre R3 es
una función F que asigna a cada punto (x, y, z) en E un vector
tridimensional F(x, y, z).
Un campo vectorial F sobre R3 se representa en la Figura 4.
Podemos expresarlo en términos de sus funciones constituyentes P ,
Q y R como

F(x, y, z) = P (x, y, z)i + Q(x, y, z)j + R(x, y, z)k.

Figura 4: Campo vectorial sobre R3 .

Campos vectoriales
I Se puede demostrar que F es continua si y sólo si sus
funciones constituyentes P , Q y R son continuas.
I Algunas veces identificamos un punto (x, y, z) con su vector
de posición x = (x, y, z) y escribimos F(x) en lugar de
F(x, y, z). Entonces F se convierte en una función que asigna
un vector F(x) a un vector x.
Ejemplo 1
Un campo vectorial sobre R2 está definido por F (x, y) = −yi + xj.
Describa F trazando algunos de sus vectores F(x, y) como en la
Figura 3.
Solución
Puesto que F(1, 0) = j, dibujamos el vector j = (0, 1) con inicio en
el punto (1, 0) en la Figura 5. Como F(0, 1) = −i, dibujamos el
vector (−1, 0) con inicio en el punto (0, 1). Al continuar de este
modo, calculamos varios valores representativos de F(x, y) en la
Tabla 1 y dibujamos los vectores correspondientes para representar
el campo vectorial en la Figura 5.
Campos vectoriales

(x, y) F(x, y) (x, y) F(x, y)

(1, 0) (0, 1) (−1, 0) (0, −1)
(2, 2) (−2, 2) (−2, −2) (2, −2)
(3, 0) (0, 3) (−3, 0) (0, −3)
(0, 1) (−1, 0) (0, −1) (1, 0)
(−2, 2) (−2, −2) (2, −2) (2, 2)
(0, 3) (−3, 0) (0, −3) (3, 0)
Cuadro 1: Valores representativos de F(x, y).
Campos vectoriales

Figura 5: F(x, y) = −yi + xj.

Al parecer, según la Figura 5, cada flecha es tangente a la

circunferencia con centro en el origen. Para confirmarlo, calculemos
el producto punto del vector de posición x = xi + yj con el vector
F(x) = F(x, y):

x · F(x) = (xi + yj) · (−yi + xj) = −xy + yx = 0.

Campos vectoriales
Esto demuestra que F(x, y) es perpendicular al vector de posición (x, y) y, por
tanto,p
es tangente a la circunferencia con centro en el origen y radio
|x| = x2 + y 2 . Observe que también
q p
|F(x, y)| = (−y)2 + x2 = x2 + y 2 = |x|

de modo que la magnitud del vector F(x, y) es igual al radio de la

circunferencia.
I Algunos sistemas algebraicos computarizados son capaces de dibujar
campos vectoriales en dos o tres dimensiones. Proporcionan una mejor
representación del campo vectorial de lo que es posible a mano, porque la
computadora puede trazar una gran cantidad de vectores representativos.
La Figura 6 muestra una gráfica por computadora del campo vectorial del
Ejemplo 1. El código en SageMath es
var(’x,y’)
plot_vector_field((-y,x), (x,-5,5), (y,-5,5),color=’blue’)
I Las Figuras 7 y 8 muestran otros dos campos vectoriales. Observe que las
computadoras dan una escala a las longitudes de los vectores de modo
que no sean demasiado grandes, pero que sean proporcionales a sus
longitudes verdaderas.
Campos vectoriales

Figura 6: Campo vectorial de F(x, y) = (−y, x).

Campos vectoriales
El código en SageMath que dibuja el campo vectorial de F(x, y) = (y, sen x) es
(ver Figura 7)
var(’x,y’)
plot_vector_field((-y,sin(x)), (x,-6,6), (y,-6,6),color=’blue’)

Figura 7: Campo vectorial de F(x, y) = (y, sen x).

Campos vectoriales
El código en SageMath que dibuja el campo vectorial de
F(x, y) = ln 1 + y 2 , ln 1 + x2 es (ver Figura 8)

var(’x,y’)
plot_vector_field((ln(1+y^2),ln(1+x^2)), (x,-5,5),
(y,-5,5),color=’blue’)

Figura 8: Campo vectorial de F(x, y) = ln 1 + y 2 , ln 1 + x2 .

Campos vectoriales

I La mayoría de los campos vectoriales tridimensionales son

virtualmente imposibles de dibujar a mano, por lo que se
necesita recurrir a un sistema algebraico de cómputo. Se
ilustran ejemplos en las figuras 9, 10 y 11 realizados con
SageMath.
I Observe que los campos vectoriales de las figuras 9 y 10 tienen
fórmulas similares, pero todos los vectores de la figura 10
apuntan en la dirección general del eje y negativo porque sus
componentes y son −2.
I Si el campo vectorial en la figura 11 representa un campo de
velocidad, entonces una partícula podría ser desplazada hacia
arriba y giraría en espiral alrededor del eje z en el sentido de
las manecillas del reloj si se ve desde arriba.
Campos vectoriales
El código en SageMath que dibuja el campo vectorial de F(x, y) = (y, z, x) es (ver Figura 9)
var(’x,y,z’)
g = plot_vector_field3d((y,z,x), (x,-2,2), (y,-2,2),(z,-2,2),colors=’blue’)
from sage.manifolds.utilities import set_axes_labels
graf = set_axes_labels(g,’x’,’y’,’z’,color=’red’)
graf

Figura 9: Campo vectorial de F(x, y) = (y, z, x).

Campos vectoriales
El código en SageMath que dibuja el campo vectorial de F(x, y) = (y, −2, x) es (ver Figura 10)
var(’x,y,z’)
g = plot_vector_field3d((y,-2,x), (x,-2,2), (y,-2,2),(z,-2,2),colors=’blue’)
from sage.manifolds.utilities import set_axes_labels
graf = set_axes_labels(g,’x’,’y’,’z’,color=’red’)
graf

Figura 10: Campo vectorial de F(x, y) = (y, −2, x).

Campos vectoriales
El código en SageMath que dibuja el campo vectorial de F(x, y) = (y/z, −x/z, z/4) es (ver Figura 11)
var(’x,y,z’)
g = plot_vector_field3d((y/z,-x/z,z/4), (x,-2,2), (y,-2,2),(z,1,5),colors=’blue’)
from sage.manifolds.utilities import set_axes_labels
graf = set_axes_labels(g,’x’,’y’,’z’,color=’red’)
graf

Figura 11: Campo vectorial de F(x, y, z) = (y/z, −x/z, z/4).

Campos vectoriales

Ejemplo 2
Imagine un fluido que corre en forma estable por una tubería, y sea ∇ (x, y, z)
el vector velocidad en un punto (x, y, z). Entonces ∇ (x, y, z) asigna un vector
a cada punto (x, y, z) en un determinado dominio E (el interior de la tubería),
de modo que ∇ (x, y, z) es un campo vectorial sobre R3 llamado campo de
velocidades. Un campo de velocidades posible se ilustra en la Figura 12. La
rapidez en cualquier punto dado se indica por la longitud de la flecha.

Figura 12: Campo de velocidades.

Campos vectoriales
Los campos de velocidades también se presentan en otras áreas de la física. Por
ejemplo, el campo vectorial del Ejemplo 1 se podría usar como campo de
velocidades para describir la rotación de una rueda en el sentido contrario al de
las manecillas del reloj. En las Figuras 1 y 2 se ven otros ejemplos de campos
de velocidad.

Ejemplo 3
La ley de la gravitación de Newton establece que la magnitud de la fuerza
gravitacional entre dos objetos con masas m y M es
mM G
|F| =
r2
donde r es la distancia entre los objetos y G es la constante gravitacional.
(Éste es un ejemplo de la ley de los cuadrados inversos.) Supongamos que el
objeto de masa M está en el origen en R3 . (Por ejemplo, M podría ser la masa
de la Tierra y el origen podría ser su centro.) Sea x = (x, y, z) el vector de
posición del objeto con masa m. Entonces, r = |x|, así que r2 = |x|2 . La
fuerza gravitacional ejercida en este segundo objeto actúa hacia el origen, y el
vector unitario en esta dirección es
x
− .
|x|
Campos vectoriales
Por lo tanto, la fuerza gravitacional que actúa sobre el objeto en x = (x, y, z) es
mM G
F(x) = − x. (1)
|x|3
[Los físicos utilizan la notación r en lugar de x para el vector de posición, de
modo que podemos encontrar la fórmula 1 escrita en la forma
mM G
F(r) = − r
r3
donde r = |r|. La función dada por la ecuación 1 es un ejemplo de un campo
vectorial, llamado campo gravitacional, porque asocia un vector [la fuerza
F(x)] con todo punto x en el espacio.
La fórmula 1 es una forma compacta de expresar el campo gravitacional, pero
también podemos escribirla en términos dep sus funciones constituyentes usando
el hecho de que x = xi + yj + zk y |x| = x2 + y 2 + z 2 :
−mM Gx −mM Gy −mM Gz
F(x, y, z) = i+ j+ k.
(x2 + y2 + z 2 )3/2 (x2 + y2 + z 2 )3/2 (x2 + y 2 + z 2 )3/2
El campo gravitacional F se representa en la Figura 13.
Campos vectoriales

Figura 13: Campo de fuerza garvitacional.

Ejemplo 4
Suponga que una carga eléctrica Q se localiza en el origen. De acuerdo con la ley de Coulomb, la fuerza
eléctrica F(x) que ejerce esta carga sobre la carga q situada en el punto (x, y, z) con vector de posición
x = (x, y, z) es
εqQ
F(x) = x (2)
|x|3

donde ε es una constante (que depende de las unidades que se utilizan). En el caso de cargas similares,
qQ > 0 y la fuerza es de repulsión; si las cargas son de signo contrario, entonces qQ < 0 y la fuerza es
de atracción. Observe la similitud entre las fórmulas 1 y 2. Ambos campos vectoriales son ejemplos de
campos de fuerza.
Campos vectoriales
En lugar de considerar la fuerza eléctrica F, los físicos toman en cuenta a
menudo la fuerza por unidad de carga:
1 εQ
E(x) = F(x) = x.
q |x|3

Entonces E es un campo vectorial sobre R3 , llamado campo eléctrico de Q.

Campos gradiente
Si f es una función escalar de dos variables, de acuerdo con el Tema 4 su
gradiente ∇f (o grad f ), se define como
∂f (x, y) ∂f (x, y)
∇f (x, y) = i+ j.
∂x ∂y

Por tanto, ∇f es realmente un campo vectorial sobre R2 y se llama campo

vectorial gradiente. Del mismo modo, si f es una función escalar de tres
variables, su gradiente es un campo vectorial sobre R3 dado por
∂f (x, y, z) ∂f (x, y, z) ∂f (x, y, z)
∇f (x, y, z) = i+ j+ k.
∂x ∂y ∂z
Campos vectoriales
Ejemplo 5
Encuentre el campo vectorial gradiente de f (x, y) = x2 y − y 3 . Dibuje el campo vectorial gradiente
junto con un mapa de contorno de f . ¿Cuál es su relación?

Solución
El campo vectorial gradiente está dado por

∂f ∂f 
2 2

∇f (x, y) = i+ j = 2xyi + x − 3y j.
∂x ∂y

En la figura 15 se muestra un mapa de contorno de f con el campo vectorial gradiente. Observe que los
vectores gradiente son perpendiculares a las curvas de nivel, como era de esperarse.

Figura 14: Relación entre el campo vectorial gradiente y el mapa de contorno de f .

Campos vectoriales
Note también que los vectores gradiente son largos donde las curvas de nivel
están cercanas entre sí, y cortos donde las curvas se separan. La razón es que la
longitud del vector gradiente es el valor de la derivada direccional de f y las
curvas de nivel cercanas indican una gráfica de fuerte pendiente.
El código en SageMath de la Figura 14 es:
var(’x,y’)
f=x^2*y-y^3
contour_plot(f,(x,-4,4),(y,-4,4),fill=False)
+plot_vector_field(f.gradient(), (x,-4,4), (y,-4,4),color=’blue’)
Un campo vectorial F se denomina campo vectorial conservativo si es el
gradiente de alguna función escalar, es decir, si existe una función f tal que
F = ∇f . En esta situación, f recibe el nombre de función potencial para F.
No todos los campos vectoriales son conservativos, pero tales campos surgen
con frecuencia en la física. Por ejemplo, el campo gravitacional F del Ejemplo 4
es conservativo porque si definimos
mM G
f (x, y, z) = p
x2 + y 2 + z 2
Campos vectoriales

entonces
∂f ∂f ∂f
∇f (x, y, z) = i+ j+ k
∂x ∂y ∂z
−mM Gx −mM Gy −mM Gz
= i+ j+ k
(x2 + y 2 + z 2 )3/2 (x2 + y 2 + z 2 )3/2 (x2 + y 2 + z 2 )3/2
= F(x, y, z).

En secciones posteriores aprenderemos la manera de determinar si un campo

vectorial dado es conservativo o no lo es.
Integrales de línea
En esta sección se define una integral que es similar a la integral simple, excepto que en lugar de
integrar sobre un intervalo [a, b], integramos sobre una curva C. Estas integrales se llaman integrales de
línea, aunque un mejor nombre es el de “integrales curvilíneas”. Fueron inventadas a principios del siglo
XIX para resolver problemas relacionados con el flujo de fluidos, fuerzas, electricidad y magnetismo.
I Iniciamos con una curva plana C dada por las ecuaciones paramétricas

x = x(t), y = y(t) (3)

o, en forma equivalente, por la ecuación vectorial r(t) = x(t)i + y(t)j, y supongamos que C es
una curva suave. Esto significa que r0 es continua y que r0 (t) 6= 0.
I  
Si dividimos el intervalo del parámetro [a, b] en n subintervalos ti−1 , ti de igual ancho y
hacemos xi = x (ti ) y yi = y (ti ), entonces los puntos correspondientes Pi (xi , yi ) dividen a
C en n subarcos de longitudes ∆s1 , ∆s2 , . . . , ∆sn (ver la Figura 15).
I Elegimos cualquier punto Pi∗ x∗
 
∗

i , yi en el i-ésimo subarco. (Esto corresponde a un punto ti
∗

en ti−1 , ti ).

I Ahora, si f es
una función de dos variables cuyo dominio incluye a la curva C, evaluamos f en el
punto x∗ ∗
i , yi , multiplicamos por la longitud ∆si del subarco, y formamos la suma

n
∗ ∗

X
f xi , yi ∆si
i=1

que es similar a la suma de Riemann.

I Luego tomamos el límite de estas sumas y establecemos la siguiente definición por analogía con
la integral simple.
Integrales de línea

Figura 15: División de C en n subarcos.

Definición 3
Si f se define sobre una curva C suave dada por las ecuaciones 3, entonces la
integral de línea de f a lo largo de C es
ˆ n
X
f (x, y)ds = lı́m f (x∗i , yi∗ ) ∆si
C n→∞
i=1

si este límite existe.

Integrales de línea
Conocemos que la longitud de C es
ˆ b  2  2
s
dx dy
L= + dt.
a dt dt

Un razonamiento similar se puede plantear para demostrar que si f es una

función continua, entonces el límite de la Definición 3 siempre existe y la
fórmula siguiente se puede usar para evaluar la integral de línea:
ˆ ˆ
s 
2  2
dx dy
f (x, y)ds = f (x(t), y(t)) + dt. (4)
C C dt dt

El valor de la integral de línea no depende de la parametrización de la curva,

siempre que ésta se recorra exactamente una vez cuando t se incrementa desde
a hasta b.
Si s(t) es la longitud de C entre r(a) y r(t), entonces,
s 
2  2
ds dx dy
= +
dt dt dt
Integrales de línea

La manera de recordar la fórmula 4 es expresar todo en términos del parámetro

t: usamos las ecuaciones paramétricas para expresar x y y en términos de t y
escribimos ds como s 
2  2
dx dy
ds = + dt.
dt dt
En el caso especial donde C es el segmento rectilíneo que une (a, 0) con (b, 0),
al usar x como parámetro, podemos escribir las ecuaciones paramétricas de C
como sigue: x = x, y = 0, a ≤ x ≤ b. La fórmula 4 se transforma en
ˆ ˆ b
f (x, y)ds = f (x, 0)dx
C a

y en este caso la integral de línea se reduce a una integral simple ordinaria.

De la misma manera que para una integral simple ordinaria, la integral de línea
de una función
´ positiva se interpreta como un área. De hecho, si f (x, y) ≥ 0,
entonces C f (x, y)ds representa el área de un lado de la “cerca” o de la
“cortina” de la Figura 16, cuya base es C y altura por arriba del punto (x, y) es
f (x, y).
Integrales de línea

Figura 16: Interpretación de una integral de línea.

Ejemplo
´
6
Evalúe C 2 + x2 y ds, donde C es la mitad superior de la circunferencia

x2 + y 2 = 1.

Solución
Con objeto de aplicar la fórmula 4 necesitamos primero las ecuaciones
paramétricas que representen a C. Recuerde que la circunferencia unitaria se
puede parametrizar por medio de las ecuaciones

x = cos t, y = sen t
Integrales de línea
y la mitad superior de la circunferencia se describe por el intervalo del
parámetro 0 ≤ t ≤ π. (Véase la Figura 17).

Figura 17: Mitad superior de la circunferencia.

Por tanto, la fórmula 4 da

ˆ ˆ
s 2 2
π 
2
2
dx dy
2 + x y ds = 2 + cos t sen t + dt
C 0 dt dt
ˆ π
p
= 2 + cos2 t sen t
sen2 t + cos2 t dt
0
ˆ π π
cos3 t

2 + cos2 t sen t dt = 2t −


=
0 3 0
2
= 2π + ≈ 6.94985.
3
Integrales de línea

El código en SageMath para resolver esta integral de línea es:

t = var(’t’)
f(x,y)=2+x^2*y
r=vector([cos(t), sin(t)])
tinicial=0
tfinal=pi
def integral_linea(integrando):
return RR(numerical_integral((integrando).subs(x=r[0], y=r[1]), tinicial, tfinal)[0])
ds=diff(r,t).norm()
dm=f(x,y)*ds
m=integral_linea(dm)
m
Supongamos que C es una curva suave por tramos; es decir, C es una unión de
una cantidad finita de curvas suaves C1 , C2 , . . . , Cn , donde, de acuerdo con
la Figura 18, el punto inicial de Ci+1 es el punto final de Ci . Entonces,
definimos la integral de f a lo largo de C como la suma de las integrales de f a
lo largo de cada una de las partes suaves de C:
ˆ ˆ ˆ
f (x, y)ds = f (x, y)ds + · · · + f (x, y)ds.
C C1 Cn
Integrales de línea

Figura 18: Una curva suave en tramos.

Ejemplo
´
7
Evalúe C 2xds, donde C consiste del arco C1 de la parábola y = x2 desde
(0, 0) hasta (1, 1) seguido por el segmento rectilíneo C2 desde (1, 1) hasta
(1, 2).

Solución
La curva C se muestra en la Figura 19. C1 es la gráfica de una función de x,
de modo que elegimos a x como el parámetro y las ecuaciones de C1 se
convierten en
x = x, y = x2 , , 0 ≤ x ≤ 1.
Integrales de línea

Figura 19: C = C1 ∪ C2 .

Por tanto,
ˆ ˆ ˆ
s 2 2
1  1
dx dy p
2xds = 2x + dx = 2x 1 + 4x2 dx
C1 0 dx dx 0
√

1

1 2 2 3/2  5 5−1
= · 1 + 4x  = ≈ 1.69672.
4 3 0 6
Integrales de línea
El código en SageMath para calcular la integral de línea para C1 es:
t = var(’t’)
f(x,y)=2*x
r=vector([t, t^2])
tinicial=0
tfinal=1
def integral_linea(integrando): return RR(numerical_integral((integrando).subs(x=r[0], y=r[1]),
tinicial, tfinal)[0])
ds=diff(r,t).norm()
dm=f(x,y)*ds
m=integral_linea(dm)
m
Sobre C2 elegimos a y como el parámetro, de modo que las ecuaciones de C2
son
x = 1, y = y, 1 ≤ y ≤ 2.
y
ˆ ˆ ˆ 2
s 
2 2  2
dx dy
2xds = 2(1) + dy = 2dy = 2
C2 1 dy dy 1

Por tanto,
ˆ ˆ ˆ √
5 5−1
2xds = 2xds + 2xds = + 2 ≈ 3.69672.
C C1 C2 6
Integrales de línea
El código en SageMath para calcular la integral de línea para C2 es:
t = var(’t’)
f(x,y)=2*x
r=vector([1, t])
tinicial=1
tfinal=2
def integral_linea(integrando): return RR(numerical_integral((integrando).subs(x=r[0], y=r[1]),
tinicial, tfinal)[0])
ds=diff(r,t).norm()
dm=f(x,y)*ds
m=integral_linea(dm)
m
´
Cualquier interpretación física de una integral de línea C f (x, y)ds depende de
la interpretación física de la función f . Suponga que ρ(x, y) representa la
densidad lineal en un punto (x, y) de un alambre delgado con forma de la curva
C. Entonces la masa de la parte del alambre desde Pi−1 hasta Pi , ver la
Figura 15, es aproximadamente ρ(x∗i , yi∗ )∆si y, así, la masa total del alambre es
P ∗ ∗
aproximadamente ρ(xi , yi )∆si . Al tomar más y más puntos sobre la curva
obtenemos la masa m del alambre como el valor límite de estas aproximaciones:
Xn ˆ
m = lı́m ρ(x∗i , yi∗ )∆si = ρ(x, y)ds.
n→∞ C
i=1

El centro de masa del alambre con función de densidad ρ se sitúa en el punto

(x̄, ȳ), donde
ˆ ˆ
1 1
x̄ = xρ(x, y)ds, ȳ = yρ(x, y)ds. (5)
m C m C
Integrales de línea
Ejemplo 8
Un alambre toma la forma de una semicircunferencia x2 + y 2 = 1, y ≥ 0, y es más grueso cerca de su
base que de la parte superior. Calcule el centro de masa del alambre si la densidad lineal en cualquier
punto es proporcional a su distancia desde la recta y = 1.

Solución
Como en el ejemplo 6, usamos la parametrización x = cos t, y = sen t, 0 ≤ t ≤ π; y encontramos que
ds = dt. La densidad lineal es
ρ(x, y) = k (1 − y)

donde k es una constante, por lo que la masa del alambre es

ˆ ˆ π
π
m= k (1 − y) ds = k (1 − sen t) dt = k [t + cos t]0 = k (π − 2) .
C 0

Según las ecuaciones 5 tenemos

ˆ ˆ
1 1
ȳ = yρ(x, y)ds = yk (1 − y) ds
m C k (π − 2) C
ˆ π π
1 1 1 1
 
2
= sen t − sen t dt = − cos t − t + sen 2t
π−2 0 π−2 2 4 0
4−π
= ≈ 0.38.
2 (π − 2)

Por simetría vemos que x̄ = 0, de modo que el centro de masa es

(0, (4 − π) / [2 (π − 2)]) ≈ (0, 0.38). Ver la Figura 20.
Integrales de línea

Figura 20: Centro de masa.

El código en SageMath para calcular el centro de masa del alambre es:

t = var(’t’)
densidad(x,y)=1-y
r=vector([cos(t), sin(t)])
tinicial=0
tfinal=pi
def integral_linea(integrando): return RR(numerical_integral((integrando).subs(x=r[0], y=r[1]),
tinicial, tfinal)[0])
ds=diff(r,t).norm()
dm=densidad(x,y)*ds
m=integral_linea(dm)
centro_masa_x=integral_linea(x*dm)/m
centro_masa_y=integral_linea(y*dm)/m
Para conocer la ordenada del centro de masa:
centro_masa_y
Obtenemos
0.375969196942054
Integrales de línea
Se pueden obtener otras dos integrales de línea reemplazando ∆si por ∆xi = xi − xi−1 o por
∆yi = yi − yi−1 en la Definición 3. Se les llama integrales de línea de f a lo largo de C respecto a x y
y:
ˆ n
∗ ∗

X
f (x, y)dx = lı́m f xi , yi ∆xi (6)
C n→∞
i=1

ˆ n
∗ ∗

X
f (x, y)dy = lı́m f xi , yi ∆yi . (7)
C n→∞
i=1
´
Cuando queramos distinguir la integral de línea original C f (x, y)ds de las integrales de línea 6 y 7, la
integral de línea original se denomina integral de línea respecto a la longitud de arco.

Las fórmulas siguientes establecen que las integrales de línea respecto a x y y se pueden también
evaluar expresando todo en términos de t: x = x(t), y = y(t), dx = x0 (t)dt y dy = y 0 (t)dt.

ˆ ˆ
0
f (x, y)dx = f (x(t), y(t))x (t)dt (8)
C C

ˆ ˆ
0
f (x, y)dy = f (x(t), y(t))y (t)dt (9)
C C

A menudo sucede que las integrales de línea respecto a x y y se presentan juntas. Cuando esto sucede,
se acostumbra abreviarlas escribiendo
ˆ ˆ ˆ
P (x, y)dx + Q(x, y)dy = P (x, y)dx + Q(x, y)dy. (10)
C C C
Integrales de línea
Algunas veces, al plantear una integral de línea lo más difícil es pensar en una
representación paramétrica de una curva cuya descripción geométrica se conoce. En
particular, con frecuencia necesitamos parametrizar un segmento rectilíneo, de modo
que es útil recordar que una representación vectorial del segmento rectilíneo que inicia
en r0 y termina en r1 está dado por
r(t) = (1 − t) r0 + tr1 , 0 ≤ t ≤ 1. (11)

Ejemplo
´
9
Evalúe C y 2 dx + xdy, donde
a) C = C1 es el segmento rectilíneo desde (−5, −3) hasta (0, 2), y
b) C = C2 es el arco de parábola x = 4 − y 2 desde (−5, −3) hasta (0, 2). Ver la
Figura 21.

Figura 21: Curva C1 ∪ C2 .

Integrales de línea
Solución
Para resolver esta integral de línea utilizaremos la forma abreviada 10.
a) Utilizando 11 con r0 = (−5, −3) y r1 = (0, 2), una representación paramétrica del segmento
rectilíneo es
x = 5t − 5, y = 5t − 3, 0 ≤ t ≤ 1.
Entonces dx = 5dt, dy = 5dt, y con las fórmulas 8 y 10 se tiene
ˆ ˆ 1
2 2
y dx + xdy = (5t − 3) (5dt) + (5t − 5) (5dt)
C1 0
ˆ 1 
2
=5 25t − 25t + 4 dt
0
#1
25t3 25t2
"
5
=5 − + 4t =− .
3 2 0
6

b) Puesto que la parábola está definida como una función de y, tomamos a y como el parámetro y
escribimos C2 como
2
x = 4 − y , y = y, −3 ≤ y ≤ 2.
Entonces dx = −2ydy y de acuerdo con las fórmulas 9 y 10 se tiene
ˆ ˆ 2  
2 2 2
y dx + xdy = y (−2y) dy + 4 − y dy
C2 −3
ˆ 2  
3 2
= −2y − y + 4 dy
−3
#2
y4 y3
"
= − − + 4y 5.
= 40 6
2 3 −3
Integrales de línea

I Observe que las respuestas de los incisos a) y b) del ejemplo 9 son

diferentes aun cuando las dos curvas tienen los mismos puntos extremos.
Por tanto, el valor de una integral de línea depende, en general, no sólo
de los puntos extremos de la curva, sino también de la trayectoria. (Véase
en la sección siguiente las condiciones en las cuales la integral es
independiente de la trayectoria.)
I Observemos también que las respuestas del ejemplo 9 dependen de la
dirección u orientación de la curva. Si −C1 denota el segmento rectilíneo
desde (0, 2) hasta (−5, −3), es posible verificar, mediante la
parametrización

x = −5t, y = 2 − 5t, 0≤t≤1

que ˆ
5
y 2 dx + xdy = .
−C1 6
Integrales de línea
En general, una parametrización dada x = x(t), y = y(t), a ≤ t ≤ b, determina
una orientación de una curva C, cuya dirección positiva corresponde a los
valores crecientes del parámetro t. (Véase la figura 22, en donde el punto inicial
A corresponde al valor del parámetro t = a y el punto terminal B corresponde
a t = b.)

Figura 22: Orientación de la integral de línea.

Si −C denota la curva que consiste de los mismos puntos que C, pero con la
orientación opuesta es decir, del punto inicial B al punto terminal A de la
Figura 22, entonces tenemos
Integrales de línea

ˆ ˆ ˆ ˆ
f (x, y) dx = − f (x, y) dx, f (x, y) dy = − f (x, y) dy.
−C C −C C

Pero si integramos respecto a la longitud de arco, el valor de la integral

de línea no cambia cuando se invierte la orientación de la curva:
ˆ ˆ
f (x, y) ds = f (x, y) ds.
−C C

La razón es que ∆si siempre es positiva, mientras que ∆xi y ∆yi

cambian de signo cuando se invierte la orientación de C.
Integrales de línea en el espacio
I Ahora supongamos que C es una curva suave en el espacio, dado por las
ecuaciones paramétricas

x = x(t), y = y(t), z = z(t), a≤t≤b

o por la ecuación vectorial r(t) = x(t)i + y(t)j + z(t)k.

I Si f es una función de tres variables que es continua en alguna región que
contiene a C, entonces definimos la integral de línea de f a lo largo de
C (respecto a la longitud de arco), de manera similar a la de las curvas
planas:
ˆ n
X
f (x, y, z)ds = lı́m f (x∗i , yi∗ , zi∗ ) ∆si .
C n→∞
i=1
I Evaluamos usando una fórmula similar a la fórmula 4:
ˆ ˆ
s 2 2 2
b  
dx dy dz
f (x, y, z)ds = f (x(t), y(t), z(t)) + + dt.
C a dt dt dt
(12)
I Observemos que las integrales en las fórmulas 4 y 12 se pueden escribir
en la forma vectorial más compacta
ˆ b
f (r(t)) r0 (t) dt.
 
a
Integrales de línea en el espacio
I En el caso especial de f (x, y, z) = 1, obtenemos
ˆ ˆ b  0 
ds = r (t) dt = L
C a

donde L es la longitud de la curva C.

I Las integrales de línea a lo largo de C respecto a x, y y z también se
pueden definir. Por ejemplo,
ˆ n
X
f (x, y, z)dz = lı́m f (x∗i , yi∗ , zi∗ ) ∆zi
C n→∞
i=1
ˆ b
= f (x(t), y(t), z(t)) z 0 (t)dt.
a

I Por tanto, como sucede con las integrales de línea en el plano, evaluamos
las integrales de la forma
ˆ
P (x, y, z)dx + Q(x, y, z)dy + R(x, y, z)dz (13)
C

expresando todo (x, y, z, dx, dy, dz) en términos del parámetro t.

Integrales de línea en el espacio

Ejemplo
´
10
Evalúe C y sen z ds, donde C es la hélice circular dada por las ecuaciones
x = cos t, y = sen t, z = t, 0 ≤ t ≤ 2π (ver la Figura 23).

Figura 23: Hélice circular.

Integrales de línea en el espacio
Solución
El resultado con la fórmula 12 es
ˆ ˆ
s 2 2 2
2π  
dx dy dz
y sen z ds = (sen t) sen t + + dt
C 0 dt dt dt
ˆ 2π √ ˆ 2π
p 1
= sen2 t sen2 t + cos2 t + 1dt = 2 (1 − cos 2t) dt
0 0 2
√  2π
2 1 √
= t − sen 2t = 2π.
2 2 0

El código en SageMath para calcular esta integral es

t = var(’t’)
f(x,y,z)= y*sin(z)
r=vector([cos(t), sin(t),t])
tinicial=0
tfinal=2*pi
def integral_linea(integrando):
return RR(numerical_integral((integrando).subs(x=r[0], y=r[1],z=r[2]), tinicial, tfinal)[0])
ds=diff(r,t).norm()
dm=f(x,y,z)*ds
m=integral_linea(dm)
m
cuyo resultado es: 4.44288293815837.
Integrales de línea en el espacio
Ejemplo
´
11
Evalúe C ydx + zdy + xdz, donde C consiste del segmento rectilíneo C1
desde (2, 0, 0) hasta (3, 4, 5) seguido por el segmento vertical C2 desde (3, 4, 5)
hasta (3, 4, 0).

Figura 24: Curva C consistente de los segmentos rectilíneos C1 y C2 .

Solución
La curva C se ilustra en la Figura 24. Utilizando la ecuación 11, expresamos a
C1 como
r(t) = (1 − t) (2, 0, 0) + t (3, 4, 5) = (2 + t, 4t, 5t)
Integrales de línea en el espacio
o bien, en forma paramétrica, como

x = 2 + t, y = 4t, z = 5t, 0 ≤ t ≤ 1.

Por tanto,

ˆ ˆ 1
ydx + zdy + xdz = (4t) dt + (5t) 4dt + (2 + t) 5dt
C1 0
ˆ 1
29 2 1
 
= (10 + 29t) dt = 10t + t = 24.5
0 2 0

De manera similar, C2 se puede expresar en la forma

r(t) = (1 − t) (3, 4, 5) + t (3, 4, 0) = (3, 4, 5 − 5t)

o bien, x = 3, y = 4, z = 5 − 5t, 0 ≤ t ≤ 1. Entonces dx = dy = 0, de modo que

ˆ ˆ 1
ydx + zdy + xdz = 3 (−5) dt = −15.
C2 0

Al sumar los valores de estas integrales

ˆ
ydx + zdy + xdz = 24.5 − 15 = 9.5.
C
Integrales de línea de campos vectoriales
I Ya sabemos que el trabajo realizado por una fuerza variable f (x) que mueve
´
una partícula desde a hasta b a lo largo del eje x es W = ab f (x)dx. También
sabemos que el trabajo que efectúa una fuerza constante F al mover un objeto
desde el punto P hasta otro punto Q en el espacio es W = F · D, donde
−−
→
D = P Q es el vector desplazamiento.
I Ahora supongamos que F = P i + Qj + Rk es un campo de fuerzas continuo
sobre R3 , tal como el campo gravitacional del Ejemplo 3 o el campo de fuerzas
eléctricas del Ejemplo 4. (Un campo de fuerzas sobre R2 se puede considerar
como un caso especial donde R = 0 y P y Q dependen sólo de x y de y.)
Deseamos calcular el trabajo que realizó esta fuerza al mover la partícula a lo
largo de la curva suave C.
I Dividimos C en subarcos Pi−1 Pi de longitudes ∆si dividiendo el intervalo del
parámetro [a, b] en subintervalos de igual ancho. (Véase en la Figura 15 el caso
bidimensional
o en la Figura 25 el caso tridimensional.) Elegimos un punto
Pi∗ x∗i , yi∗ , zi∗ sobre el i-ésimo subarco que corresponde al valor del parámetro
t∗i . Si ∆si es pequeño, entonces cuando la partícula se mueve de Pi−1 hasta Pi
a lo largo de la curva, lo hace aproximadamente en la dirección de T t∗i , el


vector unitario tangente a Pi∗ . Por tanto, el trabajo que efectúa la fuerza F al
mover la partícula desde Pi−1 hasta Pi es aproximadamente de

F (x∗i , yi∗ , zi∗ ) · (∆si T (t∗i )) = (F (x∗i , yi∗ , zi∗ ) · T (t∗i )) ∆si
Integrales de línea de campos vectoriales

Figura 25: Subdivisión de la curva C en n subarcos.

y el trabajo total realizado al mover la partícula a lo largo de C es aproximadamente

n
∗ ∗ ∗
∗


X
F xi , yi , zi · T ti ∆si (14)
i=1

donde T(x, y, z) es el vector unitario tangente en el punto (x, y, z) sobre C. Intuitivamente es posible
ver que estas aproximaciones deben ser mejores a medida que n se incrementa. Por tanto, definimos el
trabajo W realizado por el campo de fuerza F como el límite de las sumas de Riemann, a saber,

ˆ ˆ
W = F (x, y, z) · T (x, y, z) ds = F · Tds. (15)
C C

La Ecuación 15 establece que el trabajo es la integral de línea respecto a la longitud de arco de la

componente tangencial de la fuerza.
Integrales de línea de campos vectoriales
Si la curva C está dada por la ecuación vectorial r(t) = x(t)i + y(t)j + z(t)k,
entonces T(t) = r0 (t)/ |r0 (t)|, de modo que al usar la Ecuación 12 podemos volver a
expresar la Ecuación 15 en la forma
ˆ ˆ b
b r0 (t)  0 
 
W = F (r(t)) · 0 r (t) dt = F (r(t)) · r0 (t)dt.
a |r (t)| a

´
Esta integral se abrevia a menudo como C F · dr y se encuentra también en otras
especialidades de la física. Por tanto, proporcionamos la siguiente definición para la
integral de línea de cualquier campo vectorial continuo.

Definición 4
Sea F un campo vectorial continuo definido sobre una curva suave C dada por una
función vectorial r(t), a ≤ t ≤ b. Entonces la integral de línea de F a lo largo de C es
ˆ ˆ b ˆ
F · dr = F (r(t)) · r0 (t)dt = F · T ds.
C a C

Observación
Al utilizar la Definición 4, recuerde que F (r(t)) es sólo una forma de abreviar
F (x(t), y(t), z(t)), de modo que evaluamos F (r(t)) haciendo simplemente x = x(t),
y = y(t) y z = z(t) en la expresión para F(x, y, z). Observe también que podemos
escribir formalmente dr = r0 (t)dt.
Integrales de línea de campos vectoriales
Ejemplo 12
Determine el trabajo efectuado por el campo de fuerza F(x, y) = x2 i − xyj
cuando mueve una partícula a lo largo del cuarto de circunferencia
r(t) = cos ti + sen tj, 0 ≤ t ≤ π/2.

Solución
Puesto que x = cos t y y = sen t

F (r(t)) = cos2 t i − cos t sen t j

y
r0 (t) = − sen t i + cos t j.
Por tanto, el trabajo realizado es
ˆ ˆ π/2 ˆ π/2
F (r(t)) · r0 (t)dt = −2 cos2 t sen t dt


F · dr =
C 0 0
π/2
cos3 t

2
= 2 =− .
3 0 3

Ver la Figura 26.

Integrales de línea de campos vectoriales

Figura 26: Esta figura muestra el campo de fuerza y la curva del ejemplo 12. El
trabajo hecho es negativo porque el campo obstruye el movimiento a lo largo de la
curva.

Observación
´ ´
Aun cuando C F · dr = C F · T ds y las integrales respecto a la longitud de
arco permanecen sin cambio cuando se invierte la dirección, se sigue
cumpliendo que ˆ ˆ
F · dr = − F · dr
C C
porque el vector unitario tangente T es reemplazado por su negativo cuando C
es reemplazado por −C.
Integrales de línea de campos vectoriales
Ejemplo
´
13
Evalúe C F · dr, donde F(x, y, z) = xy i + yz j + zx k y C es la cúbica
alabeada dada por

x = t, y = t2 , z = t3 , 0 ≤ t ≤ 1.

Solución
Tenemos
r(t) = t i + t2 j + t3 k
r0 (t) = i + 2tj + 3t2 k
F (r(t)) = t3 i + t5 j + t4 k.
Por tanto,
ˆ ˆ 1
F · dr = F (r(t)) · r0 (t)dt
C 0
ˆ 1  4 1
t 5t7 27
t3 + 5t6 dt =


= + = .
0 4 7 0 28

Ver la Figura 27.

Integrales de línea de campos vectoriales

Figura 27: En esta figura se ilustra la cúbica torcida C del ejemplo 13 y algunos
vectores representativos que actúan en tres puntos sobre C.

Para finalizar, hacemos notar la relación entre las integrales de línea de los
campos vectoriales y las integrales de línea de los campos escalares.
Supongamos que el campo vectorial F sobre R3 está dado en la forma de
componentes mediante la ecuación F = P i + Qj + Rk. Ulilizamos la
Definición 4 para calcular su integral de línea a lo largo de C.
Integrales de línea de campos vectoriales

ˆ ˆ b
0
F · dr = F (r(t)) · r (t)dt
C a
ˆ b  
0 0 0
= (P i + Qj + Rk) · x (t)i + y (t)j + z (t)k dt
a
ˆ b h i
0 0 0
= P (x(t), y(t), z(t)) x (t) + Q (x(t), y(t), z(t)) y (t) + R (x(t), y(t), z(t)) z (t) dt.
a

Pero esta última integral es precisamente la integral de línea de 13. Por tanto,
tenemos ˆ ˆ
F · dr = P dx + Q dy + R dz
C C
donde F = P i + Qj + Rk.´
Por ejemplo,
´ la integral C y dx + z dy + x dz del ejemplo 11 se podría expresar
como C F · dr donde

F(x, y, z) = y i + z j + x k.
Teorema fundamental de las integrales de línea
Sabemos que el teorema fundamental del cálculo se puede expresar como
ˆ b
0
F (x) = F (b) − F (a) (16)
a
0
donde F es continua sobre [a, b]. A la Ecuación 16 también se le conoce con el nombre de Teorema del
cambio neto: la integral de la razón de cambio es el cambio neto.

Teorema 1
Sea C una curva suave definida por la función vectorial r(t), a ≤ t ≤ b. Sea f la función derivable de
dos o tres variables cuyo vector gradiente ∇f es continuo sobre C. Entonces

ˆ
∇f · dr = f (r(b)) − f (r(a)) .
C

Demostración.
Utilizando la Definición 4 tenemos

ˆ ˆ b
0
∇f · dr = ∇f (r(t)) · r (t)dt
C a
ˆ b
∂f dx ∂f dy ∂f dz
 
= + + dt
a ∂x dt ∂y dt ∂z dt
ˆ b
d
= f (r(t)) dt por la regla de la cadena
a dt
= f (r(b)) − f (r(a)) .

El último paso se infiere del teorema fundamental del cálculo (Ecuación 16).
Teorema fundamental de las integrales de línea
Observación
El Teorema 1 establece que se puede evaluar la integral de línea de un campo
vectorial conservativo (el campo vectorial del gradiente de la función potencial
f ) simplemente si se conoce el valor de f en los puntos extremos de C. De
hecho, el Teorema 1 establece que la integral de línea de ∇f es el cambio neto
en f . Si f es una función de dos variables y C es una curva plana cuyo punto
inicial es A (x1 , y1 ) y el terminal es B (x2 , y2 ), como en la Figura 28, entonces
el Teorema 1 se transforma en
ˆ
∇f · dr = f (x2 , y2 ) − f (x1 , y1 ) .
C

Si f es una función de tres variables y C es una curva en el espacio que une el

punto A (x1 , y1 , z1 ) con el punto B (x2 , y2 , z2 ), entonces tenemos
ˆ
∇f · dr = f (x2 , y2 , z2 ) − f (x1 , y1 , z1 ) .
C

Observación
Aunque se ha demostrado el Teorema 1 para el caso de curvas suaves, también
es válido para curvas suaves por tramos. Esto se puede ver al subdividir C en
un número finito de curvas suaves y sumar las integrales resultantes.
Teorema fundamental de las integrales de línea

Figura 28: Curva plana con punto inicial A (x1 , y1 ) y punto final B (x2 , y2 ).
Teorema fundamental de las integrales de línea
Ejemplo 14
Calcule el trabajo que realiza el campo gravitacional

mM G
F(x) = − x
|x|3

al mover una partícula de masa m desde el punto (3, 4, 12) hasta el punto (2, 2, 0) a lo
largo de la curva C suave por tramos.

Solución
Sabemos que F es un campo vectorial conservativo y, de hecho, F = ∇f , donde

mM G
f (x, y, z) = p .
x2 + y 2 + z 2

Por tanto, según el Teorema 1, el trabajo realizado es

ˆ ˆ
W = F · dr = ∇f · dr
C C
= f (2, 2, 0) − f (3, 4, 12)
 
mM G mM G 1 1
= √ −√ = mM G √ − .
2 2 + 22 32 + 42 + 122 2 2 13
Teorema fundamental de las integrales de línea
Independencia de la trayectoria
Supongamos que C1 y C2 son dos curvas suaves por tramos
(denominadas trayectorias) que tienen el mismo punto inicial A y el
punto
´ terminal ´B. Sabemos por el Ejemplo 9 que, en general,
C1 ∇f · dr 6= C2 ∇f · dr. Pero una implicación del Teorema 1 es
que ˆ ˆ
∇f · dr = ∇f · dr
C1 C2

cuando ∇f es continuo. En otras palabras, la integral de línea de

un campo vectorial conservativo depende sólo del punto inicial y del
punto terminal de la curva.
En general, si F es ´un campo vectorial continuo cuyo dominio es D,
´ integral de ´línea C F · dr es independiente de la trayectoria si
la
C1 F · dr = C2 F · dr para cualesquiera dos trayectorias C1 y C2
en D que tienen los mismos puntos iniciales y terminales. Con esta
terminología podemos decir que las integrales de línea de campos
vectoriales conservativos son independientes de la trayectoria.
Teorema fundamental de las integrales de línea
Se dice que una curva es cerrada si su punto final coincide
´ con su punto
terminal, es decir, r(b) = r(a) (véase la Figura 29). Si C F · dr es
independiente de la trayectoria en D y C es cualquier trayectoria cerrada en D,
podemos escoger dos puntos cualesquiera A y B sobre la curva C y considerar
que C está compuesta de la trayectoria C1 desde A hasta B seguida por la
trayectoria C2 desde B hasta A (véase la Figura 30)1 . Luego
˛ ˆ ˆ ˆ ˆ
F · dr = F · dr + F · dr = F · dr − F · dr = 0
C C1 C2 C1 −C2

puesto que C1 y −C2 poseen los ¸ mismos puntos inicial y terminal.

Inversamente, si es cierto que C F · dr = 0 cuando C es una trayectoria
cerrada en D, entonces se demuestra la independencia de la trayectoria como
sigue. Tomamos dos trayectorias cualesquiera C1 y C2 desde A hasta B en D y
definimos C como la curva cerrada que consiste de C1 seguida de C2 . Entonces
˛ ˆ ˆ ˆ ˆ
0= F · dr = F · dr + F · dr = F · dr − F · dr
C C1 −C2 C1 C2
´ ´
y de este modo C1 F · dr = C2
F · dr. Por tanto, hemos demostrado el
teorema siguiente.
1
Si la¸trayectoria de una integral de línea es una curva cerrada, la notación
usual es C F · dr.
Teorema fundamental de las integrales de línea

Figura 29: Curva cerrada.

Figura 30: Curva cerrada con dos puntos cualesquiera.

Teorema fundamental de las integrales de línea

Teorema
´ 2 ¸
C
F · dr es independiente de la trayectoria en D si y sólo si C
F · dr = 0 para
toda trayectoria cerrada C en D.
I Puesto que sabemos que la integral de línea de cualquier campo vectorial
conservativo
¸ F es independiente de la trayectoria, se infiere que
C
F · dr = 0 para cualquier trayectoria cerrada. La interpretación física
es que el trabajo hecho por un campo de fuerzas conservativo (como el
campo gravitacional o eléctrico estudiado en la primera sección), cuando
se desplaza un objeto alrededor de una trayectoria cerrada es 0.
I El teorema siguiente establece que los únicos campos vectoriales que son
independientes de la trayectoria son conservativos. Está establecido y
demostrado para curvas planas, pero hay una versión similar para las
curvas en el espacio. Si suponemos que D es abierto, lo que significa que
para todo punto P en D hay un disco con centro P que está totalmente
en D. (De esta manera, D no contiene ningún punto de su frontera).
Además, supongamos que D es conexa. Esto quiere decir que dos puntos
cualesquiera en D se pueden unir mediante una trayectoria que está en D.
Teorema fundamental de las integrales de línea
Teorema 3
Supongamos que F es´ un campo vectorial que es continuo sobre una región
conexa abierta D. Si C F · dr es independiente de la trayectoria en D,
entonces F es un campo vectorial conservativo sobre D, es decir, existe una
función f tal que ∇f = F.

Demostración.
Sea A(a, b) un punto fijo en D. Construimos la función potencial f deseada
definiendo ˆ (x,y)
f (x, y) = F · dr
(a,b)
´
para cualquier punto (x, y) en D. Puesto que C F · dr es independiente de la
trayectoria, no importa qué trayectoria C desde (a, b) hasta (x, y) se recorra
para evaluar f (x, y). Como D es abierto, existe un disco que está contenido en
D cuyo centro es (x, y).
Elijamos cualquier punto (x1 , y) en el disco con x1 < x, y sea C que consiste
de cualquier trayectoria C1 desde (a, b) hasta (x1 , y) seguido por el segmento
rectilíneo horizontal C2 desde (x1 , y) hasta (x, y) (véase la Figura 31).
Entonces,
Teorema fundamental de las integrales de línea
continuación de la demostración.

Figura 31: D abierto y conexo.

ˆ ˆ ˆ (x,y) ˆ
f (x, y) = F · dr + F · dr = F · dr + F · dr.
C1 C2 (a,b) C2

Observemos que la primera de estas integrales no depende de x, de modo que

ˆ
∂ ∂
f (x, y) = 0 + F · dr.
∂x ∂x C2
Teorema fundamental de las integrales de línea
continuación de la demostración.
Si escribimos F = P i + Q j, entonces
ˆ ˆ
F · dr = P dx + Qdy
C2 C2

En C2 , y es constante, de modo que dy = 0. Usando t como el parámetro,

donde x1 ≤ t ≤ x, tenemos
ˆ ˆ x
∂ ∂ ∂
f (x, y) = P dx + Qdy = P (t, y)dt = P (x, y)
∂x ∂x C2 ∂x x1

de acuerdo con el Teorema Fundamental del Cálculo. Al utilizar un segmento

vertical (véase la Figura 31), un razonamiento similar demuestra que
ˆ ˆ y
∂ ∂ ∂
f (x, y) = P dx + Qdy = Q(x, t)dt = Q(x, y)
∂y ∂y C2 ∂y y1

Por tanto,
∂f ∂
F = P i + Qj = i+ j = ∇f
∂x ∂y
lo cual establece que F es conservativo.
Teorema fundamental de las integrales de línea
La pregunta inmediata: ¿Cómo es posible determinar si un campo vectorial F
es conservativo? Supongamos que ya se sabe que F = P i + Q j es
conservativo, donde P y Q tienen derivadas parciales continuas de primer
orden. Entonces, hay una función f tal que F = ∇f , es decir,
∂f ∂f
P = y Q=
∂x ∂y
Por tanto, según el teorema de Clairaut,

∂P ∂2f ∂2f ∂Q
= = = .
∂y ∂y∂x ∂x∂y ∂x

Teorema 4
Si F(x, y) = P (x, y) i + Q(x, y) j es un campo vectorial conservativo, donde P
y Q tienen derivadas parciales continuas de primer orden sobre un dominio D,
entonces en todo D tenemos
∂P ∂Q
= .
∂y ∂x
Teorema fundamental de las integrales de línea
El inverso del Teorema 4 es válido sólo para un tipo especial de región. Para
explicarlo necesitamos primero el concepto de una curva simple, la cual es una
curva que no se corta a sí misma en ninguna parte entre sus puntos extremos.
Véase la Figura 32; r(a) = r(b) para una curva cerrada simple, pero
r(t1 ) 6= r(t2 ) cuando a < t1 < t2 < b.

Figura 32: Tipos de curvas.

Teorema fundamental de las integrales de línea
En el Teorema 3 necesitamos una región conexa abierta. En el caso del teorema
siguiente requerimos una condición más rigurosa. Una región simplemente
conexa en el plano es una región conexa D tal que toda curva simple cerrada
en D abarca sólo puntos que están en D. Observe que, según la Figura 33,
intuitivamente hablando, una región simplemente conexa no contiene agujeros
y no puede consistir de dos partes separadas.

Figura 33: Regiones simplemente conexas.

Teorema fundamental de las integrales de línea
En términos de regiones simplemente conexas, podemos ahora enunciar un
inverso parcial del Teorema 4 que proporciona un método conveniente para
comprobar que el campo vectorial sobre R2 es conservativo. La demostración se
delinea en la sección siguiente como una consecuencia del teorema de Green.

Teorema 5
Sea F = P i + Q j un campo vectorial sobre una región simplemente conexa D.
Supongamos que P y Q tienen derivadas continuas de primer orden y
∂P ∂Q
=
∂y ∂x
en toda la región D. Entonces F es conservativo.

Ejemplo 15
Determine si el campo vectorial

F(x, y) = (x − y) i + (x − 2) j

es conservativo.
Teorema fundamental de las integrales de línea
Solución
Sea P (x, y) = x − y y Q(x, y) = x − 2. Entonces
∂P ∂Q
= −1, = 1.
∂y ∂x
Como ∂P/∂y 6= ∂Q/∂x, F no es conservativo según el Teorema 4.

Ejemplo 16
Determine si el campo vectorial

F(x, y) = (3 + 2xy) i + x2 − 3y 2 j

es conservativo.

Solución
Sea P (x, y) = 3 + 2xy y Q(x, y) = x2 − 3y2 . Entonces
∂P ∂Q
= 2x = .
∂y ∂x

Asimismo, el dominio de F es todo el plano, D = R2 , el cual es abierto y

simplemente conexo. Por tanto, podemos aplicar el Teorema 5 y concluir que F
es conservativo.
Teorema fundamental de las integrales de línea

I En las Figuras 34 y 35 se muestran los campos vectoriales de los

Ejemplos 15 y 16, respectivamente. Los vectores de la Figura 34 que
inician en la curva cerrada C parecen apuntar aproximadamente en
la
¸ misma dirección que C. De este modo parece como si
C
F · dr > 0 y por consiguiente F no es conservativo. Los cálculos
del Ejemplo 15 confirma esta impresión.
I Algunos de los vectores cercanos a las curvas C1 y C2 de la
Figura 35 apuntan casi en la misma dirección que las curvas,
mientras que otros señalan la dirección opuesta. Entonces parece
posible que las integrales de línea alrededor de todas las trayectorias
cerradas son 0. En el Ejemplo 16 se demuestra que, en efecto, F es
conservativo.
Teorema fundamental de las integrales de línea

Figura 34: Campo vectorial del Ejemplo 15.

Figura 35: Campo vectorial del Ejemplo 16

Teorema fundamental de las integrales de línea
En el Ejemplo 16, el Teorema 5 señala que F es conservativo, pero no dice cómo determinar la función
potencial f tal que F = ∇f . La demostración del Teorema 3 da una pista para encontrar f . Utilizamos
la “integración parcial”, como en el ejemplo siguiente.

Ejemplo 17  
Si F(x, y) = (3 + 2xy) i + x2 − 3y 2 j,

a) determine una función f tal que F = ∇f .

´
b) Evalúe la integral de línea C F · dr, donde C es la curva definida por

t t
r(t) = e sen t i + e cos t j, 0 ≤ t ≤ π.

Solución

a) De acuerdo con el Ejemplo 16 sabemos que F es conservativo y, por tanto, existe una función f
con ∇f = F, es decir,
∂f (x, y)
= 3 + 2xy (17)
∂x
y
∂f (x, y) 2 2
= x − 3y (18)
∂y
Teorema fundamental de las integrales de línea
Al integrar 17 respecto a x, obtenemos

f (x, y) = 3x + x2 y + g(y) (19)

Observe que la constante de integración es una constante respecto a x, es decir, una

función de y, la cual se llama g(y). Luego derivamos parcialmente ambos miembros de
19 respecto a y:
∂f (x, y)
= x2 + g 0 (y). (20)
∂y
Comparando 18 y 20, vemos que

g 0 (y) = −3y 2 .

Si integramos respecto a y, tenemos

g(y) = −y 3 + K

donde K es una constante. Al sustituir en 19 tenemos

f (x, y) = 3c + x2 y − y 3 + K

como la función potencial deseada.

Teorema fundamental de las integrales de línea
b) Para utilizar el Teorema 1, todo lo que tenemos que saber es
cuáles son los puntos inicial y terminal de C, a saber,
r(0) = (0, 1) y r(π) = (0, −eπ ). En la expresión para f (x, y)
del inciso a), cualquier valor de la constante K servirá, de
modo que seleccionamos K = 0. Luego tenemos
ˆ ˆ
F·dr = ∇f ·dr = f (0, −eπ )−f (0, 1) = e3π −(−1) = e3π +1.
C C

Este método es mucho más corto que el método directo para

evaluar las integrales de línea que se tratan en la sección
anterior.
I En una sección posterior se da un criterio para determinar si el
campo vectorial F sobre R3 es conservativo o no. Mientras
tanto, el ejemplo siguiente muestra que la técnica para hallar
la función potencial es en gran medida la misma que para los
campos vectoriales sobre R2 .
Teorema fundamental de las integrales de línea
Ejemplo 18
 
Si F(x, y, z) = y 2 i + 2xy + e3z j + 3ye3z k, determine una función f tal que ∇f = F.

Solución
Si hay tal función f , entonces
∂f (x, y, z) 2
=y (21)
∂x
∂f (x, y, z) 3z
= 2xy + e (22)
∂y

∂f (x, y, z) 3z
= 3ye (23)
∂z
Al integrar 21 respecto a x obtenemos

2
f (x, y, z) = xy + g(y, z) (24)

donde g(y, z) es una constante respecto a x. Luego, al derivar parcialmente 24 respecto a y, tenemos

∂f (x, y, z) ∂g(y, z)
= 2xy +
∂y ∂y

y la comparación con 22 da
∂g(y, z) 3z
=e .
∂y
Teorema fundamental de las integrales de línea

Integrando esta última derivada parcial respecto a y, obtenemos

g(y, z) = ye3z + h(z) y reescribimos 24 como

f (x, y, z) = xy 2 + ye3z + h(z).

Para terminar, al derivar respecto a z, obtenemos

∂f (x, y, z)
= 3ye3z + h0 (z)
∂z
y comparando con 23, obtenemos h0 (z) = 0 y, por tanto, h(z) = K, una
constante. La función deseada es

f (x, y, z) = xy 2 + ye3z + K.

Se comprueba con facilidad que ∇f = F.

Teorema fundamental de las integrales de línea

Conservación de la energía
Aplicaremos las ideas de este tema a un campo de fuerzas continuo F que hace que se desplace un
objeto a lo largo de una trayectoria C definida por r(t), a ≤ t ≤ b, donde r(a) = A es el punto inicial
y r(b) = B es el punto terminal de C. De acuerdo con la segunda ley de Newton del movimiento, la
fuerza F (r(t)) en un punto sobre C se relaciona con la aceleración a(t) = r00 (t) mediante la ecuación

00
F (r(t)) = mr (t).

De modo que el trabajo que efectúa la fuerza sobre el objeto es

ˆ ˆ b ˆ b
0 00 0
W = F · dr = F (r(t)) · r (t)dt = mr (t) · r (t)dt
C a a
ˆ
m b d h 0 0
i
= r (t) · r (t) dt
2 a dt
ˆ 2 b
m b d  0 m  0

2

= r (t) dt = r (t)

2 a dt 2 a
m  0
 2  2 
 0
= r (b) − r (a) .
 
2

Por tanto,
1 2 1 2
W = m |v(b)| − m |v(a)| (25)
2 2

donde v = r0 es la velocidad.
Teorema fundamental de las integrales de línea
La cantidad 12 m |v(t)|2 es decir, la mitad de la masa por el cuadrado de la rapidez, se
llama energía cinética del objeto. Por tanto, podemos volver a escribir la Ecuación 25
como
W = K(B) − K(A) (26)
la cual establece que el trabajo hecho por el campo de fuerzas a lo largo de C es igual
al cambio de la energía cinética en los puntos extremos de C.
Supongamos además que F es un campo de fuerzas conservativo, es decir, podemos
escribir F = ∇f . En física, la energía potencial de un objeto en el punto (x, y, z) se
define como P (x, y, z) = f (x, y, z), de modo que tenemos F = −∇P . Entonces,
según el Teorema 1, tenemos
ˆ ˆ
W = F · dr = − ∇P · dr = − [P (r(b)) − P (r(a))] = P (A) − P (B).
C C

Comparando esta ecuación con la Ecuación 26, encontramos que

P (A) + K(A) = P (B) + K(B)

la cual establece que si un objeto se mueve desde un punto A hacia un punto B bajo
la influencia de un campo de fuerzas conservativo, entonces la suma de su energía
potencial y de su energía cinética permanece constante. Este enunciado recibe el
nombre de ley de la conservación de la energía, y es la razón de que el campo
vectorial se llame conservativo.
Teorema de Green
El teorema de Green establece la relación entre una integral de línea alrededor
de una curva simple cerrada C y una integral doble sobre la región plana D
acotada por C (véase la Figura 36). Supongamos que D consiste de todos los
puntos dentro de C, así como de todos los puntos sobre C. En el
planteamiento del teorema de Green se usa la convención de que la orientación
positiva de una curva simple cerrada C se refiere a un recorrido sencillo de C
en el sentido contrario al de las manecillas del reloj. Por tanto, si C está
definida por la función vectorial r(t), a ≤ t ≤ b, entonces la región D está
siempre a la izquierda cuando el punto r(t) recorre C (véase la Figura 37).

Figura 36: D consiste de todos los puntos dentro de C, así como de todos los puntos
sobre C.
Teorema de Green

Figura 37: La región D está siempre a la izquierda cuando el punto r(t) recorre C.

Teorema 6 (Teorema de Green)

Sea C una curva simple cerrada, suave por tramos con orientación positiva en el plano, y sea D la
región que delimita C. Si P y Q tienen derivadas parciales continuas sobre una región abierta que
contiene a D, entonces ‰ ¨ 
∂Q ∂P

P dx + Qdy = − dA.
C ∂x ∂y
D

Observación
Se usa la notación C P dx + Qdy para señalar que la integral de línea se calcula usando la orientación
positiva de la curva cerrada C.

Observación  
Recuerde que el primer miembro de esta ecuación C P dx + Qdy es otra forma de escribir C F · dr
donde F = P i + Qj.
Teorema de Green

I El teorema de Green se debe considerar como el equivalente del teorema

fundamental del cálculo para las integrales dobles.
I El teorema de Green no es fácil de demostrar en general.

Ejemplo

19  p 
sen x
Evalúe C
(3y − e ) dx + 7x + y 4 + 1 dy, donde C es la circunferencia

x2 + y 2 = 9.

Solución
La región D acotada por C es el disco x2 + y 2 ≤ 9, de modo que cambiemos a
coordenadas polares después de aplicar el teorema de Green:
‰ ¨ 
∂ ∂ 
 q   q  
sen x
sen x 
3y − e dx + 7x + y 4 + 1 dy = 7x + y 4 + 1 − 3y − e dA
C ∂x ∂y
D
ˆ 2π ˆ 3 ˆ 2π ˆ 3
= (7 − 3) r dr dθ = 4 dθ rdr = 36π.
0 0 0 0
Teorema de Green
Ejemplo

20
Evalúe C x4 dx + xydy, donde C es la curva triangular que consiste de los
segmentos rectilíneos de (0, 0) a (1, 0), de (1, 0) a (0, 1) y de (0, 1) a (0, 0).

Solución
Aunque la integral de línea dada se podría evaluar como se acostumbra
mediante los métodos de la segunda sección, eso significaría plantear tres
integrales separadas a lo largo de los tres lados del triángulo, de modo que en
lugar de eso aplicaremos el teorema de Green. Observe que la región D
encerrada por C es simple y C sigue una orientación positiva (véase la
Figura 38). Si hacemos P (x, y) = x4 y Q(x, y) = xy, entonces tenemos
‰ ¨  
∂ ∂  4
x4 dx + xydy = [xy] − x dA
C ∂x ∂y
D
ˆ 1 ˆ 1−x y=1−xˆ 1 
1 2
= (y − 0) dydx =
y dx
0 0 0 2 y=0
ˆ 1
1 1

1 1
= (1 − x)2 dx = − (1 − x)3 = .
2 0 6 0 6
Teorema de Green

Figura 38: C es una región triangular.

Teorema de Green

I En los Ejemplos 19 y 20 parece que es más fácil evaluar la integral doble

que la integral de línea. (¡Trate de plantear la integral de línea del
Ejemplo 19, y pronto se convencerá!) Pero algunas veces es más fácil
evaluar la integral de línea, y utilizar el teorema de Green en la dirección
inversa. Por ejemplo, si sabemos que P (x, y) = Q(x, y) = 0 sobre la
curva C, entonces el teorema de Green da
¨   ‰
∂Q ∂P
− dA = P dx + Qdy = 0
∂x ∂y C
D

sin importar qué valores tomen P y Q en la región D.

I Otra aplicación de la dirección inversa˜del teorema de Green es para
calcular áreas. Como el área de D es D 1dA, seleccionamos P y Q tales
que
∂Q ∂P
− = 1.
∂x ∂y
Teorema de Green
Hay varias posibilidades:
1
P (x, y) = 0, P (x, y) = −y,
P (x, y) = − y
2
1
Q(x, y) = x, Q(x, y) = 0, Q(x, y) = x.
2
Entonces el teorema de Green da las fórmulas siguientes para el área de D:
˛ ˛ ˛
1
A= xdy = − ydx = xdy − ydx. (27)
C C 2 C

Ejemplo 21
x2 y2
Determine el área delimitada por la elipse + 2 = 1.
a2 b
Solución
Las ecuaciones paramétricas de la elipse son x = a cos t y y = b sen t , donde
0 ≤ t ≤ 2π. Al aplicar la tercera fórmula de la Ecuación 27, tenemos
ˆ
1
A= xdy − ydx
2 C
ˆ 2π
1
= (a cos t) (b cos t) dt − (b sen t) (−a sen t) dt
2 0
ˆ 2π
ab
= dt = πab.
2 0
Teorema de Green
I La fórmula 27 se puede utilizar para explicar cómo funcionan los planímetros. Un planímetro es
un instrumento mecánico que se emplea para medir el área de una región al trazar su curva
frontera. Estos aparatos son útiles en todas las ciencias: en biología para medir el área de hojas o
alas, en medicina para medir el tamaño de secciones transversales de órganos o tumores, en
bosques para estimar el tamaño de regiones pobladas de árboles a partir de fotografías.
I La Figura 39 muestra la operación de un planímetro polar: el polo se fija y, cuando el trazador se
mueve a lo largo de la curva frontera de la región, la rueda se desliza parcialmente y rueda
también parcialmente en forma perpendicular al brazo trazador. El planímetro mide la distancia
que la rueda gira y ésta es proporcional al área encerrada. La explicación como una consecuencia
de la fórmula 27 se puede hallar en los siguientes artículos:
I R. W. Gatterman, “The planimeter as an example of Green’s Theorem”, Amer. Math.
Monthly, vol. 88 (1981), pp. 701-704.
I Tanya Leise, “As the planimeter wheel turns”, College Math, Journal, vol. 38. (2007), pp.
24-31.

Figura 39: Planímetro polar Keuffel & Esser.

Versiones extendidas del teorema de Green
Aunque hemos estudiado el teorema de Green sólo para el caso donde D es simple, ya
podemos generalizarlo al caso donde D es una unión finita de regiones simples. Por
ejemplo, si D es la región mostrada en la Figura 40, entonces podemos escribir
D = D1 ∪ D2 , donde D1 y D2 también son simples. La frontera de D1 es C1 ∪ C3 y
la frontera de D2 es C2 ∪ (−C3 ), por lo que, al aplicar el teorema de Green para D1 y
D2 por separado, obtenemos
‰ ¨  
∂Q ∂P
P dx + Qdy = − dA
C1 ∪C3 ∂x ∂y
D1

‰ ¨  
∂Q ∂P
P dx + Qdy = − dA
C2 ∪(−C3 ) ∂x ∂y
D2

Figura 40: Unión de regiones simples.

Versiones extendidas del teorema de Green
Si sumamos estas dos ecuaciones, las integrales de línea a lo largo de C3 y −C3
se cancelan, de modo que
‰ ¨  
∂Q ∂P
P dx + Qdy = − dA
C1 ∪C2 ∂x ∂y
D

que es el teorema de Green para D = D1 ∪ D2 , puesto que su frontera es

C = C1 ∪ C2 .
I El mismo tipo de razonamiento permite establecer el teorema de Green
para cualquier unión finita de regiones simples que no se traslapan (véase
la Figura 41).

Figura 41: Union finita de regiones simples.

Versiones extendidas del teorema de Green
Ejemplo

22
Evalúe C y 2 dx + 3xydy, donde C es la frontera de la región semianular D
entre las circunferencias x2 + y 2 = 1 y x2 + y 2 = 4 en el semiplano superior.

Solución
Observe que aunque D no es simple, el eje y la divide en dos regiones simples
(véase la Figura 42). En coordenadas polares podemos escribir

D = { (r, θ)| 1 ≤ r ≤ 2, 0 ≤ θ ≤ π} .

Por tanto, el teorema de Green da

‰ ¨  
∂ ∂
y 2 dx + 3xydy = y 2 dA


(3xy) −
C ∂x ∂y
D
¨ ˆ π ˆ 2
= ydA = (r sen θ) r dr dθ
0 1
D
ˆ π ˆ 2  2
1 3 14
= sen θ dθ r2 dr = [− cos θ]π0 r = .
0 1 3 1 3
Versiones extendidas del teorema de Green

Figura 42: Región D entre dos semicircunferencias.

I El teorema de Green se puede generalizar para aplicarlo a regiones con agujeros,

es decir, regiones que no son simplemente conexas. Observe que la frontera C
de la región D de la Figura 43 consiste en dos curvas simples cerradas C1 y C2 .
Supongamos que estas curvas frontera están orientadas de tal modo que la
región D siempre está a la izquierda cuando se recorre la curva C. Por tanto, la
dirección positiva es contraria a la de las manecillas del reloj en el caso de la
curva exterior C1 , pero en el sentido de las manecillas del reloj en el caso de la
curva interior C2 . Si dividimos D en dos regiones D0 y D00 por medio de las
líneas mostradas en la Figura 44 y luego aplicamos el teorema de Green a D0 y
D00 , obtenemos
Versiones extendidas del teorema de Green
¨   ¨   ¨  
∂Q ∂P ∂Q ∂P ∂Q ∂P
− dA = − dA + − dA
∂x ∂y ∂x ∂y ∂x ∂y
D D0 D 00
‰ 
= P dx + Qdy + P dx + Qdy
C1 C2

Como las integrales de línea a lo largo de las rectas frontera comunes siguen
direcciones opuestas, se cancelan y entonces
¨   ‰  ˛
∂Q ∂P
− dA = P dx + Qdy + P dx + Qdy = P dx + Qdy
∂x ∂y C1 C2 C
D

lo cual es el teorema de Green para la región D.

Figura 43: Región D consistente de dos curvas simples cerradas C1 y C2 .

Versiones extendidas del teorema de Green

Figura 44: División de D en dos regiones D0 y D00 .

Ejemplo 23   ´ ´
Si F(x, y) = (−y i + x j) / x2 + y 2 , demuestre que C F · dr = 2π para toda trayectoria simple,
cerrada, orientada positivamente y que encierra el origen.

Solución
Como C es una trayectoria cerrada arbitraria que encierra el origen, es difícil calcular en forma directa la
integral dada. De modo que consideremos una circunferencia C 0 orientada en el sentido contrario al de
las manecillas del reloj con centro en el origen y radio a, donde a se escoge de tal manera que sea tan
pequeño que C 0 quede dentro de C (véase la Figura 45). Sea
D la región acotada por C y C .
0
Entonces su frontera orientada positivamente es C ∪ −C 0 y de este modo la versión general del
teorema de Green da
Versiones extendidas del teorema de Green

Figura 45: C es una trayectoria cerrada que encierra el origen.

Versiones extendidas del teorema de Green
‰ ‰ ¨  
∂Q ∂P
P dx + Qdy + P dx + Qdy = − dA
C −C 0 ∂x ∂y
D
¨ 
y 2 − x2 y 2 − x2

= 2 − 2 dA = 0.
(x2 + y 2 ) (x2 + y 2 )
D

Por tanto, ‰ ‰
P dx + Qdy = P dx + Qdy
C C0
es decir, ˆ ˆ
F · dr = F · dr.
C C0
Ahora calculamos con facilidad esta última integral usando la parametrización
definida por r(t) = a cos t i + a sen t j, 0 ≤ t ≤ 2π. Por tanto,
ˆ ˆ ˆ 2π
F · dr = F · dr = F (r(t)) · r0 (t)dt
C C0 0
ˆ2π ˆ 2π
(−a sen t) (−a sen t) + (a cos t) (a cos t)
= dt = dt = 2π.
0 a2 cos2 t + a2 sen2 t 0
Rotacional y divergencia
En esta sección se definen dos operaciones que se pueden ejecutar sobre los campos
vectoriales y que desempeñan un papel fundamental en las aplicaciones del cálculo
vectorial al flujo de fluidos y en electricidad y magnetismo. Cada operación es similar a
la derivación, pero una genera un campo vectorial mientras que la otra proporciona un
campo escalar.

Rotacional
Si F = P i + Qj + Rk es un campo vectorial sobre R3 y existen las derivadas parciales
de P , Q y R, entonces el rotacional de F es el campo vectorial sobre R3 definido por
     
∂R ∂Q ∂P ∂R ∂Q ∂P
rot F = − i+ − j+ − k (28)
∂y ∂z ∂z ∂x ∂x ∂y

Como un auxiliar nemotécnico, escribimos la Ecuación 28 usando la notación del

operador. Introducimos el operador diferencial vectorial ∇ («nabla») como

∂ ∂ ∂
∇=i +j +k .
∂x ∂y ∂z

Tiene significado cuando opera sobre una función escalar para producir el gradiente de
f.
∂f ∂f ∂f ∂f ∂f ∂f
∇f = i +j +k = i+ j+ k.
∂x ∂y ∂z ∂x ∂y ∂z
Rotacional y divergencia

Si imaginamos que ∇ es un vector con componentes ∂/∂x, ∂/∂y y ∂/∂z,

podemos considerar el producto cruz de ∇ y el campo vectorial F como sigue:
 
 i j k
 
 
∂ ∂ ∂ 
∇ × F =  ∂x ∂y ∂z 
 
P Q R
     
∂R ∂Q ∂P ∂R ∂Q ∂P
= − i+ − j+ − k.
∂y ∂z ∂z ∂x ∂x ∂y
Por tanto, la manera más sencilla de recordar la definición 28 es por medio de
la expresión simbólica
rot F = ∇ × F. (29)

Ejemplo 24
Si F(x, y, z) = xz i + xyz j − y 2 k, determine el rotacional de F.
Rotacional y divergencia
Solución
Utilizando la ecuación 29
 
 i j k 
 
 
∂ ∂ ∂ 
rot F = ∇ × F =  ∂x ∂y ∂z 
 xz xyz −y 2 
 
   
∂ ∂ ∂ ∂
−y 2 − −y 2 −



= (xyz) i − (xz) j
∂y ∂z ∂x ∂z
 
∂ ∂
+ (xyz) − (xz) k
∂x ∂y
= (−2y − xy) i − (0 − x) j + (yz − 0) k
= −y (2 + x) i + x j + yz k.

Recuerde que el gradiente de una función f de tres variables es un campo vectorial

sobre R3 y por eso es posible calcular su rotacional. El teorema siguiente establece que
el rotacional de un campo vectorial gradiente es 0.

Teorema 7
Si f es una función de tres variables que tiene derivadas parciales continuas de
segundo orden, entonces
rot (∇f ) = 0.
Rotacional y divergencia

Demostración.
Tenemos
 
 i j k
 
 
∂ ∂ ∂ 
rot (∇f ) = ∇ × (∇f ) =  ∂x ∂y ∂z 
 ∂f ∂f ∂f 
 ∂x ∂y ∂z 
 2 2
 2
∂2f
 2
∂2f
  
∂ f ∂ f ∂ f ∂ f
= − i+ − j+ − k
∂y∂z ∂z∂y ∂z∂x ∂x∂z ∂x∂y ∂y∂x
= 0i + 0j + 0k

de acuerdo con el teorema de Clairaut.

Puesto que un campo vectorial conservativo es uno para el cual F = ∇f , el

Teorema 7 se puede volver a enunciar como sigue:
Si F es conservativo, entonces rot F = 0.
Esto proporciona una forma de verificar que un campo vectorial no es
conservativo.
Rotacional y divergencia
Ejemplo 25
Demuestre que el campo vectorial F(x, y, z) = xz i + xyz j − y 2 k no es
conservativo.

Solución
En el Ejemplo 24 demostramos que

rot F = −y (2 + x) i + x j + yz k.

Esto demuestra que rot F 6= 0 y entonces según el Teorema 7, F es no

conservativo.
El inverso del Teorema 7 no es cierto en general, pero el siguiente teorema
establece que el inverso es válido si F está definido dondequiera. (Más
generalmente, es cierto si el dominio es simplemente conexo, es decir, “no hay
agujeros”.) El Teorema 8 es la versión tridimensional del Teorema 5. Su
demostración requiere del Teorema de Stokes».

Teorema 8
Si F es un campo vectorial definido en todo R3 cuyas funciones componentes
tienen derivadas parciales continuas y rot F = 0, entonces F es un campo
vectorial conservativo.
Rotacional y divergencia
Ejemplo 26

a) Demuestre que
F(x, y, z) = y 2 z 3 i + 2xyz 3 j + 3xy 2 z 2 k
es un campo vectorial conservativo.
b) Encuentre una función f tal que F = ∇f .

Solución

a) Calculamos el rotacional de F:
 
 i j k
 

 ∂ ∂ ∂

rot F = ∇ × F =  ∂x ∂y ∂z


 2 3
2xyz 3 2 2 
y z 3xy z 
= 6xyz 2 − 6xyz 2
i − 3y 2 z 2 − 3y 2 z 2 j + 2yz 3 − 2yz 3 k

= 0.

Puesto que rot F = 0 y el dominio de F es R3 , F es un campo vectorial

conservativo de acuerdo con el Teorema 8.
Rotacional y divergencia
b) La técnica para determinar f se trató en la sección «Teorema fundamental de
las integrales de línea». Tenemos

∂f (x, y, z)
= y2 z3 (30)
∂x

∂f (x, y, z)
= 2xyz 3 (31)
∂y
∂f (x, y, z)
= 3xy 2 z 2 . (32)
∂z
Al integrar (30) respecto a x obtenemos

f (x, y, z) = xy 2 z 3 + g (y, z) . (33)

Al derivar (33) respecto a y, obtenemos

∂f (x, y, z) ∂g(y, z)
= 2xyz 3 + ,
∂y ∂y

de modo que al comparar con (31) se obtiene ∂g(y, z)/∂y = 0. Por tanto,
g(y, z) = h(z) y de (33) obtenemos

f (x, y, z) = xy 2 z 3 + h(z). (34)

Rotacional y divergencia
Derivando parcialmente (34) respecto a z, tenemos

∂f (x, y, z)
= 3xy 2 z 2 + h0 (z) (35)
∂z

Comparando (35) con (32) obtenemos h0 (z) = 0. Por tanto h(z) = K y

f (x, y, z) = xy 2 z 3 + K.

I La razón del nombre «rotacional» es que el vector rotacional se relaciona con

rotaciones.
I Hay otra cuando F representa el campo de velocidades en el flujo de fluidos (ver
el Ejemplo 2).
I Las partículas cercanas a (x, y, z) en el fluido tienden a girar alrededor del
eje que señala la dirección del rotacional de F(x, y, z), y la longitud de
este vector rotacional es una medida de qué tan rápido se desplazan las
partículas alrededor del eje (véase la Figura 46).
I Si rot F = 0, en un punto P , entonces el fluido está libre de rotaciones en
P y F se llama irrotacional en P . En otras palabras, no hay remolinos en
P.
I Si rot F = 0, entonces una pequeñísima rueda de paletas se
movería con el fluido, pero no giraría respecto a su eje.
I Si rot F 6= 0, la rueda de paletas gira alrededor de su eje.
I Hay una explicación más amplia como consecuencia del teorema de Stokes.
Rotacional y divergencia

Figura 46: Rotacional de F.

Divergencia
Si F = P i + Qj + Rk es un campo vectorial sobre R3 y existen las derivadas parciales
∂P/∂x, ∂Q/∂y y ∂R/∂z, entonces la divergencia de F es la función de tres variables
definida por
∂P ∂Q ∂R
div F = + + . (36)
∂x ∂y ∂z

Observación
Observe que el rot F es un campo vectorial, pero div F es un campo escalar. En
términos del operador gradiente ∇ = (∂/∂x) i + (∂/∂y) j + (∂/∂z) k, la divergencia
de F se puede expresar simbólicamente como el producto punto de ∇ y F:

div F = ∇ · F (37)
Rotacional y divergencia

Ejemplo 27
Si F(x, y, z) = xz i + xyz j − y 2 k, encuentre div F.

Solución
De acuerdo con la definición de divergencia (ecuación 36 o 37) tenemos
∂ ∂ ∂
−y 2 = z + xz.


div F = ∇ · F = (xz) + (xyz) +
∂x ∂y ∂z

Si F es un campo vectorial sobre R3 , entonces rot F es también un campo

vectorial sobre R3 . Si es así, podemos calcular su divergencia. El teorema
siguiente demuestra que el resultado es 0.

Teorema 9
Si F = P i + Qj + Rk es un campo vectorial sobre R3 y P , Q y R tienen
derivadas parciales de segundo orden continuas, entonces

div rot F = 0.
Rotacional y divergencia
Demostración.
Utilizando las definiciones de divergencia y rotacional, tenemos

div rot F = ∇ · (∇ × F)
     
∂ ∂R ∂Q ∂ ∂P ∂R ∂ ∂Q ∂P
= − + − + −
∂x ∂y ∂z ∂y ∂z ∂x ∂z ∂x ∂y
∂2R ∂2Q ∂2P ∂2R ∂2Q ∂2P
= − + − + −
∂x∂y ∂x∂z ∂y∂z ∂y∂x ∂z∂x ∂z∂y
=0

porque los términos se anulan en pares según el teorema de Clairaut.

Observación
Observe la analogía con el producto triple escalar a · (a × b) = 0.

Ejemplo 28
Demuestre que el campo vectorial F(x, y, z) = xz i + xyz j − y 2 k no se puede
expresar como el rotacional de otro campo vectorial, es decir, F 6= rot G.
Rotacional y divergencia

Solución
En el Ejemplo 27, demostramos que

div F = z + xz

y, por tanto, div F 6= 0. Si fuera cierto que F = rot G, entonces el Teorema 9

daría
div F = div rot G = 0
lo cual contradice div F 6= 0. Por tanto, F no es el rotacional de otro campo
vectorial.
I Otra vez, la razón del nombre de divergencia se puede entender en el
contexto del flujo de fluidos. Si F(x, y, z) es la velocidad de un fluido (o
gas), entonces div F(x, y, z) representa la razón de cambio neta (respecto
al tiempo) de la masa del fluido (o gas) que fluye desde el punto (x, y, z)
por unidad de volumen.
I En otras palabras, div F(x, y, z) mide la tendencia del fluido a divergir del
punto (x, y, z).
I Si div F = 0, entonces se dice que F es incompresible.
Rotacional y divergencia
Se presenta otro operador diferencial cuando calculamos la divergencia de un
campo vectorial gradiente ∇f . Si f es una función de tres variables tenemos

∂2f ∂2f ∂2f

div (∇f ) = ∇ · (∇f ) = + +
∂x2 ∂y 2 ∂z 2

y esta expresión se presenta con tanta frecuencia que se abrevia como ∇2 f . El

operador
∇2 = ∇ · ∇
se llama operador de Laplace debido a su relación con la ecuación de Laplace

∂2f ∂2f ∂2f

∇2 f = 2
+ 2
+ = 0.
∂x ∂y ∂z 2

También podemos aplicar el operador de Laplace ∇2 a un campo vectorial

F = P i + Qj + Rk

en términos de sus componentes:

∇2 F = ∇2 P i + ∇2 Qj + ∇2 Rk.
Rotacional y divergencia
Formas vectoriales del teorema de Green
Los operadores rotacional y divergencia permiten volver a expresar el teorema de Green
en versiones que serán útiles en trabajos posteriores. Supongamos que la región plana
D, su curva frontera C y las funciones P y Q satisfacen las hipótesis del teorema de
Green. Entonces consideremos el campo vectorial F = P i + Qj. Su integral de línea es
˛ ˛
F · dr = P dx + Qdy
C C

y, en cuanto a F como un campo vectorial sobre R3 con la tercera componente igual a

0, tenemos  
 i j k  
 
 ∂ ∂ ∂  ∂Q ∂P
rot F =  = − k.

 ∂x ∂y ∂z  ∂x ∂y
 
P (x, y) Q(x, y) 0 

Por tanto,  
∂Q ∂P ∂Q ∂P
(rot F) · k = − k·k= −
∂x ∂y ∂x ∂y
y se puede volver a expresar la ecuación del teorema de Green en la forma vectorial
˛ ¨
F · dr = (rot F) · k dA. (38)
C
D
Rotacional y divergencia

I La Ecuación 38 expresa la integral de línea de la componente tangencial

de F a lo largo de C como la integral doble de la componente vertical de
rot F sobre la región D encerrada por C. En seguida deducimos una
fórmula similar que involucra la componente normal de F .
Si C está dada por la ecuación vectorial

r(t) = x(t)i + y(t)j, a≤t≤b

entonces el vector tangente unitario es

x0 (t) y 0 (t)
T(t) = 0
i+ 0 j.
|r (t)| |r (t)|
Podemos verificar que el vector normal unitario exterior a C está dado por
y 0 (t) x0 (t)
n(t) = 0
i− 0 j
|r (t)| |r (t)|
(Véase la Figura 47.)
Rotacional y divergencia

Figura 47: Región D encerrada por la curva C.

Rotacional y divergencia
Entonces, según la Ecuación 4
˛ ˆ b
(F · n) (t) r0 (t) dt
 
F · n ds =
C a
ˆb
P (x(t), y(t)) y 0 (t) Q (x(t), y(t)) x0 (t)  0 
 
= − r (t) dt
a |r0 (t)| |r0 (t)|
ˆb
= P (x(t), y(t)) y 0 (t)dt − Q (x(t), y(t)) x0 (t)dt
a
ˆ ¨  
∂P ∂Q
= P dy − Qdx = + dA
C ∂x ∂y
D

de acuerdo con el teorema de Green. Pero el integrando en esta integral doble

es justamente la divergencia de F. Entonces una segunda forma vectorial del
teorema de Green: ˛ ¨
F · n ds = div F(x, y)dA (39)
C
D

Esta versión establece que la integral de línea de la componente normal de F a

lo largo de C es igual a la integral doble de la divergencia de F sobre la región
D encerrada por C.
Superficies paramétricas y sus áreas
Hasta ahora se han considerado tipos especiales de superficies: cilindros, superficies
cuadráticas, gráficas de funciones de dos variables y superficies de nivel de funciones
de tres variables. Aquí se usan funciones vectoriales para describir superficies más
generales, llamadas superficies paramétricas, y se calculan sus áreas. A continuación
tomamos la fórmula general del área superficial y vemos cómo se aplica a superficies
especiales.

Superficies paramétricas
Casi de la misma manera como se describió una curva en el espacio mediante una
función vectorial r(t) de un solo parámetro t, podemos describir una superficie
mediante una función vectorial r(u, v) de dos parámetros u y v. Suponemos que

r(u, v) = x(u, v)i + y(u, v)j + z(u, v)k (40)

es una función con valor vectorial definida sobre una región D en el plano uv. De este
modo, x, y y z, que son las funciones componentes de r, son funciones de dos
variables u y v con dominio D. El conjunto de todos los puntos (x, y, z) en R3 tal que

x = x(u, v), y = y(u, v), z = z(u, v) (41)

y (u, v) varía en todo el dominio D, se llama superficie paramétrica S y las

Ecuaciones 41 se llaman ecuaciones paramétricas de S. Cada elección de u y v da un
punto sobre S; luego de efectuar todas las elecciones, obtenemos todo S. En otras
palabras, la superficie S es trazada por la punta del vector de posición r(u, v) cuando
(u, v) se desplaza por toda la región D (véase la Figura 48).
Superficies paramétricas y sus áreas

Figura 48: Superficie paramétrica.

Superficies paramétricas y sus áreas

Ejemplo 29
Identifique y trace la superficie con ecuación vectorial

r(u, v) = 2 cos u i + v j + 2 sen u k.

Solución
Las ecuaciones paramétricas de esta superficie son

x = 2 cos u, y = v, z = 2 sen u.

Así, para cualquier punto (x, y, z) sobre la superficie

x2 + z 2 = 4 cos2 u + 4 sen2 u = 4.

Esto significa que las secciones transversales verticales paralelas al plano xz (es
decir, con y constante) son todas circunferencias de radio 2. Puesto que y = v
y no hay restricción sobre v, la superficie es un cilindro circular de radio 2 y su
eje es el eje de las y (véase la Figura 49).
Superficies paramétricas y sus áreas

Figura 49: Cilindro circular de radio 2 y eje el eje y.

Superficies paramétricas y sus áreas

En el Ejemplo 29 no hubo restricciones sobre los parámetros u y v, y por eso

obtuvimos todo el cilindro. Por ejemplo, si restringimos u y v expresando el
dominio del parámetro como

0 ≤ u ≤ π/2, 0≤v≤3

entonces x ≥ 0, z ≥ 0, 0 ≤ y ≤ 3, y obtenemos el cuarto de cilindro cuyo largo

es 3 y que se ilustra en la Figura 50.

Figura 50: Cilindro circular de radio 2 y eje el eje y, con u y v acotados.

Superficies paramétricas y sus áreas
I Si una superficie paramétrica S está dada por una función vectorial r(u, v),
entonces hay dos familias útiles de curvas que están sobre S; una familia con u
constante y la otra con v constante. Estas familias corresponden a las rectas
vertical y horizontal en el plano uv. Si mantenemos u constante al hacer
u = u0 , entonces r (u0 , v) se convierte en una función vectorial del parámetro v
y define una curva C1 que está sobre S (véase la Figura 51).
I En forma similar, si mantenemos constante a v haciendo v = v0 , obtenemos una
curva C2 dada por r(u, v0 ) que está en S. Estas curvas se denominan curvas
reticulares. (En el Ejemplo 29, las curvas reticulares obtenidas al hacer u
constante son rectas horizontales, en tanto que las curvas reticulares cuando v
es constante son circunferencias.)

Figura 51: Dos familias de curvas.

Superficies paramétricas y sus áreas

De hecho, cuando una computadora grafica una superficie paramétrica, por lo regular delinea una
superficie trazando estas curvas reticulares, como vemos en el ejemplo siguiente.

Ejemplo 30
Utilice un sistema algebraico computarizado para graficar la superficie

r(u, v) = ((2 + sen v) cos u, (2 + sen v) sen u, u + cos v)

¿Qué curvas tienen a u como constante? ¿Cuáles tienen a v como constante?

Solución
Graficamos la parte de la superficie correspondiente al dominio de los parámetros 0 ≤ u ≤ 4π,
0 ≤ v ≤ 2π, obteniéndose la Figura 52. Tiene el aspecto de un tubo en espiral. Para identificar las
curvas reticulares, escribimos las ecuaciones paramétricas correspondientes:

x = (2 + sen v) cos u, y = (2 + sen v) sen u, z = u + cos v.

Si v es constante, entonces sen v y cos v son constantes, por lo que las ecuaciones paramétricas son
similares a la de una hélice. Así, las curvas reticulares con v constante son las curvas en espiral de la
Figura 52. Se infiere que las curvas reticulares con u constante deben ser las curvas que se ven como
circunferencias en la Figura 52. Más evidencias para esta afirmación es que si u se mantiene constante,
u = u0 , entonces la ecuación z = u0 + cos v muestra que los valores de z varían de u0 − 1 a u0 + 1.
Superficies paramétricas y sus áreas

Figura 52: Gráfica de r(u, v) = ((2 + sen v) cos u, (2 + sen v) sen u, u + cos v) .
Superficies paramétricas y sus áreas
En los Ejemplos 29 y 30 se proporcionaba una ecuación vectorial y se pedía dibujar la superficie
paramétrica correspondiente. En los ejemplos siguientes se plantea un problema más difícil: hallar una
función vectorial que represente una superficie dada. En lo que resta del capítulo, a menudo se necesita
hacer exactamente eso.

Ejemplo 31
Formule una función vectorial que represente el plano que pasa por el punto P0 y cuyo vector de
posición es r0 y contiene dos vectores no paralelos a y b.

Solución
Si P es cualquier punto en el plano, podemos obtener desde P0 a P desplazándose una cierta distancia
en la dirección de a y otra distancia en la dirección de b. Así que hay escalares u y v tales que
−−→
P0 P = ua + vb. En la Figura 53 se ilustra cómo funciona esto mediante la ley del paralelogramo en el
caso donde u y v son positivas. Si r es el vector de posición de P , entonces

−−→ −−→
r = OP0 + P0 P = r0 + ua + vb.

Entonces, la ecuación vectorial del plano se puede expresar como

r(u, v) = r0 + ua + vb

donde u y v son números reales.

Si escribimos r = (x, y, z), r0 = (x0 , y0 , z0 ), a = (a1 , a2 , a3 ) y b = (b1 , b2 , b3 ), entonces podemos
expresar las ecuaciones paramétricas del plano que pasa por el punto (x0 , y0 , z0 ) como sigue:

x = x0 + ua1 + vb1 , y = y0 + ua2 + vb2 , z = z0 + ua3 + vb3 .

Superficies paramétricas y sus áreas

Figura 53: Plano que pasa por P0 .

Ejemplo 32
Determine una representación paramétrica de la esfera

x2 + y 2 + z 2 = a2 .

Solución
La esfera tiene una representación simple ρ = a en coordenadas esféricas, así que
escogemos los ángulos φ y θ en coordenadas esféricas como los parámetros. Luego, al
hacer ρ = a en las ecuaciones para la conversión de coordenadas esféricas a
rectangulares, obtenemos
Superficies paramétricas y sus áreas

x = a sen φ cos θ, y = a sen φ sen θ, z = a cos φ

como las ecuaciones paramétricas de la esfera. La ecuación vectorial
correspondiente es

r(φ, θ) = a sen φ cos θ i + a sen φ sen θ j + a cos φ k.

Tenemos 0 ≤ φ ≤ π y 0 ≤ θ ≤ 2π, de modo que el dominio de los parámetros

es el rectángulo D = [0, π] × [0, 2π]. Las curvas reticulares con φ constante son
las circunferencias de latitud constante (sin olvidar el ecuador). Las curvas
reticulares con θ constante son los meridianos (semicircunferencias), los cuales
conectan los polos norte y sur. (Ver la Figura 54).
Nota: Vimos en el Ejemplo 32 que las curvas reticulares para una
esfera son curvas de latitud y longitud constantes. Para una
superficie paramétrica general, realmente estamos haciendo un
mapa y las curvas reticulares son similares a las líneas de
longitud y latitud. Describimos un punto sobre una superficie
paramétrica (como la de la Figura 52) dando valores específicos
a u y v, como si se dieran la latitud y la longitud de un punto.
Superficies paramétricas y sus áreas

Figura 54: Curvas reticulares de la esfera x2 + y2 + z 2 = a2 .

Superficies paramétricas y sus áreas
Ejemplo 33
Encuentre la representación paramétrica del cilindro

2 2
x +y = 4, 0 ≤ z ≤ 1.

Solución
El cilindro tiene una representación simple r = 2 en coordenadas cilíndricas, de modo que elegimos como
parámetros a θ y z en coordenadas cilíndricas. Entonces las ecuaciones paramétricas del cilindro son

x = 2 cos θ, y = 2 sen θ, z=z

donde 0 ≤ θ ≤ 2π y 0 ≤ z ≤ 1.

Ejemplo 34
Encuentre una función vectorial que representa el paraboloide elíptico z = x2 + 2y 2 .

Solución
Si consideramos a x y y como parámetros, entonces las ecuaciones paramétricas son simplemente

2 2
x = x, y = y, z = x + 2y

y la ecuación vectorial es  
2 2
r(x, y) = x i + y j + x + 2y k.
Superficies paramétricas y sus áreas
En general, si se da una superficie como la gráfica de una función de x y de y,
es decir, con una ecuación de la forma z = f (x, y), siempre se le puede
considerar como una superficie paramétrica tomando a x y y como parámetros
y expresando las ecuaciones paramétricas como

x = x, y = y, z = f (x, y).

Las representaciones paramétricas (también conocidas como parametrizaciones)

de superficies no son únicas. El ejemplo siguiente muestra dos maneras de
parametrizar un cono.

Ejemplo 35 p
Encuentre una representación paramétrica para la superficie z = 2 x2 + y 2 , es
2 2 2
decir, la mitad superior del cono z = 4x + 4y .

Solución
Una posible representación se consigue escogiendo a x y y como parámetros:
p
x = x, y = y, z = 2 x2 + y 2 .

De modo que la ecuación vectorial es

p
r(x, y) = x i + y j + 2 x2 + y 2 k.
Superficies paramétricas y sus áreas
Solución (Segundo método)
Otra representación resulta de escoger como parámetros a las coordenadas
polares
pr y θ. Un punto (x, y, z) del cono satisface x = r cos θ, y = r sen θ, y
z = 2 x2 + y 2 = 2r. De modo que una ecuación vectorial para el cono es

r(r, θ) = r cos θ i + r sen θ j + 2r k

donde r ≥ 0 y 0 ≤ θ ≤ 2π.

Superficies de revolución
Las superficies de revolución se pueden representar en forma paramétrica y, por
tanto, se pueden graficar mediante una computadora. Por ejemplo,
consideremos la superficie S que se obtiene al hacer girar la curva y = f (x),
a ≤ x ≤ b, alrededor del eje x, donde f (x) ≥ 0. Sea θ el ángulo de rotación
como se muestra en la Figura 55. Si (x, y, z) es un punto sobre S, entonces

x = x, y = f (x) cos θ, z = f (x) sen θ. (42)

Por tanto, tomamos a x y θ como parámetros y consideramos las

Ecuaciones 42 como ecuaciones paramétricas de S. El dominio de los
parámetros están dados por a ≤ x ≤ b, 0 ≤ θ ≤ 2π.
Superficies paramétricas y sus áreas

Figura 55: Superficie de revolución.

Ejemplo 36
Determine las ecuaciones paramétricas de la superficie generada al hacer girar
la curva y = sen x, 0 ≤ x ≤ 2π, alrededor del eje x. Con estas ecuaciones,
grafique la superficie de revolución.
Superficies paramétricas y sus áreas

Según las Ecuaciones 42, las ecuaciones paramétricas son

x = x, y = sen x cos θ, z = sen x sen θ

y el dominio del parámetro es 0 ≤ x ≤ 2π, 0 ≤ θ ≤ 2π. Utilizando una

computadora para graficar estas ecuaciones y girando la imagen, obtenemos la
gráfica de la Figura 56.

Figura 56: Superficie de revolución de f (x) = sen x.

I Podemos adaptar las Ecuaciones 42 para representar una superficie

obtenida a partir de la revolución alrededor del eje y o del eje z.
Superficies paramétricas y sus áreas

Este es el código en SageMath para graficar superficies paramétricas

(Ejemplo 30):
from sage.plot.plot3d.parametric_surface import ParametricSurface;
def f(u,v): return (2+sin(v))*cos(u),(2+sin(v))*sin(u),u+cos(v);
S =
ParametricSurface(f,(srange(float(0),float(4*pi),float(0.1)),srange(float(0),float(2*pi),float(0.1))),
mesh=True,opacity=0.6)
from sage.manifolds.utilities import set_axes_labels
graf = set_axes_labels(S,’x’,’y’,’z’,color=’red’)
graf
Este es el código en SageMath para graficar superficies de revolución
(Ejemplo 36):
from sage.plot.plot3d.parametric_surface import ParametricSurface
def f(x,theta): return x,sin(x)*cos(theta),sin(x)*sin(theta)
S =
ParametricSurface(f,(srange(float(0),float(2*pi),float(0.1)),srange(float(0),float(2*pi),float(0.1))),
mesh = True)
show(S)
Superficies paramétricas y sus áreas
Plano tangente
Ahora encontraremos el plano tangente a una superficie paramétrica S trazada por una función vectorial

r(u, v) = x (u, v) i + y (u, v) j + z (u, v) k

en un punto P0 con vector de posición r(u0 , v0 ). Si mantenemos a u constante haciendo u = u0 ,

entonces r(u0 , v) se transforma en una función vectorial del parámetro único v, y define una curva
reticular C1 que está sobre S (véase la figura 12). El vector tangente a C1 en P0 se obtiene
determinando la derivada parcial de r respecto a v:

∂ ∂x  ∂y  ∂z 
  
rv = (r(u, v)) = u=u0 i +

u=u0 j +
  k (43)
∂v ∂v v=v ∂v v=v ∂v u=u 0

0 0 v=v 0

Figura 57: Vectores tangente a la superficie en el punto P0 .

Superficies paramétricas y sus áreas
De manera similar, si mantenemos a v constante haciendo v = v0 , obtenemos
una curva reticular C2 dada por r(u, v0 ) que queda sobre S, y su vector
tangente en P0 es
  
∂ ∂x  ∂y  ∂z 
ru = (r(u, v)) = i + j + k. (44)
∂u ∂u u=u0 ∂u u=u0 ∂u u=u0
v=v0 v=v0 v=v0

Si ru × rv 6= 0, entonces la superficie S se llama suave (no tiene “esquinas”).

En el caso de una superficie suave, el plano tangente es el plano que contiene
los vectores tangentes ru y rv , y el vector ru × rv , es un vector normal al plano
tangente.

Ejemplo 37
Encuentre el plano tangente a la superficie cuyas ecuaciones paramétricas son
x = u2 , y = v 2 , z = u + 2v en el punto (1, 1, 3).

Solución
Primero calculamos los vectores tangentes:
∂x ∂y ∂z
ru = i+ j+ k = 2u i + k
∂u ∂u ∂u
∂x ∂y ∂z
rv = i+ j+ k = 2v j + 2k
∂v ∂v ∂v
Superficies paramétricas y sus áreas

Por tanto, un vector normal al plano tangente es

 
i j k
 
 
ru × rv = 2u 0 1  = −2v i − 4u j + 4uv k.
 
 0 2v 2 

Observe que el punto (1, 1, 3) corresponde a los valores del parámetro u = 1 y

v = 1, así que el vector normal es

−2 i − 4 j + 4 k.

Por tanto, una ecuación del plano tangente en (1, 1, 3) es

−2 (x − 1) − 4 (y − 1) + 4 (z − 3) = 0

o
x + 2y − 2z + 3 = 0.
Superficies paramétricas y sus áreas
Área de una superficie
Ahora definimos el área de una superficie paramétrica general dada por la
Ecuación 40. Para simplificar el trabajo, primero consideramos una superficie
cuyo dominio D de los parámetros es un rectángulo, y se subdivide en más
rectángulos Rij . Escogemos u∗i , vj∗ como el vértice inferior izquierdo de Rij

(véase la Figura 58).

Figura 58: La imagen del subrectángulo Rij es el parche Sij .

Superficies paramétricas y sus áreas
La parte Sij de la superficie S que corresponde
a Rij se denomina parche y el
punto Pij cuyo vector de posición es r u∗i , vj∗ , es una de sus esquinas. Sean

r∗u = ru u∗i , vj∗ y r∗v = rv u∗i , vj∗

los vectores tangente en Pij como lo definen las Ecuaciones 44 y 43. La

Figura 59a) muestra la manera como las dos orillas del parche que se
encuentran en Pij se pueden aproximar mediante vectores. A su vez, estos
vectores se pueden aproximar por medio de los vectores ∆u r∗u y ∆v r∗v porque
las derivadas parciales se pueden aproximar por cocientes de diferencias. De
este modo, aproximamos Sij , por medio del paralelogramo definido por los
vectores ∆u r∗u y ∆v r∗v . Este paralelogramo se ilustra en la Figura 59b) y se
ubica en el plano tangente a S en Pij . El área de este paralelogramo es

|(∆u r∗u ) × (∆v r∗v )| = |r∗u × r∗v | ∆u ∆v

y de este modo una aproximación del área de S es

m X
X n
|r∗u × r∗v | ∆u ∆v.
i=1 j=1
Superficies paramétricas y sus áreas

Figura 59: Aproximación de un parche mediante un paralelogramo.

Superficies paramétricas y sus áreas
La intuición dice que esta aproximación es mejor cuando incrementamos el
número de subrectángulos, e identificamos
˜ la doble suma como una suma de
Riemann para la integral doble D |r∗u × r∗v | du dv. Esto da lugar a la definición
siguiente.

Definición 5
Si una superficie paramétrica suave S está dada por la ecuación

r(u, v) = x (u, v) i + y (u, v) j + z (u, v) k, (u, v) ∈ D

y S es cubierta sólo una vez cuando (u, v) varía en todo el dominio del
parámetro D, entonces el área de la superficie de S es
¨
A (S) = |ru × rv | dA
D

donde
∂x ∂y ∂z
ru = i+ j+ k
∂u ∂u ∂u
∂x ∂y ∂z
rv = i+ j+ k.
∂v ∂v ∂v
Superficies paramétricas y sus áreas
Ejemplo 38
Determine el área de la superficie de una esfera de radio a.

Solución
En el Ejemplo 32, encontramos que la representación paramétrica

x = a sen φ cos θ, y = a sen φ sen θ, z = a cos φ

donde el dominio del parámetro es

D = { (φ, θ)| 0 ≤ φ ≤ π, 0 ≤ θ ≤ 2π} .

Primero calculamos el producto cruz de los vectores tangentes:

 
 i j k 

 ∂x ∂y ∂z 
 
rφ × rθ =  ∂φ ∂φ ∂φ 

 ∂x ∂y ∂z 

 
∂θ ∂θ ∂θ
 

 i j k 

=  a cos φ cos θ a cos φ sen θ −a sen φ
 
 
−a sen φ sen θ a sen φ cos θ 0 
2 2 2 2 2
= a sen φ cos θ i + a sen φ sen θ j + a sen φ cos φ k.
Superficies paramétricas y sus áreas

Por consiguiente,
p
|rφ × rθ | = a4 sen4 φ cos2 θ + a4 sen4 φ sen2 θ + a4 sen2 φ cos2 φ
p p
= a4 sen4 φ + a4 sen2 φ cos2 φ = a2 sen2 φ = a2 sen φ.

Puesto que sen φ ≥ 0 para 0 ≤ φ ≤ π. Por tanto, según la Definición 5,

el área de la esfera es
¨ ˆ 2π ˆ π
A= |rφ × rθ | dA = a2 sen φ dφ dθ
0 0
D
ˆ 2π ˆ π
= a2 dθ sen φ dφ = a2 (2π) 2 = 4πa2 .
0 0
Superficies paramétricas y sus áreas
Área de la superficie de la gráfica de una función
Para el caso especial de una superficie S cuya ecuación es z = f (x, y), donde (x, y)
está en D y f tiene derivadas parciales continuas, tomamos a x y y como parámetros.
Las ecuaciones paramétricas son

x = x, y = y, z = f (x, y)

de modo que    
∂f ∂f
rx = i + k, ry = j + k
∂x ∂y
y  
i j k 

 ∂f  ∂f ∂f
rx × ry = 1 0

∂x  = − ∂x i − ∂y j + k. (45)


0 ∂f 
1
 ∂y 
Por tanto, tenemos que
s 2  2 s  2  2
∂f ∂f ∂z ∂z
|rx × ry | = + +1= 1+ + (46)
∂x ∂y ∂x ∂y

y la fórmula del área superficial de la Definición 5 se transforma en

Superficies paramétricas ¨yssus áreas
 
∂z 2 ∂z
2
A (S) = 1+ + dA (47)
∂x ∂y
D

Ejemplo 39
Encuentre el área de la parte del paraboloide z = x2 + y 2 que se ubica bajo el plano z = 9.

Solución
El plano corta el paraboloide y forma la circunferencia x2 + y 2 = 9, z = 9. Por tanto, la superficie dada
está arriba del disco D con centro en el origen y radio 3. (Véase Figura 60.) Al aplicar la Fórmula 47,
tenemos

¨ s 2 2
∂z ∂z
 
A= 1+ + dA
∂x ∂y
D
¨ q
= 1 + (2x)2 + (2y)2 dA
D
¨ q
1 + 4 x2 + y 2 dA.


=
D

Transformando a coordenadas polares

ˆ 2π ˆ 3 p ˆ 2π ˆ 3 p
A= 1 + 4r 2 r dr dθ = dθ r 1 + 4r 2 dr
0 0 0 0
3
1 2  π  √
 
2 3/2 
 
= 2π 1 + 4r  = 37 37 − 1 .
8 3 0 6
Superficies paramétricas y sus áreas

Figura 60: Paraboloide z = x2 + y2 .

Superficies paramétricas y sus áreas
La pregunta que resta es si nuestra definición de área de una superficie (5) va
de acuerdo con la fórmula del área de una superficie a partir del cálculo de una
sola variable.
Consideramos la superficie S que se obtiene al hacer girar la curva y = f (x),
a ≤ x ≤ b, alrededor del eje de las x, donde f (x) ≥ 0, y f 0 es continua. De
acuerdo con las Ecuaciones 42 sabemos que las ecuaciones paramétricas de S
son

x = x, y = f (x) cos θ, z = f (x) sen θ, a ≤ x ≤ b, 0 ≤ θ ≤ 2π.

Para calcular el área de la superficie de S, necesitamos los vectores tangentes

rx = i + f 0 (x) cos θ j + f 0 (x) sen θ k

rθ = −f (x) sen θ j + f (x) cos θ k

Por tanto,
 
i j k
 

 
rx × rθ = 1 f 0 (x) cos θ 0
f (x) sen θ
−f (x) sen θ
 
0 f (x) cos θ 
= f (x)f 0 (x) i − f (x) cos θ j − f (x) sen θ k
Superficies paramétricas y sus áreas
y entonces,
q
2 2 2 2
|rx × rθ | = [f (x)] [f 0 (x)] + [f (x)] cos2 θ + [f (x)] sen2 θ
r h i q
2 2 2
0
= [f (x)] 1 + [f (x)] = f (x) 1 + [f 0 (x)]

porque f (x) ≥ 0. Por tanto, el área de S es

¨
A= |rx × rθ | dA
D
ˆ 2π ˆ b q
2
= f (x) 1 + [f 0 (x)] dx dθ
0 a
ˆ b q
2
= 2π f (x) 1 + [f 0 (x)] dx.
a

Ésta es precisamente la fórmula que se utilizó para definir el área de una

superficie de revolución en el cálculo de una sola variable.
Integrales de superficie
La relación que existe entre las integrales de superficie y el área de una
superficie es la misma que la relación entre integrales de línea y longitud de
arco. Supongamos que f es una función de tres variables en cuyo dominio se
encuentra la superficie S. Definiremos la integral de superficie de f sobre S en
tal forma que, en el caso donde f (x, y, z) = 1, el valor de la integral de
superficie es igual al área superficial de S. Comenzamos con superficies
paramétricas y luego trabajamos con el caso especial donde S es la gráfica de
una función de dos variables.

Superficies paramétricas
Supongamos que una superficie S tiene una ecuación vectorial

r(u, v) = x (u, v) i + y (u, v) j + z (u, v) k, (u, v) ∈ D.

Primero supongamos que el dominio D del parámetro es un rectángulo y lo

dividimos en subrectángulos Rij de dimensiones ∆u y ∆v. Entonces la
superficie S se divide en los parches correspondientes Sij como en la Figura 61.
Evaluamos f en un punto Pij∗ en cada parche, multiplicamos por el área ∆Sij
del parche, y formamos la suma de Riemann.
m X
X n
f Pij∗ ∆Sij .


i=1 j=1
Integrales de superficie

Figura 61: Subdivisión de S en parches Sij .

Integrales de superficie
A continuación tomamos el límite haciendo que el número de parches se
incremente y definimos la integral de superficie de f sobre la superficie S como
¨ X m X n
f Pij∗ ∆Sij .


f (x, y, z)dS = lı́m (48)
m,n→∞
S i=1 j=1

Observe la analogía con la definición de una integral de línea y también la

analogía con la definición de una integral doble.
Para evaluar la integral de superficie en la Ecuación 48, aproximamos el área de
parche ∆Sij mediante el área de un paralelogramo de aproximación en el plano
tangente. En el análisis del área de una superficie de la sección «Superficies
paramétricas y sus áreas» hicimos la aproximación

∆Sij ≈ |ru × rv | ∆u ∆v

donde
∂x ∂y ∂z
ru = i+ j+ k
∂u ∂u ∂u
∂x ∂y ∂z
rv = i+ j+ k.
∂v ∂v ∂v
son los vectores tangente en la esquina de Sij .
Integrales de superficie
Si las componentes son continuas y ru y rv no son cero y no son paralelos en el
interior de D, se puede demostrar por la Definición 48, incluso cuando D no es
un rectángulo, que
¨ ¨
f (x, y, z) dS = f (r(u, v)) |ru × rv | dA. (49)
S D

Debemos comparar con la fórmula para una integral de línea:

ˆ ˆ b
f (r(t)) r0 (t) dt.
 
f (x, y, z) =
C a

Observe que ¨ ¨
1 dS = |ru × rv | dA = A(S).
S D

La Fórmula 49 permite evaluar una integral de superficie convirtiéndola en un

integral doble sobre el dominio D del parámetro. Al usar esta fórmula, recuerde
que f (r(u, v)) se evalúa escribiendo x = x(u, v), y = y(u, v) y z = z(u, v) en
la fórmula para f (x, y, z).
Integrales de superficie
Ejemplo 40
˜
Calcule la integral de superficie S x2 dS, donde S es esfera unitaria x2 + y 2 + z 2 = 1.

Solución
Como en el Ejemplo 32, usamos la representación paramétrica

x = sen φ cos θ, y = sen φ sen θ, z = cos φ, 0 ≤ φ ≤ π, 0 ≤ θ ≤ 2π

es decir,
r(φ, θ) = sen φ cos θ i + sen φ sen θ j + cos φ k.

Al igual que en el Ejemplo 38, se calcula que

 
rφ × rθ  = sen φ.

Por tanto, de acuerdo con la Fórmula 49,

¨ ¨
x2 dS = (sen φ cos θ)2 rφ × rθ  dA
 

S D
ˆ 2π ˆ π ˆ 2π ˆ π
= sen2 φ cos2 θ sen φ dφ dθ = cos2 θ dθ sen3 φ dφ.
0 0 0 0
Integrales de superficie
Utilizando las identidades cos2 θ = 12 (1 + cos 2θ) y sen2 φ = 1 − cos2 φ,
tenemos
¨ ˆ 2π ˆ π
x2 dS = 1
sen φ − sen φ cos2 φ dφ


2
(1 + cos 2θ) dθ
0 0
S
2π  π 4π
1 1 1
cos3 φ

= 2
θ+ 2
sen 2θ 0
− cos φ + 3 0
= .
3
Las integrales de superficie tienen aplicaciones parecidas a las de las integrales
que ya tratamos. Por ejemplo, si una hoja delgada de aluminio tiene la forma
de una superficie S y la densidad, masa por unidad de área, en el punto
(x, y, z) es ρ(x, y, z), entonces la masa total de la lámina es
¨
m= ρ(x, y, z)dS
S

y el centro de masa es (x̄, ȳ, z̄), donde

¨ ¨ ¨
1 1 1
x̄ = xρ(x, y, z)dS, ȳ = yρ(x, y, z)dS, z̄ = zρ(x, y, z)dS
m m m
S S S

Los momentos de inercia también se pueden definir como antes.

Integrales de superficie
A continuación presentamos el código en SageMath para determinar la masa, el centroide, momentos
con respecto a los planos xy, yz, xz, centro de masa, momentos de inercia, radios de giro de una
superficie S con densidad variable.

densidad(x,y,z)=x^2+z
r=vector([3*cos(t), 3*sin(t),4*t])
t_inicial=0
t_final=2*pi
t_rango=(t,t_inicial,t_final)
def integral_linea(integrando): return RR(numerical_integral((integrando).subs(x=r[0],
y=r[1],z=r[2]), t_inicial, t_final)[0])
_ = var(’x,y,z,t’)
r_prima = diff(r,t)
ds=diff(r,t).norm()
s=integral_linea(ds)
centroide_x=integral_linea(x*ds)/s
centroide_y=integral_linea(y*ds)/s
centroide_z=integral_linea(z*ds)/s
dm=densidad(x,y,z)*ds
m=integral_linea(dm) # masa
avg_densidad = m/s # Densidad promedio
momento_alrededor_plano_yz=integral_linea(x*dm)
momento_alrededor_plano_xz=integral_linea(y*dm)
momento_alrededor_plano_xy=integral_linea(z*dm)
centro_masa_x = momento_alrededor_plano_yz/m
centro_masa_y = momento_alrededor_plano_xz/m
centro_masa_z = momento_alrededor_plano_xy/m
Ix=integral_linea((y^2+z^2)*dm) # Momento de inercia respecto a x
Iy=integral_linea((x^2+z^2)*dm) # Momento de inercia respecto a y
Iz=integral_linea((x^2+y^2)*dm) # Momento de inercia respecto a z
Rx = sqrt(Ix/m) # Radio de giro de masa x
Ry = sqrt(Iy/m) # Radio de giro de masa y
Rz = sqrt(Iz/m) # Radio de giro de masa z
Integrales de superficie
Gráficas

Cualquier superficie S con ecuación z = g(x, y) se puede considerar como una

superficie paramétrica con ecuaciones paramétricas

x = x, y = y, z = g(x, y)

y así tenemos    
∂g ∂g
rx = i + k, ry = j + k.
∂x ∂y
Por tanto,
∂g ∂g
rx × ry = − i− j+k (50)
∂x ∂y
y s 2  2
∂z ∂z
|rx × ry | = + +1
∂x ∂y
Por tanto, en este caso, la Fórmula 49 se convierte en
¨ ¨
s 
2  2
∂z ∂z
f (x, y, z)dS = f (x, y, g(x, y)) + + 1dA (51)
∂x ∂y
S D
Integrales de superficie
Se aplican fórmulas similares cuando es más conveniente proyectar S en el plano yz o en el xz. Por
ejemplo, si S es una superficie cuya ecuación es y = h(x, z) y D es su proyección sobre el plano xz,
entonces
¨ ¨
s
∂y 2 ∂y 2
  
f (x, y, z)dS = f (x, h(x, z), z) + + 1dA
∂x ∂z
S D

Ejemplo 41
˜
Evalúe S y dS, donde S es la superficie z = x + y 2 , 0 ≤ x ≤ 1, 0 ≤ y ≤ 2. (Ver la Figura 62).

Figura 62: Superficie z = x + y2 .

Solución
Puesto que
∂z ∂z
=1 y = 2y
∂x ∂y
Integrales de superficie
la Fórmula 51 da

¨ ¨ s 2 2
∂z ∂z

ydS = y + + 1dA
∂x ∂y
S D
ˆ 1ˆ 2 q
= y 1 + 4y 2 + 1dydx
0 0
ˆ 1 √ ˆ 2 q
= dx 2 y 1 + 2y 2 dy
0 0
√
√ 1 2  2 13 2
  
2 3/2 

= 2 1 + 2y = .
4 3 3

0

El código en SageMath para calcular esta integral de superficie es:

var(’x,y’)
f(x,y) = y
g(x,y) = x + y^2
gx = diff(g,x)
gy = diff(g,y)
x_inicial = 0
x_final = 1
y_inicial = 0
y_final = 2
dS = sqrt(1+gx^2+gy^2)
def integral_superficie(integrando): return
RR(integral(integral(integrando,y,y_inicial,y_final),x,x_inicial,x_final))
integral_superficie(f*dS)
Integrales de superficie
Si S es una superficie suave por tramos, es decir, una unión finita de superficies suaves
S1 , S2 , . . . , Sn que se intersectan sólo a lo largo de sus fronteras, entonces la integral
de superficie de f sobre S se define mediante
¨ ¨ ¨
f (x, y, z)dS = f (x, y, z)dS + · · · + f (x, y, z)dS.
S S1 Sn

Ejemplo
˜
42
Evalúe S z dS, donde S es la superficie cuyos lados S1 los define el cilindro
x2 + y 2 = 1, cuyo fondo S2 es el disco x2 + y 2 ≤ 1 en el plano z = 0, y cuya tapa S3
es la parte del plano z = 1 + x que queda arriba de S2 .

Solución
La superficie S se ilustra en la Figura 63. (Cambiemos la posición usual de los ejes
para tener mejor visión en S). Para S1 usamos θ y z como parámetros (véase el
Ejemplo 33 de la sección «Superficies paramétricas y sus áreas») y expresamos las
ecuaciones paramétricas como

x = cos θ, y = sen θ, z=z

donde
0 ≤ θ ≤ 2π y 0 ≤ z ≤ 1 + x = 1 + cos θ.
Integrales de superficie
Por tanto,  
 i j k
 
 
rθ × rz = − sen θ cos θ 0  = cos θ i + sen θ j
 
 0 0 1
y p
|rθ × rz | = cos2 θ + sen2 θ = 1.

Figura 63: Superficie con lado S1 , fondo S2 y tapa S3 .

Integrales de superficie
Por consiguiente, la integral de superficie sobre S1 es
¨ ¨
zdS = z |rθ × rz | dA
S1 D
ˆ 2π ˆ 1+cos θ ˆ 2π
= z dz dθ = 1
2
(1 + cos θ)2 dθ
0 0 0
ˆ 2π 
1
1 + 2 cos θ + cos2 θ dθ

= 2
0
ˆ 2π
1 1
 
= 2
1 + 2 cos θ + 2
(1 + cos 2θ) dθ
0
1
3 1
2π 3π
= 2 2
θ + 2 sen θ + 4
sen 2θ 0
= .
2
Puesto que S2 queda en el plano z = 0, tenemos
¨ ¨
z dS = 0 dS = 0.
S2 S2
Integrales de superficie
La superficie de la parte superior S3 está sobre el disco unitario D y es parte del plano z = 1 + x. De
este modo, si tomamos g(x, y) = 1 + x en la Fórmula 51 y la convertimos a coordenadas polares,
tenemos

¨ ¨ s 2 2
∂z ∂z
 
z dS = (1 + x) 1+ + dA
∂x ∂y
S3 D
ˆ 2π ˆ 1
√
= (1 + r cos θ) 1 + 1 + 0 r dr dθ
0 0
√ ˆ 2π ˆ 1
2

= 2 r + r cos θ dr dθ
0 0
√ ˆ 2π  
= 2 1 + 1 cos θ dθ
2 3
0
√ θ sen θ 2π √
 
= 2 + = 2π
2 3 0

Por tanto,

¨ ¨ ¨ ¨
z dS = zdS + zdS + zdS
S S1 S2 S3

3π √ 3 √
 
= +0+ 2π = + 2 π.
2 2
Integrales de superficie
Superficies orientadas
I Para definir integrales de superficie de campos vectoriales, debemos
descartar superficies no orientables como la cinta de Möbius que se
muestra en la Figura 64. Se le dio ese nombre en honor al geómetra
alemán August Möbius (1790-1868). Usted mismo puede construirla
tomando una larga tira de papel, darle media vuelta y pegar los extremos
como se indica en la Figura 65. Si una hormiga caminara por la banda de
Möbius empezando en el punto P , finalizaría su recorrido en el otro lado
de la cinta (es decir, con su lado superior apuntando en la dirección
opuesta). Entonces si la hormiga continúa su recorrido en la misma
dirección, terminaría de nuevo en el punto P , pero arriba de él, sin haber
brincado al otro lado. (Si usted ya tiene su cinta de Möbius, dibuje una
línea por todo el centro de la cinta.) Entonces, una cinta de Möbius tiene
sólo un lado. Podemos graficar la cinta de Möbius utilizando sus
ecuaciones paramétricas.
I De aquí en adelante, sólo se consideran superficies susceptibles de ser
orientadas, es decir, que tengan dos lados. Iniciamos con una superficie S
que tiene un plano tangente en cada punto (x, y, z) sobre S (excepto en
cualquier punto de la frontera). Hay dos vectores unitarios normales n1 y
n2 = −n1 en (x, y, z) (véase la Figura 66).
Superficies orientadas

Figura 64: Cinta de Möebius.

Figura 65: Construcción de una cinta de Möebius.

Figura 66: Superficie S con dos vectores normales.

Superficies orientadas
Si es posible elegir un vector unitario normal n en todos los puntos (x, y, z) de
modo que n varíe continuamente sobre S, entonces se dice que S es una
superficie orientada y la elección dada de n proporciona a S una orientación.
Hay dos orientaciones posibles para cualquier superficie orientable (véase la
Figura 67).

Figura 67: Las dos orientaciones de una superficie orientable.

En el caso de una superficie z = g(x, y) dada como gráfica de g, usamos la

Ecuación 50 para asociar con la superficie una orientación natural dada por el
vector normal unitario
Superficies orientadas
Si S es una superficie suave y orientable dada en la forma paramétrica por
medio de una función vectorial r(u, v), entonces automáticamente adquiere la
orientación del vector unitario normal
ru × rv
n= (52)
|ru × rv |
y la orientación opuesta se consigue con −n. Por ejemplo, en el Ejemplo 32
determinamos la representación paramétrica

r (φ, θ) = a sen φ cos θ i + a sen φ sen θ j + a cos φ k

para la esfera x2 + y 2 + z 2 = a2 . Luego, en el Ejemplo 38, se determinó que

rφ × rθ = a2 sen2 φ cos θ i + a2 sen2 φ sen θ j + a2 sen φ cos φ k

y
|rφ × rθ | = a2 sen φ.
De modo que la orientación inducida por r (φ, θ) se define por medio del vector
normal unitario
rφ × rθ 1
n= = sen φ cos θ i + sen φ sen θ j + cos φ k = r (φ, θ) .
|rφ × rθ | a
Superficies orientadas
Observe que n apunta en la misma dirección que el vector de posición, es decir,
hacia fuera de la esfera (véase la Figura 68). La orientación opuesta (es decir,
hacia adentro), se obtendría si invertimos el orden de los parámetros (véase la
Figura 69) porque rθ × rφ = −rφ × rθ .

Figura 68: Esfera: orientación positiva.

Figura 69: Esfera: orientación negativa.

Superficies orientadas
I En el caso de una superficie cerrada, es decir, una superficie que es la
frontera de una región sólida E, la convención es que la orientación
positiva es aquella para la cual los vectores normales señalan hacia afuera
de E, y los vectores normales que señalan hacia el interior dan la
orientación negativa (véanse las Figuras 68 y 69).

Integrales de superficie de campos vectoriales

Suponga que S es una superficie orientada con un vector unitario normal n, e
imagine que hay un fluido de densidad ρ(x, y, z) y campo de velocidad
v(x, y, z) que circula a través de S. (Piense que S es una superficie imaginaria
que no impide el flujo de fluidos, tal como una red para pescar atravesada en
un arroyo.) Entonces, el caudal (masa por unidad de tiempo) por unidad de
área es ρv. Si dividimos S en pequeños parches Sij , como en la Figura 70
(compare con la Figura 61), entonces Sij es casi plana y podemos aproximar la
masa del fluido que atraviesa en la dirección de la normal n por unidad de
tiempo mediante la cantidad

(ρv · n) A (Sij )
Superficies orientadas

Figura 70: Subdivisión de S en parches Sij .

Superficies orientadas

donde r, v y n se evalúan en algún punto sobre Sij . (Recuerde que la

componente del vector ρv en la dirección del vector unitario n es ρv · n.)
Según la Definición 48, luego de sumar estas cantidades y obtener el límite, el
resultado es la integral de superficie de la función ρv · n sobre S:
¨ ¨
ρv · ndS = ρ(x, y, z)v(x, y, z) · n(x, y, z) dS (53)
S S

y la interpretación física es el caudal que atraviesa S.

Si escribimos F = ρv, entonces F es también un campo vectorial sobre R3 y la
integral de la Ecuación 53 se vuelve
¨
F · n dS.
S

Una integral de superficie de esta forma se presenta con frecuencia en física,

aun cuando F no sea igual a ρv, y se le llama integral de superficie (o integral
de flujo) de F sobre S.
Superficies orientadas
Definición 6
Si F es un campo vectorial continuo definido sobre una superficie orientada S con un
vector unitario normal n, entonces la integral de superficie de F sobre S es
¨ ¨
F · dS = F · n dS.
S S

Esta integral también se denomina flujo de F a través de S.

I En lenguaje común, la Definición 6 establece que la integral de superficie de un
campo vectorial en S es igual a la integral de superficie de su componente
normal en S (como se definió antes).
I Si S está definida por una función vectorial r(u, v), entonces n está dada por la
Ecuación 52, y de acuerdo con la Definición 6 y la Ecuación 49, tenemos
¨ ¨
ru × rv
F · dS = F· dS
|ru × rv |
S S
¨  
ru × rv
= F (r(u, v)) · |ru × rv | dA
|ru × rv |
D

donde D es el dominio del parámetro. Por tanto,

¨ ¨
F · dS = F · (ru × rv ) dA. (54)
S D
Superficies orientadas

Ejemplo 43
Determine el flujo del campo vectorial F(x, y, z) = z i + y j + x k a través de la
esfera unitaria x2 + y 2 + z 2 = 1.

Solución
Como en el Ejemplo 40, utilizamos la representación paramétrica

r(φ, θ) = sen φ cos θ i + sen φ sen θ j + cos φ k, 0 ≤ φ ≤ π, 0 ≤ θ ≤ 2π.

Entonces
F (r(φ, θ)) = cos φ i + sen φ sen θ j + sen φ cos θ k
y, de acuerdo con el Ejemplo 38,

rφ × rθ = sen2 φ cos θ i + sen2 φ sen θ j + sen φ cos φ k.

Por tanto,

F (r(φ, θ)) · (rφ × rθ ) = cos φ sen2 φ cos θ + sen3 φ sen2 θ + sen2 φ cos φ cos θ
y, según la fórmula 54, el flujo es
Superficies orientadas

¨ ¨
F · dS = F · (rφ × rθ ) dA
S D
ˆ 2π ˆ π
2 sen2 φ cos φ cos θ + sen3 φ sen2 θ dφ dθ


=
0 0

ˆ π *0ˆ π
ˆ
 2π ˆ 2π
2
sen3 φ dφ sen2 θ dθ
=2 sen φ cos φ dφ

cos θ dθ +

0
ˆ 0  ˆ
0 0
π 2π
=0+ sen3 φ dφ sen2 θ dθ
0 0
4π
= .
3

Observación
Si el campo vectorial del Ejemplo 43 es un campo de velocidades que describe el
flujo de un fluido cuya densidad es 1, entonces la respuesta, 4π/3, representa el
caudal a través de la esfera unitaria en unidades de masa por unidad de tiempo.
Superficies orientadas

Figura 71: Esta figura muestra el campo vectorial F del ejemplo 43 en puntos sobre la esfera
unitaria.

En el caso de una superficie S dada por una gráfica z = g(x, y), podemos considerar a x y y como
parámetros y usamos la Ecuación 50 para escribir

∂g ∂g
 


F · rx × ry = (P i + Q j + R k) . − i− j+k .
∂x ∂y

Entonces la Fórmula 54 se convierte en

¨ ¨ 
∂g ∂g

F · dS = −P −Q +R dA (55)
∂x ∂y
S D

Esta fórmula toma la orientación hacia arriba de S; para una orientación hacia abajo multiplicamos por
−1. Es posible resolver fórmulas similares si S está dada por y = h(x, z) o x = k(y, z).
Superficies orientadas
Ejemplo
˜
44
Evalúe S F · dS, donde F(x, y, z) = y i + x j + z k y S es la frontera de la región sólida E encerrada
por el paraboloide z = 1 − x2 − y 2 y el plano z = 0.

Solución
La superficie S consiste en una superficie S1 parabólica en la parte superior y una superficie S2 circular
en el fondo (véase la Figura 72). Como S es una superficie cerrada, usamos la convención de la
orientación positiva (hacia afuera). Esto significa que S1 está orientada hacia arriba y que podemos usar
la Ecuación 55 en donde D es la proyección de S1 sobre el plano xy, a saber, el disco x2 + y 2 ≤ 1.
Puesto que
2 2
P (x, y, z) = y, Q(x, y, z) = x, R(x, y, z) = z = 1 − x − y

sobre S1 y
∂g ∂g
= −2x, = −2y.
∂x ∂y

Figura 72: Región sólida encerrrada por paraboloide y el plano z=0.

Superficies
Tenemos
orientadas
¨ ¨ 
∂g ∂g

F · dS = −P −Q +R dA
∂x ∂y
S1 D
¨ h i
2 2
= −y(−2x) − x(−2y) + 1 − x − y dA
D
¨  
2 2
= 1 + 4xy − x − y dA
D
ˆ 2π ˆ 1  
2 2
= 1 + 4r cos θ sen θ − r r dr dθ
0 0
ˆ 2π ˆ 1  
3 3
= r − r + 4r cos θ sen θ dr dθ
0 0
ˆ 2π 
1

1 (2π) + 0 =
π
= 4
+ cos θ sen θ dθ = 4 .
0 2

El disco S2 está orientado hacia abajo, de modo que su vector normal unitario es n = −k y tenemos

¨ ¨ ¨ ¨
F · dS = F · (−k) dS = (−z) dA = 0dA = 0
S2 S2 D D

˜
puesto que z = 0 sobre S2 . Finalmente, calculamos por definición S F · dS como la suma de las
integrales de superficie de F sobre las partes S1 y S2 :

¨ ¨ ¨
π π
F · dS = F · dS + F · dS = +0= .
2 2
S S1 S2
Superficies orientadas
Aunque introducimos la integral de superficie de un campo vectorial usando el
ejemplo de flujo de fluidos, este concepto surge también en otras situaciones
físicas. Por ejemplo, si E es un campo eléctrico (véase Ejemplo 4), entonces la
integral de superficie ¨
E · dS
S

recibe el nombre de flujo eléctrico de E a través de la superficie S. Una de las

leyes importantes de la electrostática es la ley de Gauss, la cual establece que
la carga neta encerrada por medio de una superficie cerrada S es
¨
Q = ε0 E · dS (56)
S

donde ε0 es una constante (que se denomina permitividad del vacío), y que

depende de las unidades que se utilicen. (En el SI,
ε0 ≈ 8.8542 × 10−12 C2 /N · m2 .) Por tanto, si el campo vectorial F del
Ejemplo 43 representa un campo eléctrico, podemos concluir que la carga
encerrada por S es Q = 43 πε0 .
Superficies orientadas
Otra aplicación de las integrales de superficie se encuentra en el estudio del flujo
de calor. Supongamos que la temperatura en un punto (x, y, z) de un cuerpo es
u(x, y, z). Entonces el flujo de calor se define como el campo vectorial
F = −K∇u
donde K es una constante determinada en forma experimental que se llama
conductividad de la sustancia. El flujo de calor a través de la superficie S en el
cuerpo lo define entonces la integral de superficie
¨ ¨
F · dS = −K ∇u · dS.
S S

Ejemplo 45
La temperatura u de una bola de metal es proporcional al cuadrado de la
distancia desde el centro de la misma. Determine el flujo de calor a través de
una esfera S de radio a con centro en el centro de la bola.

Solución
Tomando el centro de la bola como el origen, tenemos

u(x, y, z) = C x2 + y 2 + z 2

donde C es la constante de proporcionalidad.

Superficies orientadas
Luego, el flujo de calor es

F(x, y, z) = −K∇u = −KC (2x i + 2y j + 2z k)

donde K es la conductividad del metal. En lugar de usar la parametrización de

la esfera como en el Ejemplo 43, observamos que la normal unitaria hacia
afuera de la esfera x2 + y 2 + z 2 = a2 en el punto (x, y, z) es
1
n= (x i + y j + z k)
a
y entonces
2KC 2
x + y2 + z2 .


F·n=−
a
Pero sobre S tenemos x2 + y 2 + z 2 = a2 , de modo que F · n = −2aKC. Por
tanto, el flujo de calor a través de S es
¨ ¨ ¨
F · dS = F · ndS = −2aKC dS
S S S

= −2aKC A(S) = −2aKC 4πa2 = −8KCπa3 .

Teorema de Stokes
I Podemos considerar que el teorema de Stokes es una versión para varias
dimensiones del teorema de Green.
I Mientras el teorema de Green relaciona una integral doble sobre una región D
plana con una integral de línea alrededor de su curva frontera plana, el teorema
de Stokes relaciona una integral de superficie sobre una superficie S con una
integral de línea alrededor de la curva frontera de S (que es una curva en el
espacio).
I En la Figura 73 se muestra una superficie orientada con vector normal unitario
n. La orientación de S induce la orientación positiva de la curva frontera C
ilustrada en la Figura 73. Esto significa que si usted camina en la dirección
positiva alrededor de C con su cabeza señalando en la dirección de n, entonces
la superficie siempre quedará a su izquierda.

Figura 73: Una superficie orientada con un vector normal unitario n.

Teorema de Stokes

Teorema 10 (Teorema de Stokes)

Sea S una superficie suave por tramos y orientada que está acotada por una
curva C suave por tramos, simple y cerrada con orientación positiva. Sea F un
campo vectorial cuyas componentes tienen derivadas parciales continuas en una
región abierta en R3 que contiene a S. Entonces,
‰ ¨
F · dr = rot F · dS. (57)
C
S

Puesto que
‰ ˆ ¨ ¨
F · dr = F · Tds y rot F · dS = rot F · n dS.
C C
S S

El teorema de Stokes establece que la integral de línea alrededor de la curva

frontera de S de la componente tangencial de F es igual a la integral de
superficie de la componente normal del rotacional de F.
Teorema de Stokes
La curva acotada orientada en forma positiva de la superficie orientada S se
escribe a menudo como ∂S, de modo que el teorema de Stokes se puede
expresar como ¨ ˆ
rot F · dS = F · dr. (58)
∂S
S

Hay una analogía entre el teorema de Stokes, el teorema de Green y el teorema

fundamental del cálculo. Como antes, hay una integral con derivadas en el lado
izquierdo de la Ecuación 58 (recuerde que rot F es una clase de derivada de F)
y el lado derecho involucra los valores de F que están solamente en la frontera
de S.
De hecho, en el caso especial donde la superficie S es plana y está en el plano
xy con orientación hacia arriba, la normal unitaria es k, la integral de superficie
se vuelve una integral doble, y el teorema de Stokes se transforma en
‰ ¨ ¨
F · dr = rot F · dS = rot F · k dA.
C
S S

Esta es precisamente la forma vectorial del teorema de Green dado en la

Ecuación 38. Por tanto, el teorema de Green es realmente un caso especial del
teorema de Stokes.
Teorema de Stokes
Ejemplo

46
Evalúe C F · dr, donde F(x, y, z) = −y 2 i + x j + z 2 k y C es la curva de intersección del plano
y + z = 2 y el cilindro x2 + y 2 = 1. (La orientación de C es en el sentido contrario al de las manecillas
del reloj cuando se le ve desde arriba).

Solución 
La curva C (una elipse) se ilustra en la Figura . Aunque C F · dr se puede evaluar en forma directa, es
más fácil aplicar el teorema de Stokes. Primero calculamos

 
 i j k 
 
 ∂ ∂ ∂  = (1 + 2y) k.

rot F =  ∂x ∂y ∂z 
−y 2 2 

x z

Figura 74: Intersección de un cilindro con un plano.

Teorema de Stokes

Aunque hay muchas superficies cuya frontera es C, la elección más conveniente

es la región elíptica S en el plano y + z = 2 que está acotada por C. Si
orientamos a S hacia arriba, entonces C tiene la orientación positiva inducida.
La proyección D de S sobre el plano xy es el disco x2 + y 2 ≤ 1, por lo que al
aplicar la Ecuación 55 con z = g(x, y) = 2 − y tenemos
‰ ¨ ¨
F · dr = rot F · dS = (1 + 2y) dA
C
S D
ˆ 2π ˆ 1
= (1 + 2r sen θ) r dr dθ
0 0
ˆ 2π 1 ˆ 2π
r2 r3

1 2


= + 2 sen θ dθ = 2
+ 3
sen θ dθ
0 2 3 0 0
1
= (2π) + 0 = π.
2
Teorema de Stokes
Ejemplo 47
˜
Mediante el teorema de Stokes, calcule la integral S
rot F · dS, donde
F(x, y, z) = xz i + yz j + xy k y S es la parte de la esfera x2 + y 2 + z 2 = 4 que
está situada en el interior del cilindro x2 + y 2 = 1 y encima del plano xy (véase
la Figura 75).

Figura 75: S es la parte de la esfera x2 + y2 + z 2 = 4.

Solución
Para determinar la curva frontera C resolvemos las ecuaciones
x2 + y 2 + z 2 = 4 y x2 + y 2 = 1. Al efectuar una diferencia, obtenemos z 2 = 3
Teorema de Stokes
√
y entonces z = 3 (porque z√ > 0). Por tanto, C es el círculo dado por las
ecuaciones x2 + y 2 = 1, z = 3. Una ecuación vectorial de C es
√
r(t) = cos t i + sen t j + 3k, 0 ≤ t ≤ 2π
así r0 (t) = − sen t i + cos t j.
Además, tenemos
√ √
F (r(t)) = 3 cos t i + 3 sen t j + cos t sen t k.
Por tanto, de acuerdo con el teorema de Stokes,
¨ ‰ ˆ 2π
rot F · dS = F · dr = F (r(t)) · r0 (t)dt
C 0
S
ˆ 2π  √ √ 
= − 3 cos t sen t + 3 sen t cos t dt
0
√ ˆ 2π
= 3 0 dt = 0.
0

Observación
Observe que en el Ejemplo 47 calculamos una integral de superficie
simplemente con el conocimiento de los valores de F sobre la curva frontera C.
Esto significa que si tenemos otra superficie orientada con la misma curva
frontera C, entonces ¡obtenemos exactamente el mismo valor para la integral
de superficie!
Teorema de Stokes
En general, si S1 y S2 son superficies orientadas con la misma curva frontera
orientada C, y ambas satisfacen las hipótesis del teorema de Stokes, entonces
¨ ‰ ¨
rot F · dS = F · dr = rot F · dS (59)
C
S1 S2

Este hecho es útil cuando es difícil integrar sobre una superficie, pero fácil
integrar en la otra.
Enseguida usaremos el teorema de Stokes para dilucidar el significado del vector
rotacional. Suponga que C es una curva cerrada orientada y v representa el
campo de velocidad en un flujo de fluidos. Consideremos la integral de línea
ˆ ˆ
v · dr = v · Tds
C C

y recordemos que v · T es la componente de v en la dirección del vector

tangente unitario T. Esto significa que, a medida que es más cercana la
´ de v a la dirección de T, es más grande el valor de v · T. Por lo
dirección
tanto, C v · dr es una medida de la tendencia del fluido a moverse alrededor
de C y se llama circulación de v alrededor de C (véase la Figura 76).
Teorema de Stokes

Figura 76: Circulación de v alrededor de C.

Ahora, sea P0 (x0 , y0 , z0 ) un punto en el fluido y sea Sa un disco pequeño con

radio a y centro P0 . Entonces (rot F) (P ) ≈ (rot F) (P0 ) para todos los puntos
P en Sa porque rot F es continua. Por lo tanto, de acuerdo con el teorema de
Stokes, obtenemos la siguiente aproximación a la circulación a lo largo de la
círculo frontera Ca :
ˆ ¨ ¨
v · dr = rot F · dS = rot v · ndS
Ca
Sa Sa
¨
= rot v(P0 ) · n(P0 )dS = rot v(P0 ) · n(P0 )πa2 .
Sa
Teorema de Stokes
Esta aproximación es mejor cuando a → 0 y entonces
ˆ
1
rot v(P0 ) · n(P0 ) = lı́m v · dr. (60)
a→0 πa2 C
a

La Ecuación 60 da la relación entre el rotacional y la circulación. Se demuestra

que rot v · n es una medida del efecto de giro del fluido respecto al eje n. El
efecto de rotación es mayor respecto al eje paralelo a rot v.
I Imagine una pequeñísima rueda de paletas en el fluido en un punto P,
como en la Figura 77; la rueda gira más rápido cuando su eje es paralelo
a rot v.

Figura 77: La rueda gira más rápido cuando su eje es paralelo a rot v.
El teorema de la divergencia

El teorema de Green en la versión vectorial se expresa según (39) como

‰ ¨
F · n ds = div F(x, y)dA
C
D

donde C es la curva frontera orientada en la dirección positiva de la región D

del plano. Si estuvieramos tratando de generalizar este teorema a los campos
vectoriales sobre R3 , podríamos plantear la conjetura de que
¨ ˚
F · n dS = div F(x, y, z)dV (61)
S E

donde S es la superficie frontera de la región sólida E. Resulta que la

Ecuación 61 es cierta, con las hipótesis adecuadas, y se llama teorema de la
divergencia. Observe su similitud con el teorema de Green y el teorema de
Stokes: este teorema relaciona la integral de una derivada de una función
(div F, en este caso) sobre una región con la integral de la función original F
en la frontera de la región.
El teorema de la divergencia
I En esta etapa debería revisar los distintos tipos de regiones en las cuales
es capaz de evaluar integrales triples.
I Plantee y demuestre el teorema de divergencia para regiones E que son
simultáneamente de los tipos 1, 2 y 3 llamadas regiones sólidas simples.
(Por ejemplo, las regiones acotadas por elipsoides o cajas rectangulares
son de este tipo.) La frontera de E es una superficie cerrada, y use la
convención según la cual la orientación positiva es hacia afuera, es decir,
el vector normal unitario n se dirige hacia afuera desde E.

Teorema 11 (Teorema de la divergencia2 )

Sea E una región sólida simple y S la superficie frontera de E, dada con
orientación positiva (hacia afuera). Sea F un campo vectorial cuyas funciones
componentes tienen derivadas parciales continuas en una región abierta que
contiene E. Entonces,
¨ ˚
F · dS = div F dV.
S E

2
El teorema de la divergencia se llama a veces teorema de Gauss en honor al matemático alemán
Karl Friedrich Gauss (1777-1855), quien descubrió este teorema durante su investigación sobre
electrostática. En Europa del Este, el teorema de la divergencia se conoce con el nombre de teorema de
Ostrogradsky, en honor al matemático ruso Mikhail Ostrogradsky (1801-1862), quien publicó este
resultado en 1826.
El teorema de la divergencia
Por lo tanto, el teorema de divergencia establece que, bajo las condiciones
dadas, el flujo de F a través de la superficie frontera de E es igual a la triple
integral de la divergencia de F sobre E.
Demostración.
Sea F = P i + Q j + R k. Entonces
∂P ∂Q ∂R
div F = + +
∂x ∂y ∂z
de modo que
˚ ˚ ˚ ˚
∂P ∂Q ∂R
div F dV = dV + dV + dV.
∂x ∂y ∂z
E E E E

Si n es el vector normal unitario hacia afuera de S, entonces la integral de

superficie en el lado izquierdo o del teorema de la divergencia es
¨ ¨ ¨
F · dS = F · n dS = (P i + Q j + R k) · n dS
S S S
¨ ¨ ¨
= P i · n dS + Q j · n dS + R k · n dS.
S S S
El teorema de la divergencia
Por tanto, para demostrar el teorema de la divergencia, es suficiente demostrar
las tres ecuaciones siguientes:
¨ ˚
∂P
P i · n dS = dV (62)
∂x
S E
¨ ˚
∂Q
Q j · n dS = dV (63)
∂y
S E
¨ ˚
∂R
R k · n dS = dV. (64)
∂z
S E

Para demostrar la Ecuación 64, recurrimos al hecho de que E es una región

tipo 1:
E = { (x, y, z)| (x, y) ∈ D, u1 (x, y) ≤ z ≤ u2 (x, y)}
donde D es la proyección de E sobre el plano xy. Como E es una región tipo 1
tenemos ˚ ¨ "ˆ u2 (x,y) #
∂R ∂R
dV = (x, y, z)dz dA
∂z u1 (x,y) ∂z
E D

y, por tanto, según el teorema fundamental del cálculo,

El teorema de la divergencia
˚ ¨
∂R
dV = [R (x, y, u2 (x, y)) − R (x, y, u1 (x, y))] dA (65)
∂z
E D

La superficie frontera S consiste en tres partes: la superficie del fondo S1 , la

superficie de la tapa S2 y posiblemente una superficie vertical S3 , la cual se
ubica encima de la curva frontera de D (véase la Figura 78. Podría ocurrir que
S3 no aparezca, como en el caso de una esfera). Observe que en S3 tenemos
k · n = 0, porque k es vertical y n es horizontal, y así
¨ ¨
R k · n dS = 0 dS = 0.
S3 S3

Por esto, sin que importe si hay una superficie vertical, podemos escribir
¨ ¨ ¨
R k · n dS = R k · n dS + R k · n dS. (66)
S S1 S2

La ecuación de S2 es z = u2 (x, y), (x, y) ∈ D, y la normal n hacia afuera

apunta hacia arriba, de modo que de la Ecuación 55 (en la que Rk reemplaza a
F), tenemos
El teorema de la divergencia

Figura 78: La superficie frontera S consiste de tres partes.

El teorema de la divergencia
¨ ¨
R k · n dS = R (x, y, u2 (x, y)) dA.
S2 D

Sobre S1 tenemos z = u1 (x, y), pero en este caso la normal n hacia afuera
apunta hacia abajo, de modo que multiplicamos por −1:
¨ ¨
R k · n dS = − R (x, y, u1 (x, y)) dA.
S1 D

Por tanto, la Ecuación 66 da

¨ ¨
R k · n dS = [R (x, y, u2 (x, y)) − R (x, y, u1 (x, y))] dA
S D

La comparación con la Ecuación 65 demuestra que

¨ ˚
∂R
R k · n dS = dV.
∂z
S E

Las Ecuaciones 62 y 63 se demuestran en forma similar usando las expresiones

para E como una región tipo 2 o 3, respectivamente.
El teorema de la divergencia

Ejemplo 48
Determine el flujo del campo vectorial F(x, y, z) = z i + y j + x k sobre la
esfera unitaria x2 + y 2 + z 2 = 1.

Solución
Primero calculamos la divergencia de F:
∂ ∂ ∂
div F = (z) + (y) + (x) = 1.
∂x ∂y ∂z
La esfera unitaria S es la frontera de la bola unitaria B definida por
x2 + y 2 + z 2 ≤ 1. En estos términos, el teorema de la divergencia da el flujo
como ¨ ˚ ˚
4π
F · dS = div F dV = div 1 dV = V (B) = .
3
S B B

Observación
La solución del Ejemplo 48 se debe comparar con la solución del Ejemplo 43.
El teorema de la divergencia
Ejemplo
˜
49
Evalúe F · dS, donde
S
 2

F(x, y, z) = xy i + y 2 + exz j + sen (xy) k

y S es la superficie de la región E acotada por el cilindro parabólico z = 1 − x2

y los planos z = 0, y = 0 y y + z = 2 (véase la Figura 79).

Figura 79: Superfice acotada por z = 1 − x2 , z = 0, y = 0 y y + z = 2.

El teorema de la divergencia
Solución
Sería difícil en extremo evaluar en forma directa la integral de superficie. (Tendríamos que evaluar
cuatro integrales de superficie correspondientes a las cuatro partes de S.) Además, la divergencia de F
es mucho menos complicada que la misma F:

∂ ∂ ∂
 
2 xz 2
div F = (xy) + y +e + sen (xy) = y + 2y = 3y.
∂x ∂y ∂z

Por tanto, usamos el teorema de la divergencia para transformar la integral de superficie dada en una
integral triple. La manera más fácil de evaluar la integral triple, es expresar E como una región tipo 3:

n o
2
E = (x, y, z)| − 1 ≤ x ≤ 1, 0 ≤ z ≤ 1 − x , 0 ≤ y ≤ 2 − z .

Entonces tenemos
¨ ˚ ˚
F · dS = div F dV = 3y dV
S E E
ˆ 1 ˆ 1−x2 ˆ 2−z ˆ 1 ˆ 1−x2
(2 − z)2
=3 y dy dz dx = 3 dz dx
−1 0 0 −1 0 2
ˆ 1 #1−x2 ˆ
(2 − z)3
"
3 1 1  2
 3 
= − dx = − x +1 dx
2 −1 3 0
2 −1
ˆ 1
6 4 2
 184
=− x + 3x + 3x − 7 dx = .
0 35
El teorema de la divergencia
I Hemos demostrado el teorema de la divergencia sólo para regiones sólidas
simples, pero se puede demostrar también para regiones que son uniones
finitas de regiones sólidas simples. (El procedimiento es similar al usado
para generalizar el teorema de Green.)
I Por ejemplo, consideremos la región E que se ubica entre las superficies
cerradas S1 y S2 , donde S1 queda en el interior de S2 . Sean n1 y n2 las
normales hacia afuera de S1 y S2 . Por tanto, la superficie frontera de E
es S = S1 ∪ S2 y su normal n está dada por n = −n1 en S1 y n = n2 en
S2 (véase la Figura 80). Al aplicar el teorema de la divergencia a S
obtenemos

Figura 80: Frontera de E es S = S1 ∪ S2 .

El teorema de la divergencia
˚ ¨ ¨
div F dV = F · dS = F · ndS
E S S
¨ ¨
= F · (−n1 ) dS + F · n2 dS
S1 S2
¨ ¨
=− F · dS + F · dS. (67)
S1 S2

Ejemplo 50
En el Ejemplo 4 consideramos al campo eléctrico:
εQ
E(x) = x
|x|3
donde la carga eléctrica Q está localizada en el origen y x = (x, y, z) es un
vector de posición. Usamos el teorema de la divergencia para demostrar que el
flujo eléctrico de E a través de cualquier superficie cerrada S2 que encierra el
origen es ¨
E · dS = 4πεQ.
S2
El teorema de la divergencia
Solución
La dificultad es que no tenemos una ecuación explícita para S2 porque es cualquier
superficie que encierra el origen. La más simple superficie como ésta es la esfera, así
que sea S1 una pequeña esfera con radio a y centro en el origen. Podemos verificar
que div E = 0. Por tanto, la Ecuación 67 da como resultado
¨ ¨ ˚ ¨
E · dS = E · dS + div F dV = E · dS.
S2 S1 E S1

Lo importante de este cálculo es que podemos evaluar la integral de superficie sobre

S1 porque S1 es una esfera. El vector normal unitario en x es x/ |x|. Por tanto,
 
εQ x εQ εQ εQ
E·n= 3
x· = x·x= =
|x| |x| |x|4 |x|2 a2

ya que la ecuación de S1 es |x| = a. Debido a esto tenemos

¨ ¨ ¨
εQ εQ εQ
E · dS = E · ndS = dS = A (S1 ) = 2 4πa2 = 4πεQ.
a2 a2 a
S2 S1 S1
El teorema de la divergencia
Esto demuestra que el flujo eléctrico de E es 4πεQ a través de cualquier
superficie cerrada S2 que contiene el origen. Es un caso especial de la ley de
Gauss (Ecuación 56) para una carga sencilla. La relación entre ε y ε0 es
ε = 1/ (4πε0 ).
Otro campo en el que se puede aplicar el teorema de la divergencia es en el
flujo de fluidos. Sea v(x, y, z) el campo de velocidad de un fluido con densidad
constante ρ. Por tanto, F = ρv es el caudal por unidad de área. Si
P0 (x0 , y0 , z0 ) es un punto en el fluido y Ba es una bola con centro P0 y radio
muy pequeño a, entonces div F (P ) ≈ div (P0 ) para todos los puntos en Ba
porque div F es continua. El valor del flujo sobre la esfera frontera Sa se
aproxima como sigue:
¨ ˚ ˚
F · dS = div F dV ≈ div F (P0 ) dV = div F (P0 ) V (Ba ) .
Sa Ba Ba

Esta aproximación mejora cuando a → 0 y sugiere

¨
1
div F (P0 ) = lı́m F · dS. (68)
a→0 V (Ba )
Sa
El teorema de la divergencia

I La Ecuación 68 establece que div F (P0 ) es el flujo neto que sale por
unidad de volumen en P0 . (Ésta es la razón del nombre divergencia). Si
div F (P0 ) > 0, el flujo neto cerca de P es hacia afuera y se dice entonces
que P es una fuente. Si div F (P0 ) < 0, el flujo neto cerca de P es hacia
adentro y se dice entonces que P es un sumidero.
I Al parecer, en el caso del campo vectorial de la Figura 81 los vectores que
terminan en las cercanías de P1 son más cortos que los vectores que
empiezan cerca de P1 . Por esto, el flujo neto es hacia afuera en la
vecindad de P1 , de modo que div F (P1 ) > 0, por lo que P1 es una
fuente. Por lo contrario, las flechas que entran cerca de P2 son más largas
que las flechas que salen. En este caso, el flujo neto es hacia adentro, de
modo que div F (P2 ) < 0 y P2 es un sumidero. Puede utilizar la fórmula
para F con el fin se confirmar esta impresión. Puesto que F = x2 i + y 2 j,
div F = 2x + 2y, la cual es positiva cuando y > −x. Por eso los puntos
por arriba de la recta y = −x son fuentes, y los que se encuentran abajo
son sumideros.
El teorema de la divergencia

Figura 81: El campo vectorial F = x2 i + y2 j.