LIVRO - DE - ..Calculo 3

Sumário
Aula 1: Integrais Duplas 11

1.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 12
1.2 Integral Dupla: Domínios Retangulares . . . . . . . . 12
1.3 Integral Dupla: Domínios Não Retangulares Limitados 14
1.4 Interpretação Geométrica . . . . . . . . . . . . . . . 15
1.5 Integrais Iteradas . . . . . . . . . . . . . . . . . . . 16
1.6 Propriedades das Integrais Duplas . . . . . . . . . . 19
1.7 Alguns Exemplos . . . . . . . . . . . . . . . . . . . 20
1.8 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 25
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 25
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 29
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 29
LEITURA COMPLEMENTAR . . . . . . . . . . . 30
Aula 2: Mudança de Variáveis em Integrais Duplas 33

2.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 34
2.2 Mudança de Variáveis em Integrais Duplas . . . . . 34
2.3 Alguns Exemplos . . . . . . . . . . . . . . . . . . . 39
2.4 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 43
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 43
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 45
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 45
Aula 3: Algumas Aplicações da Integral Dupla 47
3.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 48
3.2 Preliminares . . . . . . . . . . . . . . . . . . . . . . 48
3.3 Algumas Aplicações da Integral Dupla . . . . . . . . 52
3.4 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 58
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 58
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 59
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 59
Aula 4: Integrais triplas 63
4.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 64
4.2 Integração Tripla: Domínios Paralelepípedais . . . . 64
4.3 Integração Tripla: Domínios Não Paralelepípedais Li-

mitados . . . . . . . . . . . . . . . . . . . . . . . . 66
4.4 Interpretação Geométrica . . . . . . . . . . . . . . . 67
4.5 Integrais Iteradas . . . . . . . . . . . . . . . . . . . 67
4.6 Propriedades das Integrais Triplas . . . . . . . . . . 68
4.7 Exemplos . . . . . . . . . . . . . . . . . . . . . . . 69
4.8 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 74
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 75
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 80
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 80
Aula 5: Mudança de Variáveis em Integrais tríplas 83
5.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 84
5.2 Mudança de Variáveis em Integrais Triplas . . . . . . 84
5.3 Alguns Exemplos . . . . . . . . . . . . . . . . . . . 87
5.4 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 99
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 99
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 103
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 103
Aula 6: Algumas Aplicações das Integrais tríplas 105

6.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 106
6.2 Preliminares . . . . . . . . . . . . . . . . . . . . . . 106
6.3 Algumas Aplicações da Integral Tripla . . . . . . . . 110
6.4 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 117
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 118
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 119
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 120
Aula 7: Integrais de Funções Vetoriais sobre Curvas em R3

123
7.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 124
7.2 Curvas em R3 . . . . . . . . . . . . . . . . . . . . . 124
7.3 Massa, Momento de Massa e Momento de Inércia de
Curvas em R3 . . . . . . . . . . . . . . . . . . . . . 126
7.4 Campos Vetoriais: Trabalho, Circulação e Fluxo . . . 128
7.5 Independência do Caminho . . . . . . . . . . . . . . 130
7.6 Algumas Aplicações das Integrais de Linha . . . . . 133
7.7 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 138
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 139
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 142
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 142
Aula 8: Integrais de Superfícies 145

8.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 146
8.2 Superfícies em R3 . . . . . . . . . . . . . . . . . . . 146
8.3 Área de Superfícies em R3 . . . . . . . . . . . . . . 147
8.4 Momento de massa e Momento de Inércia de Super-
fícies de Casca Fina em R3 . . . . . . . . . . . . . . 151
8.5 Superfícies Parametrizadas . . . . . . . . . . . . . . 155
8.6 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 160
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 160
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 164
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 164
Aula 9: Teorema de Green e Teorema de Stokes 167

9.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 168
9.2 Preliminares . . . . . . . . . . . . . . . . . . . . . . 168
9.3 Teorema de Green . . . . . . . . . . . . . . . . . . 171
9.4 Estendendo o Teorema de Green para Outras Regiões 175
9.5 Verificação do Teorema de Green . . . . . . . . . . 178
9.6 Teorema de Stokes . . . . . . . . . . . . . . . . . . 181
9.7 Aplicação do Teorema de Stokes . . . . . . . . . . . 183
9.8 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 185
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 185
PRÓXIMA AULA . . . . . . . . . . . . . . . . . . . . 187
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 187
Aula 10: Teorema de Divergência 189

10.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . 190
10.2 Preliminares . . . . . . . . . . . . . . . . . . . . . . 190
10.3 Teorema da Divergência . . . . . . . . . . . . . . . 191
10.4 Estendendo o Teorema da Divergência . . . . . . . . 194
10.5 Algumas Aplicações do Teorema da Divergência . . 196
10.6 Conclusão . . . . . . . . . . . . . . . . . . . . . . . 200
RESUMO . . . . . . . . . . . . . . . . . . . . . . . . . 200
ATIVIDADES . . . . . . . . . . . . . . . . . . . . . . 202
AULA
Integrais Duplas 1
META:
Apresentar integrais duplas de funções de valores reais e domínio
em R2 .
OBJETIVOS:
Ao fim da aula os alunos deverão ser capazes de:
Definir a integral dupla de funções de valores reais e domínio em
R2 .
Calcular algumas integrais duplas de funções de valores reais e do-
mínio em R2 .
PRÉ-REQUISITOS
Os conhecimentos de integrais de funções de valores reais com do-
mínio em R, da disciplina Cálculo I.
Integrais Duplas
1.1 Introdução
Caros alunos iniciamos aqui nosso curso de Cálculo III com o

tema “Integrais Dupla”. A integração dupla, em essência, é uma
extensão natural da integral simples vista em Cálculo I e definida
como limite de somas de Riemann. Na prática, a integração dupla
é dada por duas integrações simples, cada uma efetuada sobre uma
variável e considerando as demais como constantes. É o que de-
nominamos de integrais interadas. Suas características e detalhes
próprios serão vistas ao longo do nosso curso, nas próximas duas
aulas.
1.2 Integral Dupla: Domínios Retangulares
Começamos por considerar uma função f definida em um do-

mínio retangular R = {(x, y) ∈ R2 |a ≤ x ≤ b∧c ≤ y ≤ d}. Formal-
mente f : [a, b]×[c, d] 7→ R. Usando a imaginação, pensemos em R
coberta por uma rede de retas paralelas aos eixos coordenados e que
dividem R em pequenos retângulos (Fig. 1.1) . Oficialmente, con-
sideraremos duas partições P [a, b] = {x0 = a, x1 , . . . , xj , xj+1 , . . . ,
xm = b} e P [c, d] = {y0 = c, y1 , . . . , yk , yk+1 , . . . , yn = d} onde
como visto em Cálculo I temos: x0 < x1 < · · · < xj < xj+1 <
· · · < xm e y0 < y1 < · · · < yk < yk+1 < · · · < yn .
Desta forma cada um dos Ij = [xj−1 , xj ] e Jk = [yk−1 , yk ] pe-
quenos subintervalos têm comprimentos ∆xj = xj − xj−1 e ∆yk =
yk − yk−1 , respectivamente. Definimos, agora, a uma partição para
o retângulo R por P = P [R] = P [a, b] × P [c, d], o produto carte-
siano das partições P [a, b] e P [c, d]. As retas retalham a região
R em uma série de retângulos Ajk = [xj−1 , xj ] × [yk−1 , yk ], 1 ≤
12
Cálculo III AULA
1
Figura 1.1: Partição de R = [a, b] × [c, d]
j ≤ m, 1 ≤ k ≤ n. A area de cada pequeno retângulo é dada

por ∆Ajk = ∆xj ∆yk . Como tanto ∆xj quanto ∆yk são dife-
rentes de zero, a área de cada pequeno retângulo é também di-
ferente de zero. Podemos então definir a norma da partição por:
|P | = max (∆Ajk ), que corresponde a maior área entre todos os
1≤j≤m
1≤k≤n
pequeno retângulo.
Pausa para respirar que já vamos definir a integral dupla sobre

domínios retangulares. Para isto tomamos um ponto (ξj , ζk ) ∈
BIOGRAFIA
[xj−1 , xj ] × [yk−1 , yk ] em cada pequeno retângulo e definimos a se-
Georg Friedrich
guinte soma de Riemann:
Bernhard Riemann
nasceu em Breselenz,
m X
X n Reino de Hanôver,
Smn = f (ξj , ζk )∆Ajk 17 de Setembro de
j=1 k=1
1826 e morreu em
Selasca, Itália, 20 de
A integral dupla da função f (x, y) sobre o retângulo R, denotada Junho de 1866, foi um
matemático alemão,
Z Z
f (x, y)dxdy será então definida como o seguinte limite: com contribuições
R fundamentais para a
Z Z análise e a geometria
def diferencial. Wikipedia
f (x, y)dxdy = lim Smn
R |P |→0
13
Integrais Duplas
Figura 1.2: Soma de Riemann para f (x, y) em R = [a, b] × [c, d]
1.3 Integral Dupla: Domínios Não Retangula-

res Limitados
Para definir a integral dupla de uma função f : D ⊂ R2 7→ R

onde D é não é uma região retangular, porém é limitada, co-
meçamos por considerar uma função F definida em um domí-
nio retangular R = {(x, y) 2
 ∈ R |a ≤ x ≤ b ∧ c ≤ y ≤ d} tal
 f (x, y) , (x, y) ∈ D
que D ⊂ R e F (x, y) = . Formalmente
 0 , (x, y) ∈
/D
F : [a, b] × [c, d] 7→ R é uma extensão da função f (x, y). Usando
a imaginação, pensemos em R coberta por uma rede de retas pa-
ralelas aos eixos coordenados e que dividem R em pequenos re-
tângulos e procedemos como na integral dupla sobre domínios re-
tangulares, considerando a uma partição para o retângulo R por
14
Cálculo III AULA
P = P [R] = P [a, b] × P [c, d], o produto cartesiano das partições 1
P [a, b] e P [c, d] onde P [a, b] = {x0 = a, x1 , . . . , xj , xj+1 , . . . , xm =
b} e P [c, d] = {y0 = c, y1 , . . . , yk , yk+1 , . . . , yn = d}. Do mesmo
modo definimos a norma da partição por: |P | = max (∆Ajk )
1≤j≤m
1≤k≤n
onde ∆Ajk = ∆xj ∆yk , ∆xj = xj − xj−1 e ∆yk = yk − yk−1 . To-
mamos um ponto (ξj , ζk ) ∈ [xj−1 , xj ] × [yk−1 , yk ] em cada pequeno
retângulo e definimos a seguinte soma de Riemann para a função
estendida F (x, y):
m X
X n
Smn = F (ξj , ζk )∆Ajk
j=1 k=1
A integral dupla da função f (x, y) sobre o domínio D ⊂ R2 , deno-

Z Z
tada f (x, y)dxdy será então definida como o seguinte limite:
D
Z Z
def
D |P |→0
Observem na partição (Fig. 1.3) que apenas os pequenos retân-

gulos cinza claro contribuem para a soma de Riemann os demais
têm contribuição nula visto que o ponto escolhido dentro destes
estão fora de D ⊂ R2 e portanto F (ξj , ζk ) = 0.
1.4 Interpretação Geométrica
Quando a função f (x, y) é positiva na região R, como a da

(Fig. 1.2), vemos que a soma de Riemann aproxima o volume do
prisma sólido reto limitado inferiormente por R e superiormente
pela superfície z = f (x, y) e quanto maior for o refinamento da par-
tição de R melhor será a aproximação. Podemos então, interpretar
Z Z
a integral dupla f (x, y)dxdy como o volume do prisma sólido
R
15
Integrais Duplas
Figura 1.3: Partição para F (x, y) em R = [a, b] × [c, d]
reto limitado inferiormente por R e superiormente pela superfície

z = f (x, y).
1.5 Integrais Iteradas
Do mesmo modo que para a integral simples, na integral dupla

a soma de Riemann não é um modo prático de se calcular uma
integral dupla. Vejamos agora um procedimento que facilitará o
cálculo de integrais duplas. Vamos exemplificar calculando o vo-
lume de um prisma reto de base retangular, limitado inferiormente
por [a, b] × [c, d] e superiormente pela função de valores positivos
f (x, y). para cada valor fixo de x no intervalo [a, b] consideremos
o perfil A(x) (área da seção transversal em x) (Fig. 1.4) fazemos
o produto por dx e integramos no intervalo [a, b]. Isto resulta no
16
Cálculo III AULA
1
Figura 1.4: A(x), x fixo, integramos em relação a y
volume do citado prisma.
Z b
V = A(x)dx
a
Por outro lado o perfil A(x) é dada pela área abaixo da curva
f (x, y), fixado o x, entre os valores de y no intervalo [c, d]. E como
Z d
vimos em Cálculo I A(x) = f (x, y)dy.
c
O volume do prisma pode ser então escrito como:
Z b Z d
V = f (x, y)dy dx
a c
.
Podemos alternativamente calcular o mesmo volume considerando
os perfis A(y) (área da seção transversal em y) (Fig. 1.5) fazemos
o produto por dy e integramos no intervalo [c, d]. Isto resulta no
volume do citado prisma.
17
Integrais Duplas
Figura 1.5: A(y), y fixo, integramos em relação a x
Z d
V = A(y)dy
c
Z b
Da mesma forma como vimos em Cálculo I A(y) = f (x, y)dx.
a
O volume do prisma pode ser então escrito como:
Z d Z b
V = f (x, y)dx dy
c a
.
Como o volume dado pelas duas expressões é o mesmo temos que:
Z d Z b Z b Z d
f (x, y)dx dy = f (x, y)dy dx
c a a c
ou seja a ordem em que as integrais simples são executadas não

altera o resultado final da integração dupla em domínios retangu-
lares. Este procedimento e conhecido como integrais iteradas.
18
Cálculo III AULA
1.6 Propriedades das Integrais Duplas 1
Como nosso curso é de Cálculo, apenas listaremos, sem de-
monstração, alguma das propriedades das integrais duplas. Caso
desejem conhecer a demonstração de algumas destas propriedades,
recomendo livros de Cálculo Avançado como os citados na biblio-
grafia abaixo.
Propriedade 1.1. Sejam f : D ⊂ R2 7→ R uma função de valores

reais integrável em D e c ∈ R, então vale:
Z Z Z Z
cf (x, y)dxdy = c f (x, y)dxdy
D D
Propriedade 1.2. Sejam f, g : D ⊂ R2 7→ R duas funções de

valores reais integráveis em D, então vale:
Z Z Z Z Z Z
(f + g)(x, y)dxdy = f (x, y)dxdy + g(x, y)dxdy
D D D

reais integrável em D tal que f (x, y) ≥ 0, ∀(x, y) ∈ D, então vale:
Z Z
f (x, y)dxdy ≥ 0
D
Propriedade 1.4. Sejam f, g : D ⊂ R2 7→ R duas funções de va-

lores reais integráveis em D tais que f (x, y) ≥ g(x, y), ∀(x, y) ∈ D,
então vale:
Z Z Z Z
f (x, y)dxdy ≥ g(x, y)dxdy
D D
Propriedade 1.5. Seja f : D ⊂ R2 7→ R uma função de valores

reais integrável em D onde D = A ∪ B e A ∩ B é a união de um
19
Integrais Duplas
número finito de curvas em R2 , então vale:
Z Z Z Z Z Z
f (x, y)dxdy = f (x, y)dxdy + f (x, y)dxdy
D A B
OBS 1.1. As duas primeiras propriedades diz respeito à “line-

aridade” do operador integral dupla. As terceira e quarta pro-
priedades são denominadas “dominação” enquanto que a quinta
propriedade é denominada “aditividade”.
1.7 Alguns Exemplos
Nada mais natural que ilustrar um novo conceito com exem-

plos e, vamos aqui fazer exatamente isto, ilustrar o conceito de
integral dupla com dois exemplos. Antes porém, vale observar que
a na prática uma integral dupla equivale a duas integrais simples e
neste caso uma pergunta fica no ar e não deixaremos sem resposta.
Qual das duas variáveis x ou y integraremos primeiro? Muito bem,
a resposta é dada pela propria expressão da integral dupla. Isto
Z Z
é, na integral f (x, y)dxdy primeiramente integramos na va-
R Z Z
riável x e depois na variável y. Já na integral f (x, y)dydx
R
primeiramente integramos na variável y e depois na variável x.
Vamos diretamente para o primeiro exemplo de integral dupla so-

bre domínios retangulares. A saber:
Exemplo 1.1. Considere a função f : [0, 1] × [0, 1] 7→ R (Fig.

1.6) dada por f (x, y) = exp(−x − y) e determine a integral dupla
Z Z
I = f (x, y)dxdy sobre a região R = {(x, y) ∈ R2 |0 ≤ x ≤
R
1 ∧ 0 ≤ y ≤ 1}.
20
Cálculo III AULA
1
Figura 1.6: Função f : [0, 1] × [0, 1] 7→ R: f (x, y) = exp(−x − y)
SOLUÇÃO:
Passo 1 colocaremos os limites de integração que representam a
região R dada, segundo a ordem de integração:
Z 1Z 1
I= exp(−x − y)dxdy
0 0
Lembrando que: exp(−x − y) = exp(−x) exp(−y) temos:
Z 1Z 1
I= exp(−x) exp(−y)dxdy
0 0
Passo 2 integraremos na variável x considerando a variável y
como uma constante:
Z 1 1
I= − exp(−x) exp(−y)dy

0 0
Substituindo os limites de integração temos:
Z 1
I= (− exp(−1) − (− exp(−0))) exp(−y)dy
0
Efetuando os cálculos temos:
Z 1
I= (1 − exp(−1)) exp(−y)dy
0
Passo 3 integraremos
na variável y considerando a variável:
1
I = (1 − exp(−1)) − exp(−y)

0
I = (1 − exp(−1)) (− exp(−1) − (− exp(−0)))
21
Integrais Duplas
I = (1 − exp(−1))2
OBS 1.2. Daremos aqui um método prático para determinar os

limites de integração em uma integral dupla sobre domínio não re-
tangular da forma: D.
Passo 1 Fazer um desenho da região D. (Fig. 1.7) identificando

as curvas inferior a(x) e superior b(x) que limitam a região D.
Passo 2 Atravessar toda a região D e o eixo x com um segmento
de reta paralelo e orientado na direção positiva ao eixo y (segmento
AB na Fig. 1.7)
Passo 3 Deslocar o segmento de reta AB paralelo ao eixo y na
Figura 1.7: Determinação prática dos limites para D
direção negativa do eixo x até tocar o ponto mais à esquerda de D

marcando o limite inferior de x (ponto a na Fig. 1.7).
direção positiva do eixo x até tocar o ponto mais à direita de D
marcando o limite superior de x (ponto b na Fig. 1.7).
Passo 5 Tomando um ponto qualquer x ∈ (a, b) passamos o seg-
22
Cálculo III AULA
mento de reta AB através da região D. O limite inferior para a 1
variável y será a função a(x), ponto da curva onde o segmento
entra na região D e o limite superior para a variável y será b(x),
ponto da curva onde o segmento de reta sai da região D.
Nossa integral será efetuada assim:
Z Z Z bZ b(x)
f (x, y)dxdy = f (x, y)dydx
D a a(x)
Exemplo 1.2. Considere a função f : [0, 1] × [0, 1] 7→ R (Fig.

1.8) dada por f (x, y) = y(3x − x2 − y) e determine a integral du-
Z Z
pla I = f (x, y)dxdy sobre a região D ∈ R2 interseção das
R
curvas y = 0 e y = 3x − x2 .
Figura 1.8: Função f : [0, 1] × [0, 1] 7→ R: f (x, y) = x.y
SOLUÇÃO:
Passo 1 faremos o desenho das duas curvas que determinam os
limites para a região D. A saber y = 0 e y = 3x − x2 (Fig. 1.9).
Passo 2 usando o processo prático exposto acima determinamos
os limites de integração. A saber: a = 0, b = 3, a(x) = 0 e
b(x) = 3x − x2 .
23
Integrais Duplas
Figura 1.9: Limites para o domínio D
A integral passa a ser escrita como:

Z Z Z 3 Z 3x−x2
I= f (x, y)dxdy = y(3x − x2 − y)dydx
0 0
Operando no integrando fazendo o produto por y temos:
Z 3 Z 3x−x2
I= (y(3x − x2 ) − y 2 )dydx
0 0
Passo 3 efetuando a integração em y temos:
Z 3 2
y y 3 3x−x2
I= ( (3x − x2 ) − ) dx
0 2 3 0
Substituindo os limit3es de integração temos:
Z 3
(3x − x2 )2 (3x − x2 )3
I= ( (3x − x2 ) − )dx
0 2 3
Efetuando as simplificações teremos:
Z 3
(3x − x2 )3
I= dx
0 6
Expandindo o binômio de Newton temos:
1 3
Z
I= (27x3 − 27x4 + 9x5 − x6 )dx
6 0
Passo 4 efetuando a integração em x temos:
1 x4 x5 x6 x7 3
I = (27 − 27 + 9 − )
6 4 5 6 7 0
Substituindo os limit3es de integração temos:
1 34 35 36 37
I = (27 − 27 + 9 − )
6 4 5 6 7
Efetuando os cálculos, garantido muito trabalho, temos:
24
Cálculo III AULA
I=
729
280

1
1.8 Conclusão
Na aula de hoje, vimos que a integral dupla é uma extensão

natural do conceito de integral simples visto em Cálculo I. E se
por um lado a integral simples pode ser interpretada como a área
sob a curva descrita pela função a ser integrada, a integral dupla
pode ser vista como o volume sob a superfície descrita pela função
a ser duplamente integrada.
RESUMO
No nosso resumo da Aula 01 constam os seguintes tópicos:
Integração Dupla: Domínios retangulares

Considerando uma função f : R 7→ R onde R = {(x, y) ∈ R2 |a ≤
x ≤ b ∧ c ≤ y ≤ d} é um retângulo em R2 . Podemos cobri-lo com
uma malha de retas formada pela partição: P = P [R] = P [a, b] ×
P [c, d] onde cada P [a, b] = {x0 = a, x1 , . . . , xj , xj+1 , . . . , xn = b}
e P [c, d] = {y0 = c, y1 , . . . , yk , yk+1 , . . . , ym = d} são partições
dos intervalos [a, b] em x e [c, d] em y respectivamente. A ma-
lha divide R nos retângulos Ajk = [xj−1 , xj ] × [yk−1 , yk ], 1 ≤ j ≤
n, 1 ≤ k ≤ m de área ∆Ajk = ∆xj ∆yk onde ∆xj = xj − xj−1
e ∆yk = yk − yk−1 são os comprimentos dos subintervalos Ij =
[xj−1 , xj ] e Jk = [yk−1 , yk ] respectivamente. Defini-se a norma da
partição por: |P | = max (∆Ajk ). Toma-se um ponto (ξj , ζk ) ∈
1≤j≤n
1≤k≤m
25
Integrais Duplas
[xj−1 , xj ] × [yk−1 , yk ] em cada retângulo Ajk e definimos a seguinte

soma de Riemann:
n X
X m
Snm = f (ξj , ζk )∆Ajk
j=1 k=1
A integral dupla da função f (x, ) sobre o retângulo R, denotada

Z Z
f (x, y)dxdy será então definida como o seguinte limite:
R
Z Z
def
f (x, y)dxdy = lim Snm
R |P |→0
Integração Dupla: Domínios não Retangulares

Para definir a integral dupla de uma função f : D ⊂ R2 7→ R
onde D é não é uma região retangular, porém é limitada, co-
meçamos por considerar uma função F definida em um domí-
nio retangular R = {(x, y) 2
 ∈ R |a ≤ x ≤ b ∧ c ≤ y ≤ d} tal
 f (x, y) , (x, y) ∈ D
que D ⊂ R e F (x, y) = . Formalmente
 0 , (x, y) ∈
/D
F : [a, b] × [c, d] 7→ R é uma extensão da função f (x, y). A partir
daqui todo o procedimento é semelhante ao da definição da inte-
gral dupla em domínios retangulares. Podemos definir a integral
dupla de uma função f (x, y) em um domínio não retangular D por:
Z Z
def
D |P |→0
Pm Pn
onde: Smn = j=1 k=1 F (ξj , ζk )∆Ajk . é a soma de Riemann
para F (x, y)
Integrais Iteradas
As integrais iteradas dizem que em um domínio retangular R =
26
Cálculo III AULA
[a, b] × [c, d] a ordem de execução das integrais simples não alteram 1
o valor da integral dupla, que pode ser representada por:
Z d Z b Z b Z d
f (x, y)dx dy = f (x, y)dy dx
c a a c
.
Propriedades das Integrais Duplas
As integrais duplas são de certo modo semelhantes às propriedades
das integrais simples que vimos em Cálculo I sendo quase que uma
extensão natural destas. As integrais duplas têm, entre outras, as
seguintes propriedades:
Propriedade 1 Sejam f : D ⊂ R2 7→ R uma função de valores
Z Z Z Z
cf (x, y)dxdy = c f (x, y)dxdy
D D
Propriedade 2 Sejam f, g : D ⊂ R2 7→ R duas funções de valores

reais integráveis em D, então vale:
Z Z Z Z Z Z
(f + g)(x, y)dxdy = f (x, y)dxdy + g(x, y)dxdy
D D D

reais integrável em D tal que f (x, y) ≥ 0, ∀(x, y) ∈ D, então vale:
Z Z
f (x, y)dxdy ≥ 0
D
27
Integrais Duplas
Propriedade 4 Sejam f, g : D ⊂ R2 7→ R duas funções de valo-

res reais integráveis em D tais que f (x, y) ≥ g(x, y), ∀(x, y) ∈ D,
então vale:
Z Z Z Z
f (x, y)dxdy ≥ g(x, y)dxdy
D D
Propriedade 5 Seja f : D ⊂ R2 7→ R uma função de valores

número finito de curvas em R2 , então vale:
Z Z Z Z Z Z
f (x, y)dxdy = f (x, y)dxdy + f (x, y)dxdy
D A B
Determinação dos Limites de Integração

Para determinar os limites de integração em uma integral dupla
sobre domínio não retangular da forma: D seguimos os seguintes
passos:
as curvas inferior a(x) e superior b(x) que limitam a região D.
Passo 2 Atravessar toda a região D e o eixo x com um segmento
AB na Fig. 1.7)
direção negativa do eixo x até tocar o ponto mais à esquerda de D
28
Cálculo III AULA
Passo 4 Deslocar o segmento de reta AB paralelo ao eixo y na 1
direção positiva do eixo x até tocar o ponto mais à direita de D
mento de reta AB através da região D. O limite inferior para a
variável y será a função a(x), ponto da curva onde o segmento
entra na região D e o limite superior para a variável y será b(x),
ponto da curva onde o segmento de reta sai da região D.
Z Z Z bZ b(x)
f (x, y)dxdy = f (x, y)dydx
D a a(x)
PRÓXIMA AULA
Em nossa próxima aula veremos mudança de variáveis na in-

tegração dupla. O objetivo da mudança de variáveis em uma in-
tegral dupla será a de facilitar esta integração de uma de duas
formas. A primeira será tornando o integrando mais simples. A
segunda transformando o domínio D do integrando em um domínio
de forma geométrica mais simples.
ATIVIDADES
Deixamos como atividades o cálculo de algumas integrais du-

plas.
ATIV. 1.1. Seja f : [−1, +1] × [−1, +1] 7→ R dada por f (x, y) =
29
Integrais Duplas
Z Z
x2 + y2. Determine a integral dupla f (x, y)dxdy.
R
Comentário: Volte ao texto e reveja com calma e atenção o

cálculo de integrais duplas dos exemplos acima, elas lhe servirão
de guia.
ATIV. 1.2. Seja f : D ⊂ R2 7→ R dada por f (x, y) = x2 + y 2 ,

onde D = {(x, y) ∈ R2 |x ≥ 0 ∧ 0 ≤ y ≤ 1 − x2 }.
Z Z
• Determine os limites da integral dupla f (x, y)dxdy,
D
• esboce a região de integração e

Z Z
• calcule a integral dupla f (x, y)dxdy.
D

de guia.
LEITURA COMPLEMENTAR
ÁVILA, Geraldo, Cálculo 3: Funções de Várias Variáveis, Livros

Técnicos e Científicos Editora, São Paulo, 3a edição, 1982.
LEITHOLD, Louis, O Cálculo com Geometria Analítica. Volume
2, Editora Harbra, 1994.
STEWART, James,Cálculo. Volume 3, 5a edição, Editora CEN-
GAGE Learning, 2009.
SWOKOWSKI, Earl E., Cálculo com Geometria Analítica, Volume
2, 2a edição, Makron Books do Brásil SP, 1994.
THOMAS, George B., Cálculo, Volume 2, 10a, Addilson Wesley,
30
Cálculo III AULA
2003. 1
KAPLAN, Wilfred, Cálculo Avançado Vol.1 e vol.2 Editora Edgard
Blücher 1991.// SPIEGEL, Murray R. Cálculo Avançado, Editora
McGraw-Hill do Brasil, 1971.
BOUCHARA, Jacques, Cálculo Integral Avançado, EDUSP, 2006.
31
AULA
Mudança de Variáveis em
Integrais Duplas 2
META:
Introduzir mudança de variáveis em integrais duplas de funções de
valores reais e domínio em R2 .
OBJETIVOS:
Calcular o jacobiano de aplicações de R2 em R2 .
Calcular integrais duplas de funções de valores reais e domínio em
R2 utilizando mudança de variáveis.
Calcular integrais duplas de funções de valores reais e domínio em
R2 em coordenadas polares.
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I, curvas em R2 e coordenadas
polares da disciplina Cálculo II e integrais duplas aula 01.
Mudança de Variáveis em Integrais Duplas
2.1 Introdução
Caros alunos a segunda aula do nosso curso de Cálculo III

