Docsity
Docsity

Prepara tus exámenes
Prepara tus exámenes

Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity


Consigue puntos base para descargar
Consigue puntos base para descargar

Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium


Orientación Universidad
Orientación Universidad


Integrales múltiples, Apuntes de Cálculo

Asignatura: calculo, Profesor: Pepe Aranda, Carrera: Física, Universidad: UCM

Tipo: Apuntes

2014/2015

Subido el 07/12/2015

avecont
avecont 🇪🇸

4.5

(15)

2 documentos

1 / 8

Toggle sidebar

Esta página no es visible en la vista previa

¡No te pierdas las partes importantes!

bg1
4. Integrales de línea
4.1 Integrales de campos escalares a lo largo de curvas
y
x
c
(a)c
[]
ab
t
(b)c
(t)c
(t)c'
Una función vectorial (también llamada trayectoria ocamino)
era una c:[a,b]RRn, cuya gráfica c(t)=(x1(t),...,xn(t))
es una curva Cen Rny su derivada c0(t) = (x0
1(t),...,x0
n(t))
describe el vector tangente a la curva (el vector velocidad, si c(t)
está describiendo un movimiento en Rn).
La recta tangente a Cen c(to)era: l(t)= c(to)+(tto)c0(to).
cse dice C1si es continua y c0existe y es continua t(a,b). Es C1a trozos si Ces continua
y[a,b]se puede dividir en un número finito de subintervalos en cada uno de los cuales ces C1
[su gráfica será entonces una curva continua sin recta tangente en un número finito de puntos].
y
x
(0)c
(π/2)c
c*
Ej 1. c :0,π
2R2con c(t)=(cost,sent)es una trayectoria C1pues
c0(t)= (sen t,cost)existe t0,π
2.
c:[0,1]R2con c(t)= t,1t2,c0(t)= 1,t(1t2)1/2
es otro camino C1, que describe esa misma curva en sentido opuesto.
Se dice que cycson dos parametrizaciones de la misma curva C.
Sea c:[a,b]Rnun camino C1y sea fun campo escalar en Rntal que fc(t)
es continua en [a,b]. La integral de fa lo largo de cse define:
Zc
f ds Zb
a
fc(t)kc0(t)kdt .
Si c(t)es solamente C1a trozos o si f(c(t)) es continua a trozos, definimos Rcf ds
descomponiendo [a,b]en intervalos sobre los que f(c(t))kc0(t)ksea continua y
sumando las integrales sobre cada uno de ellos.
Ej 2. Si f(x,y) = xy2yc,cson los de arriba, las integrales a lo largo de las dos trayectorias son:
Zc
f ds =Zπ/2
0costsen2t·1·d t =1
3sen3tπ/2
0=1
3.
Zc
f ds =Z1
0t(1t2)·1
1t2·dt =1
3(1t2)3/21
0=1
3.
No es casualidad que ambas integrales coincidan. Probaremos en 4.2 que:
Teor 1. Si cycdescriben la misma curva C, entonces Rcf ds =Rcf d s RCf d s .
Como la integral de línea de una fescalar no depende de la parametrización, sólo de la
curva, es lícita la notación RCf ds (integral de fsobre C) donde cno aparece por ningún lado.
Interpretemos esta integral. Sea primero f1 . Si pensamos que c(t)describe una partícula, al
ser kc0(t)kla velocidad escalar en el instante t, parece claro que ds =kc0(t)kdt (‘diferencial
de arco’) representa la distancia recorrida en un ‘diferencial de tiempo dt y por tanto:
L=RCds =Rb
akc0(t)kdt representa la longitud de la curva Cdefinida por c.
1
y
x
1
–1
Ej 3. Hallemos la longitud de la curva descrita por c(t) = t2,t3,t[1,1].
kc0(t)k=4t2+9t4L=Z1
1|t|4+9t2dt =2
27 (4+9t2)3/21
0=2(133/28)
27 2.88 .
O con otra parametrización de la misma curva (y usando su simetría):
(x,x3/2),x[0,1]L=2Z1
0q1+9x
4dx =16
27 1+9x
43/21
0=16
27 133/2
81.
[Cuando c0=0 pueden aparecer picos en las curvas, pues el vector tangente no tiene dirección definida].
[Estas integrales de campos escalares, por incluir una raíz, muchas veces llevan a primitivas no calculables].
33
pf3
pf4
pf5
pf8

