Lectura Fundamental 7
Lectura Fundamental 7
Lectura Fundamental 7
Lectura Fundamental
Integrales de lı́nea
Contenido
2 Campos vectoriales
5 Teorema de Green
6 Ejercicios
Palabras clave: campo escalar, campo vectorial, campo conservativo, integral de lı́nea.
En este escenario, primero se define el concepto de integral de lı́nea de un campo escalar y luego, el de integral
de lı́nea de un campo vectorial. Además se describen algunas de las propiedades principales de una integral de
lı́nea, ası́ como métodos que permiten calcular este tipo especial de integrales, destacando entre estos el teorema
de Green.
Recuerde que se le llama campo escalar a cualquier función f : A ⊂ Rn −→ R, es decir, a una aplicación que le
asigna a cada punto x ∈ A un número f (x) ∈ R. De igual manera, tenga en cuenta que una trayectoria es una
función r : I ⊂ R −→ Rn , donde I es un intervalo
del de números reales; se dice que una curva C ⊂ R
conjunto n
está parametrizada por una trayectoria r si C = r(t) t ∈ I . Siempre se asumirá que una curva C es suave por
partes, lo que significa que C consiste en una unión finita de curvas C1 , C2 , . . . , Ck unidas extremo por extremo,
donde cada curva Ci está parametrizada por una trayectoria diferenciable cuya derivada es continua en el interior
del intervalo donde está difinida.
Definición 1.1. Considere un campo escalar continuo f : A ⊂ Rn −→ R y una curva C parametrizada por
la trayectoria r : [a, b] −→ n
Z R tal que C está contenida en A. La integral de lı́nea de f a lo largo de la
curva C, se denota por f (x)ds, es la integral,
C
Z Z b
f r(t) |r0 (t)| dt
f (x) ds := (1)
C a
Suponga que C es la traza de la trayectoria r : [a, b] −→ Rn . Como C es una curva suave por pedazos, existe una
secuencia a = t0 < t1 < t2 < . . . < tk−1 < tk = b, tal que r es diferenciable en el intervalo abierto
Z (ti−1 , ti ), para
1 ≤ i ≤ k; ası́, aunque r0 (t) no exista para t = t0 , t1 , t2 , . . . , tk , se define la integral de lı́nea f (x) ds como la
C
Z Xk Z ti
f r(t) |r0 (t)| dt.
suma de las integrales f (x) ds =
C i=1 ti−1
2
Ejemplo 1.2. Considere la curva C parametrizada por la trayectoria r : [1, 2] −→ R3 , r(t) = t, t2 , t . Hallar la
2
Z
(x + y)z (x+y)z t2
t+ t2 t
integral de lı́nea ds. El campo escalar es f (x, y, z) = y+z , luego se tiene que f t, 2 , t = t2 = t.
C y+z √ √ 2
+t
0 0 2 2 2 2
Como r (t) = (1, t, 1), entonces |r (t)| = 1 + t + 1 = 2 + t . Por lo tanto:
Z
(x + y)z
Z 2 p 1 2 33/2 √
ds = t 2 + t2 dt = (2 + t2 )3/2 = (23/2 − 1) = 3(23/2 − 1)
C y+z 1 3 1 3
Considere una trayectoria r : [a, b] −→ Rn y una función diferenciable monótona ϕ : [c, d] −→ [a, b], luego la
n
trayectoria r̂ : [c, d] −→ R , definida por r̂(t) := r ϕ(t) traza la misma curva que la trayectoria r; la trayectoria r̂
se denomina una reparametrización de r. Si ϕ es creciente, entonces el punto inicial de r̂ es igual al punto inicial
de r, el punto final de r̂ es igual al punto inicial de r y las dos trayectorias recorren la curva C en el mismo sentido.
1
Si ϕ es decreciente, entonces el punto inicial de r̂ es igual al punto final de r, el punto final de r̂ es igual al punto
inicial de r y la trayectoria r̂ recorre a la curva C en el sentido contrario a la trayectoria r. En cualquier caso,
por la regla de la cadena, se tiene que |r̂0 (t)| = |r Zϕ(t) ||ϕ0 (t)|. Ası́, usando la fórmula del cambio de variables
para integrales, es fácil ver que la integral de lı́nea f (x)ds no depende de la parametrización de la curva C en
C
el sentido de que dicha integral tendrá el mismo valor tanto si se usa la trayectoria r como si se usa cualquier
reparametrización r̂(t) = r ϕ(t) .
Una aplicación de la integral definida consiste en determinar la masa de un alambre delgado y recto de longitud
` > 0; si dicho alambre se ubica en el intervalo cerrado [0, `] y si su densidad de masa en el punto x ∈ [0, `] es δ(x),
Z `
entonces la masa del alambre es m = δ(x)dx. De manera similar, se puede suponer que se tiene un alambre en
0
el espacio tridimensional cuya forma es una curva C, la cual está parametrizada por r : [a, b] −→ R3 . Si la densidad
3
de masa del alambre en un punto Z (x, y, z) de la curva C viene dada por la función δ : C ⊂ R −→ R, entonces
b−a
la masa m del alambre es m = δ(x, y, z)ds. En efecto, dado un entero n ≥ 2, tome ∆t = n y considere la
C
secuencia a = t0 < t1 < t2 < . . . < tn−1 < tn = b, donde ti = a + i∆t, que determina una partición de [a, b]
en n subintervalos de igual longitud ∆t. Al restrigir la trayectoria r a cada subintervalo [xi−1 , xi ] se obtiene una
curva Ci asociada a la trayectoria r : [xi−1 , xi ] −→ R3 , para 1 ≤ i ≤ n; es claro que C es la unión de los pedazos
de curvas C1 , C2 , . . . , Cn , los cuales se unen con los extremos. Suponga que la masa de cada pedazo de curva Ci
∗ ∗
es ∆mi y que su longitud es ∆`i . Tome un parámetro arbitrario
∗
ti ∈ [ti−1 , ti ], de manera que r(ti ) es un punto
de Ci , luego la masa ∆mi es aproximadamente igual a δ r(ti ) ∆`i , cuando n es suficientemente grande (lo que
hace que los pedazos de curvas Ci sean pequeños). Por otro lado, la longitud ∆`i es aproximadamente igual a la
rapidez de la trayectoria r en algún punto de Ci multiplicada por el intervalo de tiempo ∆t necesario para recorrer
la curva Ci , en particular ∆`i ≈ |r0 (t∗i )|∆t. Ası́, ∆mi ≈ δ r(t∗i ) ∆`i ≈ δ r(t∗i )|r0 (t∗i )|∆t. La suma de Riemann
Xn X n
∆mi = δ r(t∗i )|r0 (t∗i )|∆t nos da una aproximación de la masa total m del alambre; esta aproximación será
i=1 i=1
cada vez mejor a medida que el entero n sea mayor y, en consecuencia, tanto los subintervalos [xi−1 , xi ] como las
subcurvas Ci sean de menor longitud. Por lo tanto, usando la definición de integral definida, se concluye que la
Xn Z b Z
δ r(t∗i ) |r0 (t∗i )|∆t = δ r(t) |r0 (t)|dt =
masa del alambre es m = lim f (x)ds. Un argumento similar
n→+∞ a C
i=1
permite determinar el centro de masa del mismo alambre. A continuación se resumen estos resultados.
Ejemplo
√ 1.3. Considere
√ un alambre con la forma de una hélice C, la cual está parametrizada por la trayectoria
r(t) = 3 cos ti+ 3 sen tj+tk, 0 ≤√t ≤ 3π. Si √ su densidad de masa es δ(x, y, z)√= z, calcule la masa y el centro
√ de
0 0 2 2 2
masa del alambre. Como r (t) = − 3 sen ti + 3 cos tj + k, entonces |r (t)| = 3 sen t + 3 cos t + 1 = 4 = 2.
Z Z 3π √ √ Z 3π 3π
0
Luego, la masa del alambre es m = δ(x, y, z)ds = δ 3 cos t, 3 sen t, t |r (t)|dt = 2t dt = t2 = 9π 2 .
C 0 0 0
2
Las coordenadas del centro de masa (xc , yc , zc ) del alambre son:
Z 3π √ √ Z √ √
2 3 3π
Z
1 1 2 3 3π 4 3
xc = xδ(x, y, z)ds = 2 3 cos t(2t)ds = t cos t dt = t sen t + cos t = − 2
m C 9π 0 9π 2 0 9π 2 0 9π
Z 3π √ √ √ √
2 3 3π
Z Z
1 1 2 3 3π 2 3
yc = yδ(x, y, z)ds = 2 3 sen t(2t)ds = t sen t dt = −t cos t + sen t =
m C 9π 0 9π 2 0 9π 2 0 3π
Z Z 3π Z 3π 3
1 1 2 2 t 3π
zc = zδ(x, y, z)ds = 2 t(2t)ds = 2 t2 dt = 2 = 2π
m C 9π 0 9π 0 9π 3 0
√ √
Por lo tanto, el centro de masa del alambre está en el punto − 49π23 , 23π3 , 2π .
Otra aplicación interesante de las integrales de lı́nea de un campo escalar es la siguiente: considere una trayectoria
r : [a, b] −→ R2 y un campo escalar f : C ⊂ R2 −→ R, donde C es la curva de r y donde f (x, y) ≥ 0, para todo
(x, y) en la curva C. Observe que la gráfica de f está por encima de la curva C, pensando en el plano R2 como
el subespacio z = 0 de R3 . A su vez, la gráfica de f determina una curva C 0 en R3 , ver la Figura 2. Si S es la
superficie entre la curva C y la curva C 0 , es decir, S es la unión de los segmentos verticales que conectan
Z el punto
(x, y, 0) con el punto (x, y, f (x, y)), para cada (x, y) ∈ C, entonces el área de dicha superficie es A = f (x, y)ds.
C
Para deducir este resultado, una vez más, basta hacer una interpretación del problema usando sumas de Riemann;
se deja como ejercicio al lector hacer dicha deducción.
2. Campos vectoriales
En esta sección se introducen los campos vectoriales y se muestran algunos ejemplos básicos de estos. Al final de
la sección se define el rotacional y la divergencia de un campo vectorial.
3
Observe que un campo vectorial es simplemente una función vectorial cuyo dominio e imagen están en el espacio
Rn . Sin embargo, se debe interpretar un campo vectorial F como una aplicación que le asigna a cada punto
x ∈ Rn un vector F(x) ∈ Rn ; en los casos n = 2 y n = 3 se puede representar gráficamente al campo vectorial F
colocando en cada posición x del espacio al vector F(x), pero anclado en dicho punto, ver la Figura 3. Ası́, un
campo vectorial puede servir para modelizar, por ejemplo, la dinámica de un fluido que no depende del tiempo,
de manera que en un punto x, F(x) representa la velocidad del fluido en dicho punto.
Más adelante se exponen algunas propiedades de los campos vectoriales que son de este tipo (ver la definición 4.1).
Ejemplo 2.3. Suponga que una partı́cula de masa M se ubica en el punto (0, 0, 0) de R3 . Según la ley de
gravitación de Newton, la fuerza de gravedad de esta partı́cula sobre otra
GM m
de masa m en el punto (x, y, z) 6= (0, 0, 0)
viene dada por el campo vectorial F(x, y, z) = − 3/2 x, y, z , donde G es la constante de gravitación
x2 +y 2 +z 2
universal. Observe que cada vector F(x, y, z) anclado en el punto (x, y, z) apunta en la dirección del origen. Fı́jese
que F = −∇V , donde V (x, y, z) = − √ GM2
m
2 2
. El campo escalar V (x, y, z) es la energı́a potencial gravitacional
x +y +z
de las dos partı́culas de masas M y m.
Ejemplo 2.4. Sea F : R2 −→ R2 el campo vectorial F(x, y) = −yi + xj = (−y, x). Note que para todo punto
(x, y) del plano, el vector F(x, y) es perpendicular al vector (x, y). La Figura 4 muestra algunos de los vectores
F(x, y) anclados en los puntos (x, y) respectivos. Este campo vectorial describe un movimiento giratorio y en
sentido antihorario con respecto al origen, en donde los puntos giran a una velocidad mayor a medida que están
más alejados del origen.
Para finalizar esta sección definiendo el rotacional de un campo vectorial en R3 . El rotacional de un campo
vectorial F : U ⊂ R3 −→ R3 , cuyas funciones componentes tienen derivadas parciales continuas en U , es otro
campo vectorial en R3 , el cual, si se piensa en el campo vectorial F como el flujo de un fluido, describe la rotación
infinitesimal de F en cada punto de U . Esta idea se explicará con más detalle en el próximo Escenario, una vez
establecido el teorema de Stokes. Por lo pronto, se define el rotacional de un campo vectorial F de la siguiente
forma: suponga que F(x, y, z) = F1 (x, y, z)i + F2 (x, y, z)j + F3 (x, y, z)k y, como se indicó anteriormente, las
derivadas parciales de las funciones componentes F1 , F2 y F3 son continuas en U . El rotacional de F es el campo
4
vectorial rot F : U ⊂ R3 −→ R3 defindo por:
∂F3 ∂F2 ∂F1 ∂F3 ∂F2 ∂F1
rot F(x, y, z) := (x, y, z) − (x, y, z) i + (x, y, z) − (x, y, z) j + (x, y, z) − (x, y, z) k
∂y ∂z ∂z ∂x ∂x ∂y
Esta identidad es más sencilla de recordar si se piensa en el simbolo del gradiente ∇ como un “operador diferenciable
∂ ∂ ∂
vectorial” dado por el vector ∂x i + ∂y j + ∂z k, de manera que, para un campo escalar diferenciable f , el gradiente
de f es ∇f = ∂x ∂ ∂
i + ∂y ∂
j + ∂z k f = ∂f ∂f ∂f
∂x i + ∂y j + ∂z k. De esta manera, se cumple que el campo vectorial rot F
viene dado como el producto vectorial:
i j k
∂ ∂ ∂ ∂F 3 ∂F2 ∂F1 ∂F 3 ∂F 2 ∂F1
rot F = ∇ × F = ∂x ∂y ∂z = − i+ − j+ − k
F F F ∂y ∂z ∂z ∂x ∂x ∂y
1 2 3
Ejemplo 2.5. Calcular el rotacional de F(x, y, z) = y 2 z 2 i + x2 z 2 j + x2 y 2 k. Hallar rot F(2, −1, 2). Se calculará
primero rot F en un punto (x, y, z) ∈ R3 arbitrario. Ası́:
i j k
∂ ∂ ∂
rot F = ∇ × F = ∂x ∂y ∂z
y 2 z 2 x2 z 2 x2 y 2
∂ 2 2 ∂ 2 2 ∂ 2 2 ∂ 2 2 ∂ 2 2 ∂ 2 2
= (x y ) − (x z ) i + (y z ) − (x y ) j + (x z ) − (y z ) k
∂y ∂z ∂z ∂x ∂x ∂y
= 2x2 y − 2x2 z i + 2y 2 z − 2xy 2 j + 2xz 2 − 2yz 2 k = 2x2 (y − z)i − 2y 2 (x − z)j + 2z 2 (x − y)k
Se dice que un campo vectorial F en R3 es irrotacional, si rot F = 0. Intuitivamente, esto significa que si se
coloca en cada punto del dominio de F una pequeña rueda con aspas, esta se moverá por la acción de F, pero no
girará alrededor de su eje, nuevamente pensando en F como el flujo de un fluido. El siguiente resultado, de fácil
demostración, muestra que todo campo gradiente es irrotacional.
Proposición 2.6. Sea f : U ⊂ R3 −→ R un campo escalar con derivadas parciales de segundo orden continuas
en el abierto U . Entonces rot (∇f )(x, y, z) = 0, para todo (x, y, z) ∈ U .
Considere ahora un campo vectorial F(x, y) = F1 (x, y)i + F2 (x, y)j definido en un abierto U del plano R2 . Si se
piensa en R2 como el subespacio z = 0 de R3 y se extiende el campo F al campo vectorial F̂ en R 3
dado por
∂F2 ∂F1
F̂(x, y, z) := F1 (x, y), F2 (x, y), 0 , se puede comprobar fácilmente que rot F̂(x, y, z) = ∂x − ∂y k. De esta
manera, se justifica definir el rotacional de un campo vectorial F en un abierto U del plano R2 como el campo
escalar rot F(x, y) := ∂F ∂F1
∂x (x, y) − ∂y (x, y).
2
La integral de lı́nea de un campo escalar se define en esta sección para un campo vectorial en Rn . Recuerde que,
dados dos vectores x = (x1 , x2 , . . . , xn ) y y = (y1 , y2 , . . . , yn ) en Rn , el producto escalar entre x y y es el número
real x · y = x1 y1 + x2 y2 + . . . + xn yn .
5
Definición 3.1. Considere una curva C suave por partes parametrizada por la trayectoria r : [a, b] −→ Rn y
sea F : A ⊂ Rn −→ Rn un campo vectorial definido
Z en el conjunto A que contiene a la curva C. La integral
de lı́nea de F a lo largo de C, se denota por F · dr, se define como la integral:
C
Z Z b
F r(t) · r0 (t)dt
F · dr := (4)
C a
Si F(x) =Z F1 (x), F2 (x), . . . , Fn (x) , donde cada Fi (x) es la i-ésima componente de F(x), se acostumbra usar la
notación F1 dx1 + F2 dx2 + . . . + Fn dxn para indicar a la integral de lı́nea del campo vectorial F a lo largo de la
C
curva C. La razón de esta notación se debe a lo siguiente: si r(t)= x1 (t), x2 (t), . . . , xn (t) , entonces se tiene que
dr = r0 (t)dt = x01 (t)dt, x02 (t)dt, . . . , x0n (t)dt = dx1 , dx2 , . . . , dxn , por lo cual:
Z Z b Z b Z
0
F1 r(t) x01 (t)dt + . . . + Fn r(t) x0n (t)dt =
F · dr := F r(t) · r (t)dt = F1 dx1 + F2 dx2 + . . . + Fn dxn
C a a C
I
Cuando C es una curva cerrada, es decir, cuando r(a) = r(b), también se acostumbra a escribir Fdr para
C
denotar la integral de lı́nea de F a lo largo de la curva cerrada C.
√
Ejemplo 3.2. En R2 considere el campo vectorial F(x, y) = yi + (x3 + y)j y la trayectoria r(t) = t2 i + t3 j,
0 ≤ t ≤ 1. Sea C la traza de r, entonces la integral de lı́nea de F a lo largo de C es:
Z Z 1 Z 1 √
2 3 d 2 d 3
t3 , (t2 )3 + t3 · 2t, 3t2 dt =
F · dr = F(t , t ) · (t ), (t ) dt =
C 0 dt dt 0
Z 1 Z 1
t3/2 (2t) + (t6 + t3 )(3t2 ) dt = (2t5/2 + 3t8 + 3t5 dt =
=
0 0
4 7/2 1 9 1 6 1 4 1 1 59
= t + t + t = + + =
7 3 2 0 7 3 2 42
Para motivar la definición 3.1, se va a tomar en cuenta el concepto de trabajo realizado sobre un objeto que se mueve
en un campo de fuerza (fuerza gravitacional, fuerza electromagnética, etc). Considere una función F : Rn −→ Rn
que representa un campo de fuerza y suponga que una partı́cula se mueve en el espacio siguiendo una trayectoria
n
−→ R , n = 2, 3. Ası́, la fuerza ejercida sobre la partı́cula cuando se ubica en el punto r(t) es el vector
r : [a, b]
F r(t) . Si el movimiento de la partı́cula es rectilı́neo y la fuerza es constante, el trabajo realizado es el producto
escalar W = F · d, donde d es el desplazamiento recorrido por dicha partı́cula. En general, suponga que C es la
curva recorrida por la trayectoria r, luego para todo t y ∆t pequeño, el desplazamiento de la partı́cula desde r(t)
hasta r(t + ∆t) es ∆s = r(t + ∆t) − r(t) ≈ r0 (t)∆t. En consecuencia, el trabajo
realizado por la partı́cula entre los
instantes t y t + ∆t es aproximadamente igual a ∆W ≈ F r(t) · ∆s ≈ F r(t) · r0 (t)∆t. Ası́, una aproximación del
trabajo total realizado sobre todo el movimiento de la partı́cula a lo largo de la curva C es dada por la suma de
Xn
F r(ti ) · r0 (ti )∆t, donde se toma una partición a = t0 < t1 < t2 < . . . < tn−1 < tn = b del intervalo
Riemann
i=1
[a, b] en n partes iguales. A medida que se tome un entero n suficientemente grande, esta aproximación mejora
cada vez más; por tanto, al tomar el lı́mite de la suma de Riemann cuando n tiende a infinito, se obtiene que el
trabajo W realizado sobre la partı́cula cuando se mueve desde la posición r(a) hasta la posición r(b) a lo largo de
Z b Z
F r(t) · r0 (t) dt, es decir, la integral de lı́nea
C viene dado por la integral F · dr. En conclusión:
a C
6
El trabajo W realizado por el campo de fuerza F sobre Z una partı́cula cuyo movimiento es descrito por la
trayectoria r(t) con traza C es la integral de lı́nea W = F · dr.
C
Ejemplo 3.3. Calcular el trabajo realizado por el campo de fuerza F(x, y) = x2 i − xyj sobre una partı́cula que se
mueve a lo largo del cuarto de cı́rculo r(t) = cos ti + sen tj, 0 ≤ t ≤ π2 .
Z
x2 dx − xydy. Como r0 (t) = − sen ti + cos tj y F r(t) = F(cos t, sen t) = cos2 ti − cos t sen tj,
Se busca calcular
C
entonces el trabajo realizado es:
Z Z π
2
2
cos2 ti − cos t sen tj · − sen ti + cos tj dt =
W = x dx − xy dy =
C 0
Z π Z π π
2
2
2 2 2 2
(−2 cos2 t sen t)dt = cos3 t = −
= (cos t)(− sen t) + (− cos t sen t)(cos t) dt =
0 0 3 0 3
Si en cambio ϕ es una función monótona decreciente, entonces en particular ϕ(c) = b y ϕ(d) = a, por lo cual:
Z Z d Z a Z b Z
0 0
F · dr̂ = F r ϕ(t) · r ϕ(t) ϕ (t)dt = F r(u) · r(u)du = − F r(u) · r(u)du = − F · dr
C c b a C
Proposición 3.4. Sea r̂(t) = r ϕ(t) una reparametrización de la trayectoria r cuya traza es la curva C. Si
F es un campo vectorial cuyo dominio contiene a la curva C, entonces:
Z Z
i) F · dr̂ = F · dr, si ϕ es una función creciente.
C C
Z Z
ii) F · dr̂ = − F · dr, si ϕ es una función decreciente.
C C
Aquı́ el campo vectorial F viene dado por F(x, y) = y 2 i−x2 j. Observe que la curva cerrada C se puede escribir como
la concatenación de las tres curvas C1 , C2 y C3 , donde cada una de las curvas está parametrizada respectivamente
7
Figura 5. Curva del ejemplo 3.3 Figura 6. Curva cerrada del ejemplo 3.5
Fuente: Elaboración propia Fuente: Elaboración propia
por las trayectorias r1 (t) = ti, 0 ≤ t ≤ 2, r2 (t) = 2i + tj, 0 ≤ t ≤ 4 y r3 (t) = ti + t2 j, 1 ≤ t ≤ 2. En este caso, las
trayectorias r1 y r2 recorren las curvas C1 y C2 en el mismo sentido que se indica en la Figura 6, mientras que la
trayectoria r3 recorre la curva C3 en el sentido inverso que se indica en la Figura 6. Por la proposición anterior
se sabe que: Z Z Z Z
y 2 dx − x2 dy = F · dr1 + F · dr2 − F · dr3
C C1 C2 C3
De esta manera, basta calcular por separado cada una de estas integrales de lı́nea y luego hacer la suma mostrada
arriba. Ası́ se tiene que: Z Z 2 Z 2
2 2
F · dr1 = 0 i − t j · i dt = 0dt = 0
C1 0 0
Z Z 4 Z 4 4
t2 i − 22 j · j dt = −
F · dr2 = 4dt = −4t = −16
C2 0 0 0
2 2 t5 t4 2
Z Z Z
2 2
2
4 3 8
F · dr3 = (t ) i − (t) j · i + 2tj dt = (t − 2t )dt = − =−
C3 0 0 5 2 0 5
Z 8 72
Por lo tanto, y 2 dx − x2 dy = 0 + (−16) − − =−
C 5 5
Definición 4.1. Un campo vectorial F : U ⊂ Rn −→ Rn definido sobre un conjunto abierto U se dice que es
un campo conservativo si existe una función diferenciable f : U ⊂ Rn −→ R tal que F = ∇f .
Los campos conservativos resultan muy interesantes pues la integral de lı́nea de este tipo de campo vectorial a lo
largo de una trayectoria solo depende del punto inicial y final de la trayectoria y no de los puntos intermedios.
8
Esto es consecuencia directa del siguiente resultado, conocido como el teorema fundamental de las integrales de
lı́nea.
Teorema 4.2. Sea f : U ⊂ Rn −→ R una función diferenciable con derivadas parciales continuas en U . Si
r : [a, b] −→ Rn es una trayectoria cuya traza es la curva C ⊂ U , entonces:
Z
∇f · dr = f r(b) − f r(a)
C
Observe la semejanza de este resultado con el segundo teorema fundamental del cálculo para las integrales de
funciones reales de variable real. Para
la demostración de
este teorema basta aplicar la regla de la cadena a
d
f r(t) = ∇f r(t) · r0 (t). La palabra conservativo proviene de la fı́sica,
la función compuesta f r(t) pues dt
principalmente de la mecánica, en donde hay campos vectoriales F que representan la fuerza de un sistema fı́sico
que satisface el principio de conservación de la energı́a; en dicho caso, el campo F viene dado como el campo
gradiente F = −∇U, siendo U la energı́a potencial del sistema. Si F I es un campo conservativo definido en un
abierto U , entonces para toda curva cerrada C contenida en U se tiene F · dr = 0.
C
Z
Ejemplo 4.3. Calcular la integral de lı́nea F · dr, donde F(x, y, z) = 2xyi + (x2 + z 2 )j + 2yzk y C es la curva
C
que se muestra en la Figura 7, cuyo punto inicial es (1, 1, 0) y su punto final es (1, 2, 3). Observe que no se cuenta
con información sobre la curva C para escoger una parametrización adecuada de tal curva. Sin embargo, el campo
vectorial F es conservativo pues si se define f (x, y, z) = x2 y + z 2 y, entonces ∂f ∂f 2 2 ∂f
∂x = 2xy, ∂y = x + z y ∂z = 2zy,
por lo cual ∇f (x, y, z) = 2xyi + (x2 + z 2 )j + 2yzk = F(x, y, z). Ası́, por el teorema fundamental de las integrales
de lı́nea, indiferentemente de como sea la curva C siempre que su punto inicial sea (1, 1, 0) y su punto final sea
(1, 2, 3), se tiene que:
Z Z
2 2
∇f ·dr = f (1, 2, 3)−f (1, 1, 0) = 12 2+32 2 − 12 1+02 1 = 2+18−1−0 = 19
2xydx+(x +z )dy +2yzdz =
C C
Figura 7. Trayectoria del ejemplo 4.3 Figura 8. Dos regiones de R2 ; U1 es simplemente conexa
Fuente: Elaboración propia Fuente: Elaboración propia
No es fácil determinar a simple vista si un campo vectorial dado es conservativo o no. En el caso de un campo
vectorial en el plano R2 existe un criterio que permite identificar los campos vectoriales conservativos, siempre
que este campo esté definido sobre una región simplemente conexa: se dice que U ⊂ R2 es simplemente conexo si
toda curva cerrada contenida en U se puede “contraer” en U hasta reducirla a un punto, sin salirse de la región
U . Intuitivamente una región simplemente conexa no tiene “huecos”. La Figura 8 muestra dos regiones de R2 ,
la de la izquierda es una región simplemente conexa, mientras que la de la derecha no lo és; observe que la curva
cerrada naranja no se puede contraer a un punto sin salirse de la región.
9
Proposición 4.4. Sea U ⊂ R2 una región simplemente conexa y sea F : U −→ R2 , F(x, y) = P (x, y)i +
Q(x, y)j un campo vectorial con derivadas parciales continuas en U . El campo vectorial F es conservativo en
∂Q
U si, y solamanete si, ∂P
∂y = ∂x para todo punto en U .
Ejemplo 4.5. Determinar si el campo vectorial F(x, y) = (x2 − 2y 3 )i + (x + 5y)j es conservativo en el plano.
∂Q
En este caso, P (x, y) = x2 − 2y 3 y Q(x, y) = x + 5y, por lo cual, ∂P 2
∂y = −6y 6= 1 = ∂x . Por lo tanto, F no es
conservativo.
Ejemplo 4.6. Calcular el trabajo realizado por el campo de fuerza F(x, y) = −ye−xy i − xe−xy j sobre una partı́cula
2 2
cuando da un giro en sentido anti-horario sobre la elipse (x−1)
4 + (y+3)
25 = 1. En este caso, las funciones compo-
∂
−ye−xy = −e−xy − xye−xy = ∂x ∂
−xe−xy , luego el campo F es conservativo en R2 .
nentes de F satisfacen ∂y
(x−1)2 (y+3)2
Si C denota
I a la elipse 4 + 25 = 1, siendo C una curva cerrada, se concluye que el trabajo realizado es
W= F · dr = 0.
C
Observe que en la proposición 4.4, decir que el campo vectorial F(x, y) = P (x, y)i + Q(x, y)j en R2 es conservativo
equivale a afirmar que el rotacional de F es igual a cero, pues rotF(x, y) = ∂Q ∂P
∂x (x, y) − ∂y (x, y). Si F(x, y, z) =
P (x, y, z)i + Q(x, y, z)j + R(x, y, z)k es un campo vectorial en R3 , se tiene de la proposición 2.6 que si F es
conservativo, entonces rotF = 0. Resulta que el recı́proco de esta proposición es cierto, por lo cual se cuenta con
el siguiente criterio para identificar campos conservativos en el espacio R3 .
Proposición 4.7. Sea U ⊂ R3 una región simplemente conexa y sea F : U −→ R3 , F(x, y, z) = P (x, y, z)i +
Q(x, y, z)j + R(x, y, z)k, un campo vectorial con derivadas parciales continuas en U . El campo vectorial F es
conservativo en U si, y solamanete si, rot F (x, y, z) = 0 para todo punto (x, y, z) ∈ U . Esto equivale a las
∂Q ∂P ∂R ∂Q
identidades ∂P ∂R
∂y = ∂x , ∂z = ∂x , ∂z = ∂y .
5. Teorema de Green
Sea C ⊂ R2 una curva cerrada y simple, es decir, C está parametrizada por una trayectoria r : [a, b] −→ R2 , donde
r(a) = r(b) y r(t1 ) 6= r(t2 ), para todo a ≤ t1 < t2 < b. Ası́, C limita una región del plano D, como se indica
en la Figura 9. Se asumirá siempre que la curva C tiene orientación positiva, lo cual significa que se recorre a la
curva en sentido antihorario o, equivalentemente, que si una persona camina sobre la curva C, la región D estará
siempre a su izquierda. El teorema de Green permite convertir el cálculo de una integral de lı́nea de un campo
vectorial en el plano, a lo largo de una curva cerrada simple, en una integral doble sobre la región encerrada por
dicha curva.
10
Figura 9. Región D limitada por la curva cerrada Figura 10. Región D del ejemplo 5.2 limitada por
C positivamente orientada las curvas y = x2 y y = x3
Fuente: Elaboración propia Fuente: Elaboración propia
Teorema 5.1 (Teorema de Green). Sea C ⊂ R2 una curva cerrada simple con orientación positiva y suponga
que C es el borde de la región simplemente conexa D ⊂ R2 . Si F(x, y) = P (x, y)i + Q(x, y)j es un campo
vectorial definido sobre una región abierta U que contiene a D y cuyas derivadas parciales son continuas en
U , entonces: I ZZ
∂Q ∂P
F · dr = − dxdy (5)
C D ∂x ∂y
∂Q ∂P
Recuerde que ∂x − ∂y es el rotacional del campo vectorial F, luego otra forma de escribir la fórmula (5) del
I ZZ
teorema de Green es F · dr = rot F(x, y) dA. Esta igualdad es un caso particular del teorema de Stokes
C D
que se verá en el próximo Escenario.
I
Ejemplo 5.2. Calcular la integral de lı́nea (x2 − y 2 )dx + (2y − x)dy, donde C consiste en la frontera de la
C
región D limitada por las curvas y = x2 y y = x3 , ver la Figura 10. En este ejemplo el campo vectorial viene dado
por P (x, y) = x2 − y 2 y Q(x, y) = 2y − x. Observe que D es una región del tipo I; de hecho, todo punto (x, y) en
D cumple las desigualdades 0 ≤ x ≤ 1, x3 ≤ y ≤ y 2 . Como ∂Q ∂P
∂x = −1 y ∂y = −2y, entonces por el teorema de
Green se concluye que:
I ZZ Z 1 Z x2
(x2 − y 2 )dx + (2y − x)dy =
−1 − (−2y) dy dx = (2y − 1)dy dx
C D 0 x3
Z 1 x2 Z 1
2
= (y − y) 3 dx = (−x6 + x4 + x3 − x2 )dx
0 x 0
x7 x5 x4 x3 1 11
= − + + − =−
7 5 4 3 0 420
Ejemplo 5.3. Una partı́cula recorre en sentido antihorario el cı́rculo de radio 3 y centro (0, 0)y está sometida al
campo de fuerza F(x, y) = y 3 i + (x3 + 3xy 2 )j. Hallar el trabajo realizado por F sobre la partı́cula. Por el teorema
de Green, se tiene que el trabajo realizado por F es:
I ZZ ZZ ZZ
3 3 2 ∂ 3 2 ∂ 3 2 2 2
3x2 dA,
W= y dx + (x + 3xy )dy = x − 3xy − y dA = (3x + 3y ) − 3y dA =
C D ∂x ∂y D D
donde C es el cı́rculo x2 + y 2 = 9 y D es el disco limitado por dicho cı́rculo. Usando la fórmula del cambio de
11
variables para coordenadas polares se obtiene que:
ZZ Z 2π Z 3 Z 2π Z 3 Z 2π 4
2 2 3 2 r 2
3
3x dxdy = 3 (r cos θ) (r dr dθ) = 3 r cos θ dr dθ = 3 cos θ dθ =
D 0 0 0 0 0 4 0
Z 2π Z 2π
81 243 243 sen 2θ 2π
243π
= 3 cos2 θ dθ = (1 + cos 2θ)dθ == θ+ =
0 4 8 0 8 2 0 4
243π
Por lo tanto, el trabajo realizado es W = 4 .
Considere el campo vectorial en R2 dado por F(x, y) = P (x, y)i + Q(x, y)j, donde P (x, y) = − y2 y Q(x, y) = x2 . En
este caso, se cumple ∂Q ∂P 1 1
∂x − ∂y = 2 + 2 = 1, luego si D es la región del plano limitada por la curva cerrada simple
I ZZ ZZ
y x ∂Q ∂P
C, se tiene por el teorema de Green que − dx + dy = − dA = 1 dA = área de D.
C 2 2 D ∂x ∂y D
Por lo tanto, una consecuencia del teorema de Green es la siguiente proposición que expresa el área de una región
compacta del plano como una integral de lı́nea.
Proposición 5.4. Sea D ⊂ R2 una región simplemente conexa limitada por la curva cerrada simple C.
Entonces el área de D viene dada por:
I
1
A(D) = xdy − ydx (6)
2 C
x2 y2
Ejemplo 5.5. Hallar el área de la región D limitada por la elipse dada por la ecuación a2
+ b2
= 1.
Para usar la fórmula (6) se debe conseguir primero una parametrización de la elipse. Es fácil ver que la trayectoria
r : [0, 2π] −→ R2 definida por r(t) = a cos ti + b sen tj sirve como parametrización de la elipse. Ası́, si A denota el
área de D, entonces:
1 2π ab 2π ab 2π
I Z Z Z
1 d d
[cos2 t+sen2 t]dt =
A= xdy−ydx = (a cos t) (b sen t)−(b sen t) (a cos t) dt = dt = abπ
2 C 2 0 dt dt 2 0 2 0
Recuerde que los valores a y b corresponden a las longitudes del semieje mayor y menor de la elipse, por lo que
se concluye que el área de la región encerrada por una elipse es igual a π por el producto de las longitudes de sus
semiejes.
El teorema de Green se puede generalizar para regiones compactas D ⊂ R2 que están limitadas por un número
finito de curvas, como la que se muestra en la Figura 11. En este caso, suponga que el borde de D es la unión
disjunta de las curvas cerradas simples C1 , C2 , . . . , Ck , las cuales están orientadas de manera que una persona que
camina sobre la curva, siguiendo la orientación de cada curva Ci , 1 ≤ i ≤ k, tiene a la región D siempre a su
izquierda; se afirma, entonces, que la unión finita de las curvas tienen orientación positiva. Si se denota por C a
la unión C = C1 ∪ C2 ∪ . . . ∪ Ck , la integral de lı́nea de un campo vectorial F a lo largo de C se define como la
suma de las integrales de lı́nea de F a lo largo de las curvas C1 , C2 , . . . , Ck , es decir,
I I I I
F · dr = F · dr1 + F · dr2 + . . . + F · drk ,
C C1 C2 Ck
donde ri denota una parametrización de la curva Ci , para cada 1 ≤ i ≤ k. Si F(x, y) = P (x, y)i + Q(x, y)j y las
funciones P y Q tienen derivadas parciales continuas en un abierto que contiene a la región D, entonces:
I ZZ
∂Q ∂P
F · dr = − dA
C1 ∪C2 ∪...∪Ck D ∂x ∂y
12
Figura 11. Región D cuyos bordes son las curvas Figura 12. Curva C del ejemplo 5.6. C 0 es un
cerradas C1 , C2 y C3 cı́rculo de radio r pequeño
Fuente: Elaboración propia Fuente: Elaboración propia
∂ −y i
I I I I I ZZ h
∂ x
F · dr − F · dr = F · dr + F · dr = F · dr = 2 2
− dA = 0
C C0 C −C 0 C∪−C 0 D ∂x x + y ∂y x2 + y 2
I I
Por lo tanto, F · dr = F · dr. Si se parametriza al cı́rculo C 0 por la trayectoria r cos ti + r sen tj, 0 ≤ t ≤ 2π,
C C0
entonces:
2π
−r sen t
I Z
d r cos t d
F · dr = (r cos t) + (r sen t) dt
C0 0 (r cos t)2 + (r sen t)2 dt (r cos t)2 + (r sen t)2 dt
Z 2π
(−r sen t)(−r sen t) (r cos t)(r cos t)
= + dt
0 r2 r2
Z 2π Z 2π
2 2
= (sen t + cos t)dt = dt = 2π
0 0
I
Ası́ se concluye que F · dr = 2π, para toda curva C con orientación positiva que rodea al origen.
C
6. Ejercicios
13
x2
Z
c) ds, C es el gráfico de la función y = 32 x2/3 , 1 ≤ x ≤ 8.
C y3
3. Hallar la masa y el centro de masa del alambre con densidad de masa δ(x, y) y parametrizada por la
trayectoria r descrita.
a) δ(x, y, z) = x2 + y 2 y r(t) = et cos ti + et sen tj + tk , 0 ≤ t ≤ π.
√ √
b) δ(x, y) = x2 y r(t) = ti + log tj , 3 ≤ t ≤ 8.
√
c) δ(x, y) = y y r(t) = 4(t − sen t)i + 4(1 − cos t)j , 0 ≤ t ≤ 2π.
d) δ(x, y, z) es constante y r(t) = a cos ti + a sen tj + btk , 0 ≤ t ≤ 3π, a, b > 0 son constantes.
5. Calcular el rotacional de los siguientes campos vectoriales. Determinar cuáles son conservativos.
a) F(x, y, z) = (2y−z)i+ez j+xyzk b) F(x, y, z) = x sen yi−y cos xj+yz 2 c) F(x, y, z) = y 2 zi+2xyzj+xy 2 k
2 +y 2 2 +z 2 xy
d) F(x, y, z) = ex i + ey j + xyzk e) F(x, y, z) = y log zi − x log zj + z k
7. Calcule el trabajo total realizado al mover una partı́cula a lo largo de la curva C si el movimiento es
ocasionado por el campo de fuerza F.
a) F(x, y) = 2xyi + (x2 + y 2 )j, C es el segmento de recta que va desde el origen al punto (1, 1).
2 x2 y+y
b) F(x, y) = y x+1
3 i − x2
, C es cualquier curva diferenciable que va del punto (2, 3) hasta el punto (1, 0),
sin cruzar el eje x.
c) F(x, y) = (x − y)i + xyj, C es el arco del cı́rculo x2 + y 2 = 4, desde el punto (2, 0) hasta el punto (0, −2),
en sentido antihorario.
d) F(x, y, z) = (xyz + x)i + (x2 z + y)j + (x2 y + z)k, C es la traza de r(t) = ti + t2 j + t3 k, 0 ≤ t ≤ 2.
7π
e) F(x, y, z) = sen zi + cos zj + (xy)1/3 k, C es la traza de r(t) = cos3 ti + sen3 tj + tk, 0 ≤ t ≤ 2 .
8. Evaluar las siguientes integrales de lı́nea mediante dos métodos: a) De manera directa. b) Usando el teorema
de Green. Suponga siempre que la orientación de la curva C es positiva.
I
a) (2y −x)dx+(5x−y)dy, C es la frontera de la región limitada por la recta y = x y la parábola y = x2 −x.
C
14
I
b) (ex − 2x2 y)dx + (ey + 2xy 2 )dy, C es el cı́rculo x2 + y 2 = 9.
C
I
c) (2x + ex sen y)dx + (x + ex cos y)dy, C es el trapecio de vértices (0, −2), (4, −1), (−2, 3) y (0, 2).
C
15
INFORMACIÓN TÉCNICA
16