tem com o tema “Mudança de Variáveis em Integrais Duplas”. As
Z Z
vezes, na integral dupla f (x, y)dxdy, dada a natureza ou de
D
HISTÓRIA f (x, y) ao do seu domínio D, fica mais fácil integrar se fizermos
O teorema de mu- uma mudança nas variáveis de integração, como quando D é uma
dança de variáveis
em integrais duplas disco, um semi-disco, um setor circular ou mesmo uma faixa de
foi primeiro proposto
por Euler quando ele disco, usando-se o sistema de coordenadas polares de modo geral
desenvolveu a noção
de integral dupla em a integral dupla é mais fácil de se determinar que em coordenadas
1769. Usado por
Legendre, Laplace e cartesianas.
Gauss, foi primeira-
mente generalizado
para n variáveis por
Mikhail Ostrogradski 2.2 Mudança de Variáveis em Integrais Duplas
em 1836, resistiu a
uma demonstração
Caros alunos começaremos revendo mudança de variáveis em
mais rigorosa por longo
tempo (cerca de 125 integrais simples. Considere uma função f : [a, b] 7→ R. A idéia
anos). E foi satisfató-
riamente demonstrado é mudar a variável inicial x para uma nova variável ξ relaciona-
por Elie Cartan em
uma série de artigos das por x = g(ξ), onde g(ξ) é uma função biunívoca estritamente
nos anos 1890.
crescente ou estritamente decrescente em [a, b]. Isto garante que
podemos inverter a mudança de variáveis.
Seja F (x) uma anti-derivada de f (x) tal que F 0 (x) = f (x). Então,
da regra da cadeia temos:
d
F (g(ξ)) = F 0 (g(ξ))g 0 (ξ) = f (g(ξ))g 0 (ξ).
dξ
Integrando com respeito a ξ temos:
Z Z
d
F (g(ξ))dξ = f (g(ξ))g 0 (ξ)dξ
dξ
Das propriedades da integral temos:
Z
F (g(ξ)) + C = f (g(ξ))g 0 (ξ)dξ
Como x = g(ξ) temos:
Z
F (x) + C = f (g(ξ))g 0 (ξ)dξ
34
Cálculo III AULA
Como F (x) é uma primitiva de f (x) a primeira expressão é a in- 2
tegral indefinida de f (x) com respeito a x e temos:
Z Z
f (x)dx = f (g(ξ))g 0 (ξ)dξ
Que representa a mudança de variáveis em uma integral simples.
Para integrais definidas, se c = g(a) e d = g(b) então:
Z b Z d
f (x)dx = f (g(ξ))g 0 (ξ)dξ
a c
A expressão acima funciona bem quando g(ξ) é crescente neste
caso a < b e c < d. Porém, no caso de g(ξ) decrescente (g 0 (ξ) < 0)
pois neste caso a < b e d < c e portanto o limite inferior da se-
gunda integral não conhecide com o limite inferior do intervalo da
imagem de g(ξ) o mesmo acontecendo com o limite superior. Neste
caso, usando as propriedades da integral simples temos:
Z b Z c
f (x)dx = − f (g(ξ))g 0 (ξ)dξ
a d
De outra forma escrevemos:
Z b Z c
f (x)dx = f (g(ξ))|g 0 (ξ)|dξ.
a d
e operaremos os limites inferiores e superiors das integrais como os
limites inferiores e superiores dos domínios (intervalos) e a expres-
são acima vale tanto pra g(ξ) crescente quanto decrescente. Vamos
agora diretamente ao assunto dando uma argumentação heurística OBSERVAÇÃO
para a expressão da mudança de variáveis em integrais duplas. heurística heu.rís.ti.ca

Z Z sf (gr heuristiké) 1
Para isto, consideremos a integral dupla f (x, y)dxdy sobre Ciência ou arte do pro-
D cedimento heurístico.
uma região D ∈ R2 do plano (x, y) e a transformação (x, y) = 2 Método de ensino
que consiste em que
T (u, v) tal que o domínio D do plano (x, y) seja a imagem do o educando chegue à
verdade por seus pró-
domínio D0 do plano (u, v) (podemos expressar este fato como prios meios. 3 Ramo
da ciência histórica que
D = T (D0 )). Mais especificamente podemos escrever: x = x̂(u, v) consiste na pesquisa
dos documentos do
e y = ŷ(u, v) tomando uma partição para o domínio D0 no plano
passado.
(u, v) cobrindo-o com pequenos retângulos e usando a transforma-
ção T podemos levar o pequeno retângulo A0jk na pequena figura
35
plana Ajk = T (A0jk ) (ver Fig 2.1 e Fig 2.2). A área do pequeno
retângulo no plano (u, v) é ∆A0jk a área da pequena figura Ajk no
plano (x, j), e ai é que reside a argumentação heurística, será apro-
∂T
ximada pela área do paralelogramo formado pelos vetores ∆vk
∂v
∂T
e ∆uj e pelas linhas tracejadas (paralelas aos respectivos veto-
∂u
res). Do calculo vetorial temos:
∂T ∂ x̂ ∂ ŷ
∆uj = ∆uj~i + ∆uj~j + 0~k.
∂u ∂u ∂u
∂T ∂ x̂ ∂ ŷ
∆vk = ∆vk~i + ∆vk~j + 0~k
∂v ∂v ∂v
Vistos como vetores de R3 e a área do paralelogramo (ver Vetores
e Geometria Analítica) dada pelo módulo do seguinte produto ve-
torial:
∂T ∂T
∆Ajk = ∆uj × ∆vk .

∂u ∂v
Fazendo o cálculo do produto  vetorial temos: 
~i ~j ~k
 
∂T ∂T  ∂ x̂ ∂ ŷ 
∆uj × ∆vk = det  ∆uj ∆uj 0 
∂u ∂v  ∂u ∂u
 
∂ x̂ ∂ ŷ 
∆vk ∆vk 0
∂v ∂v
Fazendo os cálculos temos:
∂T ∂T ∂ x̂ ∂ ŷ ∂ x̂ ∂ ŷ
∆uj × ∆vk = − ∆uj ∆vk~k.
∂u ∂v ∂u ∂v ∂v ∂u
Tomando o módulo da expressão acima, para a área de Ajk , temos:
∂ x̂ ∂ ŷ ∂ x̂ ∂ ŷ
∆Ajk ≈ − ∆uj ∆vk .

∂u ∂v ∂v ∂u
A expressão dentro do módulo é o determinante de uma matrix
2 × 2 conhecida como jacobiano da transformação x = x̂(u, v) e
y = ŷ(u, v) e denotado:
 
∂ x̂ ∂ ŷ
∂(x, y) ∂ x̂ ∂ ŷ ∂ x̂ ∂ ŷ
= det  ∂u ∂u  = − .

∂(u, v) ∂ x̂ ∂ ŷ ∂u ∂v ∂v ∂u
∂v ∂v
Como a área do pequeno retângulo A0jk é dada por ∆A0jk = ∆uj ∆vk
temos:
∂(x, y)
∆Ajk ≈ ∆A0jk .

∂(u, v)
36
Cálculo III AULA
2
Figura 2.1: Plano (u, v) Figura 2.2: Plano (x, y)
O que nos leva a considerar a seguinte fórmula para a mudança de

variáveis em integrais duplas:
Z Z Z Z ∂(x, y)
f (x, y)dxdy = f (x̂(u, v), ŷ(u, v)) dudv.

D D0 ∂(u, v)
Que representa a mudança de variáveis na integral dupla pela
transformação (x, y) = T (u, v).
OBS 2.1. Para o caso particular da mudança de variáveis do sis-

tema de coordenadas cartesianas para o sistema de coordenadas
polares (x, y) = T (r, ϑ) = (r cos(ϑ), r sin(ϑ)) onde x = x̂(r, ϑ) =
r cos(ϑ) e y =ŷ(r, ϑ) = r sin(ϑ),
 o jacobiano é dado por:
∂ x̂ ∂ ŷ  
∂(x, y)  ∂r ∂r  cos(ϑ) sin(ϑ)
= det  ∂ x̂ ∂ ŷ  = det   = r.
∂(r, ϑ) −r sin(ϑ) r cos(ϑ)
∂ϑ ∂ϑ
∂(x, y)
Portanto o jacobiano da transformação = r a mudança de
∂(r, ϑ)
variáveis na integral dupla toma a forma:
Z Z Z Z
f (x, y)dxdy = f (r cos(ϑ), r sin(ϑ))rdrdϑ.
D D0

limites de integração em uma integral dupla sobre domínio não re-
tangular da forma: D em coordenadas polares.
37
Passo 1 Fazer um desenho da região D (Fig. 2.3), identificando

as curvas que limitam a região D.
Passo 2 Atravessar toda a região D com um raio ~r (ϑ) orientado
na direção positiva (Fig. 2.3)
Passo 3 Deslocar o raio ~r (ϑ) na direção negativa do ângulo ϑ (di-
NOTA
reção horária) até tocar o ponto mais à negativa de D marcando
Por convenção a
medida de ângulo tem
sinal positivo quando o
deslocamento é feito na
direção anti-horária,
direção contrária ao
movimento dos pon-
teiros do relógio e tem
sinal negativo quando
o deslocamento é feito
na direção horária,
direção do movimento
dos ponteiros do
relógio.
o limite inferior de ϑ (ângulo α na Fig. 2.3).

Passo 4 Deslocar o raio ~r (ϑ) na direção positiva do ângulo ϑ
(direção anti-horária) até tocar o ponto mais à positiva de D mar-
cando o limite inferior de ϑ (ângulo β na Fig. 2.3).
Passo 5 Tomando um ponto qualquer ϑ ∈ (α, β) passamos o raio
~r (ϑ) através de D o limite inferior para a variável r será a função
α(ϑ), ponto da curva onde o raio ~r (ϑ) entra na região D e o limite
superior para a variável r será β(ϑ), ponto da curva onde o raio
~r (ϑ) sai da região D.
Z Z Z β Z β(ϑ)
f (x, y)dxdy = f (r cos(ϑ), r sin(ϑ))rdrdϑ
D α α(ϑ)
38
Cálculo III AULA
2.3 Alguns Exemplos 2
Caros alunos, nesta seção ilustraremos, com dois exemplos, a
mudança de variáveis em integrais duplas. A rigor, trataremos
apenas de exemplos em coordenadas polares.
Z Z
Exemplo 2.1. Determinar a integral dupla f (x, y)dxdy onde
D
D = {(x, y) ∈ R2 |x ≥ 0 ∧ y ≥ 0 ∧ x2 + y 2 ≤ 1} e f (x, y) =
exp(−x2 − y 2 ). O domínio da função representa um quarto de
disco (Fig 2.4).
Figura 2.4: Gráfico do exemplo 1
SOLUÇÃO:
Passo 1 Como o domínio D é um quarto de disco, o mais ade-
quado é utilizar o sistema de coordenadas polares. Podemos usar
o método prático de determinação dos limites da integral dupla
π
em coordenadas polares (Fig 2.5) e verificar que: α = 0, β = ,
2
α(ϑ) = 0 e β(ϑ) = 1.
Neste caso podemos descrever o domínio como: D0 = {(r, ϑ) ∈
R2 |0 ≤ r ≤ 1 ∧ 0 ≤ ϑ ≤ π/2}. E como x = r cos(ϑ) e y = r sin(ϑ) e
39

∂(x, y)
o módulo do jacobiano da transformação é dado por: = r.
∂(r, ϑ)
Quanto a variável r varia no intervalo [0, 1] independentemente de
ϑ e a variável ϑ varia no intervalo [0, π/2] ( a variação de ângulo no
primeiro quadrante). Podemos reescrever a integral dupla como:
Z Z Z 1 Z π/2
I= f (x, y)dxdy = f (r cos(ϑ), r sin(ϑ))rdϑdr Subs-
D 0 0
tituindo f (x, y) temos:
Z 1 Z π/2
I= exp(−(r cos(ϑ))2 − (r sin(ϑ))2 )rdϑdr
0 0
Efetuando as simplificações temos:
Z 1 Z π/2
I= exp(−r2 )rdϑdr
0 0
Passo 2 Integrando primeiramente na variável ϑ e como o inte-
grando não depende de ϑ temos:
Z 1 π/2
I= exp(−r2 )ϑ rdr

0 0
Z 1
I = π/2 exp(−r2 )rdr
0
Passo 3 A última integral (variável r) podemos efetuar por mu-
dança de variáveis pondo ξ = r2   temos: dξ = 2rdr
deste modo
1  1  1
ou seja rdr = − dξ e os limites r eξ . Daí, a integral
2  0  0
40
Cálculo III AULA
passará a forma:
Z 1
2
I = π/4 exp(−ξ)dξ
0
Cuja integração é fácil e da forma:
1
I = π/4 − exp(−ξ)

0
Efetuando os cálculo temos:
π
I = (1 − exp(−1))
4
Vamos agora, diretamente ao nosso segundo exemplo. Trata-se de

uma curva já conhecida de vocês (Cálculo II) a lemniscata.
Exemplo 2.2. Determinar a área da região D, a parte da lemnis-

p
cata, r = cos(2ϑ), que situa-se no primeiro quadrante. ver parte
cinza da (Fig 2.6).
SOLUÇÃO:
Passo 1 Como o domínio D é um quarto de uma lemniscata, o
mais adequado é utilizar o sistema de coordenadas polares. Pode-
mos usar o método prático de determinação dos limites da integral
dupla em coordenadas polares (Fig 2.7) e verificar que: α = 0,
41
π p
β= , α(ϑ) = 0 e β(ϑ) = cos(2ϑ).
4
Neste caso podemos descrever o domínio como: D0 = {(r, ϑ) ∈

p
R2 |0 ≤ ϑ ≤ π/4 ∧ 0 ≤ r ≤ cos(2ϑ)}. E como, neste exemplo,
queremos calcular área temos que f (x, y) = 1 e em coordenadas
polares podemos escrever na forma da seguinte integral dupla:
Z Z Z π/4 Z √cos(2ϑ)
A= dxdy = rdrdϑ
D 0 0
Integrando em√r temos:
Z π/4 2
r cos(2ϑ)
A= dϑ
2 0

0
Substituindopos limites de integração temos:
Z π/4 2
cos(2ϑ)
A= dϑ
0 2
Simplificando o integrando temos:
Z π/4
cos(2ϑ)
A= dϑ
0 2
Integrando na variável ϑ temos:
sin(2ϑ) π/4
A=
4

0
sin(π/2) − sin(0)
A=
4
Portanto:
1
A=
4
42
Cálculo III AULA
2
OBS 2.3. Caros alunos, é muito importante neste ponto uma revi-
são cuidadosa e detalhada dos dois exemplos dados acima. Efetuar
uma mudança de varáveis em integrais duplas não é tão simples
quanto efetuar uma mudança de variáveis em integrais simples.
2.4 Conclusão
Na aula de hoje, vimos que a mudança de variáveis em integrais

dupla, nos permite, facilitar o cálculo das ditas integrais quando
trabalhamos com domínios de integração de geometrias específicas,
como a induzida pelas coordenadas polares.
RESUMO
No nosso resumo da Aula 02 constam os seguintes tópicos:

Consideramos a transformação (x, y) = T (u, v) tal que o domínio
D do plano (x, y) seja transformado no domínio D0 do plano (u, v)
(D = T (D0 )) e mais especificamente x = x̂(u, v) e y = ŷ(u, v).
∂(x, y)
Definindo o jacobiano da transformação, denotado , por:
∂(u, v)
 
∂ x̂ ∂ ŷ
∂(x, y) ∂ x̂ ∂ ŷ ∂ x̂ ∂ ŷ
= det  ∂u ∂u  = − .

∂(u, v) ∂ x̂ ∂ ŷ ∂u ∂v ∂v ∂u
∂v ∂v
Vale então,a seguinte fórmula para a mudança de variáveis en in-
tegrais duplas:
43
Z Z Z Z ∂(x, y)
f (x, y)dxdy = f (x̂(u, v), ŷ(u, v)) dudv.

D D 0 ∂(u, v)
Sistema de Coordenadas Polares

Para o caso particular da mudança de variáveis do sistema de co-
ordenadas cartesianas para o sistema de coordenadas polares no
cálculo de integrais duplas temos:
(x, y) = T (r, ϑ) = (r cos(ϑ), r sin(ϑ)) onde x = x̂(r, ϑ) = r cos(ϑ) e
y = ŷ(r, ϑ) = r sin(ϑ).
Vale a seguinte transformação de variáveis:
Z Z Z Z
f (x, y)dxdy = f (r cos(ϑ), r sin(ϑ))rdrdϑ.
D D0
Determinação dos Limites de Integração em Coordena-

das Polares
Daremos aqui um método prático para determinar os limites de
integração em uma integral dupla sobre domínio não retangular
da forma: D em coordenadas polares.
Passo 1 Fazer um desenho da região D (Fig. 2.3), identificando
as curvas que limitam a região D.
Passo 2 Atravessar toda a região D com um raio ~r (ϑ) orientado
na direção positiva (Fig. 2.3)
Passo 3 Deslocar o raio ~r (ϑ) na direção negativa do ângulo ϑ (di-
reção horária) até tocar o ponto mais à negativa de D marcando
o limite inferior de ϑ (ângulo α na Fig. 2.3).
Passo 4 Deslocar o raio ~r (ϑ) na direção positiva do ângulo ϑ
(direção anti-horária) até tocar o ponto mais à positiva de D mar-
cando o limite inferior de ϑ (ângulo β na Fig. 2.3).
Passo 5 Tomando um ponto qualquer ϑ ∈ (α, β) passamos o raio
44
Cálculo III AULA
~r (ϑ) através de D o limite inferior para a variável r será a função 2
α(ϑ), ponto da curva onde o raio ~r (ϑ) entra na região D e o limite
superior para a variável r será β(ϑ), ponto da curva onde o raio
~r (ϑ) sai da região D.
Z Z Z β Z β(ϑ)
f (x, y)dxdy = f (r cos(ϑ), r sin(ϑ))rdrdϑ
D α α(ϑ)
PRÓXIMA AULA
Em nossa próxima aula veremos algumas das inúmeras aplica-

ções da integral dupla. Nossa atenção estará voltada para o cálculo
do centro de massa de perfis planos bem como no cálculo de seus
momentos de inércia.
ATIVIDADES
Deixamos como atividades as seguintes questões.
ATIV. 2.1. Determine a área da parte da cardióide r(ϑ) = 1 +

cos(ϑ) que fica acima do eixo dos x (Fig 2.8) que está em cinza.
Comentário: Volte ao texto e reveja com calma e atenção os
exemplos acima, eles lhe servirão de guia.
ATIV. 2.2. Determine a área entre a cardióide r(ϑ) = 1+cos(ϑ) e

o círculo r(ϑ) = 1 acima do eixo do x (Fig 2.9) que está em cinza.
exemplos acima, eles lhe servirão de guia.
45
Figura 2.8: Atividade 1 Figura 2.9: Atividade 2

2003.
46
AULA
Algumas Aplicações da
Integral Dupla 3
META:
Apresentar algumas aplicações das integrais duplas de funções de
OBJETIVOS:
Determinar área, massa, centro de massa, momento de massa e
momento de inércia de figuras planas usando integrais duplas de
funções de valores reais e domínio em R2 .
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I, curvas em R2 e coordenadas
polares da disciplina Cálculo II e integrais duplas aula 01 e aula 02.
Algumas Aplicações da Integral Dupla
3.1 Introdução
Caros alunos nesta terceira aula do nosso curso de Cálculo III

com o tema “Algumas Aplicações das Integrais Duplas”. Dentre as
inúmeras aplicações da integral dupla, veremos apenas duas pelo
pouco tempo que dispomos. Veremos apenas como usar as inte-
grais duplas para calcular a massa de uma região plana dada sua
distribuição de densidade e como calcular seu centro de gravidade.
Para outras aplicações recomendo uma busca na INTERNET
3.2 Preliminares
Consideraremos uma região D ⊂ R2 finita, com uma distri-

buição de densidade mássica superficial (massa por unidade de
superfície) %(x, y), ∀(x, y) ∈ D.
Determinação da massa
Para determinar a massa consideremos uma função Φ definida
em um domínio retangular R = {(x, y)  ∈ R2 |a ≤ x ≤ b ∧ c ≤
 %(x, y) , (x, y) ∈ D
y ≤ d} tal que D ⊂ R e Φ(x, y) = .
 0 , (x, y) ∈
/D
Considerando a uma partição para o retângulo R dada por P =
P [R] = P [a, b] × P [c, d], o produto cartesiano das partições P [a, b]
e P [c, d] onde P [a, b] = {x0 = a, x1 , . . . , xj , xj+1 , . . . , xm = b}
e P [c, d] = {y0 = c, y1 , . . . , yk , yk+1 , . . . , yn = d}. Tomamos um
ponto (ξj , ζk ) ∈ [xj−1 , xj ] × [yk−1 , yk ] em cada pequeno retângulo
e definimos a seguinte soma de Riemann:
m X
X n
Smn = Φ(ξj , ζk )∆Ajk .
j=1 k=1
48
Cálculo III AULA
A massa da região D, denotada m(D), será a integral dupla da fun-
2
Z Z
3
ção %(x, y) sobre o domínio D ⊂ R , denotada %(x, y)dxdy
D
será então definida como o seguinte limite:
Z Z
def
m(D) = %(x, y)dxdy = lim Smn
D |P |→0
OBS 3.1. Para a determinação do peso da região D toma-se a

seguinte soma de Riemann:
m X
X n
Smn = g(ξj , ζk )Φ(ξj , ζk )∆Ajk
j=1 k=1
onde g(ξj , ζk ) é a aceleração da gravidade no ponto (ξj , ζk ). E o

peso da região D, denotado p(D), será dado pela integral dupla:
Z Z
def
p(D) = g(x, y)%(x, y)dxdy = lim Smn
D |P |→0
Determinação do Momento de Massa

Usando as mesmas considerações acima para o cálculo da massa
de uma região D limitada com distribuição de densidade %(x, y).
Para calcular o momento de massa de um pequeno retângulo com
relação ao eixo y tomamos o seguinte produto ξj Φ(ξj , ζk )∆Ajk . O
momento de massa total em relação ao eixo y para a região D será
aproximado pelo limite da soma de Riemann:
m X
X n
Smn = ξj Φ(ξj , ζk )∆Ajk
j=1 k=1
.
O momento de massa da região D em relação ao eixo y será dada
49
pelo limite:
Z Z
def
My (D) = x%(x, y)dxdy = lim Smn
D |P |→0
De forma semelhante chega-se ao momento de massa da região D

em relação ao eixo x tomando-se a seguinte soma de Riemann:
m X
X n
Smn = ζk Φ(ξj , ζk )∆Ajk
j=1 k=1
.
O momento de massa da região D em relação ao eixo x será dada
pelo limite:
Z Z
def
Mx (D) = y%(x, y)dxdy = lim Smn
D |P |→0
Determinação do Centro de Massa

O centro de massa de uma região plana D ⊂ R2 finita, com uma
distribuição de densidade mássica superficial %(x, y), ∀(x, y) ∈ D,
é o ponto (x̄, ȳ)Zdefinido
Z por:
x%(x, y)dxdy
My (D)
x̄ = = Z ZD
m(d)
%(x, y)dxdy
Z ZD
y%(x, y)dxdy
Mx (D) D
ȳ = = Z Z
m(d)
%(x, y)dxdy
D
Determinação do Momento de Inércia

50
Cálculo III AULA
de uma região D limitada com distribuição de densidade %(x, y). 3
Para calcular o momento de inércia de um pequeno retângulo com
relação ao eixo y tomamos o seguinte produto ξj2 Φ(ξj , ζk )∆Ajk . O
momento de inércia total em relação ao eixo y para a região D será
aproximado pelo limite da soma de Riemann:
m X
X n
Smn = ξj2 Φ(ξj , ζk )∆Ajk
j=1 k=1
.
O momento de inércia da região D em relação ao eixo y será dada
pelo limite:
Z Z
def
Iy (D) = x2 %(x, y)dxdy = lim Smn
D |P |→0
De forma semelhante chega-se ao momento de inércia da região D

em relação ao eixo x tomando-se a seguinte soma de Riemann:
m X
X n
Smn = ζk2 Φ(ξj , ζk )∆Ajk
j=1 k=1
.
O momento da região D em relação ao eixo x será dada pelo limite:
Z Z
def
Ix (D) = y 2 %(x, y)dxdy = lim Smn
D |P |→0
O momento de inércia em relação a origem é dado pela seguinte

integral dupla:
51
Z Z
I0 (D) = (x2 + y 2 )%(x, y)dxdy
D
3.3 Algumas Aplicações da Integral Dupla
Faremos duas aplicações da integral dupla ao cálculo do centro

de massa de duas figuras planas. Na primeira usaremos o sistema
de coordenadas cartesiano. Na segunda usaremos uma mudança
de variáveis para o sistema de coordenadas polares.
Vamos aos nossos exemplos.
Exemplo 3.1. Para o primeiro exemplo desejamos determinar o

centro de massa de uma região triangular D dada pela interseção
das retas x = 0, y = 0 e a reta que passa pelos pontos (0, a) e
(b, 0) com a, b > 0 (Fig 3.1), cuja densidade superficial de massa
é constante %(x, y) = %.
52
Cálculo III AULA
SOLUÇÃO: 3
Começaremos por determinar os limites de integração inspecio-
x
nando a (Fig 3.1) e verificando que 0 ≤ x ≤ a e 0 ≤ y ≤ b 1 − .
a
Em segundo calcularemos a massa da região D, m(D) e os respec-
tivos momentos de massa com relação ao eixo x e ao eixo y Mx (D)
e My (D) respectivamente.
Passo 1 determinar a massa m(D), dada pela integral dupla:
Z Z Z a Z b(1−x/a)
m(D) = %(x, y)dxdy = %dydx
D 0 0
Integrando em y temos:
Z a
b(1−x/a)
m(D) = % y dx
0 0
Z a
x
m(D) = % b 1 − dx
0 a
Integrando em x temos:
x2 a
m(D) = %b x −
2a 0

a2
m(D) = %b a −
2a
Simplificando temos:
ab
m(D) = %
2
Passo 2 calcular o momento de massa Mx (D) dado pela integral
dupla:
Z Z
Mx (D) = %(x, y)ydxdy
D
Substituindo os limites temos:
Mx (D) = %(x, y)ydxdy = %ydydx
D 0 0
Integrando em y teremos:
Z a 2
y b(1−x/a)
Mx (D) = % dx
0 2 0
Z a
(b(1 − x/a))2
Mx (D) = dx
0 2
Simplificando o integrando temos:
53
a
b2 b2 x b2 x2
Z
Mx (D) = % − + dx
0 2 a 2a2
Integrando em x teremos:
b2 x b2 x2 b2 x3 a
Mx (D) = % − +
2 2a 6a2 0

b2 a b2 a2 b2 a3
Mx (D) = % − +
2 2a 6a2
Simplificando as frações temos:
b2 a
Mx (D) = %
6
Passo 3 calcular o momento de massa My (D) dado pela integral
dupla:
Z Z
My (D) = %(x, y)xdxdy
Z ZD
D
Substituindo os limites temos:
My (D) = %(x, y)xdxdy = %xdydx
D 0 0
Integrando em y teremos:
Z a b(1−x/a)
My (D) = %xy dx

0 0
Z a
x
My (D) = %bx 1 − dx
0 a
Integrando em x teremos:
x2 x3 a
My (D) = %b −
2 3a 0

a2 a3
My (D) = %b −
2 3a
Simplificando as frações temos:
ba2
My (D) = %
6
Passo 4 Determinar o centro de massa de D pelas fórmulas:
My (D) Mx (D)
x̄ = e ȳ = .
m(D) m(D)
Usando os resultados anteriores temos:
54
Cálculo III AULA
%
ba2
%
b2 a 3
x̄ = 6 e ȳ = 6
ab ab
% %
2 2
a b
x̄ = e ȳ =
3 3
Como segundo exemplo usaremos uma região em que o sistema de
coordenadas polares facilita os cálculos.
Exemplo 3.2. Para o segundo exemplo desejamos determinar o

centro de massa de uma região D dada pelo quarto da coroa cir-
cular de raio interno a e raio externo b que situa-se no primeiro
quadrante (Fig 3.2), cuja densidade superficial de massa é cons-
tante %(x, y) = %.
SOLUÇÃO:
nando a (Fig 3.2) e verificando que 0 ≤ ϑ ≤ π/2 e a ≤ r ≤ b.
Em segundo calcularemos a massa da região D, m(D) e os respec-
tivos momentos de massa com relação ao eixo x e ao eixo y Mx (D)
e My (D) respectivamente.
55
Passo 1 determinar a massa m(D), dada pela integral dupla:

Z Z Z π/2 Z b
m(D) = %(x, y)dxdy = %rdrdϑ
D 0 a
Integrando em r temos:
Z π/2 2
r b
m(D) = % dϑ
0 2 a
Z π/2 2
b a2
m(D) = % − dϑ
0 2 2
Integrando em ϑ temos:
b2 a2 π/2
m(D) = % − ϑ
2 2 0
1
m(D) = %π(b2 − a2 )
4
Passo 2 calcular o momento de massa Mx (D) dado pela integral
dupla:
Z Z
D
Substituindo os limites em coordenadas polares e sabendo que
y = r sin(ϑ) temos:
Z Z Z π/2 Z b
Mx (D) = %(x, y)ydxdy = %r sin(ϑ)rdrdϑ
D 0 a
Z π/2
r3 b
Mx (D) = % sin(ϑ) dϑ
0 3 a
Z π/2
b3 a3
Mx (D) = % sin(ϑ) − dϑ
0 3 3
b3 a3 π/2
Mx (D) = % − (− cos(ϑ))

3 3 0
b3 a3
Mx (D) = % − (− cos(π/2) − − cos(0))
3 3
1
Mx (D) = %(b3 − a3 )
3
Passo 3 calcular o momento de massa My (D) dado pela integral
dupla:
56
Cálculo III AULA
My (D) =
Z Z
Z ZD
%(x, y)xdxdy
3
D
Substituindo os limites em coordenadas polares e sabendo que
x = r cos(ϑ) temos:
Z Z Z π/2 Z b
Mx (D) = %(x, y)ydxdy = %r cos(ϑ)rdrdϑ
D 0 a
Z π/2
r3 b
Mx (D) = % cos(ϑ) dϑ
0 3 a
Z π/2
b3 a3
Mx (D) = % cos(ϑ) − dϑ
0 3 3
b3 a3 π/2
Mx (D) = % − (sin(ϑ))

3 3 0
b3 a3
Mx (D) = % − (sin(π/2) − sin(0))
3 3
1
Mx (D) = %(b3 − a3 )
3
Passo 4 Determinar o centro de massa de D pelas fórmulas:
My (D) Mx (D)
x̄ = e ȳ = .
m(D) m(D)
Usando os resultados anteriores temos:
1 3
%(b − a3 )
x̄ = ȳ = 3
1
%π(b2 − a2 )
4
Levando em conta que b3 − a3 = (b − a)(b2 + ba + a2 ) e b2 − a2 =
(b − a)(b + a) temos:
1
%(b − a)(b2 + ba + a2 )
x̄ = ȳ = 3
1
%π(b − a)(b + a)
4
4 b2 + ba + a2
x̄ = ȳ = .
3π b+a
57
3.4 Conclusão
Na aula de hoje, vimos que dentre as inúmeras aplicações da

integral dupla, dentro da área da física destacamos, entre outras,
algumas das mais importantes que são: a determinação da massa
de uma região plana limitada por curvas, dada sua distribuição de
densidade, o cálculo do momento de massa de uma região plana
limitada por curvas, dada sua distribuição de densidade, o mo-
mento de inércia de uma região plana limitada por curvas, dada
sua distribuição de densidade e o cálculo do centro de massa de
uma região plana limitada por curvas, dada sua distribuição de
densidade.
RESUMO
Massa, Momento de Massa e Momento de Inércia