Vista previa parcial del texto

¡Descarga Integrales múltiples y más Apuntes en PDF de Cálculo solo en Docsity!

4. Integrales de línea

4.1 Integrales de campos escalares a lo largo de curvas

y

x

c

c (a)

[ ] a b

t

c (b)

c (t)

Una función vectorial (también llamada trayectoria o camino) c' (t)

era una c : [a, b] ⊂ R → R

n

, cuya gráfica c(t) = (x 1 (t),... , xn(t))

es una curva C en R

n

y su derivada c

(t) = (x

1 (t),... ,^ x

n(t))

describe el vector tangente a la curva (el vector velocidad, si c(t)

está describiendo un movimiento en R

n

[

La recta tangente a C en c(to) era: l(t) = c(to)+(t−to) c

(to)

]

c se dice C

1

si es continua y c

existe y es continua ∀t ∈(a, b). Es C

1

a trozos si C es continua

y [a, b] se puede dividir en un número finito de subintervalos en cada uno de los cuales c es C

1

[su gráfica será entonces una curva continua sin recta tangente en un número finito de puntos].

y

x c (0)

c (π/2)

c *

Ej 1. c :

[

0 , π 2

]

R^2 con c(t) = (cost, sent) es una trayectoria C^1 pues

c′(t) = (− sent, cost) existe ∀t ∈

π 2

c∗ : [ 0 , 1 ] → R 2 con c∗(t) =

t,

1 −t^2

, c∗′(t) =

1 , −t( 1 −t^2 )−^1 /^2

es otro camino C^1 , que describe esa misma curva en sentido opuesto. Se dice que c y c∗ son dos parametrizaciones de la misma curva C.

Sea c : [a, b] → R

n

un camino C

1

y sea f un campo escalar en R

n

tal que f

c(t)

es continua en [a, b]. La integral de f a lo largo de c se define:

c

f ds ≡

∫ (^) b

a

f

c(t)

‖ c

(t)‖ dt.

Si c(t) es solamente C

1

a trozos o si f (c(t)) es continua a trozos, definimos

c f ds

descomponiendo [a, b] en intervalos sobre los que f (c(t)) ‖ c

(t)‖ sea continua y

sumando las integrales sobre cada uno de ellos.

Ej 2. Si f (x, y) = xy 2 y c , c∗ son los de arriba, las integrales a lo largo de las dos trayectorias son: ∫

c

f ds =

∫ (^) π/ 2

0

cost sen^2 t · 1 · dt = 1 3 sen^3 t

]π/ 2 0

1 3

c∗

f ds =

∫ (^1)

0

t ( 1 −t 2 ) ·

√^1

1 −t^2

· dt = − 1 3 ( 1 −t 2 ) 3 / 2

] 1

0

1 3

No es casualidad que ambas integrales coincidan. Probaremos en 4.2 que:

Teor 1. Si c y c∗ describen la misma curva C , entonces

c

f ds =

c∗

f ds ≡

C

f ds.

Como la integral de línea de una f escalar no depende de la parametrización, sólo de la

curva, es lícita la notación

C f ds^ (integral de^ f^ sobre^ C^ ) donde^ c^ no aparece por ningún lado.

Interpretemos esta integral. Sea primero f ≡1. Si pensamos que c(t) describe una partícula, al

ser ‖c

(t)‖ la velocidad escalar en el instante t , parece claro que ds = ‖ c

(t)‖ dt (‘diferencial

de arco’) representa la distancia recorrida en un ‘diferencial de tiempo dt ’ y por tanto:

L =

C ds^ =^

∫ (^) b

a ‖^ c

(t)‖ dt representa la longitud de la curva C definida por c.

1

y

x

1

Ej 3. Hallemos la longitud de la curva descrita por c(t) =