Dada uma região D ∈ R2 plana limitada com distribuição de densi-
dade superficial %(x, y) podemos calcular a massa de D, o momento
de massa em relação ao eixo x, o momento de massa relativo ao
eixo y, o momento de inércia em relação ao eixo x, o momento de
inércia relativo ao eixo y e momento de inércia relativo a origem,
denotados respectivamente m(D), Mx (D), My (D), Ix (D), Iy (D)
e I0 (D), pelas integrais duplas:
Z Z
m(D) = %(x, y)dxdy
Z DZ
Z ZD
Z Z D
Ix (D) = %(x, y)y 2 dxdy
D
58
Cálculo III AULA
Iy (D) =
Z Z
Z ZD
%(x, y)x2 dxdy e
3
I0 (D) = %(x, y)(x2 + y 2 )dxdy
D
Centro de Massa
Podemos também calcular o centro de massa, denotado (x̄, ȳ) usando
as seguintes fórmulas:
Z Z
x%(x, y)dxdy
My (D) D
x̄ = = Z Z
m(d)
%(x, y)dxdy
Z ZD
y%(x, y)dxdy
Mx (D)
ȳ = = Z ZD
m(d)
%(x, y)dxdy
D
PRÓXIMA AULA
Em nossa próxima aula veremos as integrais triplas. Primeira-

mente definindo-as para funções de domínios retangulares através
do limite de somas de riemann estendendo a definição para funções
definidas em domínios não retangulares porém limitados.
ATIVIDADES
Deixamos como atividades dois problemas de determinação do

centro de massa.
ATIV. 3.1. Determine o centro de massa da região D dada pela

interseção das retas y = 0, x = 1 e y = ax2 (Fig 3.3) região em
cinza.
Comentário: Volte ao texto e reveja com calma e atenção as
59
Figura 3.3: Atividade 1 Figura 3.4: Atividade 2
demonstrações acima, elas lhe servirão de guia. Use para este caso
coordenadas cartesianas.
ATIV. 3.2. Determine o centro de massa da região D dada pelo

semi-círculo superior x2 + y 2 = a2 (Fig 3.4) região em cinza.
demonstrações acima, elas lhe servirão de guia. Use para este caso
coordenadas polares.

60
Cálculo III AULA
THOMAS, George B., Cálculo, Volume 2, 10a, Addilson Wesley, 3
2003.
61
AULA
Integrais triplas
4
META:
Apresentar integrais triplas de funções de valores reais e domínio
em R3 .
OBJETIVOS:
Definir integral tripla e calcular algumas integrais triplas de fun-
ções de valores reais e domínio em R3 .
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I.
Integrais triplas
4.1 Introdução
Caros alunos a quarta aula do nosso curso de Cálculo III com o

tema “Integrais Triplas”. Bem como a integral dupla, vista na nossa
primeira aula, a integração tripla, em essência, é uma extensão
natural da integral simples vista em Cálculo I e definida como
limite de somas de Riemann. Na prática, a integração tripla é dada
HISTÓRIA por três integrações simples, cada uma efetuada sobre uma variável
A primeira técnica e considerando as demais como constantes. É o que denominamos
sistemática documen-
tada para o cálculo de integrais interadas. As características e detalhes próprios das
de integrais triplas no
cálculo de volume foi integrais triplas serão vistas ao longo do nosso curso, nas próximas
o método da exaustão
de Eudoxus cerca três aulas.
de 370AC. O maior
avanço no cálculo
de integrais triplas
veio do Iraque, no
século 11, na figura 4.2 Integração Tripla: Domínios Paralelepípe-
de Ibn AL-Haythan
(conhecido na Europa dais
por Alhazen ). En-
quanto resolvia o que
ficou conhecido como Começamos por considerar uma função φ definida em um do-
“Problema de Alhazen”
(um problema de ótica) mínio paralelepipedal R = {(x, y, z) ∈ R3 |a ≤ x ≤ b ∧ c ≤ y ≤
ele calculou o volume
de um parabolsóide d ∧ e ≤ z ≤ f }. Formalmente φ : [a, b] × [c, d] × [e, f ] 7→ R.
usando um método de
indução. Wikipédia. Usando a imaginação, pensemos em R retalhada por uma rede
de planos paralelos aos planos coordenados e que dividem R em
pequenos paralelepípedos. Oficialmente, consideraremos três par-
tições P [a, b] = {x0 = a, x1 , . . . , xi , xi+1 , . . . , xl = b}, P [c, d] =
{y0 = c, y1 , . . . , yj , yj+1 , . . . , ym = d} e P [e, f ] = {z0 = e, z1 ,
. . . , zk , zk+1 , . . . , zn = f } onde como visto em Cálculo I temos:
x0 < x1 < · · · < xi < xi+1 < · · · < xl , y0 < y1 < · · · < yj <
yj+1 < · · · < ym e z0 < z1 < · · · < zk < zk+1 < · · · < zn . Desta
forma cada um dos pequenos subintervalos Ii = [xi−1 , xi ], Jj =
64
Cálculo III AULA
[yj−1 , yj ] e Kk = [zk−1 , zk ] têm comprimentos ∆xi = xi − xi−1 , 4
∆yj = yj − yj−1 e ∆zk = zk − zk−1 , respectivamente. Definimos,
agora, a uma partição para o paralelepípedo R por P = P [R] =
P [a, b]×P [c, d]×P [e, f ], o produto cartesiano das partições P [a, b],
P [c, d] e P [e, f ]. Os planos retalham a região R em uma série de pe-
quenos paralelepípedos Vijk = [xi−1 , xi ] × [yj−1 , yj ] × [zk−1 , zk ], 1 ≤
i ≤ l, 1 ≤ j ≤ m, 1 ≤ k ≤ n. O volume de cada pequeno paralele-
pípedo é dado por ∆Vijk = ∆xi ∆yj ∆zk . Como tanto ∆xi quanto
∆yj quanto ∆zk são diferentes de zero, o volume de cada pequeno
paralelepípedo é também diferente de zero. Podemos então definir
a norma da partição por: |P | = max (∆Vijk ), que corresponde
1≤i≤l
1≤j≤m
1≤k≤n
ao maior volume entre todos os pequenos paralelepípedos.
Pausa para respirar que já vamos definir a integral tripla sobre do-
mínios paralelepípedais. Para isto tomemos um ponto (ξi , ζj , ηk ) ∈
[xi−1 , xi ] × [yj−1 , yj ] × [zk−1 , zk ] em cada pequeno paralelepípedo
e definimos a seguinte soma de Riemann:
l X
X m X
n
Slmn = φ(ξi , ζj , ηk )∆Vijk
i=1 j=1 k=1
A integral tripla da função φ(x, y, z) sobre o paralelepípedo R,

Z Z Z
denotada φ(x, y, z)dxdydz será então definida como o se-
R
guinte limite:
Z Z Z
def
φ(x, y, z)dxdydz = lim Slmn
R |P |→0
65
Integrais triplas
4.3 Integração Tripla: Domínios Não Paralele-

pípedais Limitados
Para definir a integral tripla de uma função φ : D ⊂ R3 7→

R onde D é limitado não paralelepipedal, começamos por con-
siderar uma função Φ definida em um domínio paralelepipedal
R = {(x, y, z) ∈ R3 |a ≤ x≤ b ∧ c ≤ y ≤ d ∧ e ≤ z ≤ f } tal
 φ(x, y, z) , (x, y, z) ∈ D
que D ⊂ R e Φ(x, y, z) = . Formal-
 0 , (x, y, z) ∈
/D
mente Φ : [a, b] × [c, d] × [e, f ] 7→ R é uma extensão da função
φ(x, y, z). Usando a imaginação, pensemos em R coberta por uma
rede de planos paralelos aos planos coordenados e que dividem
R em pequenos paralelepípedos e procedemos como na integral
tripla sobre domínios paralelepípedais, considerando a uma parti-
ção para o paralelepípedo R por P = P [R] = P [a, b] × P [c, d] ×
P [e, f ], o produto cartesiano das partições P [a, b], P [c, d] e P [e, f ]
onde P [a, b] = {x0 = a, x1 , . . . , xi , xi+1 , . . . , xl = b}, P [c, d] =
{y0 = c, y1 , . . . , yj , yj+1 , . . . , ym = d} e P [e, f ] = {z0 = e, z1 ,
. . . , zk , zk+1 , . . . , zn = f }. Do mesmo modo definimos a norma
da partição por: |P | = max (∆Vijk ) onde ∆Vijk = ∆xi ∆yj ∆zk ,
1≤i≤l
1≤j≤m
1≤k≤n
∆xi = xi − xi−1 , ∆yj = yj − yj−1 e ∆zk = zk − zk−1 . Tomamos
um ponto (ξi , ζj , ηk ) ∈ [xi−1 , xi ] × [yj−1 , yj ] × [zk−1 , zk ] em cada
pequeno paralelepípedo e definimos a seguinte soma de Riemann
para a função estendida Φ(x, y, z):
l X
X m X
n
Slmn = Φ(ξi , ζj , ηk )∆Vijk
i=1 j=1 k=1
66
Cálculo III AULA
A integral tripla da função φ(x, y, z) sobre o domínio D ⊂ R3 ,
Z Z Z
4
D
guinte limite:
Z Z Z
def
φ(x, y, z)dxdydz = lim Slmn .
D |P |→0
Observem que, semelhante ao caso das integrais duplas, apenas os

pequenos paralelepípedos cujo ponto escolhido pertence ao domínio
D ⊂ R3 , contribuem para a soma de Riemann os demais têm
contribuição nula visto que o ponto escolhido dentro destes estão
fora de D ⊂ R2 e portanto Φ(ξi , ζj , ηk ) = 0.
4.4 Interpretação Geométrica
Quando a função φ : D ⊂ R3 7→ R é constante e igual a um

(φ(x, y, z) = 1, ∀(x, y, z) ∈ D) e a região domínio D é limitada,
vemos que a soma de Riemann aproxima o volume da região D e
quanto maior for o refinamento da partição de R3 ⊃ R ⊃ D melhor
será a aproximação. Podemos então, interpretar a integral tripla
Z Z Z
dxdydz como o volume da região D ⊂ R3 .
D
4.5 Integrais Iteradas
Dada uma função φ : R 7→ R onde R = [a, b] × [c, d] × [e, f ], do

mesmo modo que na integral dupla, valem as integrais interadas:
Z Z Z Z bhZ dhZ f i i
1. φ(x, y, z)dxdydz = φ(x, y, z)dz dy dx
R a c e
Z Z Z Z bhZ f hZ d i i
2. φ(x, y, z)dxdydz = φ(x, y, z)dy dz dx
R a e c
67
Integrais triplas
Z Z Z Z dhZ bhZ f i i
3. φ(x, y, z)dxdydz = φ(x, y, z)dz dx dy
R c a e
Z Z Z Z dhZ f hZ b i i
4. φ(x, y, z)dxdydz = φ(x, y, z)dx dz dy
R c e a
Z Z Z Z f hZ dhZ b i i
5. φ(x, y, z)dxdydz = φ(x, y, z)dx dy dz
R e c a
Z Z Z Z f hZ bhZ d i i
6. φ(x, y, z)dxdydz = φ(x, y, z)dy dx dz
R e a c
Em outras palavras, quando o domínio da integral tripla é parale-

lepipedal a ordem de integração não importa.
4.6 Propriedades das Integrais Triplas
Como nosso curso é de Cálculo, apenas listaremos, sem de-

monstração, alguma das propriedades das integrais triplas. Caso
desejem conhecer a demonstração de algumas destas propriedades,
recomendo livros de Cálculo Avançado como os citados na biblio-
grafia abaixo.

Z Z Z Z Z Z
cf (x, y, z)dxdydz = c f (x, y, z)dxdydz
D D
Propriedade 4.2. Sejam f, g : D ⊂ R3 7→ R duas funções de

valores reais integráveis em D, então vale:
Z Z Z Z Z Z
(f + g)(x, y, z)dxdydz = f (x, y, z)dxdydz
D Z Z ZD
+ g(x, y, z)dxdydz
D
68
Cálculo III AULA
Propriedade 4.3. Sejam f : D ⊂ R3 7→ R uma função de valores 4
reais integrável em D tal que f (x, y, z) ≥ 0, ∀(x, y, z) ∈ D, então
vale:
Z Z Z
f (x, y, z)dxdydz ≥ 0
D
Propriedade 4.4. Sejam f, g : D ⊂ R3 7→ R duas funções de valo-

res reais integráveis em D tais que f (x, y, z) ≥ g(x, y, z), ∀(x, y, z) ∈
D, então vale:
Z Z Z Z Z Z
f (x, y, z)dxdydz ≥ g(x, y, z)dxdydz
D D
Propriedade 4.5. Seja f : D ⊂ R3 7→ R uma função de valores

número finito de superfícies em R3 , então vale:
Z Z Z Z Z Z
f (x, y, z)dxdydz = f (x, y, z)dxdydz
D Z Z ZA
+ f (x, y, z)dxdydz
B
OBS 4.1. As duas primeiras propriedades diz respeito à “line-

aridade” do operador integral tripla. As terceira e quarta pro-
priedades são denominadas “dominação” enquanto que a quinta
propriedade é denominada “aditividade”.
4.7 Exemplos
Nada mais natural que ilustrar um novo conceito com exem-

plos e, vamos aqui fazer exatamente isto. Ilustrar o conceito de
integral tripla com dois exemplos. Antes porém, vale observar
69
Integrais triplas
que a na prática uma integral tripla equivale a três integrais sim-

ples e neste caso uma pergunta fica no ar. Qual das duas variá-
veis x, y ou z integraremos primeiro? Muito bem, a resposta é
dada pela propria expressão da integral tripla. Isto é, na integral
Z Z Z
f (x, y, z)dxdydz primeiramente integramos na variável x,
R
depois na variável y e por último na variável z. Já na integral
Z Z Z
f (x, y, z)dzdydx primeiramente integramos na variável z,
R
depois na variável y e por último na variável x.
Exemplo 4.1. Considere a função f : [0, 1] × [0, 1] × [0, 1] 7→ R

dada por f (x, y) = x2 + y 2 + z 2 e determine a integral tripla
Z Z
I= f (x, y, z)dxdydz sobre a região R = {(x, y, z) ∈ R3 |0 ≤
R
x ≤ 1 ∧ 0 ≤ y ≤ 1 ∧ 0 ≤ z ≤ 1}.
SOLUÇÃO:
Passo 1 colocaremos os limites de integração que representam a
região R dada, segundo a ordem de integração:
Z 1Z 1Z 1
I= (x2 + y 2 + z 2 )dxdydz
0 0 0
Passo 2 integraremos na variável x considerando as variáveis y e

z como constantes:
Z 1Z 1 3
x
I= + y 2 x + z 2 x dydz
0 0 3
Z 1Z 1
13 03
I= − + y 2 (1 − 0) + z 2 (1 − 0) dydz
0 0 3 3
Z 1Z 1
1 2 2
I= + y + z dy
0 0 3
Passo 3 integraremos na variável y considerando a variável x

como constante:
70
Cálculo III AULA
I=
Z
0
1
1
3
y+
y3
3
2

+ z y dz
1
0
4
Z 1
13 03

1 2
I= (1 − 0) + − + z (1 − 0) dz
0 3 3 3
Efetuando os cálculostemos:
Z 1
1 1
I= + + z 2 dz
0 3 3
Passo
4 último passo, integraremos na variável z:
1 1 z 3 1
I= z+ z+
3 3 3

0
Substituindo os limites de integraçãotemos:
13 03

1 1
I= (1 − 0) + (1 − 0) + −
3 3 3 3
1 1 1
I= + + =1
3 3 3

limites de integração em uma integral tripla sobre domínio não re-
tangular da forma: D.

as superfícies inferior a(x, y) e superior b(x, y) que limitam a região
71
Integrais triplas
D, bem como a sombra projetada no plano xy por D, denotada

D∗ e identificar as curvas limites da região D∗ a(x) curva inferior
e b(x) curva superior, como na AULA01.
Passo 2 Atravessar toda a região D∗ e o eixo x com um segmento
r na Fig. 4.1)
Passo 3 Deslocar o segmento de reta r paralelo ao eixo y na di-
reção negativa do eixo x até tocar o ponto mais à esquerda de D∗
Passo 4 Deslocar o segmento de reta r paralelo ao eixo y na di-
reção positiva do eixo x até tocar o ponto mais à direita de D∗
mento de reta r através da região D∗ paralelo ao eixo y na direção
positiva do eixo x. O limite inferior para a variável y será a fun-
ção a(x), ponto da curva onde o segmento entra na região D∗ e o
limite superior para a variável y será b(x), ponto da curva onde o
segmento de reta sai da região D∗ .
Passo 6 Tomando um ponto qualquer (x, y) ∈ D∗ passamos o
segmento de reta s através da região D, paralelo ao eixo z orien-
tado na direção positiva de z. O limite inferior para a variável z
será a função a(x, y), ponto da superfície onde o segmento entra
na região D e o limite superior para a variável z será b(x, y), ponto
da superfície onde o segmento de reta sai da região D.
Z Z Z Z bZ b(x) Z b(x,y)
f (x, y)dxdy = f (x, y, z)dzdydx
D a a(x) a(x,y)
72
Cálculo III AULA
4
Vamos diretamente para um segundo exemplo de integral dupla
sobre domínios não retangulares. A saber:
Exemplo 4.2. Considere a função f : D ⊂ R3 7→ R dada por

Z Z
f (x, y) = xyz e determine a integral dupla f (x, y, z)dxdydz
D
sobre a região D = {(x, y, z) ∈ R3 |0 ≤ x ≤ 1 ∧ 0 ≤ y ≤ x2 ∧ 0 ≤
z ≤ 1}, (Fig. 4.2).
Figura 4.2: Domínio D para o exemplo 2
SOLUÇÃO:
Passo 1 faremos o desenho das superfícies que determinam os
limites para a região D. A saber x = 0, x = 1, y = x2 , x = 0 e
z = 1 (Fig. 4.2).
Usando o processo prático exposto acima determinamos os limi-
tes de integração. A saber: a = 0, b = 1, a(x) = 0, b(x) = x2 ,
a(x, y) = 0 e b(x, y) = 1.
Z 1 Z x2 Z 1
I= xyzdzdydx
0 0 0
Passo 2 integraremos na variável z considerando a variável y
73
Integrais triplas
e x como uma constante:

Z 1 Z x2
z 2 1

I= xy dydx
0 0 2 0
Z 1 Z x2
12 02

I= xy − xy ) dydx
0 0 2 2
Z Z 2
1 1 x
I= xydydx
2 0 0
Passo 3 integraremos na variável y considerando a variável x

constante temos:
Z 1 2 2
y x
I= x dx
0 2 0
Substituindo os limites deintegração temos:
1 1 (x2 )2 02
Z
I= x −x dx
2 0 2 2
1 1 5
Z
I= x dx
4 0
Integrando
6 , finalmente , na variável x temos:
1 x 1
I=
4 6 0

Substituindo os
limites de integração temos:
1 16 06

I= −
4 6 6
1
Efetuando os cálculos temos: I =
24
4.8 Conclusão
Na aula de hoje, vimos que a integral tripla é uma extensão

natural do conceito de integral simples visto em Cálculo I e também
uma extensão natural do conceito de integral dupla, vista em nossa
primeira aula do curso de Cálculo III. E se por um lado a integral
simples pode ser interpretada como a área sob a curva descrita por
74
Cálculo III AULA
função positiva f (x) em um domínio [a, b] e a integral dupla pode 4
ser vista como o volume de um prisma reto limitado superiormente
pela a superfície descrita por uma função positiva f (x, y) e limitado
inferiormente pelo domínio [a, b] × [c, d], a integral tripla só tem
interpretação geométrica no caso simples em que f (x, y, z) = 1.
Neste caso a integral tripla representa o volume da região limitada
D ⊂ R3 .
RESUMO
Integração Tripla: Domínios Paralelepipedais

Considerando uma função φ definida em um domínio paralelepi-
pedal R = {(x, y, z) ∈ R3 |a ≤ x ≤ b ∧ c ≤ y ≤ d ∧ e ≤ z ≤ f }.
Podemos dividir R em pequenos paralelepípedos considerando os
planos paralelos ao planos cartesianos gerados pela partição P =
P [R] = P [a, b]×P [c, d]×P [e, f ], o produto cartesiano das partições
P [a, b], P [c, d] e P [e, f ] onde P [a, b] = {x0 = a, x1 , . . . , xi , xi+1 , . . . ,
xl = b}, P [c, d] = {y0 = c, y1 , . . . , yj , yj+1 , . . . , ym = d} e P [e, f ] =
{z0 = e, z1 , . . . , zk , zk+1 , . . . , zn = f }. Os planos retalham a região
R em uma série de pequenos paralelepípedos Vijk = [xi−1 , xi ] ×
[yj−1 , yj ] × [zk−1 , zk ], 1 ≤ i ≤ l, 1 ≤ j ≤ m, 1 ≤ k ≤ n. O volume
de cada pequeno paralelepípedo é dado por ∆Vijk = ∆xi ∆yj ∆zk .
A norma da partição fica estabelecida como: |P | = max (∆Vijk ).
1≤i≤l
1≤j≤m
1≤k≤n
Toma-se um ponto (ξi , ζj , ηk ) ∈ [xi−1 , xi ] × [yj−1 , yj ] × [zk−1 , zk ]
em cada pequeno paralelepípedo e definimos a seguinte soma de
Riemann:
75
Integrais triplas
l X
X m X
n
Slmn = φ(ξi , ζj , ηk )∆Vijk
i=1 j=1 k=1
A integral tripla da função φ(x, y, z) sobre o paralelepípedo R,

Z Z Z
R
guinte limite:
Z Z Z
def
R |P |→0
Integração Tripla: Domínios Não Paralelepípedais Limi-

tados
Para definir a integral tripla de uma função φ : D ⊂ R3 7→ R
onde D é não paralelepipedal limitado, começamos por conside-
rar uma função Φ definida em um domínio paralelepipedal R =
{(x, y, z) ∈ R3 |a ≤ x≤ b ∧ c ≤ y ≤ d ∧ e ≤ z ≤ f } tal que
 φ(x, y, z) , (x, y, z) ∈ D
D ⊂ R e Φ(x, y, z) = . Formalmente
 0 , (x, y, z) ∈
/D
Φ : [a, b] × [c, d] × [e, f ] 7→ R é uma extensão da função φ(x, y, z).
A partir daqui todo o procedimento é semelhante ao da definição
da integral tripla em domínios paralelepipedais. Podemos definir
a integral tripla de uma função φ(x, y, z) em um domínio não re-
tangular D por:
Z Z Z
def
D |P |→0
Pl Pm Pn
onde Slmn = i=1 j=1 k=1 Φ(ξi , ζj , ηk )∆Vijk é a soma de Rie-
mann para Φ(x, y, z.
Integrais Iteradas
As integrais iteradas dizem que em um domínio retangular R =
[a, b] × [c, d] × [e, f ] a ordem de execução das integrais simples não
76
Cálculo III AULA
alteram o valor da integral tripla, que pode ser representada por: 4
Z Z Z Z bhZ dhZ f i i
1. φ(x, y, z)dxdydz = φ(x, y, z)dz dy dx
R a c e
Z Z Z Z bhZ f hZ d i i
2. φ(x, y, z)dxdydz = φ(x, y, z)dy dz dx
R a e c
Z Z Z Z dhZ bhZ f i i
3. φ(x, y, z)dxdydz = φ(x, y, z)dz dx dy
R c a e
Z Z Z Z dhZ f hZ b i i
4. φ(x, y, z)dxdydz = φ(x, y, z)dx dz dy
R c e a
Z Z Z Z f hZ dhZ b i i
5. φ(x, y, z)dxdydz = φ(x, y, z)dx dy dz
R e c a
Z Z Z Z f hZ bhZ d i i
6. φ(x, y, z)dxdydz = φ(x, y, z)dy dx dz
R e a c
Propriedades das Integrais triplas

As integrais triplas são, de certo modo, semelhantes às proprie-
dades das integrais simples que vimos em Cálculo I sendo quase
que uma extensão natural destas. As integrais triplas têm, entre
outras, as seguintes propriedades:
Z Z Z Z Z Z
cf (x, y, z)dxdydz = c f (x, y, z)dxdydz
D D

reais integráveis em D, então vale:
Z Z Z Z Z Z
(f + g)(x, y, z)dxdydz = f (x, y, z)dxdydz
D Z Z ZD
+ g(x, y, z)dxdydz
D
77
Integrais triplas

reais integrável em D tal que f (x, y, z) ≥ 0, ∀(x, y, z) ∈ D, então
vale:
Z Z Z
f (x, y, z)dxdydz ≥ 0
D

reais integráveis em D tais que f (x, y, z) ≥ g(x, y, z), ∀(x, y, z) ∈
D, então vale:
Z Z Z Z Z Z
f (x, y, z)dxdydz ≥ g(x, y, z)dxdydz
D D
Propriedade 5 Seja f : D ⊂ R3 7→ R uma função de valores

número finito de superfícies em R3 , então vale:
Z Z Z Z Z Z
f (x, y, z)dxdydz = f (x, y, z)dxdydz
D Z Z ZA
+ f (x, y, z)dxdydz
B
Determinação dos Limites de Integração para Integrais

Triplas
Daremos aqui um método prático para determinar os limites de
integração em uma integral tripla sobre domínio não retangular da
forma: D.

as superfícies inferior a(x, y) e superior b(x, y) que limitam a região
D, bem como a sombra projetada no plano xy por D, denotada
78
Cálculo III AULA
D∗ e identificar as curvas limites da região D∗ a(x) curva inferior 4
e b(x) curva superior, como na AULA01.
Passo 2 Atravessar toda a região D∗ e o eixo x com um seg-

mento de reta paralelo e orientado na direção positiva ao eixo y
(segmento r na Fig. 4.1)
Passo 3 Deslocar o segmento de reta r paralelo ao eixo y na

direção negativa do eixo x até tocar o ponto mais à esquerda de
D∗ marcando o limite inferior de x (ponto a na Fig. 4.1).
Passo 4 Deslocar o segmento de reta r paralelo ao eixo y na

direção positiva do eixo x até tocar o ponto mais à direita de D∗
Passo 5 Tomando um ponto qualquer x ∈ (a, b) passamos o

segmento de reta r através da região D∗ paralelo ao eixo y na di-
reção positiva do eixo x. O limite inferior para a variável y será a
função a(x), ponto da curva onde o segmento entra na região D∗ e
o limite superior para a variável y será b(x), ponto da curva onde
o segmento de reta sai da região D∗ .
Passo 6 Tomando um ponto qualquer (x, y) ∈ D∗ passamos o

segmento de reta s através da região D, paralelo ao eixo z orien-
tado na direção positiva de z. O limite inferior para a variável z
será a função a(x, y), ponto da superfície onde o segmento entra
na região D e o limite superior para a variável z será b(x, y), ponto
da superfície onde o segmento de reta sai da região D.
79
Integrais triplas
Z Z Z Z bZ b(x) Z b(x,y)
f (x, y)dxdy = f (x, y, z)dzdydx
D a a(x) a(x,y)
PRÓXIMA AULA
Em nossa próxima aula veremos mudança de variáveis na in-

tegração tripla. O objetivo da mudança de variáveis em uma in-
tegral tripla será a de facilitar esta integração de uma de duas
formas. A primeira será tornando o integrando mais simples. A
segunda transformando o domínio D do integrando em um domínio
de forma geométrica mais simples.
ATIVIDADES
Deixamos como atividades o cálculo de algumas integrais trí-

plas.
ATIV. 4.1. Seja f : [−1, +1] × [−1, +1] × [−1, +1] 7→ R dada por
f (x, y, z) = x2 + y 2 + z 2 . Determine a integral tripla:
Z Z Z
f (x, y, z)dxdydz.
R

de guia.
ATIV. 4.2. Seja f : D ⊂ R3 7→ R dada por f (x, y, z) = 1, onde

D = {(x, y, z) ∈ R3 |x ≥ 0 ∧ 0 ≤ y ≤ 1 − x2 ∧ 0 ≤ z ≤ 1 − x2 }.
80
Cálculo III AULA
• Esboce a região de integração 4
• Determine os limites da integral tripla:
Z Z Z
f (x, y, z)dxdydz
D
Z Z Z
• Calcule a integral tripla f (x, y, z)dxdydz.
D

de guia.

2003.
81
AULA
Mudança de Variáveis
em Integrais tríplas 5
META:
Introduzir mudança de variáveis em integrais triplas de funções de
OBJETIVOS:
Calcular integrais triplas de funções de valores reais e domínio em
R3 utilizando mudança de variáveis.
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I, superfícies em R3 , de coorde-
nadas polares da disciplina Cálculo II e integrais triplas aula 04.
Mudança de Variáveis em Integrais tríplas
5.1 Introdução
Caros alunos o problema da mudança de variáveis em integrais

triplas é inteiramente análogo ao problema de mudança de variá-
veis em integrais duplas. Analogias a parte, o fato de do espaço R3
ter uma dimensão a mais que o R2 , traz um esforço algébrico adi-
cional ao tratamento geral da mudança de variáveis em integrais
HISTÓRIA triplas. Veremos dois casos particulares de mudança de variáveis
O teorema de mu- em integrais tripla que correspondem aos: sistemas de coordenadas
dança de variáveis em
integrais tríplas foi cilíndricos e sistema de coordenadas esféricos.
primeiro proposto por
Lagrange em 1773 e
usado por Legendre,
Laplace e Gauss, e
primeiramente ge- 5.2 Mudança de Variáveis em Integrais Triplas
neralizado para n
variáveis por Mikhail
Ostrogradski em 1836, Vamos considerar a integração de uma função f : D ⊂ R3 7→
resistiu a uma demons-
tração mais rigorosa R onde (x, y, z) ∈ D e conseqüentemente, ∀(x, y, z) ∈ D temos
por longo tempo (cerca
de 125 anos). E f (x, y, z) ∈ R. Consideraremos também, uma transformação T :
foi satisfatóriamente
demonstrado por Elie D ⊂ R3 7→ D0 ⊂ R3 , biunívoca de modo que D = T −1 (D0 ),
Cartan em uma série
∀(u, v, w) ∈ D0 , (x, y, z) = T −1 (u, v, w) ∈ D. Trocando em miú-
de artigos nos anos
1890. dos: x = x̂(u, v, w), y = ŷ(u, v, w) e z = ẑ(u, v, w). E suponhamos
as funções contínuas
ederiváveis e seu jacobiano, denotado J, de-
x, y, z ∂(x, y, z)
finido por: J ou :
u, v, w ∂(u, v, w)
∂ x̂ ∂ ŷ ∂ ẑ
 
 ∂u ∂u ∂u 
x, y, z ∂(x, y, z)  ∂ x̂ ∂ ŷ ∂ ẑ 
J = = det  .
u, v, w ∂(u, v, w)  ∂v
 ∂ x̂ ∂v ∂v 
∂ ŷ ∂ ẑ 
∂w ∂w ∂w
Suponhamos uma partição de D0 feita partindo de planos para-
lelos aos planos coordenados vw (u constante), uw (v constante) e
uv (w constante). Denotando ui+1 = ui + ∆ui , vj+1 = vj + ∆vj e
84
Cálculo III AULA
wk+1 = wk + ∆wk , destacamos o pequeno paralelepípedo indexado 5
por ijk, (Fig 5.1). Suponhamos que
Figura 5.1: Elemento de vo- Figura 5.2: Elemento de vo-

lume em D0 lume em D
0 = ∆u ∆v ∆w seja
este pequeno paralelepípedo de volume ∆Vijk i j k
mapeado por T −1 em um subdomínio em D de volume ∆Vijk (Fig

5.2). Seja: P = P (u, v, w) = (x̂(u, v, w), ŷ(u, v, w), ẑ(u, v, w)).
Os segmentos de reta (u, vj , wk ), (ui , v, wk ) e (ui , vj , w) come-
çando no ponto (ui , vj , wk ) são mapeados por T −1 em P (u, vj , wk )
P (ui , v, wk ) P (ui , vj , w) ver (Fig 5.2).
∂P
No subdomínio Vijk ⊂ D traçamos os vetores tangentes ∆ui ,
∂u
∂P ∂P
∆vj e ∆wk , ver (Fig 5.3). Em seguida traçamos segmentos
∂v ∂w
de reta paralelos aos vetores tangentes completando um paralele-
pípedo em D, ver (Fig 5.4), cujo volume admitiremos aproxima-
damente igual ao ∆Vijk (esta é a argumentação heurística). Este
volume é dado por:

∂P ∂P ∂P
∆Vijk ≈
∆ui × ∆vj • ∆wk .
∂u ∂v ∂w
85
Levando em conta que:

∂P ∂ x̂ ∂ ŷ ∂ ẑ ~
= ~i + ~j + k
∂u ∂u ∂u ∂u
∂P ∂ x̂ ∂ ŷ ∂ ẑ
= ~i + ~j + ~k
∂v ∂v ∂v ∂v
∂P ∂ x̂~ ∂ ŷ ~ ∂ ẑ ~
= i+ j+ k
∂w ∂w ∂w ∂w
e calculando o produto vetorial mixto teremos:
Figura 5.3: Elemento de vo- Figura 5.4: Elemento de vo-

lume em D0 lume em D
∂ x̂ ∂ ŷ ∂ ẑ
 
 ∂u ∂u ∂u 
∂P ∂P ∂P  ∂ x̂ ∂ ŷ ∂ ẑ 
∆ui × ∆vj • ∆wk = det   ∆ui ∆vj ∆wk
∂u ∂v ∂w  ∂v ∂v ∂v 
 ∂ x̂ ∂ ŷ ∂ ẑ 
∂w ∂w ∂w
Daí, levando em conta a expressão do jacobiano em R3 , dada
acima, temos:

∂(x̂, ŷ, ẑ)
∆Vijk ≈
∆ui ∆vj ∆wk
∂(u, v, w)
O que nos leva à seguinte expressão para a mudança de variáveis
em integrais triplas:
Z Z Z Z Z Z
f (x, y, z)dxdydz = F (u, v, w) |J| dudvdw
D D0
86
Cálculo III AULA
onde: F (u, v,w) = f (x̂(u,
x̂, ŷ, ẑ
v, w), ŷ(u, v, w), ẑ(u, v, w)) e J é o jaco-
5
biano J = J .
u, v, w
5.3 Alguns Exemplos
Nesta seção veremos dois exemplos de integrais triplas com

mudança de variáveis. No primeiro aplicaremos a mudança de va-
riáveis dada pelo sistema de coordenadas cilíndricas e no segundo
o sistema de coordenadas esféricas
Primeiramente veremos um exemplo em coordenadas cilíndricas.

Antes porém, veremos como determinar os limites de integração
em coordenadas cilíndricas.
Figura 5.5: Coordenadas ci- Figura 5.6: Coordenadas ci-

líndricas 1 líndricas 2
Passo 1 Esboçar o domínio D bem como sua projeção D∗ no

plano xy (ver Fig. 5.5).
Passo 2 Identificar as curvas que limitam a região D∗ . Atravessar
a região D∗ com uma reta r começando na origem (ver Fig. 5.6).
87
À medida em que a reta r percorre a região D∗ o ângulo ϑ que ela

forma com o eixo x positivo varia do mínimo α que será o limite
inferior da variável ϑ ao máximo β que será o limite superior da
variável ϑ. Daí, a variável ϑ ∈ [α, β].
Passo 3 Para cada valor fixo da variável ϑ ∈ [α, β] atravessar a
região D∗ com a reta r (ver Fig. 5.7). O ponto onde a reta r
entra na região D∗ é o limite inferior α(ϑ) para a variável r e o
ponto onde a reta r sai da região D∗ é o limite inferior β(ϑ) para
a variável r. Daí, r ∈ [α(ϑ), β(ϑ)].
Passo 4 Para cada valor fixo da variável ϑ ∈ [α, β] e da variável
r ∈ [α(ϑ), β(ϑ)] tomar o ponto (r, ϑ) ∈ D∗ em coordenadas pola-
res e levantar a reta s atravessando a região D (ver Fig. 5.8). O
ponto onde a reta s entra na região D é o limite inferior α(r, ϑ)
para a variável z e o ponto onde a reta s sai da região D é o limite
superior β(r, ϑ) para a variável z. Daí, z ∈ [α(r, ϑ), β(r, ϑ)].
Figura 5.7: Coordenadas ci- Figura 5.8: Coordenadas ci-

líndricas 3 líndricas 4
Podemos agora encarar o nosso primeiro exemplo onde colocaremos
88
Cálculo III AULA
em prática a determinação dos limites de integração em coordena- 5
das cilíndricas.
Exemplo 5.1. Considere o sólido gerado pela intersecção das su-

perfícies: z = y + a, (plano) x2 + y 2 − 2ay = 0, (cilindro) e z = 0,
(plano) (Fig 5.9) e determine seu volume.
Figura 5.9: Superfícies do Figura 5.10: Interseção das

exemplo 1 superfícies do exemplo 1
SOLUÇÃO: Para uma melhor compreensão mostramos na (Fig

5.10) o sólido gerado pela interseção das superfícies dadas e na
(Fig 5.11) as superfícies que compõem o sólido separadas no es-
paço.
Usaremos para o caso o sistema de coordenadas cilíndricas, dada

pela transformação: (x, y, z) 7→ (r, ϑ, z) onde x = r cos(ϑ), y =
r sin(ϑ) e z = z. O jacobiano da transformação é dado por:
89
∂ x̂ ∂ ŷ ∂ ẑ
 
 ∂r ∂r ∂r 
x, y, z  ∂ x̂ ∂ ŷ ∂ ẑ 
J =J = det  
r, ϑ, z  ∂ϑ
 ∂ x̂ ∂ϑ ∂ϑ 
∂ ŷ ∂ ẑ 
∂z ∂z ∂z
Efetuando as derivadas parciais temos:
 
cos(ϑ) sin(ϑ) 0
 
J = det  −r sin(ϑ) r cos(ϑ) 0 
 
 
0 0 1
Fazendo as contas do determinante temos:
J =r
Aproveitaremos o exemplo para aplicar os passos, vistos acima,

para determinação dos limites de integração de uma integral tripla
no sistema de coordenadas cilíndricas.
Passo 1: Esboçar a interseção das superfícies (sólido D), bem
como sua projeção sobre o plano xy (superfície D∗ ), (ver Fig 5.10).
A projeção sobre o plano xy (superfície D∗ ), conhecide com a su-
perfície inferior do sólido, sendo o disco dado por x2 +y 2 −2ay ≤ 0.
90
Cálculo III AULA
Passo 2: Os limites para r e ϑ são determinados em D∗ do mesmo 5
modo que para coordenadas polares em R2 . Neste caso 0 ≤ ϑ ≤ 2π
e para r temos: que r vai de zero até a borda de D∗ que é dada
por x2 + y 2 − 2ay = (r cos(ϑ))2 + (r sin(ϑ))2 − 2ar sin(ϑ) = 0.
Daí, r2 − 2ar sin(ϑ) = r(r − 2a sin(ϑ)) = 0 Simplificando temos:
0 ≤ r ≤ 2a sin(ϑ).
Passo 3: Para determinar os limites para z. Por cada par
(r, ϑ) ∈ D∗ , traçamos uma reta paralela ao eixo z orientada no
sentido positivo do eixo z atravessando o sólido. O limite inferior
de z é o ponto onde a reta entra no sólido e o limite superior o
ponto onde a reta sai do sólido. Neste caso: 0 ≤ z ≤ a + x ou
como x = r cos(ϑ) temos: 0 ≤ z ≤ a + r cos(ϑ).
Daí, o cálculo do volume de D será dado pela integral:
Z 2π Z 2a sin(ϑ) Z a+r sin(ϑ)
V ol(D) = rdzdrdϑ
0 0 0
Integrando na variável z temos:
Z 2π Z 2a sin(ϑ) a+r sin(ϑ)

V ol(D) = rz drdϑ

0 0 0

Z 2π Z 2a sin(ϑ)
V ol(D) = r(a + r sin(ϑ) − 0)drdϑ
0 0
Fazendo as contas temos:

Z 2π Z 2a sin(ϑ)
V ol(D) = (ar + r2 sin(ϑ))drdϑ
0 0
Integrando em na variável r temos:
2π
r2 r3
Z 2a sin(ϑ)
V ol(D) = (a + sin(ϑ)) dϑ

0 2 3 0
91
Substituindo o limite superior pois, o limite inferior por ser r = 0

não contribui, temos:
Z 2π
(2a sin(ϑ))2 (2a sin(ϑ))3
V ol(D) = (a + sin(ϑ))dϑ
0 2 3
2π
8a3 sin4 (ϑ)
Z
V ol(D) = (2a3 sin2 (ϑ) + )dϑ
0 3
Reescrevendo temos:
2π 2π
8a3
Z Z
3 2
V ol(D) = 2a sin (ϑ)dϑ + sin4 (ϑ)dϑ
0 3 0
Das tabelas de integrais temos:

α sinn−1 (αu) cos(αu)
Z
sinn (αu)du = −
an
Z
n−1
+ sinn−2 (αu)du
n
Z
u sin(2αu)
sin2 (αu)du = −
2 4α
Dai, temos:
Z
ϑ sin(2ϑ)
sin2 (ϑ)dϑ = −
2 4
Z 3 Z
4 sin (ϑ) cos(ϑ) 3
sin (ϑ)dϑ = − + sin2 (ϑ)dϑ
4 4
sin3 (ϑ) cos(ϑ) 3 ϑ sin(2ϑ)

= − + −
4 4 2 4
Podemos agora calcular as integrais. Para a integral de sin2 (ϑ)
temos:
Z 2π
ϑ sin(2ϑ) 2π
sin2 (ϑ)dϑ = −
2 4

0 0
2π sin(4π)
= + −
2 4
0 sin(0)
− +
2 4
= π
92
Cálculo III AULA
Para a integral de sin4 (ϑ) temos:
Z 2π
5
sin3 (ϑ) cos(ϑ) 3 ϑ sin(2ϑ)

4 2π
sin (ϑ)dϑ = − + −
4 4 2 4

0 0
3
sin (2π) cos(2π) 3 2π sin(4π)
= + − + −
4 4 2 4
3
sin (0) cos(0) 3 0 sin(0)
− − + −
4 4 2 4
3π
=
4
Substituindo no cálculo de V ol(D) temos:
8a3 3π
V ol(D) = 2πa3 +
3 4
= 4πa3
Em nosso segundo exemplo utilizaremos coordenadas esféricas, An-

tes porém, veremos como determinar os limites de integração em
coordenadas esféricas.
Figura 5.12: Coordenadas es- Figura 5.13: Coordenadas es-

féricas 1 féricas 2

93
Passo 2 Identificar as curvas que limitam a região D∗ . Atravessar

a região D∗ com uma reta r começando na origem (ver Fig. 5.13).
Á medida em que a reta r percorre a região D∗ o ângulo ϑ que ela
forma com o eixo x positivo varia do mínimo α que será o limite
inferior da variável ϑ ao máximo β que será o limite superior da
variável ϑ. Daí, a variável ϑ ∈ [α, β].
Passo 3 Para cada valor fixo da variável ϑ ∈ [α, β] atravessar
Figura 5.14: Coordenadas es- Figura 5.15: Coordenadas es-

féricas 3 féricas 4
a região D com o plano P que contem o eixo z e forma ângulo ϑ

com o eixo x positivo (ver Fig. 5.14). Traçamos uma reta r que
começa na origem e está contida no plano que corta D. À medida
em que a reta r percorre a região D o ângulo ϕ que ela forma com
o eixo z positivo varia do mínimo α(ϑ) que será o limite inferior da
variável ϕ ao máximo β(ϑ) que será o limite superior da variável
ϕ. Daí, a variável ϕ ∈ [α(ϑ), β(ϑ)].
Passo 4 Para cada valor fixo da variável ϑ ∈ [α, β] e da variável
r ∈ [α(ϑ), β(ϑ)] plano P que contem o eixo z e forma ângulo ϑ com
o eixo x positivo. No plano P traçar a reta s que forma ângulo ϕ
94
Cálculo III AULA
com o eixo z positivo atravessando a região D (ver Fig. 5.15). O 5
ponto onde a reta s entra na região D é o limite inferior α(ϑ, ϕ)
para a variável r e o ponto onde a reta s sai da região D é o limite
superior β(ϑ, ϕ) para a variável r. Daí, r ∈ [α(ϑ, ϕ), β(ϑ, ϕ)].
Podemos agora encarar o nosso segundo exemplo onde colocaremos

em prática a determinação dos limites de integração em coordena-
das esféricas.
Exemplo 5.2. Considere o sólido gerado pela interseção das su-

p p
perfícies: z = x2 + y 2 (cone), z = a2 − x2 − y 2 (esfera)(Fig
5.16) e determine seu volume.
Figura 5.16: Superfícies de Figura 5.17: Interseção das

exemplo 2 superfícies de exemplo 2
SOLUÇÃO: Para uma melhor compreensão mostramos na (Fig

5.17) o sólido gerado pela interseção das superfícies dadas e na
(Fig 5.18) as superfícies que compõem o sólido separadas no es-
paço.
95
Usaremos para o caso o sistema de coordenadas esféricas, dada

pela transformação: (x, y, z) 7→ (r, ϑ, ϕ) onde x = r cos(ϑ) cos(ϕ),
y = r sin(ϑ) cos(ϕ) e z = r sin(ϕ). O jacobiano da transformação
é dado por:
 ∂ x̂ ∂ ŷ ∂ ẑ 
 ∂r ∂r ∂r 
x, y, z  ∂ x̂ ∂ ŷ ∂ ẑ 
J =J = det 
 ∂ϑ

r, ϑ, z  ∂ x̂ ∂ϑ ∂ϑ 
∂ ŷ ∂ ẑ 
∂ϕ ∂ϕ ∂ϕ
Efetuando as derivadas parciais temos:
 
cos(ϑ) cos(ϕ) sin(ϑ) cos(ϕ) sin(ϕ)
 
J = det  −r sin(ϑ) cos(ϕ) r cos(ϑ) cos(ϕ)
 
0 
 
−r cos(ϑ) sin(ϕ) −r sin(ϑ) sin(ϕ) r cos(ϕ)
Fazendo as contas do determinante temos:
J = r2 sin(ϕ)
Aproveitaremos o exemplo para aplicar os passos na determina-

ção dos limites de integração de uma integral tripla no sistema de
96
Cálculo III AULA
5
coordenadas esféricas expostos acima.

Passo 1: Esboçar a interseção das superfícies (sólido D), bem
como sua projeção sobre o plano xy (superfície D∗ ), ver (Fig
5.19). A projeção sobre o plano xy (superfície D∗ ), é dada por
x2 + y 2 ≤ a2 .
Passo 2: Os limites para a variável ϑ são determinados em D∗
∗
√ de coordenadas polares. No caso como D é
como em um sistema
a
um disco de raio temos que: 0 ≤ ϑ ≤ 2π.
2
Passo 3: Os limites para a variável ϕ são determinados em D
do seguinte modo: para cada valor fixo de ϑ, em D∗ , cortamos o
domínio D por um plano que passa no eixo z e forma ângulo ϑ com
o eixo x. Traçamos uma reta M que passa na origem, pertence ao
plano ϑ e atravessa o domínio D. O ângulo ϕ é o ângulo formado
por M e o eixo z positivo. Para o caso o menor valor é ϕ = 0,
quando M conhecide com o eixo Z e o maior valor de ϕ em D
p
é quando M conhecide com a geratriz do cone z = x2 + y 2 e
π
ϕ= .
4
Passo 4: Os limites para a variável r são determinados em D
97
do seguinte modo: para cada par fixo ϑ, ϕ) percorremos a reta