t 2 , t 3

, t ∈ [− 1 , 1 ].

‖c′(t)‖ =

4 t^2 + 9 t^4 ⇒ L =

∫ (^1)

− 1

|t|

4 + 9 t^2 dt = 2 27 ( 4 + 9 t^2 )

3 / 2 ] 1

0

2 ( 133 /^2 − 8 ) 27

O con otra parametrización de la misma curva (y usando su simetría):

(x, x 3 / 2 ) , x ∈ [ 0 , 1 ] ⇒ L = 2

∫ (^1)

0

9 x 4 dx^ =^

16 27

9 x 4

) 3 / 2 ] 1

0 =^

16 27

[

133 /^2 8 −^1

]

[Cuando c′^ = 0 pueden aparecer picos en las curvas, pues el vector tangente no tiene dirección definida].

[Estas integrales de campos escalares, por incluir una raíz, muchas veces llevan a primitivas no calculables].

y

z f ( x (t), y (t))

x c (t)=( x (t), y (t))

Si ahora f es cualquier campo con f (c(t))≥ 0 ∀t ∈[a, b] ,

c f ds^ representa para^ n=2 el^ área de la valla^ de altura

f (x, y) en cada (x, y) de la curva C , pues un ‘diferencial

de valla’ tiene área f (c(t)) ds = f (c(t)) ‖ c

(t)‖ dt.

Tiene otra posible interpretación cuando n = 2 o n = 3 :

Si c(t) describe un alambre de densidad variable dada por ρ(x) , la masa del alambre será

M =

C ρ^ ds^ (tanto en el plano como en el espacio).

El centro de gravedad (centroide si ρ constante) del alambre, tendrá por coordenadas:

x =

1 M

C x^ ρ^ ds^ ,^ y^ =^

1 M

C y^ ρ^ ds^ ;^ si^ n^ =^ 3 , además^ z^ =^

1 M

C z^ ρ^ ds^.

Estas integrales también sirven para hallar, el valor medio de una f sobre C : f =

1 L

C f ds^.

y x

z

Ej 4. Sea el alambre en forma de hélice: c : [ 0 , 2 π] → R 3 , t → (cost, sent,t)

y de densidad ρ(x, y, z) = x^2 +y^2 +z^2. Como ‖ c′‖ =

2 , se tiene que:

Su longitud es M =

∫ (^2) π 0

2 dt = 2 π

Su masa es M =

∫ (^2) π 0 (cos

2 t+ sen 2 t+t 2 ) dt =

2 π + 8 3 π

3

≈ 125. [ Su densidad media es, por tanto: 1+ 4 3 π^2 ≈ 14.

] .

Su centro de gravedad: x = 1 M