M partindo da origem. O limite inferior de r é o ponto onde a
reta entra em D e o limite superior o ponto onde M sai de D.
Para o nosso caso: 0 ≤ r ≤ a (a reta sai na superfície da esfera
p
z = a2 − x2 − y 2 .
Podemos determinar o volume de D pela integral tripla:
Z 2π Z π/4 Z a
V ol(D) = r2 sin(ϕ)drdϕdϑ
0 0 0
Integrando primeiramente na variável r temos:

Z 2π Z π/4 3
r a
V ol(D) = sin(ϕ)dϕdϑ
0 0 3 0

Z 2π Z π/4 3
03

a
V ol(D) = − sin(ϕ)dϕdϑ
0 0 3 3
2π π/4
a3
Z Z
V ol(D) = sin(ϕ)dϕdϑ
3 0 0
Integrando na variável ϕ temos:
a3 2π
Z π/4
V ol(D) = − cos(ϕ) dϑ

3 0 0
a3 2π
Z
V ol(D) = (− cos(π/4) + cos(0)) dϑ
3 0
√ Z
a3 2 − 2 2π
V ol(D) = dϑ
3 2 0
Integrando na variável ϑ temos:

√
a3 2 − 2 2π
V ol(D) = ϑ
3 2 0
98
Cálculo III AULA
Substituindo os limites de integração temos: 5
√
a3 2 − 2
V ol(D) = (2π − 0)
3 2
Finalmente, simplificando temos:

√
πa3 (2 − 2)
V ol(D) =
3
5.4 Conclusão
Na aula de hoje, vimos que algumas vezes é conveniente fazer

uma mudança nas variáveis de integração em uma integral tripla,
para facilitar o cálculo da mesma. Vimos em particularmente duas
mudanças de variáveis são muito importantes e correspondem aos:
sistema de coordenadas cilíndricas e sistema de coordenadas esfé-
ricas.
RESUMO
Para o nosso resumo da aula 05 necessitamos algumas consi-

derações iniciais para tratar da mudança de variáveis em integrais
triplas. A saber:
Consideramos a transformação (x, y, z) = T (u, v, w) tal que o do-
mínio um ponto do domínio D ⊂ R3 , (x, y, z) seja transformado
no ponto (u, v, w) do domínio D0 ⊂ R3 , (D = T (D0 )) e mais
especificamente x = x̂(u, v, w), y = ŷ(u, v, w) e z = ẑ(u,
v, w).De-
x, y, z
finindo o jacobiano da transformação, denotado J, J ou
u, v, w
∂(x, y, z)
, por:
∂(u, v, w)
99
∂ x̂ ∂ ŷ ∂ ẑ
 
 ∂u ∂u ∂u 
x, y, z ∂(x, y, z)  ∂ x̂ ∂ ŷ ∂ ẑ 
J = = det  
u, v, w ∂(u, v, w)  ∂v ∂v ∂v 
 ∂ x̂ ∂ ŷ ∂ ẑ 
∂w ∂w ∂w
Mudança de Variáveis em Integrais Triplas

Vale então,a seguinte fórmula para a mudança de variáveis en in-
tegrais duplas:
Z Z Z Z Z Z
f (x, y, z)dxdydz = F (u, v, w) |J| dudvdw
D D0
onde: F (u, v,w) = f (x̂(u,

v, w), ŷ(u, v, w), ẑ(u, v, w)) e J é o jaco-
x̂, ŷ, ẑ
biano J = J .
u, v, w
Sistema de Coordenadas Cilíndricas

O sistema de coordenadas cilíndricas, dado pela transformação:
(x, y, z) 7→ (r, ϑ, z) onde x = r cos(ϑ), y = r sin(ϑ) e z = z. O
jacobiano da transformação é dado por: J = r e a integral tripla
pela expressão:
Z Z Z Z Z Z
f (x, y, z)dxdydz = F (r, ϑ, z)rdzdrdϑ
D D0
onde: F (r, ϑ, z) = f (r cos(ϑ), r sin(ϑ)z)
Sistema de Coordenadas esféricas

O sistema de coordenadas esféricas, que é dado pela transformação:
(x, y, z) 7→ (r, ϑ, ϕ) onde x = r cos(ϑ) cos(ϕ), y = r sin(ϑ) cos(ϕ)
e z = r sin(ϕ). O jacobiano da transformação é dado por: J =
r2 sin(ϕ) e a integral tripla pela expressão:
Z Z Z Z Z Z
f (x, y, z)dxdydz = F (r, ϑ, ϕ)r2 sin(ϕ)drdϕdϑ
D D0
onde: F (r, ϑ, ϕ) = f (r cos(ϑ) cos(ϕ), r sin(ϑ) cos(ϕ), r sin(ϕ))
100
Cálculo III AULA
5
Determinação dos Limites para Integração em Coorde-
nadas Cilíndricas
Para determinação dos limites de integração tripla em coordenadas
cilíndricas utiliza-se os seguintes passos:

Passo 2 Identificar as curvas que limitam a região D∗ . Atra-

vessar a região D∗ com uma reta r começando na origem (ver Fig.
5.6). Á medida em que a reta r percorre a região D∗ o ângulo ϑ
que ela forma com o eixo x positivo varia do mínimo α que será o
limite inferior da variável ϑ ao máximo β que será o limite superior
da variável ϑ. Daí, a variável ϑ ∈ [α, β].

a região D∗ com a reta r (ver Fig. 5.7). O ponto onde a reta r
entra na região D∗ é o limite inferior α(ϑ) para a variável r e o
ponto onde a reta r sai da região D∗ é o limite inferior β(ϑ) para
a variável r. Daí, r ∈ [α(ϑ), β(ϑ)].
Passo 4 Para cada valor fixo da variável ϑ ∈ [α, β] e da va-

riável r ∈ [α(ϑ), β(ϑ)] tomar o ponto (r, ϑ) ∈ D∗ em coordenadas
polares e levantar a reta s atravessando a região D (ver Fig. 5.8).
O ponto onde a reta s entra na região D é o limite inferior α(r, ϑ)
para a variável z e o ponto onde a reta s sai da região D é o limite
superior β(r, ϑ) para a variável z. Daí, z ∈ [α(r, ϑ), β(r, ϑ)].
101
Determinação dos Limites para Integração em Coordena-

das Esféricas
Para determinação dos limites de integração tripla em coordenadas
esféricas utiliza-se os seguintes passos:

Passo 2 Identificar as curvas que limitam a região D∗ . Atra-

vessar a região D∗ com uma reta r começando na origem (ver Fig.
5.13). Á medida em que a reta r percorre a região D∗ o ângulo ϑ
que ela forma com o eixo x positivo varia do mínimo α que será o
limite inferior da variável ϑ ao máximo β que será o limite superior
da variável ϑ. Daí, a variável ϑ ∈ [α, β].

a região D com o plano P que contem o eixo z e forma ângulo ϑ
com o eixo x positivo (ver Fig. 5.14). Traçamos uma reta r que
começa na origem e está contida no plano que corta D. À medida
em que a reta r percorre a região D o ângulo ϕ que ela forma com
o eixo z positivo varia do mínimo α(ϑ) que será o limite inferior da
variável ϕ ao máximo β(ϑ) que será o limite superior da variável
ϕ. Daí, a variável ϕ ∈ [α(ϑ), β(ϑ)].
Passo 4 Para cada valor fixo da variável ϑ ∈ [α, β] e da va-
102
Cálculo III AULA
riável r ∈ [α(ϑ), β(ϑ)] plano P que contem o eixo z e forma ângulo 5
ϑ com o eixo x positivo. No plano P traçar a reta s que forma
ângulo ϕ com o eixo z positivo atravessando a região D (ver Fig.
5.15). O ponto onde a reta s entra na região D é o limite inferior
α(ϑ, ϕ) para a variável r e o ponto onde a reta s sai da região D é o
limite superior β(ϑ, ϕ) para a variável r. Daí, r ∈ [α(ϑ, ϕ), β(ϑ, ϕ)].
PRÓXIMA AULA
Em nossa próxima aula veremos algumas das inúmeras aplica-

ções da integral tripla. Nossa atenção estará voltada para o cálculo
do centro de massa e momentos de inércia de sólidos gerados por
intersecções de superfícies em R3 .
ATIVIDADES
Deixamos como atividades dois problemas envolvendo mudança

de variáveis em integrais triplas.
ATIV. 5.1. Determine o volume do sólido formado pela intersec-

ção das superfícies z = 0, z = 1 + x2 + 3y 2 e x2 + y 2 = 1.
primeiro exemplo e use o sistema de coordenadas cilíndricas.
ATIV. 5.2. Seja D ⊂ R3 a região formada pela intersecção das

superfícies z = 0 e x2 + yh2 + z 2 = 1 e f : R3 7→ R dada por
103
Z Z Z
f (x, y, z) = z. Determine a integral tripla f (x, y, z)dxdydz.
D
Comentário: Volte ao texto e reveja com calma e atenção o se-
gundo exemplo e use o sistema de coordenadas esféricas.

2003.
104
AULA
Algumas Aplicações
das Integrais tríplas 6
META:
Apresentar algumas aplicações das integrais triplas de funções de
OBJETIVOS:
Determinar o volume, o centro de massa momento de massa e o
momento de inércia de alguns sólidos em R3 .
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I, superfícies em R3 , coordenadas
polares da disciplina Cálculo II, coordenadas cilíndricas, coordena-
das esféricas e integrais duplas aula 04 e aula 05.
Algumas Aplicações das Integrais tríplas
6.1 Introdução
Caros alunos, nossa sexta aula tem como objetivo introduzir

algumas aplicações da integral tripla. Em particular veremos como
calcular a massa de uma região D ⊂ R3 dada sua distribuição
de densidade, bem como calcular, para a mesma, seu centro de
gravidade e momentos de massa. É um bocado de cálculo, mais
chegaremos lá.
6.2 Preliminares
Consideraremos uma região D ⊂ R3 finita, com uma distribui-

ção de densidade (massa por unidade de volume) % : D 7→ R+ i.e.
%(x, y, z) > 0, ∀(x, y, z) ∈ D.
Determinação da massa
Para determinar a massa consideremos uma função Φ definida
em um domínio paralelepipedal R = {(x, y, z) ∈ R3 |a ≤ x ≤
b∧ c ≤ y ≤ d ∧ e ≤ z ≤ f } tal que D ⊂ R e Φ(x, y, z) =
 %(x, y, z) , (x, y, z) ∈ D
. Considerando a uma partição para
 0 , (x, y, z) ∈
/D
o retângulo R dada por P = P [R] = P [a, b] × P [c, d] × P [e, f ],
o produto cartesiano de partições P [a, b], P [c, d] e P [e, f ] onde
P [a, b] = {x0 = a, x1 , . . . , xi , xi+1 , . . . , xl = b},x0 < x1 < · · · <
xi < xi+1 < · · · < xl , P [c, d] = {y0 = c, y1 , . . . , yj , yj+1 , . . . , ym =
d}, y0 < y1 < · · · < yj < yj+1 < · · · < ym e P [e, f ] = {z0 =
e, z1 , . . . , zk , zk+1 , . . . , zn = f }, z0 < z1 < · · · < zk < zk+1 <
· · · < zn . Tomamos um ponto (ξi , ζj , ηk ) ∈ [xi−1 , xi ] × [yj−1 , yj ] ×
[zk−1 , zk ] em cada pequeno paralelepípedo e definimos a seguinte
106
Cálculo III AULA
soma de Riemann:
Xl Xm Xn
6
Slmn = Φ(ξi , ζj , ηk )∆Vijk .
i=1 j=1 k=1
A massa da região D, denotada m(D), será a integral integral
tripla da função %(x, y, z) sobre o domínio D ⊂ R3 , dada por
Z Z Z
%(x, y, z)dxdydz e definida como o seguinte limite:
D Z Z Z
def
m(D) = %(x, y, z)dxdydz = lim Slmn .
D |P |→0
OBS 6.1. Para a determinação do peso da região D toma-se a

seguinte soma de Riemann:
Xl Xm Xn
Slmn = g(ξi , ζj , ηk )Φ(ξi , ζj , ηk )∆Vijk ,
i=1 j=1 k=1
onde g(ξi , ζj , ηk ) é a aceleração da gravidade no ponto (ξi , ζj , ηk ).
E o peso da região D, denotado p(D), será dado pela integral tri-
pla:
Z Z Z
def
p(D) = g(x, y, z)%(x, y, z)dxdydz = lim Slmn .
D |P |→0
Determinação dos Momentos de Massa

de uma região D ⊂ R3 limitada com distribuição de densidade
%(x, y, z), ∀(x, y, z) ∈ D. Para calcular o momento de massa de
um pequeno paralelepípedo com relação ao plano coordenado yz
tomamos o seguinte produto ξi Φ(ξi , ζj , ηk )∆Vijk . Aqui ξi repre-
senta uma aproximação da distância do pequeno paralelepípedo
∆ξi ∆ζj ∆ηk ao plano coordenado yz. O momento total em relação
ao plano yz para a região D será aproximado pelo limite da soma
de Riemann:
107
l X
X m X
n
Slmn = ξi Φ(ξi , ζj , ηk )∆Vijk .
i=1 j=1 k=1
O momento de massa da região D em relação ao plano yz, denotado
Z Z Z
Myz (D), será dado pela integral tripla x%(x, y, z)dxdydz
D
definida pelo limite:
Z Z Z
def
Myz (D) = x%(x, y, z)dxdydz = lim Slmn .
D |P |→0
De forma semelhante chega-se ao momento de massa da região D

em relação ao plano xz tomando-se a seguinte soma de Riemann:
X l Xm X n
Slmn = ζj Φ(ξi , ζj , ηk )∆Vijk .
i=1 j=1 k=1
O momento de massa da região D em relação ao plano xz, denotado
Z Z Z
Mxz (D), será dado pela integral tripla y%(x, y, z)dxdydz
D
Z Z Z
def
Mxz (D) = y%(x, y, z)dxdydz = lim Slmn .
D |P |→0
E o momento de massa da região D em relação ao plano xy tomando-

se a seguinte soma de Riemann:
Xl X m Xn
Slmn = ηk Φ(ξi , ζj , ηk )∆Vijk .
i=1 j=1 k=1
O momento de massa da região D em relação ao plano xy, deno-
Z Z Z
tado Mxy (D), será dado pela integral tripla z%(x, y, z)dxdydz
D
Z Z Z
def
Mxy (D) = z%(x, y, z)dxdydz = lim Slmn .
D |P |→0
Determinação dos Momentos de Inércia

de uma região D ⊂ R3 limitada com distribuição de densidade
%(x, y, z), ∀(x, y, z) ∈ D. Para calcular, aproximadamente, o mo-
mento de inércia de um pequeno paralelepípedo com relação a uma
108
Cálculo III AULA
reta r, tomamos o seguinte produto d2 (ξi , ζj , ηk )Φ(ξi , ζj , ηk )∆Vijk , 6
onde d(ξi , ζj , ηk ) representa a distância do ponto (ξi , ζj , ηk ) à reta
r. Em particular a distância do ponto (ξi , ζj , ηk ) ao eixo x é
q
dada por: d(ξi , ζj , ηk ) = ζj2 + ηk2 e o momento de inércia do
pequeno paralelepípedo em relação ao eixo x será aproximado por:
(ζj2 + ηk2 )Φ(ξi , ζj , ηk )∆Vijk . O momento de inércia total em relação
ao eixo x para a região D será aproximado pelo limite da soma de
Riemann:
l X
X m X
n
Slmn = (ζj2 + ηk2 )Φ(ξi , ζj , ηk )∆Vijk .
i=1 j=1 k=1
O momento de inércia da região D em relação ao eixo x, denotado
Z Z Z
Ix é dado pela integral (y 2 + z 2 )%(x, y, z)dxdydz calculada
D
pelo limite:
Z Z Z
def
Ix (D) = (y 2 + z 2 )%(x, y, z)dxdydz = lim Slmn .
D |P |→0
De forma semelhante chega-se ao momento de inércia da região D

em relação ao eixo y tomando-se a seguinte soma de Riemann:
X l Xm X n
Slmn = (ξi2 + ηk2 )Φ(ξi , ζj , ηk )∆Vijk .
i=1 j=1 k=1
O momento de inércia da região D em relação ao eixo y, denotado
Z Z Z
Iy é dado pela integral (x2 + z 2 )%(x, y, z)dxdydz calculada
D
pelo limite:
Z Z Z
def
Iy (D) = (x2 + z 2 )%(x, y, z)dxdydz = lim Slmn .
D |P |→0
Também chega-se ao momento de inércia da região D em rela-

ção ao eixo z tomando-se a seguinte soma de Riemann: Slmn =
Xl Xm X n
(ξi2 + ζj2 )Φ(ξi , ζj , ηk )∆Vijk .
i=1 j=1 k=1
O momento de inércia da região D em relação ao eixo z é dada pela
Z Z Z
integral (x2 + y 2 )%(x, y, z)dxdydz calculada pelo limite:
D
109
Z Z Z
def
Iz (D) = (x2 + y 2 )%(x, y, z)dxdydz = lim Slmn .
D |P |→0
Determinação do Centro de Massa

O centro de massa de uma região D ⊂ R3 finita, com uma dis-
tribuição de densidade mássica %(x, y, z), ∀(x, y, z) ∈ D, é o ponto
(x̄, ȳ, z̄) definidoZpor:
Z Z
x%(x, y, z)dxdydz
Myz (D)
x̄ = = Z Z ZD ,
m(d)
%(x, y, z)dxdydz
Z Z Z D
y%(x, y, z)dxdydz
Mxz (D)
ȳ = = Z Z ZD e
m(d)
%(x, y, z)dxdydz
Z Z ZD
z%(x, y, z)dxdydz
Mxy (D)
z̄ = = Z Z ZD .
m(d)
%(x, y, z)dxdydz
D
6.3 Algumas Aplicações da Integral Tripla
Faremos duas aplicações da integral tripla. A primeira refere-se

ao cálculo do centro de massa de de um sólido gerado pela intersec-
ção de superfícies, usando o sistema de coordenadas cartesiano. A
segunda trata-se da determinação da massa e do momento de inér-
cia Iz de um sólido gerado pela intersecção de superfícies, usando
o sistema de coordenadas cilíndricas. Vamos aos nossos exemplos.
Exemplo 6.1. Considerando a intersecção das superfícies: x = 0,

x = a, y = 0, y = x2 , z = 0 e z = x2 , (Fig 6.1), determinar sua
massa e seu centro de massa levando en conta uma distribuição de
densidade constante %(x, y, z) = %.
110
Cálculo III AULA
6
SOLUÇÃO:
nando a (Fig 6.1) e verificando que 0 ≤ x ≤ a, 0 ≤ y ≤ x2 e
0 ≤ z ≤ x2 .
Em segundo calcularemos a massa da região D, m(D) e os res-
pectivos momentos de massa com relação ao planos yz, xz e xy,
respectivamente.
Passo 1 determinar a massa m(D), dada pela integral tripla:

Z Z Z Z a Z x2 Z x2
m(D) = %(x, y, z)dxdydz = %dzdydx
D 0 0 0
Integrando em z temos:
Z a Z x2 2
x
m(D) = %z dydx
0 0 0
Z a Z x2
m(D) = %(x2 − 0)dydx
0 0
Z a Z x2
m(D) = % x2 dydx
0 0
Z a x2
m(D) = % x2 y dx

0 0
111

Z a
m(D) = % x2 (x2 − 0)dx
0
Z a
m(D) = % x4 dx
0
Finalmente, integrando em x temos:
x5 a
m(D) = %
5 0
a5 05
m(D) = % −
5 5
a5
m(D) = %
5
Passo 2 determinar o momento de massa relativo ao plano yz

Myz (D), dada pela integral tripla:
Z Z Z Z a Z x2 Z x2
Myz (D) = %(x, y, z)xdxdydz = %xdzdydx
D 0 0 0
Z a Z x2 x2
Myz (D) = %xz dydx

0 0 0
Z a Z x2
Myz (D) = %x(x2 − 0)dydx
0 0
Z a Z x2
Myz (D) = % x3 dydx
0 0
Z a x2
Myz (D) = % x3 y dx

0 0
Z a
Myz (D) = % x3 (x2 − 0)dx
0
Z a
Myz (D) = % x5 dx
0
112
Cálculo III AULA
Myz (D) = %
x6 a
6 0

6
a6 06
Myz (D) = % −
6 6
a6
Myz (D) = %
6
Passo 3 determinar o momento de massa relativo ao plano xz

Mxz (D), dada pela integral tripla:
Z Z Z Z a Z x2 Z x2
Mxz (D) = %(x, y, z)ydxdydz = %ydzdydx
D 0 0 0
Z a Z x2 x2
Mxz (D) = %yz dydx

0 0 0
Z a Z x2
Mxz (D) = %y(x2 − 0)dydx
0 0
Z a Z x2
Mxz (D) = % x2 ydydx
0 0
Z a
y 2 x2
Mxz (D) = % x2 dx
0 2 0
Z a
(x2 )2 02
Mxz (D) = % x2 − dx
0 2 2
Z a 6
x
Mxz (D) = % dx
0 2
x7 a
Mxz (D) = %
14 0
a7 07
Mxz (D) = % −
14 14
a7
Mxz (D) = %
14
113
Passo 4 determinar o momento de massa relativo ao plano xy.

Como a região D tem simetria com relação às variáveis y e z, e a
distribuição de densidade também (por ser constante) temos que
Mxz (D) = Mxy (D). De qualquer forma vamos verificar:
Z Z Z Z a Z x2 Z x2
Mxy (D) = %(x, y, z)zdxdydz = %zdzdydx
D 0 0 0
Z a Z x2 2 2
z x
Mxy (D) = % dydx
0 0 2 0
Z a Z x2
(x2 )2 02
Mxz (D) = % − dydx
0 0 2 2
Z a Z x2 4
x
Mxz (D) = % dydx
0 0 2
Z a 4 2
x x
Mxz (D) = % y dx
0 2 0
Z a 4
x
Mxz (D) = % (x2 − 0)dx
0 2
Z a 6
x
Mxz (D) = % dx
0 2
x7 a
Mxz (D) = %
14 0
a7 07
Mxz (D) = % −
14 14
a7
Mxz (D) = %
14
Passo 5 determinar o centro de massa (x̄, ȳ, z̄) da região D, A

saber:
114
Cálculo III AULA
Myz (D)
a6
% 5a
6
x̄ = = 65 = ,
m(d) a 6
%
5
a7
Mxz (D) % 5a2
ȳ = = 145 = e
m(d) a 14
%
5
a7
Mxy (D) % 5a2
z̄ = = 145 = .
m(d) a 14
%
5
Vamos rapidinho ao nosso segundo exemplo.
Exemplo 6.2. Considerando a interseção das superfícies: x = 0,

x2 + y 2 = b2 , z = 0 e z = a, (Fig 6.2), determinar sua massa e
seu momento de inércia Iz (D), relativo ao eixo z, levando en conta
uma distribuição de densidade constante %(x, y, z) = %.
SOLUÇÃO:
nando a (Fig 6.2) e verificando que −b ≤ x ≤ +b, 0 ≤ y ≤
√
+ b2 − x2 e 0 ≤ z ≤ a. Observemos que para este caso é mais ade-
115
quado usar o sistema de coordenadas cilíndrico, dado pela trans-

formação (x, y, z) 7→ (r, ϑ, z) onde: x = r cos(ϑ), y = r sin(ϑ),
z = z e os limites de integração passam a: 0 ≤ r ≤ b, 0 ≤ ϑ ≤ π e
0 ≤ z ≤ a.
Em segundo, calcularemos a massa da região D, m(D) e o mo-
mento de inércia Iz (D), relativo ao eixo z, respectivamente.
Passo 1 determinar a massa m(D), dada pela integral tripla:

Z Z Z Z πZ bZ a
m(D) = %(x, y, z)dxdydz = %rdzdrdϑ
D 0 0 0
Z πZ b
a
m(D) = %z rdrdϑ
0 0 0
Z a Z x2
m(D) = %(a − 0)rdrdϑ
0 0
Z πZ b
m(D) = %a rdrdϑ
0 0
Z π 2
r b
m(D) = %a dϑ
0 2 0
Z π 2
b 02
m(D) = %a − dϑ
0 2 2
b2 π
Z
m(D) = %a dϑ
2 0
Finalmente, integrando em ϑ temos:
b2 π
m(D) = %a ϑ
2 0
b2
m(D) = %a (π − 0)
2
ab2
m(D) = %π
2
116
Cálculo III AULA
Passo 2 Levando em conta que: x2 +y 2 = (r cos(ϑ))2 +(r sin(ϑ))2 = 6
r2 , determinar o momento de inércia Iz (D), relativo ao eixo z, dada
pela integral tripla:
Z Z Z Z π Z bZ a
Iz (D) = %(x, y, z)(x2 +y 2 )dxdydz = %r2 rdzdrdϑ
D 0 0 0
Z πZ b
a
Iz (D) = %z r3 drdϑ
0 0 0
Z a Z x2
Iz (D) = %(a − 0)r3 drdϑ
0 0
Z πZ b
Iz (D) = %a r3 drdϑ
0 0
Z π 4
r b
Iz (D) = %a dϑ
0 4 0
Z π 4
b 04
Iz (D) = %a − dϑ
0 4 4
b4 π
Z
Iz (D) = %a dϑ
4 0
Finalmente, integrando em ϑ temos:
b4 π
Iz (D) = %a ϑ
4 0
b4
Iz (D) = %a (π − 0)
4
ab4
Iz (D) = %π
4
6.4 Conclusão
Na aula de hoje, vimos que dentre as inúmeras aplicações da

integral tripla, algumas das mais importantes são: dada uma região
117
D ⊂ R3 e sua distribuição de densidade volumétrica de massa

%(x, y, z), ∀(x, y, z) ∈ D, as determinação da massa m(d), dos seus
momentos de massa Myz relativo ao plano yz, Mxz relativo ao
plano xz e Mxy relativo ao plano xy, dos momentos de inércia Ix
relativo ao eixo x, Iy relativo ao eixo y e Iz relativo ao eixo z.
RESUMO
O nosso resumo de hoje constará de uma série de fórmulas

para os cálculo da massa, momento de massa, momento de inér-
cia e centro de gravidade de regiões D ∈ R3 limitadas no es-
paço dada sua distribuição de densidade volumétrica de massa
%(x, y, z), ∀(x, y, z) ∈ D. No corpo do texto temos umas pequenas
argumentações heurísticas de como chegar a tais fórmulas, usando
partições e somas de Riemman.
Dada uma região D ∈ R3 limitada com distribuição de densi-

dade volumétrica de massa %(x, y, z) podemos calcular:
A Massa da região D
Z Z Z
m(D) = %(x, y, z)dxdydz.
D
O Momento de Massa de D em Relação ao Plano yz

Z Z Z
Myz (D) = x%(x, y, z)dxdydz.
D
O Momento de Massa de D em Relação ao Plano xz

Z Z Z
Mxz (D) = y%(x, y, z)dxdydz.
D
118
Cálculo III AULA
O Momento de Massa em Relação ao Plano xy
Z Z Z
6
Mxy (D) = z%(x, y, z)dxdydz.
D
O Momento de inércia de D em Relação ao eixo x

Z Z Z
Ix (D) = (y 2 + z 2 )%(x, y, z)dxdydz.
D
O Momento de inércia de D em Relação ao eixo y

Z Z Z
Iy (D) = (x2 + z 2 )%(x, y, z)dxdydz.
D
O Momento de inércia de D em Relação ao eixo z

Z Z Z
Iz (D) = (x2 + y 2 )%(x, y, z)dxdydz.
D
Z Z Z (x̄, ȳ, z̄) de D

O Centro de Massa
x%(x, y, z)dxdydz
Myz (D)
x̄ = = Z Z ZD ,
m(d)
%(x, y, z)dxdydz
Z Z ZD
y%(x, y, z)dxdydz
Mxz (D) D
ȳ = = Z Z Z e
m(d)
%(x, y, z)dxdydz
Z Z ZD
z%(x, y, z)dxdydz
Mxy (D)
z̄ = = Z Z ZD .
m(d)
%(x, y, z)dxdydz
D
PRÓXIMA AULA
Em nossa próxima aula passaremos a estudar funções vetoriais

f : C ⊂ R3 7→ R3 onde C é uma curva no espaço R3 , dada para-
metricamente por: x = x̂(t), y = ŷ(t) e z = ẑ(t), t ∈ [a, b]. Não
estaremos, como no Cálculo II, interessados na geometria intrín-
119
seca das curvas e sim na contribuição de sua geometria no cálculo

de integrais de campos de vetores definidos sobre tais curvas.
ATIVIDADES
Deixamos como atividades os seguintes problemas:
ATIV. 6.1. Considerando a interseção das superfícies: x = 0,

x = a, y = 0, y = x2 , z = 0 e z = x2 , (Fig 6.1), determinar seu
momento de inércia Iz relativo ao eixo z, levando em conta uma
distribuição de densidade constante %(x, y, z) = %.
primeiro exemplo, ele lhe servirá de guia.
ATIV. 6.2. Considerando a intersecção das superfícies: x = 0,

x2 +y 2 = b2 , z = 0 e z = a, (Fig 6.2), determinar seu momento de
massa Myz relativo ao plano yz, levando em conta uma distribuição
de densidade constante %(x, y, z) = %.
segundo exemplo, ele lhe servirá de guia.

120
Cálculo III AULA
SWOKOWSKI, Earl E., Cálculo com Geometria Analítica, Volume 6
2003.
121
AULA
Integrais de Funções
Vetoriais sobre Curvas em R3 7
META:
Apresentar integrais de funções vetoriais definidas sobre curvas em
R3 .
OBJETIVOS:
Definir integrais de funções vetoriais definidas sobre curvas em R3
e calcular integrais de algumas funções vetoriais definidas sobre
curvas em R3 .
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I, vetores da disciplina Vetores
e Geometria analítica e curvas em R3 da disciplina Cálculo II.
Integrais de Funções Vetoriais sobre Curvas em R3
7.1 Introdução
Caros alunos nossa aula de hoje “Integrais de Funções Vetori-

ais sobre Curvas” tem um forte sabor de física pois, veremos coisas
como: calculo do trabalho de uma força (função vetorial) ao longo
de uma trajetória (curva) ou fluxo de um campo de vetores atra-
vés de uma curva (o termo fluxo é tipicamente da física). Isto,
não quer dizer que vocês tenham que se empenhar nos aspectos
físicos, devendo apenas ater-se aos aspectos matemáticos que são
os objetivos de nossa aula.
7.2 Curvas em R3
Nesta seção faremos uma pequena recapitulação sobre curvas

em R3 , que vocês já viram em Cálculo II. Será um breve resumo
onde estaremos recapitulando as definições e principais resultados.
Consideremos uma curva C ⊂ R3 dada parametricamente por:

x = x̂(t), y = ŷ(t) e z = ẑ(t), t ∈ [a, b] ou em sua forma vetorial
~r (t) = x̂(t)~i + ŷ(t)~j + ẑ(t)~kk.
O vetor tangente unitário à C é dado por:
d~~r (t) −1 d~~r (t)

~
T (t) =

dt dt
A velocidade e a aceleração de uma partícula seguindo a curva C,

com movimento dado por ~r (t), no instante t são dadas por:
124
Cálculo III AULA
7
d~~r (t) dx̂(t)~ dŷ(t)~ ẑ(t)~
~v (t) = = i+ j+ k
dt dt dt dt
d2~r (t) d2 x̂(t)~ d2 ŷ(t)~ d2 ẑ(t)~
~a (t) = = i+ j+ k
dt2 dt2 dt2 dt2
O comprimento de arco da curva C ⊂ R3 parametrizada por x =

x̂(t), y = ŷ(t) e z = ẑ(t), no intervalo [a, t] é dado por:
s 2 2 2
Z t
dx̂(t) dŷ(t) dẑ(t)
ŝ(t) = + + dt
a dt dt dt
Podemos inverter s = ŝ(t) como t = t̂(s) e descrever a curva

C ⊂ R3 parametrizada por comprimento de arco x = x̂(t̂(s)),
y = ŷ(t̂(s)) e z = ẑ(t̂(s)).
A curvatura de C é definida por:

dT~ (s)
k(s) =

ds

e pode ser calculada usando-se a fórmula:

1 dT~ (t)
k(t) =
|~~v (t)| dt

O vetor normal unitário é definido por:
dT~ (t) −1 dT~ (t)

~ (t) =
N

dt dt

O vetor binormal à curva C ⊂ R3 é definido por:
~
B
B(t) = T~ (t) × N
~ (t)
Finalmente a torção da curva C ⊂ R3 é definida por:
125
~
dB
B(s)
τ (s) = − ~
•N
ds
7.3 Massa, Momento de Massa e Momento de

Inércia de Curvas em R3
Muito embora o cálculo da massa, momento de massa e centro

de massa de uma curva C ⊂ R3 não envolva integração de funções
vetoriais, começaremos por aqui.
Seja C ⊂ R3 , uma curva contínua e lisa, parametrizada por com-

primento de arco e % : C ⊂ R3 7→ R+ , a densidade linear de massa
de C onde: ∀(x, y, z) ∈ C, %(x, y, z) > 0.
Definição 7.1. A massa de C ⊂ R3 , denotada m(C), é definida

por:
Z
m(C) = %(x, y, z)ds
C
Definição 7.2. O momento de massa de C ⊂ R3 relativo ao plano

yz, denotada Myz (C), é definido por:
Z
Myz (C) = %(x, y, z)xds
C

xz, denotada Mxz (C), é definido por:
Z
Mxz (C) = %(x, y, z)yds
C

xy, denotada Mxy (C), é definido por:
Z
Mxy (C) = %(x, y, z)zds
C
126
Cálculo III AULA
Definição 7.5. O centro de Massa de C ⊂ R3 é dado por (x̄, ȳ, z̄), 7
onde:
Z
%(x, y, z)xds
Myz (C)
x̄ = = ZC
m(C)
%(x, y, z)ds
ZC
%(x, y, z)yds
Mxz (C)
ȳ = = ZC
m(C)
%(x, y, z)ds
ZC
%(x, y, z)zds
Mxy (C)
z̄ = = ZC
m(C)
%(x, y, z)ds
C
Definição 7.6. O momento de inércia de C ⊂ R3 relativo ao eixo

x, denotada Ix (C), é definido por:
Z
Ix (C) = %(x, y, z)(y 2 + z 2 )ds
C

y, denotada Iy (C), é definido por:
Z
Iy (C) = %(x, y, z)(x2 + z 2 )ds
C

z, denotada Iz (C), é definido por:
Z
Iz (C) = %(x, y, z)(x2 + y 2 )ds
C
OBS 7.1. Se a curva C ⊂ R3 é dada parametricamente por: x =

x̂(t), y = ŷ(t) e z = ẑ(t), t ∈ [a, b], a massa, momento de massa
relativo aos planos yz, xz e xy, momento de inércia relativo aos
eixos x, y e z respectivamente pode ser calculados por:
127
Z b
m(C) = %(x̂(t), ŷ(t), ẑ(t))|~~v (t)|dt
a
Z b
Myz (C) = %(x̂(t), ŷ(t), ẑ(t))x̂(t)|~~v (t)|dt
a
Z b
Mxz (C) = %(x̂(t), ŷ(t), ẑ(t))ŷ(t)|~~v (t)|dt
a
Z b
Mxy (C) = %(x̂(t), ŷ(t), ẑ(t))ẑ(t)|~~v (t)|dt
a
Z b
Ix (C) = %(x̂(t), ŷ(t), ẑ(t))(ŷ 2 (t) + ẑ 2 (t))|~~v (t)|dt
a
Z b
Iy (C) = %(x̂(t), ŷ(t), ẑ(t))(x̂2 (t) + ẑ 2 (t))|~~v (t)|dt
a
Z b
Iz (C) = %(x̂(t), ŷ(t), ẑ(t))(x̂2 (t) + ŷ 2 (t))|~~v (t)|dt
a
7.4 Campos Vetoriais: Trabalho, Circulação e

Fluxo
Consideremos um campo de vetores F~ : D ⊂ R3 7→ R3 e uma

curva C ⊂ D contínua e suave.
Definição 7.9. Definimos o fluxo integral de escoamento do campo

vetorial F~ ao longo da curva C ⊂ R3 por:
Z
Φ(F, C) = F~ • T~ ds
C
OBS 7.2. Quando a curva é simples e fechada, o fluxo integral de

escoamento é denominado de circulação e escrevemos:
I
Φ(F, C) = F~ • T~ ds
C
OBS 7.3. Se a curva C ⊂ D ⊂ R3 é parametrizada por: x = x̂(t),

y = ŷ(t) e z = ẑ(t), t ∈ [a, b]. Podemos interpretar o campo
128
Cálculo III AULA
vetorial F~ : D ⊂ R3 7→ R3 como um campo de força no espaço, 7
a curva C ⊂ D ⊂ R3 como uma trajetória, a parametrização da
curva C ⊂ D ⊂ R3 dada por x = x̂(t), y = ŷ(t) e z = ẑ(t), t ∈ [a, b]
como o movimento de uma partícula seguindo a trajetória C e o
fluxo integral de escoamento como o trabalho executado pela força
F~ ao longo de C e dado por:
Z b
T (F, C) = F~ (x̂(t), ŷ(t), ẑ(t)) • T~ (t)dt
a
OBS 7.4. Se a curva C ⊂ R3 é representada vetorialmente por:

~r = x̂(t)~ii+ŷ(t)~jj+ẑ(t)~kk, t ∈ [a, b], e o campo vetorial F~ : D ⊂ R3 7→
R3 representado por: F~ = f1 (x, y, z)~i + f2 (x, y, z)~j + f3 (x, y, z)~kk,
onde f1 , f2 , f3 : D ⊂ R3 7→ R são funções de valores reais, o fluxo
integral de escoamento pode ser escrito como uma das três formas:
Z
T (F, C) = F~ • d~~r
ZC
T (F, C) = (f1 dx + f2 dy + f3 dz)
C
Z b
T (F, C) = f1 + f2 + f3 dt
a dt dt dt
Consideraremos, agora o caso particular de uma curva plana C ⊂

D ⊂ R2 simples e fechada e de um campo vetorial F~ : D ⊂ R2 7→
R2 . Interpretaremos o campo vetorial F~ como o campo de veloci-
dade de um fluido que atravessa a região D ⊂ R2 .
Definição 7.10. Definimos o fluxo de F através de C por:

I
def
φ(F, C) = F~ • N
~ ds
C
~ é a normal unitária exterior a C.

onde: N
129
7.5 Independência do Caminho
Consideremos um campo vetorial F~ : D ⊂ R3 7→ R3 , dois

pontos A, B ∈ D e um caminho C ⊂ D ligando o ponto A ao ponto
B. O trabalho realizado para mover uma partícula do ponto A ao
Z B
ponto B ao longo da trajetória C, dado por ~ de modo
F~ • dr
A
geral depende do caminho C que liga os dois pontos. Porém, para
alguns campos vetoriais este trabalho depende apenas dos pontos
A e B.
Definição 7.11. Sejam F~ : D ⊂ R3 7→ R3 um campo vetorial e

Z B
dois pontos A, B ∈ D. Se ~ é a mesma ∀C ⊂ D para-
F~ • dr
A
metrizada por: x = x̂(t), y = ŷ(t) e z = ẑ(t), t ∈ [a, b] tal que
A = C(a) e B = C(b) dizemos que F~ é um campo conservativo
em D.
Vamos em seguida definir um operador diferencial vetorial muito

importante denominado gradiente, A saber:
Definição 7.12. Sejam f : D ⊂ R3 7→ R uma função derivável de

valores reais. Definimos o gradiente de f , denotado ∇f , como o
campo vetorial ∇f : D ⊂ R3 7→ R3 dado por:
def ∂f~ ∂f ~ ∂f ~
∇f = i+ j+ k
∂x ∂y ∂z
Quando um campo vetorial pode ser dado pelo gradiente de um

campo escalar, dizemos que o campo escalar é uma função poten-
cial para o campo vetorial. A saber:
Definição 7.13. Sejam F~ : D ⊂ R3 7→ R3 um campo vetorial e

f : D ⊂ R3 7→ R uma função derivável de valores reais tais que
130
Cálculo III AULA
F~ = ∇f então f é dita uma função potencial para o campo vetorial 7
F~ em D
Daqui por diante consideraremos C uma curva lisa i.e. constituída

por um número finito de curvas simples unidas pelas extremidades
e D um conjunto aberto e conexo i.e. dado qualquer ponto de D
existe uma bola de centro no ponto inteiramente contida em D e
dado dois pontos quaisquer de D o segmento de reta que os une está
inteiramente contido em D. Consideraremos o campo vetorial F~ :
D ⊂ R3 7→ R3 dado por F~ = f1 (x, y, z)~i + f2 (x, y, z)~j + f3 (x, y, z)~k
onde f1 , f2 , f3 : D ⊂ R3 7→ R são funções de valores reais contínuas
e com derivadas de primeira ordem contínuas.
~ : D ⊂ R3 7→ R3 dado por F
Teorema 7.1. Sejam F ~ = f1 (x, y, z)~ii+
f2 (x, y, z)~j + f3 (x, y, z)~k onde f1 , f2 , f3 : D ⊂ R3 7→ R são funções

de valores reais contínuas e com derivadas de primeira ordem con-
tínuas em uma região D ⊂ R3 aberta e conexa do espaço. Então
existe uma função f : D ⊂ R3 7→ R contínua e diferenciável em
∂f ∂f ∂f
D ⊂ R3 tal que F~ = ∇f = ~i + ~j + ~k se somente se F~ for
∂x ∂y ∂z
um campo conservativo.
PROVA: Provaremos aqui a suficiência do teorema i.e. Se existe

uma função f : D ⊂ R3 7→ R contínua e diferenciável em D ⊂ R3
~ = ∇f = ∂f~i + ∂f ~j + ∂f ~k então F
tal que F ~ é um campo conser-
∂x ∂y ∂z
vativo.
Suponha dois pontos A, B ∈ D e uma curva C ⊂ D parametrizada
por ~r (t) = x̂(t)~i + ŷ(t)~j + ẑ(t)~kk, t ∈ [a, b] tal que A = C(a) e
B = C(b). Ao longo da curva C f é uma função f (x̂(t), ŷ(t), ẑ(t))
derivável com relação a t e levando em conta a regra da cadeia
temos:
131
df ∂f dx̂ ∂f dŷ ∂f dẑ

= + +
dt ∂x dt ∂y dt ∂z dt
Por outro lado o gradiente de f e a derivada do vetor posição ~r
com relação a t são dados por:
∂f ∂f ∂f ~
∇f = ~i + ~j + k
∂x ∂y ∂z
d~~r dx̂ dŷ dẑ
= ~i + ~j + ~k
dt dt dt dt
d~~r
Fazendo o produto escalar de ∇f por ao longo de C temos:
dt
d~~r ∂f dx̂ ∂f dŷ ∂f dẑ
∇f • = + +
Como F~ = ∇f ao longo de C temos:
d~~r ∂f dx̂ ∂f dŷ ∂f dẑ
F~ • = + +
~
O trabalho realizado por F ao longo da curva C do ponto A até o
ponto B é dado por:
Z Z
~ d~~r
F • d~~r = F~ • dt
C C dt
Aproveitando asequações acima podemosescrever:
Z Z
~ ∂f dx̂ ∂f dŷ ∂f dẑ
F • d~~r = + + dt
∂x dt ∂y dt ∂z dt
ZC ZCb
df
F~ • d~~r = dt
ZC a dt
F~ • d~~r = f (x̂(b), ŷ(b), ẑ(b)) − f (x̂(a), ŷ(a), ẑ(a))
C Z
Portanto F~ • d~~r é independente do caminho C pois, depende
C
apenas dos valores de f nos pontos A e B, provando assim que F~
é um campo conservativo. Caros alunos deixamos como desafio a
prova da necessidade. Novamente vocês podem recorrer aos livros
de Cálculo Avançado.
Temos outro teorema que caracteriza campos vetoriais conservati-

vos. A saber:
Teorema 7.2. Seja F~ : D ⊂ R3 7→ R3 um campo vetorial dado

por: F~ = f1~i + f2~j + f3~kk, cujas funções componentes f1 , f2 , f3 :
132
Cálculo III AULA
D ⊂ R3 7→ R tem derivadas parciais de primeira ordem contínuas.
∂f1 ∂f2 ∂f1 ∂f3
7
Então F~ é conservativo se, somente se = , = e
∂y ∂x ∂z ∂x
∂f2 ∂f3
=
∂z ∂y
PROVA: Novamente vamos provar a suficiência. Se F~ é conser-

vativo, existe f : D ⊂ R3 7→ R tal que:
∂f ∂f ∂f ~
F~ = ~i + ~j + kk.
∂x ∂y ∂z
∂f ∂f ∂f
Em outras palavras: f1 = , f2 = e f3 = .
∂x ∂y ∂z
Daí temos:
∂f1 ∂2f ∂f2 ∂2f
= e = .
∂y ∂y∂x ∂x ∂x∂y
Da continuidade das derivadas parciais de f1 e f2 temos:
∂f1 ∂f2
= . De forma semelhante temos:
∂y ∂x
∂f1 ∂2f ∂f3 ∂2f
= e = .
∂z ∂z∂x ∂x ∂x∂z
∂f1 ∂f3
= .
∂z ∂x
E finalmente:
∂f2 ∂2f ∂f3 ∂2f
= e = .
∂z ∂z∂y ∂y ∂y∂z
∂f2 ∂f3
= .
∂z ∂y
Deixamos a demonstração da necessidade para vocês. Novamente
consultem livros de Cálculo Avançado.
Na proxima seção, veremos como determinar o campo potencial
quando ele existe, utilizando um exemplo.
7.6 Algumas Aplicações das Integrais de Linha
Veremos agora três aplicações das integrais de linha de campos

vetoriais sobre curvas no espaço.
133
Exemplo 7.1. Calcular o trabalho realizado pelo campo de força

F : R3 7→ R3 dado por F~ (t) = z~i + 0~j + xy~k ao longo da hélice
C ⊂ R3 dada por ~r = a cos(t)~i + a sin(t)~j + bt~kk, t ∈ [0, 4π] (ver

Fig. 7.1 ).
SOLUÇÃO: Derivando o vetor posição ~r (t) com relação a t te-
mos:
d~~r (t)
= −a sin(t)~i + a cos(t)~j + b~k
dt
O campo de força F~ ao longo da curva C ⊂ R3 é dado por:
F~ (t) = bt~i + 0~j + a2 sin(t) cos(t)~k
d~~r (t)
Fazendo o produto escalar de F~ (t) por temos:
dt
d~~r (t)
F~ (t) • = −bat sin(t) + ba2 sin(t) cos(t)
dt
Calculando o trabalho realizado pela força F~ (t) ao longo da curva
134
Cálculo III AULA
C temos:
Z
7
T (F~ , C) = F~ (t) • d~~r (t)
C
Z
d~~r (t)
= F~ (t) • dt
C dt
Z 4π
= (−bat sin(t) + ba2 sin(t) cos(t))dt
0
ba2 4π
= ba(− sin(t) + t cos(t)) + sin2 (t)

2 0
ba2
= ba(− sin(4π) + 4π cos(4π)) + sin2 (4π) −
2
ba2
−(ba(− sin(0) + 0 cos(0)) + sin2 (0))
2
= 4πba
E agora sem demora o segundo exemplo.
Exemplo 7.2. Calcular o trabalho realizado pelo campo de força

constante F : R3 7→ R3 dado por F~ = K~i + Ky~j + K~k ao longo da
curva C ⊂ R3 da intersecção da esfera (x−a)2 +(y−a)2 +(z−a)2 =

a2 e do plano x − z = 0 (ver Fig. 7.2 ).
SOLUÇÃO: Primeira coisa a fazer é obter uma parametriza-
 (x − a)2 + (y − a)2 + (z − a)2 = a2
ção para a curva . Como a
 x−z =0
135
curva C ⊂ R3 pertence a reta podemos eliminar z = y na equação

da esfera e temos:
2(x−a)2 +(y −a)2 = a2 , podemos propor como parametrização
√ sa-
2
tisfazendo a equação acima: y = a + a sin(t) e x = a + a cos(t).
2
Como z = √ x temos:
2
z =a+ a cos(t).
2
Resumindo temos a seguinte parametrização para a intersecção da
esfera
 como √
plano dados:
2

 x=a+ a cos(t)
2



y = a + a sin(t) ∀t ∈ [−π, +π] .
 √
 z = a + 2 a cos(t)



2
Podemos √escrever o vetor posição ~r como: √

2 2
~r = (a + a cos(t))~i + (a + a sin(t))~j + (a + a cos(t))~kk. Deri-
2 2
vando o vetor posição ~r (t) com relação a t temos:
√ √
d~~r (t) 2 2
=− ~ ~
a sin(t)i + a cos(t)j − a sin(t)~k
dt 2 2
Ao longo da curva C ⊂ R3 o campo de força é dado por:

F~ (t) = K~i + Ka(1 + sin(t))~j + K~kk.
d~~r (t)
Fazendo o produto escalar de F~ (t) por temos:
dt
√
d~~r (t) 2
F~ (t) • = − Ka sin(t) + Ka2 (1 + sin(t)) cos(t) −
dt √2
2
− Ka sin(t)
2
√
= − 2Ka sin(t) + Ka2 (1 + sin(t)) cos(t)
Calculando o trabalho realizado pela força F~ (t) ao longo da curva
136
Cálculo III AULA
C temos:
Z
7
T (F~ , C) = F~ (t) • d~~r (t)
C
Z
d~~r (t)
= F~ (t) • dt
C dt
Z +π √
= (− 2Ka sin(t) + Ka2 (1 + sin(t)) cos(t))dt
−π
√
√ 2 +π
= ( 2Ka cos(t) − Ka2 (1 + sin(t))2 )

2 −π
= 0
Vejamos mais um exemplo. Desta vez veremos como determinar a

função potencial para um campo conservativo.
Exemplo 7.3. Seja F~ : R3 7→ R3 um campo vetorial conservativo

dado por: F~ = yz~i + (xz + 1)~j + xy~kk. Determine sua função po-
tencial.
SOLUÇÃO: Primeiramente testaremos se o campo vetorial F~ é
um campo conservativo, usando a condição necessária e suficiente
dada por:
∂f1 ∂f2 ∂f1 ∂f3 ∂f2 ∂f3
= , = e = .
∂y ∂x ∂z ∂x ∂z ∂y
Como para o F~ dado f1 = yz, f2 = xz + 1 e f3 = xy temos:
∂f1 ∂f2 ∂f1 ∂f3 ∂f2 ∂f3
=z= , =y= e =x= .
∂y ∂x ∂z ∂x ∂z ∂y
A condição está satisfeita e F~ é um campo vetorial conservativo e
podemos procurar o f : R3 7→ R tal que:
∂f ∂f ∂f ~
F~ = ∇f = ~i + ~j + kk.
∂x ∂y ∂z
De onde tiramos:
∂f


 = f1 = yz

 ∂x
 ∂f
= f2 = xz + 1
 ∂y
 ∂f


 = f3 = xy
∂z
137
∂f
Integrando a primeira equação = yz com relação a x temos:
∂x
∂f
f = xyz + g(y, z) pois daí tiramos = yz.
∂x
Temos agora que determinar o g(y, z) de modo que a segunda equa-
ção sejam satisfeita.
Derivando f = xyz + g(y, z) com relação a y temos:
∂f ∂g
= xz + .
∂y ∂y
∂f
Comparando com a segunda equação = xz + 1 temos:
∂y
∂g
= 1.
∂y
Integrando com relação a y temos:
∂g
g(y, z) = y + h(z) pois daí tiramos = 1.
∂y
Podemos reescrever f comos:
f = xyz + y + h(z).
Temos agora que determinar o h(z) de modo que a terceira equa-
ção sejam satisfeita.
Derivando f = xyz + y + h(z) com relação a z temos:
∂f dh
= xy + .
∂y dz
∂f
Comparando com a terceira equação = xy temos:
∂z
dh
= 0.
dz
Logo h(z) = K é uma constante que podemos sem perda de gene-
ralidade fazer igual a zero e f passa a ter a forma final:
f (x, y, z) = xyz + y
7.7 Conclusão
Na aula de hoje, vimos como integrar campos vetoriais (funções

vetoriais) ao longo de curvas no espaço e no plano. Que, essenci-
almente, os conceitos por trás da integração de campos vetoriais
138
Cálculo III AULA
como circulação e fluxo sobre curvas estão intimamente ligados à 7
Física.
RESUMO
Seja C ⊂ R3 , uma curva contínua e lisa, parametrizada por

comprimento de arco e % : C ⊂ R3 7→ R+ , a densidade linear de
massa de C onde: ∀(x, y, z) ∈ C, %(x, y, z) > 0.
Massa
A massa de C ⊂ R3 , denotada m(C), é definida por:
Z
m(C) = %(x, y, z)ds
C
Momento de Massa relativo aos planos yz

yz, xz e xy
xy.
O momento de massa de C ⊂ R3 relativo ao plano yz, xz e xy de-
notados Myz (C), Mxz (C) e Mxy (C) são definidos respectivamente
por:
Z
Myz (C) = %(x, y, z)xds
ZC
Mxz (C) = %(x, y, z)yds
ZC
Mxy (C) = %(x, y, z)zds
C
Se a curva C ⊂ R3 é dada parametricamente por: x = x̂(t), y =

ŷ(t) e z = ẑ(t), t ∈ [a, b], a massa, momento de massa relativo aos
planos yz, xz e xy, respectivamente pode ser calculados por:
139
Z b
m(C) = %(x̂(t), ŷ(t), ẑ(t))|~~v (t)|dt
a
Z b
Myz (C) = %(x̂(t), ŷ(t), ẑ(t))x̂(t)|~~v (t)|dt
a
Z b
Mxz (C) = %(x̂(t), ŷ(t), ẑ(t))ŷ(t)|~~v (t)|dt
a
Z b
Mxy (C) = %(x̂(t), ŷ(t), ẑ(t))ẑ(t)|~~v (t)|dt
a
Centro de Massa.
O centro de Massa de C ⊂ R3 é dado por (x̄, ȳ, z̄), onde:
Z
%(x, y, z)xds
Myz (C) C
x̄ = = Z
m(C)
%(x, y, z)ds
Z C
%(x, y, z)yds
Mxz (C)
ȳ = = ZC
m(C)
%(x, y, z)ds
Z C
%(x, y, z)zds
Mxy (C)
z̄ = = ZC
m(C)
%(x, y, z)ds
C
Momento de Inércia relativo aos eixos x , y e z .

O momento de inércia de C ⊂ R3 relativo ao eixos x, y e z deno-
tados Ix (C), Iy (C) e Iz (C) são definidos respectivamente por:
Z
Ix (C) = %(x, y, z)(y 2 + z 2 )ds
ZC
Iy (C) = %(x, y, z)(x2 + z 2 )ds
ZC
Iz (C) = %(x, y, z)(x2 + y 2 )ds
C
140
Cálculo III AULA
Se a curva C ⊂ R3 é dada parametricamente por: x = x̂(t), y = 7
ŷ(t) e z = ẑ(t), t ∈ [a, b], o momento de inércia relativo aos eixos
x, y e z, respectivamente pode ser calculados por:
Z b
Ix (C) = %(x̂(t), ŷ(t), ẑ(t))(ŷ 2 (t) + ẑ 2 (t))|~~v (t)|dt
a
Z b
Iy (C) = %(x̂(t), ŷ(t), ẑ(t))(x̂2 (t) + ẑ 2 (t))|~~v (t)|dt
a
Z b
Iz (C) = %(x̂(t), ŷ(t), ẑ(t))(x̂2 (t) + ŷ 2 (t))|~~v (t)|dt
a
Fluxo Integral de Escoamento.

Seja um campo de vetores F~ : D ⊂ R3 7→ R3 e uma curva C ⊂ D
contínua e suave. Definimos o fluxo integral de escoamento do
~ ao longo da curva C ⊂ R3 por:
campo vetorial F
Z
Φ(F, C) = F~ • T~ ds
C
Alternativamente se a curva C ⊂ R3 é representada vetorialmente

por: ~r = x̂(t)~i + ŷ(t)~j + ẑ(t)~kk, t ∈ [a, b] e o campo vetorial F~ :
D ⊂ R3 7→ R3 por: F~ = f1 (x, y, z)~i + f2 (x, y, z)~j + f3 (x, y, z)~kk,
com f1 , f2 , f3 : D ⊂ R3 7→ R são funções de valores reais, o fluxo
integral de escoamento pode ser escrito como:
Z
T (F, C) = F~ • d~~r
ZC
T (F, C) = (f1 dx + f2 dy + f3 dz)

C
Z b
T (F, C) = f1 + f2 + f3 dt
a dt dt dt
Campo Conservativo.
Sejam F~ : D ⊂ R3 7→ R3 um campo vetorial, dois pontos A, B ∈ D.
141
Z B
Se ~((F ) • dr
~ é constante ∀C ⊂ D parametrizada por: x = x̂(t),
A
y = ŷ(t) e z = ẑ(t), t ∈ [a, b] tal que A = C(a) e B = C(b) dizemos
que F~ é um campo conservativo em D.
Gradiente.
Sejam f : D ⊂ R3 7→ R uma função derivável de valores reais.
Definimos o gradiente de f , denotado ∇f , como o campo vetorial
∇f : D ⊂ R3 7→ R3 dado por:
def ∂f~ ∂f ~ ∂f ~
∇f = i+ j+ k
∂x ∂y ∂z
Função Potencial.
Sejam F~ : D ⊂ R3 7→ R3 um campo vetorial e f : D ⊂ R3 7→ R
uma função derivável de valores reais tais que F~ = ∇f então, f é
dita uma função potencial para o campo vetorial F~ em D
PRÓXIMA AULA
Em nossa próxima aula veremos, essencialmente, integração

de funções reais e campos vetoriais (funções vetoriais) sobre su-
perfícies S ⊂ R3 . Veremos também como calcular área, massa,
momento de massa e centro de massa de superfícies.
ATIVIDADES
Deixamos como atividades os seguintes problemas envolvendo

integração de campos vetoriais sobre curvas no espaço.
142
Cálculo III AULA
ATIV. 7.1. Seja F~ : R3 7→ R3 o campo vetorial dado por: F~ (x, y, z) = 7
y(2xyz 2 + exy )~i + x(2xyz 2 + exy )~j + 2x2 y 2 z~kk:
• Mostre que campo vetorial F~ é conservativo.
• Determine uma função potencial f : R3 7→ R tal que F~ =

∇f .

problemas resolvidos acima, eles lhe servirão de guia.
ATIV. 7.2. Sejam F~ : R3 7→ R3 o campo vetorial dado por:

F~ (x, y, z) = y~i + z~j + b~kk, b 6= 0 e C ⊂ R3 a curva no espaço dada
por ~r (t) = a cos(t)~i + a sin(t)~j + c~kk, ∀t ∈ [0, 2π], a, c > 0. Deter-
mine o trabalho realizado por F~ ao longo da curva C.

problemas resolvidos acima, eles lhe servirão de guia.

143

2003.
144
AULA
Integrais de Superfícies
8
META:
Apresentar integrais de funções definidas sobre superfícies em R3 .
OBJETIVOS:
Definir integrais de funções definidas sobre superfícies em R3 e
calcular algumas integrais de funções vetoriais definidas sobre su-
perfícies em R3 .
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I, vetores da disciplina Vetores
e Geometria analítica e superfícies em R3 .
8.1 Introdução
Caros alunos nossa aula de hoje “Integrais de Superfícies” tem,

como a nossa aula anterior “Integrais de Funções Vetoriais sobre
Curvas em R3 ”, um sabor de física. Desde a determinação da
massa, momento de massa e centro de massa de uma superfície
até a determinação do fluxo de um campo vetorial através de uma
superfície. Da mesma forma que na aula anterior, vocês devem
ater-se apenas aos aspectos matemáticos da matéria abordada.
8.2 Superfícies em R3
Bom, vamos começar, bem do começo, com algumas formas de

representação de superfícies. A primeira forma de representação
de uma superfície é considerar uma função f : D ⊂ R3 7→ R e
tomar um ponto c ∈ Img(f ) da imagem de f . Desta forma, de
modo geral, f (x, y, z) = c representa uma superfície S ⊂ R3 .
Exemplo 8.1. Sejam a, b, c > 0 e f : R3 7→ R dada por: f (x, y, z) =

x2 y 2 z 2
+ + . Desta forma f (x, y, z) = d representa elipsóides para
a2 b2 c2
valores positivos de d.
Outra forma de representação de uma superfície é através de uma

parametrização. Representar S ⊂ R3 por: x = x̂(u, v), y = ŷ(u, v)
e z = ẑ(u, v), ∀(u, v) ∈ [a, b] × [c, d].
Exemplo 8.2. Tomando o exemplo anterior podemos parametri-

√ √
zar os elipsóides por: x = a d cos(u) cos(v), y = b d sin(u) cos(v)
√
e z = c d sin(v), ∀(u, v) ∈ [−π, +π] × [−π, +π].
146
Cálculo III AULA
8.3 Área de Superfícies em R3 8
Vamos usar nesta seção uma argumentação heurística objeti-
vando encontrar uma fórmula para determinação da área de uma
superfície S ⊂ R3 . A argumentação baseia-se na possibilidade
(ver, Fig. 8.1 ) de determinar a área de uma superfície através
de uma integração dupla sobre sua projeção (sombra) no plano
D ⊂ xy. Suponhamos que a projeção da superfície S ⊂ R3 dada
Figura 8.1: Superfície S ⊂ R3 e sua projeção D
por: z = f (x, y) sobre o plano xy seja a região D ⊂ xy e seja

D ⊂ R ⊂ xy um retângulo do plano xy paralelo aos eixos coor-
denados e que contenha a região D (ver, Fig. 8.2 ). Podemos
subdividir R em pequenos retângulos (através de partições como
vimos em nossa primeira aula) ∆ij de área ∆xi ∆yj . Podemos
aproximar (ai está a argumentação heurística) a pequena área da
superfície S, denotada ∆σij cuja projeção é o pequeno retângulo
∆ij pela parte do plano tangente a S no ponto (xi , yj , f (xi , yj )),
denotada ∆Pij , que tem a forma de um paralelogramo, (ver, Fig.
8.2 ) cuja projeção no plano xy é também o pequeno retângulo
∆ij . A área de ∆Pij é, de modo geral, maior que a área de ∆ij
147
Figura 8.2: Detalhe do elemento de área ∆σij
(a área da sombra é sempre menor ou igual à área do objeto). Da

geometria vetorial |~u
~ i × ~v j • p~ | é a área da projeção do paralelo-
gramo ∆Pij onde p~ é a normal a ∆ij (no caso para projeções no
plano xy p~ = ~k mas, deixaremos p~ nas fórmulas caso seja escolhido
outro plano de projeção).
|~u
~ i × ~v j • p~ | = ∆ij
Também da geometria vetorial temos:
|~u
~ i × ~v j • p~ | = |~u
~ i × ~v j |.|~
p~|.| cos(ϕij )| = ∆ij
onde ϕij é o ângulo formado pelo vetor normal p~ (xi , yj ) e o vetor

~u i × ~v j .
Como, da geometria vetorial, (ver em livros de Cálculo Avançado)
|~u
~ i × ~v j | = ∆Pij e |~
p~| = 1, temos:
∆ij
∆Pij =
| cos(ϕ)|
Como cada pedaço ∆Pij aproxima o pedaço da superfície ∆σij
então a soma:
n−1
X m−1 n−1
X m−1
X X ∆ij
∆Pij =
| cos(ϕij )|
i=0 j=0 i=0 j=0
148
Cálculo III AULA
aproxima a área de S. Um refinamento da partição de D ⊂ xy me- 8
lhora a aproximação e podemos então (argumentação heurística)
escrever:
Z Z
1
Are(S) = dxdy
D | cos(ϕ)|
Para uma superfície dada por f (x, y, z) = c, temos |∇f • p~ | =
|∇f |.|~
p~|.| cos(ϕ)| e como |~
p~| = 1 portanto:
|∇f |
Z Z Z Z
Are(S) = dσ = dA
S D |∇f • p~ |
Por outro lado podemos estender a argumentação e determinar a

integral de uma função g : D ⊂ R3 7→ R definida sobre a superfície
S ⊂ R3 na forma:
|∇f |
Z Z Z Z
g(x, y, z)dσ = g(x, y, z)) dA
S D |∇f • p~ |
Vamos a um exemplo para ilustrar os conceitos acima expostos.
Exemplo 8.3. Considere a superfície S ⊂ R3 do espaço dada por

z = a + x2 + y 2 , cuja projeção no plano xy é a região D ⊂ xy dada
por x2 + y 2 ≤ b2 e determine sua área (ver Fig. 8.3).
SOLUÇÃO: Deixamos como a primeira atividade mostrar que:
Figura 8.3: Parabolóide z = a + x2 + y 2
149
se a superfície é dada por z = f (x, y) projetada no plano xy em

D ⊂ xy sua área é dada por:
s
∂f 2
Z Z 2
∂f
Are(S) = + + 1dxdy
D ∂x ∂y
Como z = f (x, y) = a + x2 + y 2 temos:

∂f
= 2x
∂x
∂f
= 2y
∂y
Daí, substituindo na expressão da área temos:
s
∂f 2
Z Z 2
∂f
Are(S) = + + 1dxdy
D ∂x ∂y
Z Z p
= (2x)2 + (2y)2 + 1dxdy
x2 +y 2 ≤b2
Usando so sistema de coordenadas polares, x = r cos(ϑ) e r sin(ϑ),

para o cálculo da integraldupla temospara a projeção x2 +y 2 ≤ b2
 b  2π
os seguintes limites r = eϑ= e podemos reescrever
 0  0
a integral dupla como:
ZZ bp2π
Are(S) = (2r cos(ϑ))2 + (2r sin(ϑ))2 + 1rdrdϑ
0 0
Z 2π Z b q
= 4r2 cos2 (ϑ) + 4r2 sin2 (ϑ) + 1rdrdϑ
0 0
Z 2π Z b p
= 4r2 + 1rdrdϑ
0 0
dz
Fazendo a mudança de variáveis z = 4r2 + 1 temos: = 8r e os
  dr
 b  4b2 + 1
limites de integração r = passa a z = e podemos
 0  1
150
Cálculo III AULA
reescrever a integral dupla como: 8
Z 2π Z 4b2 +1 √
Are(D) = zdzdϑ
0 1
Integrando primeiro em z depois em ϑ temos:
Z 2π Z 4b2 +1 √
Are(D) = zdzdϑ
0 1
Z 2π √ 4b2 +1
= z3 dϑ
0 1
Z 2π
2 p
= ( (4b2 + 1)3 − 1)dϑ
3 0
2 p 2 2π
= ( (4b + 1)3 − 1)ϑ

3 0
4π p 2
= ( (4b + 1)3 − 1)
3
8.4 Momento de massa e Momento de Inércia

de Superfícies de Casca Fina em R3
Seja uma superfície S ⊂ R3 de casca fina dada por f (x, y, z) = c

e com densidade superficial % : S ⊂ R3 7→ R, a massa, o momento
de massa em relação aos planos yz, xz e xy são dados, respectiva-
mente, por:
|∇f |
Z Z Z Z
m(S) = %(x, y, z)dσ = %(x, y, z)
dA
|∇f • p~ |
Z ZS Z D
|∇f |
Z
Myz (S) = %(x, y, z)xdσ = %(x, y, z)x dA
S D |∇f • p~ |
|∇f |
Z Z Z Z
Mxz (S) = %(x, y, z)ydσ = %(x, y, z)y dA
S D |∇f • p~ |
|∇f |
Z Z Z Z
Mxy (S) = %(x, y, z)zdσ = %(x, y, z)z dA
S D |∇f • p~ |
151
O centro de massa, denotado (x̄, ȳ, z̄), é dado por:
|∇f |
Z Z
%(x, y, z)x dA
Myz (S) D |∇f • p~ |
x̄ = = Z Z
m(S) |∇f |
%(x, y, z) dA
|∇f • p~|
Z ZD
|∇f |
%(x, y, z)y dA
Mxz (S) D |∇f • p~ |
ȳ = = Z Z
m(S) |∇f |
%(x, y, z) dA
D |∇f • p~ |
|∇f |
Z Z
%(x, y, z)z dA
Mxy (S) |∇f • p~ |
z̄ = = Z ZD
m(S) |∇f |
%(x, y, z) dA
D |∇f • p~ |
Os momentos de inércia com relação aos eixos x, y e z são dados,
respectivamente, por:
|∇f |
Z Z
Ix (S) = %(x, y, z)(y 2 + z 2 ) dA
D |∇f • p~ |
|∇f |
Z Z
Iy (S) = %(x, y, z)(x2 + z 2 ) dA
D |∇f • p~ |
|∇f |
Z Z
Iz (S) = %(x, y, z)(x2 + y 2 ) dA
D |∇f • p~ |
Vamos ilustrar com um exemplo.
Exemplo 8.4. Considere a casca fina descrita pela superfície S ⊂

R3 dada pela parte do cone x2 + y 2 − z 2 = 0 que situa-se acima do
plano z = 0 e e abaixo do plano z = a, cuja densidade é constante
e igual a % e determine seu centro de massa (ver Fig. 8.4).
SOLUÇÃO: Em primeiro lugar determinaremos a massa da casca
fina, levando em conta que a projeção de S ⊂ R3 no plano xy é a
região circular D ⊂ xy dada por: x2 + y 2 ≤ a2 .
Para o caso f (x, y, z) = x2 +y 2 −z 2 = 0 e p~ = ~kk. daí, seu gradiente
será:
152
Cálculo III AULA
8
Figura 8.4: Cone x2 + y 2 − z 2 = 0
∇f = −2x~i − 2y~j − 2z~k
E temos:
|∇f • p~ | = 2z
p p
|∇f | = x2 + y 2 + z 2 = 2x2 + 2y 2
A massa da casca fina será:
|∇f |
Z Z
m(S) = % dA
D |∇f • p~ |
p
2x2 + 2y 2
Z Z
= % dxdy
x2 +y 2 ≤a2 2z
√ Z Z √
2 z2
= % dxdy
2 x2 +y 2 ≤a2 z
√ Z Z
2
= % dxdy
2 x2 +y 2 ≤a2

para o cálculo da integraldupla temospara a projeção x2 +y 2 ≤ a2
 a  2π
 0  0
153
√ Z 2π Z a
2
m66(S) = % rdrdϑ
2 0
√ Z 2π 02
2 r a
= % dϑ
2 0 2 0
√ Z
a2 2 2π
= % dϑ
4 0
√
2 2π
= % ϑ
4 √0
a2 2
= π%
2
Para determinar o centro de massa temos que determinar apenas

Mxy pois, pela simetria da superfície e como % é constante temos
que x̄ = ȳ = 0.
O momento de massa Mxy da casca fina será:
|∇f |
Z Z
Mxy (S) = %z dA
D |∇f • p~ |
p
2x2 + 2y 2
Z Z
= %z dxdy
x2 +y 2 ≤a2 2z
√ Z Z
2
= % zdxdy
2 x2 +y 2 ≤a2
√ Z Z
2 p
= % x2 + y 2 dxdy
2 x2 +y 2 ≤a2

para o cálculo da integraldupla temospara a projeção x2 +y 2 ≤ a2
 a  2π
 0  0
154
Cálculo III AULA
√ Z 2π Z a
8
2 p
Mxy (S) = % (r cos(ϑ))2 + (r sin(ϑ))2 rdrdϑ
2 0
√ Z 2π Z0 a q
2
= % r2 cos2 (ϑ) + r2 sin2 (ϑ)rdrdϑ
2 0 0
√ Z 2π Z a
2 √
= % r2 rdrdϑ
2 0 0
√ Z 2π Z a
2
= % r2 drdϑ
2 0
√ Z 2π 03
2 r a
= % dϑ
2 0 3 0
√ Z
a2 2 2π
= % dϑ
6 0
√
2 2π
= % ϑ
6 √0
a3 2
= π%
3
O valor de z̄ será dado por:
√
a3 2
Mxy (S) π%
z̄ = = 3√ = 2a
m(S) 2
a 2 3
π%
2
8.5 Superfícies Parametrizadas
Nesta seção veremos como calcular integrais de superfícies para

superfícies parametrizadas.
Seja uma superfície lisa S ⊂ R3 parametrizada por: x = x̂(u, v),
y = ŷ(u, v) e z = ẑ(u, v), ∀(u, v) ∈ [a, b] × [c, d] onde x̂, ŷ e ẑ
possuem derivadas contínuas com relação a u e a v. Podemos
representar a superfície pelo vetor posição ~r (u, v) = x̂(u, v)~i +
ŷ(u, v)~j + ẑ(u, v)~kk. Representaremos as derivadas do vetor r com
relação a u e a v respectivamente por ~r u , ~r v . Consideraremos em
R = [a, b] × [c, d] as quatro retas u = u0 , u = u0 + ∆u, v = v0 e
155
v = v0 + ∆v u = u0 + ∆u e denotamos ∆uv o pequeno retângulo

formado pela intersecção das quatro retas (ver Fig. 8.5). O pe-
Figura 8.5: Domínio da parametrização
queno retângulo ∆uv é mapeado pela parametrização no pequeno

elemento de área ∆σuv sobre a superfície S. O paralelogramo for-
mado pelos vetores ∆u~~r u e ∆v~~r v aproximam (por falta) o elemento
de área ∆σuv (ver Fig. 8.6). A área do paralelogramo é dada por:
Figura 8.6: Elemento de área ∆σuv em S.
|∆u~~r u × ∆v~~r v | = |~~r u × ~r v |∆u∆v.
A suposição de que S é uma superfície lisa garante que o produto

vetorial ~r u × ~r v não é o vetor nulo e portanto a área do pequeno
156
Cálculo III AULA
paralelogramo também não é nula. Podemos então fazer um par- 8
tição da região R do plano uv e mapeando-a pela parametrização
sobre a superfície S. Aproximando cada ∆σuv pela área do para-
lelogramo associado podemos aproximar a área de S pela soma de
Riemann:
XX
|~~r u × ~r v |∆u∆v.
u v
Fazendo ∆u e ∆v tenderem a zero independentemente, a conti-
nuidade das derivadas ~r v do vetor posição garante que a soma de
Riemann aproxime-se da integral dupla que dá a área Are(S) da
superfície S i.e.
Z bZ d
Are(S) = |~~r u × ~r v |dudv.
a c
Esta argumentação heurística nos permite estender os conceitos

acima desenvolvidos para definir a integral de uma função f : S ⊂
R3 7→ R definida sobre a superfície S da seguinte forma:
Definição 8.1. Sejam S ⊂ R3 uma superfície lisa definida para-

metricamente por ~r (u, v) = x̂(u, v)~i + ŷ(u, v)~j + ẑ(u, v)~kk, ∀(u, v) ∈
[a, b] × [c, d] e f : S ⊂ R3 7→ R uma função de valores reais definida
sobre S então, a integral de f sobre S será:
Z Z Z bZ d
def
f (x, y, z)dσ = f (x̂(u, v), ŷ(u, v), ẑ(u, v))|~~r u ×~~r v |dudv.
S a c
Um conceito, vindo da Física, muito importante é o do fluxo de

um campo vetorial através de uma superfície no espaço. Como
exemplo temos o fluxo de massa (massa por unidade de tempo por
unidade de área) de um fluido que é calculado através do seu campo
de velocidade e da sua densidade de massa. Vamos à definição:
Definição 8.2. Sejam S ⊂ R3 uma superfície lisa no espaço e

F~ : S ⊂ R3 7→ R3 uma função de valores vetoriais. Definimos o
157
fluxo de F~ através de S, denotado φ(F~ ), por:

Z Z
~ def
φ(F ) = F~ (x, y, z) • ~n dσ.
S
Onde n é a normal unitária em S.
OBS 8.1. Alternativamente, se S ⊂ R3 é lisa e definida parame-

tricamente por ~r (u, v) = x̂(u, v)~i + ŷ(u, v)~j + ẑ(u, v)~kk, ∀(u, v) ∈
[a, b] × [c, d] e F~ : S ⊂ R3 7→ R3 uma função de valores vetoriais.
O fluxo de F~ através de S, é dado por:
Z bZ d
~
φ(F ) = F~ (x̂(u, v), ŷ(u, v), ẑ(u, v)) • ~n |~~r u × ~r v |dudv.
a c
1
Como podemos calcular o vetor normal por ~n = ·(~~r u ×~~r v )
|~ u × ~r v |
~
r
a integral para o fluxo do campo vetorial F~ através da superfície
S pode ser reescrita como:
Z bZ d
φ(F~ ) = F~ (x̂(u, v), ŷ(u, v), ẑ(u, v)) • (~~r u × ~r v )dudv.
a c
Vejamos um exemplo envolvendo a determinação do fluxo de um

campo vetorial através de uma superfície no espaço.
Exemplo 8.5. Determine o fluxo do campo vetorial F~ : R3 7→ R3

dado por F~ (x, y, z) = z~i + z~j + xy~k através da superfície do para-
bolóide z = a2 − x2 − y 2 , que fica acima do plano z = 0 (ver Fig.
8.7).
SOLUÇÃO: Começaremos por parametrizar a superfície do pa-
rabolóide fazendo x̂ = v cos(u), ŷ = v sin(u) e z = a2 − v 2 . Desta
forma o vetor posição para a superfície fica expresso por:
~r (u, v) = v cos(u)~i + v sin(u)~j + (a2 − v 2 )~kk.
158
Cálculo III AULA
8
Figura 8.7: Parabolóide z = a2 − x2 − y 2 .
As derivadas parciais do vetor posição ~r com relação a u e a v são:

 ~r u = −v sin(u)~i + v cos(u)~j + 0k
 ~r = cos(u)~i + sin(u)~j − 2v~k
v
Podemos calcular ~r u × ~r v operacionalmente por:

~i ~j ~k

op
~r u × ~r v = det −v sin(u) v cos(u)

0

cos(u) sin(u) −2v
Calculando o determinante temos:
~r u × ~r v = −2v 2 sin(u)~i + 2v 2 sin(u)~j − v~k
O campo vetorial F~ sobre a superfície pode ser escrito como:
F~ (u, v) = (a2 − v 2 )~i + (a2 − v 2 )~j + v 2 sin(u) cos(u)~kk.
Fazendo o produto escalar do campo vetorial F~ por ~r u ×~~r v temos:
F~ •(~~r u ×~~r v ) = −2v 2 (a2 −v 2 ) sin(u)+2v 2 (a2 −v 2 ) sin(u)−v 3 sin(u) cos(u)
159
Simplificando e calculando o fluxo do campo vetorial F~ sobre a

superfície S temos:
Z a Z +π
φ(F~ ) = −v 3 sin(u) cos(u)dudv
0 −π
a
sin2 (u) +π
Z
= −v 3 dv
0 2 −π
Z a
1
= −v 3 (sin2 (+π) − sin2 (−π))dv
2 0
= 0
8.6 Conclusão
Na aula de hoje, vimos como integrar funções definidas so-

bre uma superfície no espaço. Funções escalares de valores reais
ao longo de superfícies no espaço com as quais podemos determi-
nar área, massa, momento de massa, centro de massa momento
de inércia de uma superfície representando uma casca fina. Vi-
mos também como calcular integrais de campos vetoriais (funções
vetoriais) definidos sobre uma superfície no espaço, que essenci-
almente, os conceitos por trás da integração de campos vetoriais
como o fluxo através de superfícies estão intimamente ligados à
Física.
RESUMO
Caros alunos, em nossa aula de hoje, sobre integrais de fun-

ções definidas sobre superfícies no espaço, tanto funções escalares
quanto campos vetoriais o conteúdo visto pode ser resumido como:
160
Cálculo III AULA
Área de uma Superfície S ⊂ R3 8
Sejam S ⊂ R3 uma superfície dada por f (x, y, z) = c cuja projeção
em um dos planos coordenados seja D e p~ a normal a D. A área
da superfície S é dada por:
|∇f |
Z Z Z Z
Are(S) = dσ = dA
S D |∇f • p~ |
Integral de Superfície S ⊂ R3
Sejam S ⊂ R3 uma superfície dada por f (x, y, z) = c cuja projeção
em um dos planos coordenados seja D e p~ a normal a D e g : D ⊂
R3 7→ R uma função de valores reais cujo domínio é a superfície S.
A a integral de g sobre a superfície S é dada por:
|∇f |
Z Z Z Z
g(x, y, z)dσ = g(x, y, z) dA
S D |∇f • p~ |
Massa e Momento de Massa de uma Superfície S ⊂ R3

Se a superfície S ⊂ R3 representa uma casca fina de densidade
superficial % : S ⊂ R3 7→ R+ então, a massa, momento de massa
relativo aos planos coordenados yz, xz e xy, são dados respectiva-
mente por:
|∇f |
Z Z Z Z
m(S) = %(x, y, z)dσ = %(x, y, z)
dA
S D |∇f • p~ |
|∇f |
Z Z Z Z
Myz (S) = %(x, y, z)xdσ = %(x, y, z)x dA
S D |∇f • p~ |
|∇f |
Z Z Z Z
Mxz (S) = %(x, y, z)ydσ = %(x, y, z)y dA
S D |∇f • p~ |
|∇f |
Z Z Z Z
Mxy (S) = %(x, y, z)zdσ = %(x, y, z)z dA
S D |∇f • p~ |
Centro de Massa de uma Superfície S ⊂ R3

161
superficial % : S ⊂ R3 7→ R+ então, O centro de massa da casca

fina, denotado (x̄, ȳ, z̄), é dado por:
|∇f |
Z Z
%(x, y, z)x dA
Myz (S) D |∇f • p~ |
x̄ = = Z Z
m(S) |∇f |
%(x, y, z) dA
D |∇f • p~ |
|∇f |
Z Z
%(x, y, z)y dA
Mxz (S) |∇f • p~ |
ȳ = = Z ZD
m(S) |∇f |
%(x, y, z) dA
D |∇f • p~ |
|∇f |
Z Z
%(x, y, z)z dA
Mxy (S) |∇f • p~ |
z̄ = = Z ZD
m(S) |∇f |
%(x, y, z) dA
D |∇f • p~ |
Momento de Inércia de uma Superfície S ⊂ R3

superficial % : S ⊂ R3 7→ R+ então, os momentos de inércia com
relação aos eixos x, y e z são dados, respectivamente, por:
|∇f |
Z Z
Ix (S) = %(x, y, z)(y 2 + z 2 ) dA
D |∇f • p~ |
|∇f |
Z Z
Iy (S) = %(x, y, z)(x2 + z 2 ) dA
|∇f • p~|
Z ZD
|∇f |
Iz (S) = %(x, y, z)(x2 + y 2 ) dA
D |∇f • p~ |
Fluxo de um Campo Vetorial Através de uma Superfície

S ⊂ R3
Sejam S ⊂ R3 uma superfície lisa no espaço e F~ : S ⊂ R3 7→ R3
~ através
uma função de valores vetoriais. Definimos o fluxo de F
de S, denotado φ(F~ ), por:
Z Z
def
φ(F~ ) = F~ (x, y, z) • ~n dσ.
S
162
Cálculo III AULA
Onde n é a normal unitária em S. 8
Área de uma Superfície S ⊂ R3 Parametrizada
Sejam S ⊂ R3 uma superfície lisa definida parametricamente por
~r (u, v) = x̂(u, v)~i + ŷ(u, v)~j + ẑ(u, v)~kk, ∀(u, v) ∈ [a, b] × [c, d] então,
a área de S será:
Z bZ d
Are(S) = |~~r u × ~r v |dudv.
a c
Integral de Superfície S ⊂ R3 Parametrizada

Sejam S ⊂ R3 uma superfície lisa definida parametricamente por
~r (u, v) = x̂(u, v)~i + ŷ(u, v)~j + ẑ(u, v)~kk, ∀(u, v) ∈ [a, b] × [c, d] e
f : S ⊂ R3 7→ R uma função de valores reais definida sobre S
então, a integral de f sobre S será:
Z Z Z bZ d
def
f (x, y, z)dσ = f (x̂(u, v), ŷ(u, v), ẑ(u, v))|~~r u ×~~r v |dudv.
S a c
Fluxo de um Campo Vetorial Através de uma Superfície

S ⊂ R3 Parametrizada
Se S ⊂ R3 é lisa e definida parametricamente por ~r (u, v) = x̂(u, v)~ii+
ŷ(u, v)~j + ẑ(u, v)~kk, ∀(u, v) ∈ [a, b] × [c, d] e F~ : S ⊂ R3 7→ R3 uma
função de valores vetoriais. O fluxo de F~ através de S, é dado por:
Z bZ d
φ(F~ ) = F~ (x̂(u, v), ŷ(u, v), ẑ(u, v)) • ~n |~~r u × ~r v |dudv.
a c
1
Como podemos calcular o vetor normal por ~n = ·(~~r u ×~~r v )
|~~r u × ~r v |
a integral para o fluxo do campo vetorial F~ através da superfície
163
S pode ser reescrita como:
Z bZ d
φ(F~ ) = F~ (x̂(u, v), ŷ(u, v), ẑ(u, v)) • (~~r u × ~r v )dudv.
a c
PRÓXIMA AULA
Em nossa próxima aula veremos dois importantíssimos teore-

mas do Cálculo. São eles o “Teorema de Green” e “Teorema de
Stokes”. Dizem respeito a integração de campos vetoriais ao longo
de curvas fechadas no plano (caso do teorema de Green) e de curvas
fechadas no espaço (caso do teorema de Stokes).
ATIVIDADES
Deixamos como atividades as seguintes questões:
ATIV. 8.1. Seja S ⊂ R3 uma superfície dada por z = f (x, y)

cuja projeção no plano xy é D ⊂ xy. Mostre que sua área pode
ser dada por:
s 2 2
Z Z
∂f ∂f
+ + 1dxdy
D ∂x ∂y

demonstrações acima, elas lhe servirão de guia.
ATIV. 8.2. Seja uma casca fina dada pela superfície S ⊂ R3 des-
crita por f (x, y, z) = a2 − x2 − y 2 − z 2 = 0, y < 0 e z > 0 (ver
Fig. 8.8) e determine seu centro de gravidade.
164
Cálculo III AULA
8
Figura 8.8: 1/4 da esfera x2 + y 2 + z 2 = a2
Comentário: Observe que a superfície tem simetria e como a

densidade é constante temos: x̄ = 0 e ȳ = −z̄.

2003.
165
166
AULA
Teorema de Green
e Teorema de Stokes 9
META:
Apresentar o teorema de Green e o teorema de Stokes e algumas
de suas aplicações.
OBJETIVOS:
Enunciar o teorema de Green e o teorema de Stokes e determinar
o fluxo rotacional de um dado campo vetorial.
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I, aula 07 e aula 08.
Teorema de Green e Teorema de Stokes
9.1 Introdução
Caros alunos, nosso tema de hoje: “Teorema de Green e Te-

orema de Stokes", visa apresentar dois dos mais importantes te-
oremas do Cálculo, envolvendo integrais de linha. O teorema de
Green, converte integrais de linha de curvas fechadas no plano em
integrais duplas e o teorema de Stokes, que é uma generalização do
teorema de Green, converte integrais d linha de campos vetoriais
sobre curvas fechadas no espaço em integrais de superfície.
9.2 Preliminares
Antes de partirmos para a demonstração do teorema de Green

é necessária a introdução de alguns conceitos, que estabelecerão
BIOGRAFIA
a definição de dois novos operadores diferenciais vetoriais. O pri-
George Green nasceu
meiro deles é o conceito de “densidade de fluxo em um ponto”.
em Sneinton, condado
de Nottinghamshire 14 Seja F~ : D ⊂ R2 7→ R2 um campo vetorial (velocidade de es-
de Julho de 1793 e mo-
erreu em Nottingham, ~ (x, y) = f1 (x, y)~i + f2 (x, y)~jj,
coamento de um fluido) onde: F
31 de Maio de 1841,
foi um matemático e
físico inglês. Na sua
obra Essay on the
Application of Mathe-
matical Analysis to the
Theory of Electricity
and Magnetism (1828)
introduziu a noção de
função potencial no
estudo dos campos
magnéticos. O teorema
de Green, que demons-
trou em 1828 facilitou
bastante o estudo das
funções. Wikipedia Figura 9.1: Fluxo de F~ através do Retângulo
f1 , f2 : D ⊂ R2 7→ R são contínuas e deriváveis em D. Consi-

deremos um pequeno retângulo contido em D com vértices em:
168
Cálculo III AULA
(x, y), (x + ∆x, y), (x, y + ∆y) e (x + ∆x, y + ∆y) (ver Fig. 9.1 9
). A taxa com que F~ atravessa as arestas do retângulo são:
~ (x, y + ∆y) • ∆x~j = f2 (x, y + ∆y)∆x

em cima : F
em baixo : F~ (x, y) • (−∆x)~j = −f2 (x, y)∆x
à direita : F~ (x + ∆x, y) • ∆y~i = f1 (x + ∆x, y)∆y
à esquerda : F~ (x, y) • (−∆y)~i = −f1 (x, y)∆y
O fluxo total através das arestas do retângulo é:
(f1 (x + ∆x, y) − f1 (x, y))∆y + (f2 (x, y + ∆y) − f2 (x, y))∆x
Como podemos escrever as seguintes aproximações: f1 (x+∆x, y)−

∂f1 ∂f2
f1 (x, y) ≈ ∆x e f2 (x, y + ∆y) − f1 (x, y) ≈ ∆y temos a
∂x ∂x
seguinte aproximação para o fluxo de F~ através das arestas do
retângulo:
∂f1 ∂f2
∆x∆y + ∆y∆x
∂x ∂y
Dividindo pela área do retângulo, que é dada por: ∆x∆y teremos
uma aproximação para a densidade de fluxo de F~ no retângulo.
∂f1 ∂f2
+
∂x ∂y
Agora vamos à argumentação heurística em que fazemos ∆x e ∆y

tenderem a zero e podemos definir a densidade de fluxo do campo
vetorial F~ no ponto (x, y).
Definição 9.1. Seja F~ : D ⊂ R2 7→ R2 um campo vetorial onde:

F~ (x, y) = f1 (x, y)~ii+f2 (x, y)~jj, f1 , f2 : D ⊂ R2 7→ R são contínuas e
deriváveis em D. Definimos a densidade de fluxo do campo vetorial
169
F~ no ponto (x, y), denominado divergente de F~ , denotado DivF~

ou ∇ • F~ por:
def ∂f1 ∂f2
∇ • F~ = +
∂x ∂y
Completaremos com “densidade de circulação em um ponto” os
dois novos conceitos necessários ao estabelecimento do teorema de
Green. Seja F~ : D ⊂ R2 7→ R2 um campo vetorial (velocidade
~ (x, y) = f1 (x, y)~i + f2 (x, y)~jj,
de escoamento de um fluido) onde: F
Figura 9.2: Circulação de F~ ao longo do Retângulo
f1 , f2 : D ⊂ R2 7→ R são contínuas e deriváveis em D. Considere-

mos um pequeno retângulo contido em D com vértices em: (x, y),
(x + ∆x, y), (x, y + ∆y) e (x + ∆x, y + ∆y) (ver Fig. 9.2 ). A
soma das taxas de escoamento dos quatro lados do retângulo mede
a circulação de F~ no sentido anti-horário. As taxas de escoamento
ao longo de cada lado do retângulo são dadas por:
~ (x, y + ∆y) • (−∆x~ii) = −f1 (x, y + ∆y)∆x

em cima : F
em baixo : F~ (x, y) • ∆x~i = f1 (x, y)∆x
à direita : F~ (x + ∆x, y) • ∆y~j = f2 (x + ∆x, y)∆y
à esquerda : F~ (x, y) • (−∆y)~j = −f2 (x, y)∆y
170
Cálculo III AULA
O escoamento total ao longo das arestas do retângulo é: 9
(f2 (x + ∆x, y) − f2 (x, y))∆y + (−f2 (x, y + ∆y) + f2 (x, y))∆x
Como podemos escrever as seguintes aproximações: f1 (x+∆x, y)−

∂f1 ∂f2
f1 (x, y) ≈ ∆x e f2 (x, y + ∆y) − f1 (x, y) ≈ ∆y temos a
∂x ∂x
seguinte aproximação para o escoamento de F~ ao longo das arestas
do retângulo:
∂f2 ∂f1
∆x∆y − ∆y∆x
∂x ∂y
Dividindo pela área do retângulo, que é dada por: ∆x∆y teremos
uma aproximação para a densidade da circulação de F~ no retân-
gulo.
∂f2 ∂f1
−
∂x ∂y
Fazendo ∆x e ∆y tenderem independentemente a zero podemos
definir a densidade de circulação no ponto (x, y) por:
Definição 9.2. Seja F~ : D ⊂ R2 7→ R2 um campo vetorial onde:

F~ (x, y) = f1 (x, y)~i + f2 (x, y)~jj, f1 , f2 : D ⊂ R2 7→ R são contí-
nuas e deriváveis em D. Definimos a componente ~k da densidade
de circulação do campo vetorial F~ no ponto (x, y), denominado
rotacional de F~ , denotado RotF~ ou ∇ × F~ por:
def ∂f2 ∂f1

(∇ × F~ ) • ~k = −
∂x ∂y
9.3 Teorema de Green
O teorema de Green pode ser exposto de duas formas na pri-

meira diz que o fluxo exterior de um campo vetorial através de
171
uma curva fechada simples no plano é igual a integral dupla do di-

vergente do campo vetorial sobre a região D limitada pela curva,
i.e.
I Z Z
∂f1 ∂f2
F~ • ~n ds = + dxdy
C D ∂x ∂y
Na outra forma diz que a circulação no sentido anti-horário de um
campo de força ao longo de uma curva simples e fechada no plano
é igual a integral dupla da componente ~k do rotacional do campo
vetorial sobre a região D limitada pela curva.
I Z Z
∂f2 ∂f1
F~ • ~tds = − dxdy
C D ∂x ∂y
Vamos agora e sem muita demora à demonstração do teorema de
Green.
Teorema 9.1. Sejam F~ : D ⊂ R2 7→ R2 um campo vetorial dado

por: F~ (x, y) = f1 (x, y)~i + f2 (x, y)~jj, f1 , f2 : D ⊂ R2 7→ R são
contínuas e deriváveis em D e C ⊂ D ⊂ R2 uma curva simple,
fechada tal que retas paralelas aos eixos coordenados não a cortem
Figura 9.3: Teorema de Green
172
Cálculo III AULA
em mais que dois pontos e Seja R a região limitada por C então: 9
I Z Z
∂f2 ∂f1
f1 dx + f2 dy = − ) dxdy
C R ∂x ∂y
PROVA: Observando a Fig. 9.3 vemos C formada por duas
parte orientadas dadas por:

 C1 : y = g1 (x) a ≤ x ≤ b
 C : y = g (x) b ≥ x ≥ a
2 2
∂f1
Tomando um ponto arbitrário x ∈ (a, b) podemos integrar em
∂y
relação a y nos limites y = g1 (x) até y = g2 (x). A saber:
Z g2 (x)
∂f1 g2 (x)
dy = f1 (x, y) = f1 (x, g2 (x)) − f1 (x, g1 (x))

g1 (x) ∂y g1 (x)
Podemos então integrar este resultado na variável x nos limites

x = a até x = b e temos:
Z bZ g2 (x) Z b
∂f1
dydx = (f1 (x, g2 (x)) − f1 (x, g1 (x)))dx
a g1 (x) ∂y a
Z a Z b
= − (f1 (x, g2 (x))dx − f1 (x, g1 (x)))dx
Zb Z a
= − f1 dx − f1 dx
C2 C1
I
= − f1 dx
C
Que podemos reescrever como:
I Z Z
∂f1
f1 dx = − dxdy
C R ∂y
Por outro lado, observando a Fig. 9.4 a curva C pode ser dada
por:
173

 C1 : x = h1 (y) c ≤ y ≤ d
 C : x = h (y) d ≥ y ≥ c
2 2
Figura 9.4: Teorema de Green
∂f2
Tomando um ponto arbitrário y ∈ (c, d) podemos integrar em
∂x
relação a x nos limites x = h1 (y) até x = h2 (y). A saber:
Z h2 (y)
∂f2 h2 (y)
dx = f2 (x, y) = f2 (h2 (y), y) − f2 (h1 (y), y)