∫ (^2) π 0 cost^ (^1 +t

2 dt = 6 3 + 4 π^2

y = 1 M

∫ (^2) π 0 sent^ (^1 +t

2 )

2 dt = − 6 π 3 + 4 π^2 ≈ − 0.44 , z = 1 M

∫ (^2) π 0 t^ (^1 +t

2 )

2 dt = 3 π 1 + 2 π^2 3 + 4 π^2

Su centroide: x = 1 L

∫ (^2) π 0 cost^

2 dt = 0 , y = 1 L

∫ (^2) π 0 sent^

2 dt = 0 , z = 1 L

∫ (^2) π 0 t^

2 dt = π.

a (^) b

Otra aplicación de este tipo de integrales. Si una curva C viene ds

dada por x=x(t) , y=y(t)≥0 , t ∈[a, b] , el área de la superficie

de revolución obtenida al hacer girar C en torno al eje y = 0 es:

A = 2 π

∫ (^) b a

y ds = 2 π

∫ (^) b

a

y(t)

x′(t)

y′(t)

dt

[

pues el área de la banda de anchura ds es ≈ 2 πy ds

]

En el caso particular de que C sea la gráfica de y = y(x) queda: A = 2 π

∫ (^) b

a

y

1 +(y′)^2 dx.

y=√a –x^2^2

Ej 5. Hallemos el área de la superficie esférica de radio a.

Se puede ver como el giro de y =

a^2 −x^2 en torno a y = 0 , y por tanto:

A = 2 π

∫ (^) a

−a

a^2 −x^2

x^2 a^2 −x^2 dx = 2 π

∫ (^) a −a a dx^ =^4 πa

2 .

O escribiendo c(t) = (a cost, a sent) , t ∈ [ 0 , π] , y con la primera fórmula:

ds = a dt , A = 2 π

∫ (^) π

0

a sent a dt = 4 πa 2 .

Veamos lo que sucede con la integral al hacer un cambio en el parámetro que describe la curva.

Sea c : [a, b]→ R

n

y h : [a 1 , b 1 ]→[a, b] una biyección C

1

. Si llamamos p = c◦h : [a 1 , b 1 ]→ R

n

h

a (^) b

h

a (^) b

a a

b b

(^1 1 1 )

las trayectorias p(u) , u ∈ [a 1 , b 1 ] y c(t) , t ∈ [a, b] ,

describen la misma curva, en el mismo sentido o en

sentido opuesto según sea, respectivamente,

h(a 1 ) = a y h(b 1 ) = b o h(a 1 ) = b y h(b 1 ) = a

[se dice que la reparametrización conserva o invierte la orientación de la curva].

Teor 1.

Si c y p describen la misma curva C , entonces según c y p lo hagan en el

mismo sentido o en el opuesto se tiene, respectivamente:

c

f · ds =

p

f · ds , o bien

c

f · ds = −

p

f · ds

Si C es C

1

(si no, dividimos y sumamos las integrales), como p

(u) = c

(h(u)) h

(u) es:

p f^ ·^ ds^ =^

∫ (^) b 1

a 1 f

c(h(u))

· c

(h(u)) h

(u) du e

c f^ ·^ ds^ =^

∫ (^) b

a f

c(t)

· c

(t) dt.

Haciendo en la integral de la izquierda h(u) =t , se tiene + o − la integral de la derecha,

según se conserve o no la orientación, pues no cambian o sí los límites de integración.

Si estuviésemos integrando un campo escalar f , se tendría siempre que:

p f ds^ =^

∫ (^) b 1

a 1 f^

c(h(u))

‖c

(h(u))‖ |h

(u)| du =

∫ (^) b

a f^

c(t)

‖c

(t)‖ dt =

c f ds^ ,

y con esto tenemos probado el teorema de la sección anterior.

La integral de línea de un campo vectorial sólo depende de la curva C y el sentido en que

se recorre (la un campo escalar sólo de C ). Podemos elegir las c más sencillas.

c

Ej 5. Tiene un sentido preciso hablar de la integral de f(x, y) = (y, 0 ) a lo largo de la circunferencia unidad recorrida en el sentido de las agujas del reloj. [ Las integrales sobre curvas cerradas suelen representarse con el símbolo

∮ ]

Eligiendo c(t) = (cost, − sent) , t ∈ [ 0 , 2 π]

o [−π, π] , o...

) , ∮

c

f · ds =

∫ (^2) π 0 (−^ sent,^0 )^ ·^ (−^ sent,^ −^ cost)^ dt^ =^ π^.

Con cualquier parametrización que proporcione el mismo sentido se llega a lo mismo.

Ej 6. Calculemos la integral de línea del campo vectorial f(x, y, z) =

z , ey, xy

desde ( 1 , 0 , 0 ) hasta ( 1 , 2 , 3 ) a lo largo del segmento que une los puntos.

2

3

y

z

a b

a + (ba) t = c ( t )

Hay muchas formas de parametrizar un segmento en el espacio. Para este,

con x = 1 constante, una c casi salta a la vista: c(t) = ( 1 , 2 t, 3 t) , t ∈ [ 0 , 1 ]. [ O también, con el dibujo de la derecha: c∗(y) =

1 , y , 3 2 y

, y ∈ [ 0 , 2 ]

]

En general, vimos en 1.1 que c(t) = a +(b−a)t , t ∈ [ 0 , 1 ] da el segmento

que une a y b [es una recta, si t = 0 estamos en a y si t = 1 en b ].