h1 (y) ∂x h1 (y)
Podemos então integrar este resultado na variável y nos limites

y = c até y = d e temos:
Z d Z h2 (y) Z d
∂f2
dxdy = (f2 (h2 (y), y) − f2 (h1 (y), y))dy
c h1 (y) ∂x c
Z d Z c
= f2 (h2 (y), y)dy + f2 (h1 (y), y)dy
cZ Z d
= − f2 dy − f2 dy
C2 C1
I
= f2 dy
C
Que podemos reescrever como:
174
Cálculo III AULA
I Z Z
∂f2
9
f2 dy = dxdy
C R ∂x
Adicionando os dois resultados temos:
I Z Z
∂f2 ∂f1
C R ∂x ∂y
9.4 Estendendo o Teorema de Green para Ou-

tras Regiões
Caros alunos, muito embora a demonstração acima do teorema

de Green comportem um grande número de curvas, muitas outras
não se enquadram na categoria que goza da propriedade de que
toda reta
Figura 9.5: Curva Figura 9.6: Partição
paralela aos eixos coordenados a cortem em no máximo dois pon-

tos É o caso da curva dada na Fig. 9.5 visto que retas paralelas
ao eixo x podem corta-la em até 4 pontos na parte de cima da
curva. Nestes casos podemos estender o teorema de Green com o
seguinte procedimento: Traçamos retas unindo pontos específicos
da curva C como na Fig. 9.6 os pontos A, B e D de modo que
175
as curvas formadas pela união de C1 com a reta AB bem como

a união da curva C2 com a reta BD e a curva C3 com as retas
DB e BA orientadas todas no sentido anti-horário contornando as
regiões R1 R2 e R3 , respectivamente. As curvas acima descritas
todas satisfazem a condição da demonstração do teorema de Green
acima i.e. toda reta paralela aos eixos coordenados as cortam em
apenas dois pontos. Desta forma aplicando o teorema de Green a
cada uma das curvas temos:
I Z Z Z
∂f2 ∂f1
f1 dx + f2 dy + f1 dx + f2 dy = − ) dxdy
IC1 ZAB Z ZR1 ∂x ∂y
∂f2 ∂f1
IC2 ZBD Z R2 ∂x ∂y
f1 dx + f2 dy + f1 dx + f2 dy + f1 dx + f2 dy =
ZC3Z BA
DB
∂f2 ∂f1
− ) dxdy
R3 ∂x ∂y
Adicionando todas as equações e levando em conta que:
Z Z
f1 dx + f2 dy = − f1 dx + f2 dy
Z AB ZBA
f1 dx + f2 dy = − f1 dx + f2 dy,
BD DB
temos:
Z Z Z
∂f2 ∂f1
C1 ∪C2 ∪C3 R1 ∪R2 ∪R3 ∂x ∂y
E como C = C1 ∪ C2 ∪ C3 e R = R1 ∪ R2 ∪ R3 podemos finalmente

escrever que para a curva C da Fig. 9.5 vale o teorema de Green:
Z Z Z
∂f2 ∂f1
f1 dx + f2 dy = − ) dxdy.
C R ∂x ∂y
176
Cálculo III AULA
Uma outra categoria de regiões que não satisfaz a exigência da 9
demonstração do teorema de Green são regiões R com buracos
Figura 9.7: Curva Figura 9.8: Partição
como a da Fig. 9.7 onde a curva C ∗ que define a região do buraco

está orientada no sentido horário. Uma solução está na partição
dada na Fig. 9.8 onde repartimos R em oito sub-regiões R1 a
R8 da seguinte forma: R1 limitada pelas curvas C1 , A2 A3 e A3 A1 ,
R2 limitada pelas curvas C2 , A4 A6 , C9 e A3 A2 , R3 limitada pelas
curvas C3 , A5 A6 e A6 A4 , R4 limitada pelas curvas C4 , A7 A9 , C10
e A6 A5 , R5 limitada pelas curvas C5 , A8 A9 e A9 A7 , R6 limitada
pelas curvas C6 , A10 A12 , C11 e A9 A8 , R7 limitada pelas curvas
C7 , A11 A12 e A12 A10 , R8 limitada pelas curvas C8 , A1 A3 , C12 e
A12 A11 . Em cada uma das oito regiões o teorema de Green pode
ser aplicado individualmente e levando em conta que as integrais
ao longo dos segmentos de reta se anulam mutuamente podemos
escrever para a região R que:
Z Z Z Z
∂f2 ∂f1
C C∗ R ∂x ∂y
Podemos então relaxar um pouco as exigências do teorema de
177
Green e reformula-lo como:
Teorema 9.2. Sejam F~ : D ⊂ R2 7→ R2 um campo vetorial dado

~ (x, y) = f1 (x, y)~i + f2 (x, y)~jj, f1 , f2 : D ⊂ R2 7→ R são
por: F
contínuas e deriváveis em D e C ⊂ D ⊂ R2 uma curva simple,
fechada e Seja R a região limitada por C então:
I Z Z
∂f2 ∂f1
C R ∂x ∂y
9.5 Verificação do Teorema de Green
Caros alunos, nesta seção faremos uma verificação do teorema

de Green. A saber:
Exemplo 9.1. Considere a campo vetorial F~ : R2 7→ R2 dado por

y ~ x ~
F~ (x, y) = − 2 2
i+ 2 j e a região com buraco limitadas
x +y x + y2
pelos círculos (ver Fig. 9.9):
Figura 9.9: Verificação do teorema de Green
178
Cálculo III AULA
C1 : x = b cos(t) y = b sin(t) t ∈ [0, 2π]

9
C2 : x = a cos(t) y = −a sin(t) t ∈ [0, 2π]
onde a < b, C1 é percorrida no sentido anti-horário e C2 no sen-

tido horário. Verificar o teorema de Green para os dados campo
vetorial e região.
SOLUÇÃO: Como o campo vetorial F~ (x, y) = −x2 y~i + xy 2~j te-
mos suas componentes dadas por f1 (x, y) = −x2 y e f2 (x, y) = xy 2
∂f1 ∂f2
e temos as derivadas e dadas por:
∂y ∂x
∂f1 (x2 + y 2 )(−1) − (−y)(2y)

=
∂y (x2 + y 2 )2
x2 − y 2
=
(x2 + y 2 )2
∂f2 (x2 + y 2 )(+1) − (−x)(2x)

=
∂x (x2 + y 2 )2
x2 − y 2
=
(x2 + y 2 )2
∂f1 ∂f2
Portanto = e temos:
∂y ∂x
Z Z
∂f2 ∂f1
− dxdy = 0
R ∂x ∂y
Calculando a integral de linha em C1 temos:
179
I I
1
f1 dx + f2 dy = (xdx − ydy)
C1 C1 x2 + y 2
2π
b2 cos2 (t) + b2 sin2 (t)
Z
= dt +
0 b2 cos2 (t) + b2 sin2 (t)
Z 2π
= dt
0
2π
= t

0
= 2π
E calculando a integral de linha em C2 temos:
I I
1
f1 dx + f2 dy = (xdx − ydy)
C2 C1 x2 + y2
2π
a2 cos2 (t) + a2 sin2 (t)
Z
= − dt +
0 a2 cos2 (t) + a2 sin2 (t)
Z 2π
= − dt
0
2π
= −t

0
= −2π
Temos então que:
I I
f1 dx + f2 dy + f1 dx + f2 dy = 0
C1 C2
E o teorema de Green é verificado:
I I Z Z
∂f2 ∂f1
f1 dx + f2 dy + f1 dx + f2 dy = − dxdy
C1 C2 R ∂x ∂y
180
Cálculo III AULA
9.6 Teorema de Stokes 9
Caros alunos, nesta seção trataremos do Teorema de Stokes,
que é uma generalização do Teorema de Green. Porém, devido a
dificuldade das técnicas usada em sua demonstração e que escapam
ao escopo deste curso, nos limitaremos a apresenta-lo e fazer uma
BIOGRAFIA
aplicação do mesmo.
Sir George Gabriel
Começamos por estender o conceito de densidade de circulação
Stokes nasceu em
(rotacional). Como vimos anteriormente a componente ~k da densi- Skreen, Condando de
Sligo, 13 de Agosto
dade de circulação de um campo vetorial bi-dimensional F~ (x, y) = de 1819 e morreu
∂f2 ∂f1 em Cambridge, 1 de
f1 (x, y)~i + f2 (x, y)~j é dada por: − . Em dimensão três, a Fevereiro de 1903,
∂x ∂y foi um matemático e
densidade de circulação (rotacional) é um vetor normal ao plano físico irlandês que se
distinguiu pelas suas
de circulação cuja direção satisfaz a regra da mão direita. A taxa contribuições na dinâ-
mica de fluidos (por
de rotação do fluido é medida pelo módulo do vetor rotacional. exemplo, as equações
de Navier-Stokes),
Vamos à definição:
na óptica e física
matemática (Teorema
Definição 9.3. Seja F~ : D ⊂ R3 7→ R3 uma campo vetorial de Stokes). Wikipedia
tridimensional dado por: F~ (x, y, z) = f1 (x, y, z)~i + f2 (x, y, z)~j +