Hallemos ya la integral pedida (con la primera de las parametrizaciones dadas): ∫

c

f · ds =

∫ (^1)

0

3 t, e^2 t^ , 2 t

· ( 0 , 2 , 3 ) dt =

∫ (^1)

0

2e^2 t^ + 6 t

dt = 2e^2 + 2.

Ej 7. Hallar el trabajo realizado por el campo de fuerzas f(x, y, z) = (y−z, z−x, x−y) a lo largo de la curva C intersección de la esfera x 2 +y 2 +z 2 = 4 y el plano z = y , si C es recorrida de modo que, vista desde las z positivas, el sentido es contrario a las agujas del reloj.

y x

z (^) z=y

c

Sobre C , como z = y , es x 2

  • 2 y 2 = 4 (elipse).

c(t) =

2 cost,

2 sent,

2 sent

, t ∈ [ 0 , 2 π].

y

x

2

T =

∫ (^2) π 0 (^0 ,^

√ 2 sent−2 cost, 2 cost−

√ 2 sent) · (−2 sent,

√ 2 cost,

√ 2 cost) dt

=

∫ (^2) π 0 0 dt^ =^0 (el campo es perpendicular a la trayectoria).

[Se verá más adelante que serán 0 las integrales sobre un camino cerrado de los campos que sean gradientes de un campo escalar; pero, a pesar de ser

∮ C =^ 0 , este campo no lo será, pues rot^ f^ = (−^2 ,^ −^2 ,^ −^2 )^ ,^0 ].

4.3. Integrales de gradientes y teorema de Green

Generalizamos el 2

teorema fundamental del cálculo infinitesimal

∫ (^) b

a g

(x) dx = g(b)−g(a).

Teor 1.

Sea U : R

n

→ R un campo escalar C

1

y c : [a, b] → R

n

un camino C

1

a trozos.

Entonces:

c

∇U · ds = U

c(b)

−U

c(a)

Si c ∈C

1

(si no, dividimos),

∫ (^) b

a ∇U(c(t))·^ c

(t) dt =

∫ (^) b a

U ◦ c)

(t) dt = U(c(b))−U(c(a)).

regla de la cadena

a

b

Por tanto, la integral de línea de un gradiente no depende del camino,

sólo del punto inicial y final. Si sabemos que un campo es un gradiente, su

integral es inmediata. Además, si c(a) = c(b) , la curva descrita por c es

cerrada e

c ∇U^ ·^ ds^ =^ 0 :^ la integral de línea de un gradiente sobre una curva cerrada es^ 0.

Si un campo vectorial f es gradiente de alguna función U, a U se le llama

función potencial para f , y el campo f se dice conservativo.

¿Cómo saber si f es conservativo? Una condición necesaria sencilla para n = 2 y n = 3 es:

Teor 2.

Si f(x, y) =

f (x, y), g(x, y)

de C

1

es conservativo ⇒ fy ≡ gx.

Si f(x, y, z) =

f (x, y, z), g(x, y, z), h(x, y, z)

de C

1

es conservativo ⇒ rot f = 0.

Si f = ( f , g) = (Ux,Uy) = ∇U , con U∈C

2

, por la igualdad de Schwartz de las derivadas

cruzadas [ Uxy = Uyx ] debe ser fy ≡ gx. Y el mismo argumento se aplica si n = 3.

Si las derivadas cruzadas no coinciden, no puede f ser gradiente. Si son iguales, muchas veces

es sencillo hallar una U tal que ∇U = f [aunque la implicación ⇐ no sea cierta en general].

Ej 1. Sea f(x, y)=(y^2 , 2 xy). Hallemos la integral entre ( 0 , 0 ) y ( 1 , 1 ) a lo largo de diferentes curvas:

a) la recta que une los puntos, b) la parábola y = x^2 , c) la circunferencia x^2 +y^2 = 2 x.

Posibles parametrizaciones: a) ca =(t,t) , b) cb =(t,t^2 ) , c) cc =

t,

2 t−t^2

, con t ∈[ 0 , 1 ] todas.

y

x 1

1 c c

c (^) b

a

c

Las integrales en cada caso son:

a)

∫ (^1) 0 (t

2 , 2 t 2 ) · ( 1 , 1 ) dt =

∫ (^1) 0 3 t

2 dt = 1 ,

b)

∫ (^1) 0 (t

(^4) , 2 t (^3) ) · ( 1 , 2 t) dt = ∫ (^1) 0 5 t

(^4) dt = 1 ,

c)

∫ (^1)

0

( 2 t−t 2 , 2 t

2 t−t^2 ) ·

√^1 −t 2 t−t^2

dt =

∫ (^1) 0 (^4 t−^3 t

2 ) dt = 1.

Como se cumple ∂ ∂ y (y

2 ) = 2 y = ∂ ∂ x (^2 xy)^ , esto nos hace sospechar que^ f^ es campo conservativo. Es fácil en este caso identificar una función potencial:

Si Ux = y 2 , debe ser U = xy 2 +p(y) para alguna función p

Si Uy = 2 xy , debe ser U = xy 2 +q(x) para alguna función q

⇒ U(x, y) = xy 2 .

[ U = xy 2 +C para cualquier constante C es también potencial, desde luego

]

Por tanto, las parametrizaciones y cálculos de integrales anteriores han sido inútiles, puesto que la integral a lo largo de cualquier trayectoria debía valer U( 1 , 1 ) −U( 0 , 0 ) = 1 − 0 = 1.

Por no depender del camino, hay otras formas de calcular una U : eligiendo caminos sencillos que unan el origen con el punto (x, y) y evaluando la integral de línea. Por ejemplo, para el C 1 a trozos:

c(t) =

(t, 0 ) , t ∈ [ 0 , x] (x,t) , t ∈ [ 0 , y]

c

f · ds =

∫ (^) x

0

( 0 , 0 )·( 1 , 0 ) dt +

∫ (^) y

0

t^2 , 2 xt

·( 0 , 1 ) dt = xy^2.

Haciendo lo mismo para un f = ( f , g) general, obtendríamos:

∫ (^) x

0

f (t, 0 ) dt +

∫ (^) y

0

g(x,t) dt = U(x, y).

Veamos qué campos de los ejemplos de la sección anterior son conservativos. El (x 2 , 2 y+x) del 1 no lo es por ser fy = 0. 1 = gx. El (y, x− 4 y) del 2 sí lo es: 1 ≡ 1 y U = xy− 2 y 2 es su potencial. No lo es el del 3, (zx 2 , xy, y 3 ) por no ser rot f ≡ 0. El f = d del 4 lo es

[

con U = (d 1 x, d 2 y, d 3 z)

]

El (y, 0 ) del 5 no: no podía serlo pues vimos que su integral sobre un camino cerrado era no nula y también lo asegura 1.0. Para 6 y 7 el rot f es, respectivamente, (x, 1 −y, 0 ) y (− 2 ,− 2 ,− 2 ) con lo que no pueden derivar de una función potencial U.

Teoremas de Green y de la divergencia

Son teoremas que relacionan integrales dobles e integrales de línea sobre curvas cerradas

en el plano. Veremos otros sobre integrales de campos vectoriales en el espacio (estos se pueden

considerar casos particulares de aquellos), tratando las integrales de superficie.

curva no simple

curva simple

curva cerrada simple

Una curva simple es la imagen de un camino C

1

a trozos inyectivo.

Una curva cerrada simple en R

2

será la imagen de un c : [a, b]→ R

2

inyectivo en [a, b] y tal que c(a) = c(b). Una curva de estas puede

recorrerse en dos sentidos diferentes; para indicar que una integral se

recorre en sentido antihorario escribiremos.

Teorema de Green:

Sea D⊂ R

2

limitado por ∂ D curva cerrada simple y el campo f =( f , g)∈C

1

(D). Entonces:

"

D

[ gx − fy ] dx dy =

∂ D

f dx + g dy ≡

∂ D

f · ds.

[

Obsérvese que si f es conservativo, el teorema de Green dice 0 =

∮ f · ds como debía ser

]

a

( x )

c ( x )

b

y

x

d

D

El el caso de que D fuese del primer tipo de recintos que consideramos en las integrales dobles: D =

(x, y) : a ≤ x ≤ b, c(x) ≤ y ≤ d(x)

se tendría:

D fy dx dy = −

∫ (^) b a

∫ (^) d(x) c(x) fy(x, y) dy dx =

∫ (^) b a

[

f (x, c(x))− f (x, d(x))

]

dx.

Parametrizando los cuatro tramos de la frontera:

(x, c(x)) , x ∈ [a, b] ; (b, y) , y ∈ [c(b), d(b)] , (x, d(x)) , x ∈ [a, b] ; (a, y) ; y ∈ [c(a), d(a)] →

∂ D

f dx =

∫ (^) b a f^ (x,^ c(x))·^1 dx^ +^

∫ (^) d(b) c(b) f (y, b)· 0 dy −

∫ (^) b a f^ (x,^ d(x))·^1 dx^ −

∫ (^) c(b) c(a) f (y, a)· 0 dy = −

"

D

fy.

Análogamente, para recintos D =

(x, y) : c ≤ y ≤ b, a(y) ≤ x ≤ b(y)

se ve que ∂ D g dy =

D gx.

Se cumplirá, pues ∂ D f dx + g dy =

D [gx − fy] en un recinto D que sea de los dos tipos anteriores. y

x

D

Dividiendo cualquier recinto general D en otros de estos últimos y teniendo en cuenta que se cancelan las integrales de línea sobre segmentos comunes, por ser recorridos en sentido opuesto, se llega al resultado.

c

Ej 5. Encontremos, usando Green, el valor π de la integral de f(x, y) = (y, 0 ) sobre la circunferencia unidad calculado ya en el ejemplo 5 de 4.2.

Como el sentido de recorrido es opuesto al que pide Green y gx − fy = −1 es:

∂ D

f · ds = −

"

D

(− 1 ) dx dy =

"

D

dx dy = área de D = π · 12 = π.

y

2 x

π/

(0, t ) (^) (2, t )

( t ,π/2)

( t ,0)

Ej 6. Hallemos ∂ D

e x sen y dx+ e 2 x cos y dy siendo D el rectángulo [ 0 , 2 ]×[ 0 , π 2

].

∫ (^2)

0

( 0 · 1 + 0 ) dt +

∫ (^) π/ 2

0

( 0 +e^4 cost · 1 ) dt −

∫ (^2)

0

(et^ · 1 + 0 ) dt −

∫ (^) π/ 2

0

( 0 +cost · 1 ) dt

Con Green:

∫ (^) π/ 2

0

∫ (^2)

0

(2e 2 x cos y− e x cos y) dx dy = [sen y]

π/ 2 0

[

e 2 x − e x

] 2

0

|| e 4 − e 2 .

Normalmente utilizaremos Green para reducir integrales de línea a integrales dobles, que suelen ser más

sencillas. Pero en el siguiente ejemplo procedemos al contrario.

Ej 7. Calculemos el área encerrada por la ‘hipocicloide’ x 2 / 3

  • y 2 / 3 = a 2 / 3 .

a

a

–a

–a

Green dice que: A =

D dxdy = 1 (^2) ∂ D x dy^ −^ y dx^ , pues^ gx^ −^ fy^ =^ 2.

c(θ ) = (a cos 3 θ , a sen 3 θ ) , θ ∈ [ 0 , 2 π] es la mejor parametrización.

A =

1 2

∫ (^2) π

0

[

(a cos 3 θ )( 3 a sen 2 θ cos θ ) − (a sen 3 θ )(− 3 a cos 2 θ sen θ )

]

3 2 a

2 ∫ (^2) π 0 sen

2 θ cos 2 θ dθ = 3 16 a

2 ∫ (^2) π 0 [^1 −cos 4θ^ ]^ dθ^ =^

3 8 πa

2 .

[los picos de la curva aparecen para c′^ =0 ].

Del teorema de Green se obtiene fácilmente:

Teorema de la divergencia en el plano:

Sean D⊂ R

2

limitado por ∂ D curva cerrada simple, f : D → R

2

campo vectorial C

1

y n el vector normal unitario exterior a ∂ D. Entonces

"

D

div f dx dy =

∂ D

f · n ds.

n

( x ', y ')

D

D

Si ∂ D viene dada por c(t) =

x(t), y(t)

, la normal es n =

(y′(t),−x′(t)) ‖c′(t)‖

Si f = ( f , g) ,

∂ D

f · n ds =

∫ (^) b a

[

f (x(t), y(t)) y′(t) − g(x(t), y(t)) x′(t)

]

dt

∂ D

f dy−gdx

Green

D

( fx +gy) dx dy.

D

C n

n F F

[Imaginemos una curva cerrada C sobre la superficie de un fluido y sea F = f v , donde f es la densidad del fluido y v su velocidad. Entonces ∮

C

F · n ds mide el ritmo con el que el fluido entra o sale de D.

Si la cantidad de fluido en D disminuye (aumenta) será

∮ C <^0

C >^0

. La integral coincide con !

D div F. Por tanto, la div F describe la tendencia del fluido a acumularse o dispersarse].

Ej 8. Comprobemos este teorema para f(x, y) = ( 7 , y 2 − 1 ) en el semicírculo r ≤ 3 , 0 ≤ θ ≤ π :

C

C

D

div f = 2 y ,

"

D

2 y dx dy =

∫ (^) π

0

∫ (^3)

0

2 r 2 sen θ dr dθ = 36.

Para C 1 , si c(x) = (x, 0 ) , x ∈ [− 3 , 3 ] , n = ( 0 ,− 1 ) ,

C 1

( 1 −y^2 ) ds =

∫ (^3) − 3 dx^ =^ 6.

Para C 2 , si c(t) = (3 cost, 3 sent) , t ∈ [ 0 ,π] , ‖c′(t)‖ = 3.

Como n = (cost, sent) ,

C 2

f · n ds = 3

∫ (^) π

0

7 cost+9 sen 3 t−sent

dt = 30. 36 = 6 +.

Ej 9. Comprobemos los teoremas de Green y de la divergencia para el campo f(x, y)=( x^3 , x^2 y ) y el recinto D del primer cuadrante acotada por y = 2 x e y = x^2. y (1,2)

x

y =2 x

c (^) 1

c (^) 2

y = x^2

n 1

n 2 (1,2 x )

Green:

"

D

[gx − fy] =

∫ (^2)

0

∫ (^2) x

x^2

2 xy dy dx =

∫ (^2)

0

( 4 x 3 −x 5 ) dx =

[

x 4 − x^6 6

] 2

0

16

Posibles parametrizaciones de los dos tramos de ∂ D :

c 1 (x) = (x, x 2 ) , x ∈ [ 0 , 2 ] , c ′ 1 = (^1 ,^2 x)^ ,^ f^ (c^1 ) = (x

3 , x 4 )

c 2 (x) = (x, 2 x) , x ∈ [ 2 , 0 ] , c′ 2 = ( 1 , 2 ) , f (c 2 ) = (x^3 , 2 x^3 )

∂ D

f · ds =

∫ (^2)

0

(x 3

  • 2 x 5 ) dx −

∫ (^2)

0

(x 3

  • 4 x 3 ) dx = 26 3 −^2

4

16

Divergencia: div f = 4 x^2. n 1 =

( 2 x ,− 1 ) √ 1 + 4 x^2

, n 2 =

(− 2 , 1 ) √ 5

∫ (^2)

0

∫ (^2) x

x^2

4 x^2 dy dx =

∫ (^2)

0

( 8 x^3 − 4 x^2 ) dx = 32 5

∂ D

f · n ds =

∫ (^2)

0

(x 3 , x 4 ) ·

( 2 x ,− 1 ) √ 1 + 4 x^2

1 + 4 x^2 dx −

∫ (^2)

0

(x 3 , 2 x 3 ) ·

( 2 ,− 1 ) √ 5

5 dx =

∫ (^2)

0

x 4 dx = 32