f3 (x, y, z)~k onde f1 , f2 , f3 : D ⊂ R3 7→ R são funções contínuas e
com derivadas parciais de primeira ordem contínuas. O rotacional
de F~ , denotado RotF~ ou ∇ × F~ , é definido por:

def ∂f3 ∂f2 ~ ∂f1 ∂f3 ~ ∂f2 ∂f1 ~
∇ × F~ = − i+ − j+ − k
∂y ∂z ∂z ∂x ∂x ∂y
Vamos agora enunciar o Teorema de Stokes que nos diz que sob
certas condições a circulação de um campo vetorial ao longo da
borda de uma superfície S ⊂ R3 no sentido anti-horário com rela-
ção ao campo de vetores normais a S é igual a integral de superfície
do componente normal do rotacional.
181
Teorema 9.3. Sejam F~ : D ⊂ R3 7→ R3 um campo vetorial tridi-

mensional cujas componentes são contínuas e têm derivadas parci-
ais de primeira ordem contínuas e S ⊂ D ⊂ R3 uma superfície lisa
com borda C ⊂ D ⊂ R3 lisa então, a circulação do campo vetorial
F~ ao longo da borda C de S, no sentido anti-horário em relação
aos vetores normais unitários de S é igual a integral de superfície
da componente normal do rotacional i.e.:
I Z
F~ • d~~r = (∇ × f~) • ~n dσ
C S
OBS 9.1. Se duas superfícies orientadas diferentes S1 e S2 tem a

mesma borda C as integrais de superfície da componente normal
do rotacional de um campo vetorial que atravessa ambas são iguais:
Z Z
(∇ × f~) • ~n dσ = (∇ × f~) • ~n dσ
S1 S2
OBS 9.2. Se C é uma curva lisa do plano xy, orientada no sentido

anti-horário e R a região de xy delimitada por C o vetor normal a
R é ~n = ~kk. Daí, temos:
∂f2 ∂f1
∇ × F~ • ~n = ∇ × F~ • ~k = −
∂x ∂y
E o teorema de Stokes pode ser escrito como:
I Z Z
∂f2 ∂f1
F~ • d~~r = − dxdy
C R ∂x ∂y
Que é o teorema de Green. Daí, se diz que o teorema de Green é

um caso particular do teorema de Stokes.
182
Cálculo III AULA
9.7 Aplicação do Teorema de Stokes 9
Caros alunos, nesta seção faremos uma aplicação do teorema
de Stokes. A saber:
Figura 9.10: Verificação do teorema de Stokes
Exemplo 9.2. Considere a campo vetorial F~ : R3 7→ R3 dado por

F~ (x, y) = (yz + xz)~i + (xz + x)~j + (xy − y 2 /2)~k e a Curva C na
p
qual o plano z = a > 0 corta o cone z = x2 + y 2 (ver Fig. 9.10)
e determine a circulação de F~ ao longo de C.
~ ao longo de C
SOLUÇÃO: Para determinação da circulação de F
usaremos o teorema de Stokes. Percorrer C no sentido anti-horário
visto de cima corresponde a tomar a normal ~n ao cone apontando
para dentro.
O cone pode ser parametrizado por:
~r (r, ϑ) = ar cos(ϑ)~i + ar sin(ϑ)~j + ar~kk, ∀r ∈ [0, 1], ∀ϑ ∈ [0, 2π]
As derivadas parciais ~r r e ~r ϑ são dadas por:
183
~r r = a cos(ϑ)~i + a sin(ϑ)~j + a~k
~r ϑ = −ar sin(ϑ)~i + ar cos(ϑ)~j + 0~k
Daí, o produto vetorial ~r r × ~r ϑ pode ser calculado como:

~i ~j ~k

~r r × ~r ϑ = det a cos(ϑ)

a sin(ϑ) a

−ar sin(ϑ) ar cos(ϑ) 0
Daí, temos:
~r r × ~r ϑ = −ar cos(ϑ)~i − ar sin(ϑ)~j + ar~k
Como, para uma superfície parametrizada temos:
Z Z Z bZ d
F~ • ~n dσ = (∇ × F~ ) • (~~r u × ~r v )dudv
S a c
Calculando o rotacional de F~ temos:

~i ~j ~k

∂ ∂ ∂

∇ × F~ = det

∂x ∂t ∂z

yz + xz xz + x xy − y 2 /2
Daí, temos:
∇ × F~ = −y~i + x~j + ~k
E como ∇ × F~ (x, y, z) = −y~i + x~j + ~k temos:
184
Cálculo III AULA
9
∇ × F~ • (~~r u × ~r v ) = +a2 r2 cos(ϑ) sin(ϑ) − a2 r2 cos(ϑ) sin(ϑ) + ar
= ar
E podemos calcular a integral de superfície como:
Z Z Z 2π Z 1
F~ • ~n dσ = ardrdϑ = aπ
S 0 0
Usando o teorema de Stokes para calcular a circulação temos:
I Z Z
F~ • d~~r = (∇ × F~ ) • ~n dσ = aπ
C S
9.8 Conclusão
Na aula de hoje, vimos dois grandes teoremas do Cálculo. O

teorema de Green que relaciona o fluxo exterior através de uma
curva lisa do plano de um campo vetorial com a integral dupla
do divergente do campo vetorial sobre a região delimitada pela
curva. O teorema de Stokes, que relaciona a circulação de um
campo de vetores tridimensional ao longo da borda de uma super-
fície no espaço com a integral de superfície da componente normal
do rotacional do campo vetorial.
RESUMO
No nosso resumo de hoje constam as seguintes definições e te-

oremas:
185
Densidade de Fluxo de um Campo Vetorial Bi-dimensional
Seja F~ : D ⊂ R2 7→ R2 um campo vetorial onde: F~ (x, y) =

f1 (x, y)~ii+f2 (x, y)~jj, f1 , f2 : D ⊂ R2 7→ R são contínuas e deriváveis
em D. Definimos a densidade de fluxo do campo vetorial F~ no
ponto (x, y), denominado divergente de F~ , denotado DivF~ ou ∇•F~
por:
def ∂f1 ∂f2
∇ • F~ = +
∂x ∂y
Densidade de Circulação de um Campo Vetorial Bi-di-
mensional
Seja F~ : D ⊂ R2 7→ R2 um campo vetorial onde: F~ (x, y) =

f1 (x, y)~ii+f2 (x, y)~jj, f1 , f2 : D ⊂ R2 7→ R são contínuas e deriváveis
em D. Definimos a componente ~k da densidade de circulação do
~ no ponto (x, y), denominado rotacional de F
campo vetorial F ~,
denotado RotF~ ou ∇ × F~ por:
def ∂f2 ∂f1

∇ × F~ = −
∂x ∂y
Teorema de Green
Sejam F~ : D ⊂ R2 7→ R2 um campo vetorial dado por: F~ (x, y) =
f1 (x, y)~i + f2 (x, y)~jj, f1 , f2 : D ⊂ R2 7→ R são contínuas e derivá-
veis em D e C ⊂ D ⊂ R2 uma curva simple, fechada e Seja R a
região limitada por C então:
I Z Z
∂f2 ∂f1
C R ∂x ∂y
Rotacional
Seja F~ : D ⊂ R3 7→ R3 uma campo vetorial tridimensional dado
por: F~ (x, y, z) = f1 (x, y, z)~ii+f2 (x, y, z)~jj+f3 (x, y, z)~k onde f1 , f2 , f3 :
186
Cálculo III AULA
D ⊂ R3 7→ R são funções contínuas e com derivadas parciais de 9
primeira ordem contínuas. O rotacional de F~ , denotado RotF~ ou
∇ × F~ , é definido por:

def ∂f3 ∂f2 ~ ∂f1 ∂f3 ~ ∂f2 ∂f1 ~
∇ × F~ = − i+ − j+ − k
∂y ∂z ∂z ∂x ∂x ∂y
Teorema de Stokes
Sejam F~ : D ⊂ R3 7→ R3 um campo vetorial tridimensional cujas
componentes são contínuas e têm derivadas parciais de primeira
ordem contínuas e S ⊂ D ⊂ R3 uma superfície lisa com borda
C ⊂ D ⊂ R3 lisa então, a circulação do campo vetorial F~ ao
l.ongo da borda C de S, no sentido anti-horário em relação aos
vetores normais unitários de S é igual a integral de superfície da
componente normal do rotacional i.e.:
I Z
F~ • d~~r = (∇ × f~) • ~n dσ
C S
PRÓXIMA AULA
Na próxima aula estudaremos outro importante teorema do

Cálculo atribuído ao Matemático alemão Johann Carl Friedrich
Gauss resumidamente denominado de “Teorema da Divergência”.
ATIVIDADES
ATIV. 9.1. Sejam F~ : R2 7→ R2 o campo vetorial bi-dimensional

dado por: F~ (x, y) = y~ii+0~j e R ⊂ R2 a região limitada pelo círculo
187
C ⊂ R2 dado por: x2 + y 2 = a2 . Verifique o teorema de Green

para este campo vetorial e esta região.
texto e as verificações acima, elas lhe servirão de guia.
ATIV. 9.2. Sejam F~ : R3 7→ R3 o campo vetorial tridimensional

dado por: F~ (x, y) = x2~i + x~j + z~k e R ⊂ R2 a região limitada
pela elipse C ⊂ R2 dado por: a2 x2 + y 2 = a2 . Use o teorema de
Stokes e determine a circulação do campo F~ ao longo da curva C
no sentido anti-horário quando vista de cima.
texto e as verificações acima, elas lhe servirão de guia.

2003.
188
AULA
Teorema de Divergência
10
META:
Apresentar o teorema de Gauss e algumas de suas aplicações.
OBJETIVOS:
Enunciar o teorema de Gauss.
Determinar o divergente de um campo vetorial e determinar o fluxo
de um campo vetorial através de uma superfície fechada em R3 .
PRÉ-REQUISITOS
mínio em R, da disciplina Cálculo I e aula 08.
10.1 Introdução
Caros alunos terminamos aqui nosso curso de Cálculo III com o

tema “Teorema da Divergência”, atribuído ao Matemático alemão
BIOGRAFIA
Johann Carl Friedrich Gauss e mais tarde atribuido também ao
Johann Carl Frie-
Matemático russo Mikhail Vasilievich Ostrogradsky. O teorema
drich Gauss nasceu em
Braunschweig, 30 de de Gauss, ou teorema da divergência, relaciona uma integral tripla
Abril de 1777 e morreu
em Göttingen, 23 de num sólido D ⊂ R3 com a integral sobre a superfície S ⊂ R3 que
Fevereiro de 1855,
foi um matemático, é fronteira deste sólido.
astrônomo e físico
alemão. Conhecido
como o príncipe dos
matemáticos. Muitos
consideram Gauss o 10.2 Preliminares
maior gênio da história
da matemática. Seu QI
foi estimado em cerca Como preliminares precisaremos apenas estender a definição
de 240. Wikipedia
de divergente de um campo vetorial bi-dimensional, visto na aula
anterior, para o divergente de um campo vetorial tridimensional.
Vamos logo à tarefa.
Definição 10.1. Seja F~ : D ⊂ R3 7→ R3 um campo vetorial

tridimensional dado por F~ (x, y, z) = f1 (x, y, z)~i + f2 (x, y, z)~j +
f3 (x, y, z)~k tal que suas componentes f1 , f2 , f3 : D ⊂ R3 7→ R
sejam contínuas e tenham derivadas parciais de primeira ordem
contínuas em D. Definimos o divergente de F~ , denotado DivF~ ou
∇ • F~ , por:
def ∂f1 ∂f2 ∂f3

∇ • F~ = + +
∂x ∂y ∂z
Exemplo 10.1. Seja F~ : D ⊂ R3 7→ R3 o campo vetorial tridi-

mensional dado por F~ (x, y, z) = x2 yz~i + xy 2 z~j + xyz 2~kk. O seu
divergente será:
190
Cálculo III AULA
10
∂ 2 ∂ ∂
∇ • F~ = (x yz) + (xy 2 z) + (xyz 2 )
∂x ∂y ∂z
= 2xyz + 2xyz + 2xyz
= 6xyz
10.3 Teorema da Divergência
Caros alunos, nesta seção vamos estudar a transformação de

certas integrais de volume em integrais de superfícies que é o aná-
logo no espaço do teorema de Green (em sua forma divergente)
estudado na aula anterior. Como condições impostas primeira-
mente ao campo vetorial F~ é que ele tenha componentes contínuas
e com derivadas parciais de primeira ordem contínuas o que já
basta para o nosso propósito. Quanto a região D ⊂ R3 , deseja-
mos que ela seja regular e suave, tenha fronteira S ⊂ R3 regular
e suave, que suas projeções Sxy no plano xy, Syz no plano yz e
Sxz no plano xz sejam regiões fechadas de R2 com fronteira suave
e que retas paralelas aos eixos coordenados que atravessem suas
projeções cortem S em no máximo dois pontos Uma tal região
será aqui chamada de região simples. Um exemplo de tal região é
a limitada pela esfera x2 + y 2 + z 2 = a2 . Vamos ao enunciado e a
demonstração do teorema da divergência.
Teorema 10.1. Seja F~ : D ⊂ R3 7→ R3 um campo vetorial

191
contínuas em D e D ⊂ R3 uma região do espaço regular e suave

tal que sua fronteira fronteira S ⊂ R3 seja regular e suave, que
suas projeções Sxy no plano xy, Syz no plano yz e Sxz no plano
xz sejam regiões fechadas de R2 com fronteira suave e que retas
paralelas aos eixos coordenados que atravessem suas projeções cor-
tem S em no máximo dois pontos (ver Fig. 10.1) então:
Figura 10.1: Teorema da divergência
Z Z Z Z Z
∇ • F~ dxdydz = F~ • ~n dσ
D S
PROVA: Quando projetamos uma região regular e simples D no

plano xy, sua fronteira S consiste de duas partes S1 e S2 dadas
pelas funções: z1 (x, y) e z2 (x, y) respectivamente (ver Fig. 10.1).
Seja f3 (x, y, z) uma função contínua com derivadas parciais de pri-
meira ordem contínuas em D. Então:
Z Z Z Z Z Z z2 (x,y)
∂f3 ∂f3
dV ol = dxdy dz
D ∂z Sxy z1 (x,y) ∂z
192
Cálculo III AULA
Integrando e substituindo os limites temos: 10
Z Z Z Z Z
∂f3
dV ol = f3 (x, y, z2 (x, y))dxdy −
D ∂z Sxy
Z Z
− f3 (x, y, z1 (x, y))dxdy
Sxy
Seja ~n = nx~i + ny~j + nz~k a normal a S em cada ponto (x, y, z) ∈ S

apontando para fora de D e seja γ o ângulo entre ~n e ~k desta
forma temos: ~n • ~k = nz = cos(γ).|~~n |.|~kk| i.e. cos(γ) = nz em S2 e
cos(γ) = −nz em S1 . Daí, e do fato de que dxdy = cos(γ)dσ onde
dσ é o elemento de área em S, temos:
Z Z Z Z
f3 (x, y, z2 (x, y))dxdy = f3 nz dσ
Sxy S2
e também:
Z Z Z Z
f3 (x, y, z1 (x, y))dxdy = − f3 nz dσ
Sxy S1
Daí, e da expressão anterior temos:
Z Z Z Z Z Z Z
∂f3
dV ol = f3 nz dσ + f3 nz dσ
D ∂z S2 S1
Z Z
= f3 nz dσ
S2 ∪S1
Como S = S2 ∪ S1 temos:
Z Z Z Z Z
∂f3
dV ol = f3 nz dσ
D ∂z S
De forma análoga, usando as projeções de D sobre os planos coor-
denados yz e xz podemos deduzir que:
193
Z Z Z Z Z
∂f1
dV ol = f1 nx dσ
D ∂x S
e também:
Z Z Z Z Z
∂f2
dV ol = f2 ny dσ
D ∂y S
Somando as três equações acima temos:
Z Z Z Z Z
∂f1 ∂f2 ∂f3
+ + dV ol = (f1 nx + f2 ny + f3 nz )dσ
D ∂x ∂y ∂z S
Como F~ = f1~i + f2~j + f3~k temos:
F~ • ~n = f1 nx + f2 ny + f3 nz
∂f1 ∂f2 ∂f3
∇ • F~ = + +
∂x ∂y ∂z
E temos então:
Z Z Z Z Z
∇ • F~ dV ol = F~ • ~n dσ
D S
10.4 Estendendo o Teorema da Divergência
Caros alunos, muito embora a demonstração acima do teorema

da divergência comportem um grande número de regiões D, muitas
outras não se enquadram na categoria que goza da propriedade de
ser uma região do espaço regular e suave tal que sua fronteira fron-
teira S ⊂ R3 seja regular e suave, que suas projeções Sxy no plano
xy, Syz no plano yz e Sxz no plano xz sejam regiões fechadas de
R2 com fronteira suave e que retas paralelas aos eixos coordenados
194
Cálculo III AULA
10
Figura 10.2: Teorema da divergência
que atravessem suas projeções cortem S em no máximo dois pon-

tos. È o caso da região dada pela Fig. 10.2. Observando porém,
que se a região D puder ser decomposta em um número finito de
sub-regiões simples, podemos escrever o teorema da divergência
em cada uma das sub-regiões e somar o resultado, de forma que
para a região D o teorema continue válido. Desta forma podemos
enunciar uma forma mais ampla do teorema da divergência. A
saber:
Teorema 10.2. Seja F~ : D ⊂ R3 7→ R3 um campo vetorial

contínuas em D e D ⊂ R3 uma região do espaço regular e suave
tal que possa ser subdividida em um número finito de regiões sim-
ples e sua fronteira fronteira S ⊂ R3 seja regular e suave, então:
Z Z Z Z Z
D S
195
10.5 Algumas Aplicações do Teorema da Di-

vergência
Caros alunos, nesta seção faremos algumas aplicações do teo-

rema da divergência. A primeira é no calculo do fluxo exterior de
um campo vetorial tridimensional, a segunda no cálculo do fluxo
exterior do campo elétrico gerado por uma carga pontual através
de uma superfície fechada em cujo interior encontra-se a carga elé-
trica e a terceira na redução da forma integral da lei de balanço de
massa à sua forma diferencial pontual.
~ : R3 7→ R3 o campo vetorial tridimensi-

Exemplo 10.2. Sejam F
onal dado por: F~ (x, y) = x~i + y~j + z~k e D a região delimitada pela
esfera x2 + y 2 + z 2 ≤ a2 . Determine o fluxo exterior de F~ dado por
Z Z
F~ • ~n dσ através da superfície da esfera.
S
Vamos ao calculo do fluxo exterior de um campo vetorial tridimen-
sional.
SOLUÇÃO: Do teorema da divergência temos:
Z Z Z Z Z
F~ • ~n dσ = ∇ • F~ dV ol
S D
Portanto basta calcular a integral tripla sobre a região da esfera

do divergente de F~ .
Como F~ (x, y) = x~i + y~j + z~k seu divergente será:
∂ ∂ ∂
∇ • F~ = (x) + (y) + (z)
∂x ∂y ∂z
= 3
196
Cálculo III AULA
Daí, temos: 10
Z Z Z Z Z
F~ • ~n dσ = 3dV ol
S Z Z D
Z
= 3 dV ol
D
= 4πa3
Vamos ao cálculo do fluxo exterior do campo elétrico gerado por

uma carga pontual através de uma superfície S ⊂ R3 em cujo in-
terior encontra-se a carga elétrica.
Exemplo 10.3. O campo elétrico gerado por uma carga elétrica

pontual q localizada na origem é dado por:
~ = q
E ~r
4π0 |~~r |3
Como ~r = x~i + y~j + z~kk, colocando φ = |~~r | =
p
x2 + y 2 + z 2 temos:
~ = q ~
E F
4π0
1
onde F~ = 3 ~r .
φ
Podemos escrever o vetor F~ em suas componentes como:
x y z
F~ = 3~i + 3~j + 3 ~k
φ φ φ
Calculando as derivadas parciais das componentes de F~ temos:
∂φ
∂

x
φ3 − 3xφ2
= ∂x
∂x φ3 φ6
197
p
Como φ = |~~r | = x2 + y 2 + z 2 temos:
∂φ x x
=p =
∂x 2 2
x +y +z 2 φ
Logo:
φ3 − 3x2 φ

∂ x
=
∂x φ3 φ6
De modo análogo temos:
φ3 − 3y 2 φ

∂ y
=
∂y φ3 φ6
e também:
φ3 − 3z 2 φ

∂ z
=
∂z φ3 φ6
Somando as três equações temos:
3φ3 − 3(x2 + y 2 + z 2 )φ

∂ x ∂ y ∂ z
+ + =
∂x φ3 ∂y φ3 ∂z φ3 φ6
3φ3 − 3φ2 φ
=
φ6
= 0
Logo como:

∂ ∂ y ∂ x
z
∇ • F~ =
+ + =0
∂x
∂y φ 3 φ3
∂z φ3
~ = q ∇ • F~ temos:
E também como ∇ • E
4π0
~ =0
∇•E
198
Cálculo III AULA
Tomando D∗ como a região entre S e uma esfera centrada na 10
origem Sa e cujo raio a seja suficiente para que S permaneça no
interior da esfera eaplicando o teorema da divergência temos:
Z Z Z Z Z
~
E • ~n dσ = ∇•E ~
EdV ol = 0
S∪Sa D∗
~ através de S no sentido que se afasta da origem

Logo o fluxo de E
~ através de Sa no sentido que se afasta
é o mesmo que o fluxo de E
da origem.
~ através de Sa no sentido que se afasta da origem
Como fluxo de E
q
é temos:
0
Z Z
~ • ~n dσ = q
E
S 0
Como última aplicação veremos como reduzir da forma integral
para a forma diferencial pontual a equação de balanço de massa
conhecida como Lei de Lavoisier.
Sejam D ⊂ R3 uma região regular e suave do espaço e ~v : D ⊂
R3 7→ R3 um campo de velocidade de um fluido cuja densidade de
massa é dada por % : D ⊂ R3 7→ R+ e que preenche D. Em sua
forma integral o balanço de massa estabelece que ∀D∗ ⊂ D ⊂ R3
regular e suave com fronteira S ∗ regular e suave vale:
Z Z Z Z Z
d
%dV ol + %~~v • ~n dσ = 0
dt D∗ S∗
onde a primeira integral representa a variação total da massa den-
tro da região D∗ e a segunda representa a variação de massa que
penetra em D∗ pela superfície S ∗ .
Usando o teorema da divergência temos:
Z Z Z Z Z Z
d
%dV ol + ∇ • (%~~v )dV ol = 0
dt D∗ D∗
199
Daí, tomando regiões D∗ que não variem com o tempo
Z Z Z Z Z Z
d ∂%
%dV ol = dV ol
dt D∗ D∗ ∂t
E podemos reformular a equação de balanço de massa para a forma:
Z Z Z
∂%
+ ∇ • (%~~v ) dV ol = 0
D∗ ∂t
Como a integral acima vale ∀D∗ ⊂ D ⊂ R3 podemos dividi-la
por V ol(D∗ ), usar o teorema do valor médio para integrais, fazer
V ol(D∗ ) tender a zero e concluir que:
∂%
+ ∇ • (%~~v ) = 0
∂t
Que a forma diferencial pontual da equação de balanço de massa.
10.6 Conclusão
Na aula de hoje, vimos um importante teorema do Cálculo

denominado “Teorema da Divergência” atribuído aos Matemáticos
Gauss e Ostrogradsky. Tem forte conotação física e é utilizado para
reduzir as leis de conservação de sua forma integral para forma
diferencial pontual.
RESUMO
No nosso resumo de hoje constam as seguintes definições e te-

oremas:
200
Cálculo III AULA
Divergente de um Campo Vetorial Tridimensional 10
Seja F~ : D ⊂ R3 7→ R3 um campo vetorial tridimensional dado
por F~ (x, y, z) = f1 (x, y, z)~i + f2 (x, y, z)~j + f3 (x, y, z)~k tal que suas
componentes f1 , f2 , f3 : D ⊂ R3 7→ R sejam contínuas e tenham
derivadas parciais de primeira ordem contínuas em D. Definimos
o divergente de F~ , denotado DivF~ ou ∇ • F~ , por:
def ∂f1 ∂f2 ∂f3

DivF~ = + +
∂x ∂y ∂z
Teorema da Divergência: Forma Restritiva

~ (x, y, z) = f1 (x, y, z)~i + f2 (x, y, z)~j + f3 (x, y, z)~k tal que suas
por F
derivadas parciais de primeira ordem contínuas em D e D ⊂ R3
uma região do espaço regular e suave tal que sua fronteira fronteira
S ⊂ R3 seja regular e suave, que suas projeções Sxy no plano xy,
Syz no plano yz e Sxz no plano xz sejam regiões fechadas de R2
com fronteira suave e que retas paralelas aos eixos coordenados
que atravessem suas projeções cortem S em no máximo dois pon-
tos então:
Z Z Z Z Z
D S
Teorema da Divergência: Forma mais Ampla

por F~ (x, y, z) = f1 (x, y, z)~i + f2 (x, y, z)~j + f3 (x, y, z)~k tal que suas
201

derivadas parciais de primeira ordem contínuas em D e D ⊂ R3
uma região do espaço regular e suave tal que possa ser subdividida
em um número finito de regiões simples e sua fronteira fronteira
S ⊂ R3 seja regular e suave, então:
Z Z Z Z Z
D S
ATIVIDADES
ATIV. 10.1. Sejam F~ : R3 7→ R3 o campo vetorial tridimensional

dado por: F~ (x, y) = x2~i − 2xy~j + 3xz~k e D ⊂ R3 a região limitada
pela esfera x2 + y 2 + z 2 ≤ a2 e z ≥ 0 (acima do plano z = 0). Use
o teorema da divergência e determine o fluxo exterior através da
fronteira da região D.
texto e as aplicações acima, elas lhe servirão de guia.
ATIV. 10.2. Sejam F~ : R3 7→ R3 o campo vetorial tridimensi-

onal dado por: F~ (x, y) = x~i + y~j + z~k e suponha que a região
D ⊂ R3 seja uma região do espaço regular e suave tal que possa
ser subdividida em um número finito de regiões simples e sua fron-
teira fronteira S ⊂ R3 seja regular e suave,. Mostre que o volume
V ol(D) da região D é dado pela fórmula:
Z Z
1
V ol(D) = F~ • ~n dσ
3 S
202
Cálculo III AULA
. 10
texto e as aplicações acima, elas lhe servirão de guia.

2003.
203

LIVRO - DE - ..Calculo 3

Enviado por

Dados do documentoclique para ver informações do documento

Direitos autorais:

Formatos disponíveis

LIVRO - DE - ..Calculo 3

Enviado por

Dados do documento

Descrição original:

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

LIVRO - DE - ..Calculo 3

Enviado por

Direitos autorais:

Formatos disponíveis

Sumário

Aula 1: Integrais Duplas 11

Aula 2: Mudança de Variáveis em Integrais Duplas 33

Aula 3: Algumas Aplicações da Integral Dupla 47

3.3 Algumas Aplicações da Integral Dupla . . . . . . . . 52

Aula 4: Integrais triplas 63

4.2 Integração Tripla: Domínios Paralelepípedais . . . . 64

4.3 Integração Tripla: Domínios Não Paralelepípedais Li-

4.4 Interpretação Geométrica . . . . . . . . . . . . . . . 67

4.5 Integrais Iteradas . . . . . . . . . . . . . . . . . . . 67

4.6 Propriedades das Integrais Triplas . . . . . . . . . . 68

Aula 6: Algumas Aplicações das Integrais tríplas 105

Aula 7: Integrais de Funções Vetoriais sobre Curvas em R3

Aula 8: Integrais de Superfícies 145

Aula 9: Teorema de Green e Teorema de Stokes 167

Aula 10: Teorema de Divergência 189

Caros alunos iniciamos aqui nosso curso de Cálculo III com o

1.2 Integral Dupla: Domínios Retangulares

Começamos por considerar uma função f definida em um do-

Figura 1.1: Partição de R = [a, b] × [c, d]

j ≤ m, 1 ≤ k ≤ n. A area de cada pequeno retângulo é dada

Pausa para respirar que já vamos definir a integral dupla sobre

Figura 1.2: Soma de Riemann para f (x, y) em R = [a, b] × [c, d]

1.3 Integral Dupla: Domínios Não Retangula-

Para definir a integral dupla de uma função f : D ⊂ R2 7→ R

A integral dupla da função f (x, y) sobre o domínio D ⊂ R2 , deno-

Observem na partição (Fig. 1.3) que apenas os pequenos retân-

1.4 Interpretação Geométrica

Quando a função f (x, y) é positiva na região R, como a da

Figura 1.3: Partição para F (x, y) em R = [a, b] × [c, d]

reto limitado inferiormente por R e superiormente pela superfície

1.5 Integrais Iteradas

Do mesmo modo que para a integral simples, na integral dupla

Figura 1.4: A(x), x fixo, integramos em relação a y

volume do citado prisma.

Figura 1.5: A(y), y fixo, integramos em relação a x

ou seja a ordem em que as integrais simples são executadas não

Propriedade 1.1. Sejam f : D ⊂ R2 7→ R uma função de valores

Propriedade 1.2. Sejam f, g : D ⊂ R2 7→ R duas funções de

Propriedade 1.3. Sejam f : D ⊂ R2 7→ R uma função de valores

Propriedade 1.4. Sejam f, g : D ⊂ R2 7→ R duas funções de va-

Propriedade 1.5. Seja f : D ⊂ R2 7→ R uma função de valores

número finito de curvas em R2 , então vale:

OBS 1.1. As duas primeiras propriedades diz respeito à “line-

1.7 Alguns Exemplos

Nada mais natural que ilustrar um novo conceito com exem-

Vamos diretamente para o primeiro exemplo de integral dupla so-

Exemplo 1.1. Considere a função f : [0, 1] × [0, 1] 7→ R (Fig.

Figura 1.6: Função f : [0, 1] × [0, 1] 7→ R: f (x, y) = exp(−x − y)

OBS 1.2. Daremos aqui um método prático para determinar os

Passo 1 Fazer um desenho da região D. (Fig. 1.7) identificando

Figura 1.7: Determinação prática dos limites para D

direção negativa do eixo x até tocar o ponto mais à esquerda de D

Exemplo 1.2. Considere a função f : [0, 1] × [0, 1] 7→ R (Fig.

Figura 1.8: Função f : [0, 1] × [0, 1] 7→ R: f (x, y) = x.y

Figura 1.9: Limites para o domínio D

A integral passa a ser escrita como: