Apontamentos Análise Matemática-III 2008-2009

Apontamentos de Análise Matemática III
Hermenegildo Borges de Oliveira
Fevereiro de 2009
Conteúdo
1 Funções reais de várias variáveis 1

1.1 Introdução . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.1 Noções algébricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.2 Noções topológicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.1.3 Sucessões de vectores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.4 Funções reais de várias variáveis . . . . . . . . . . . . . . . . . . . . . . . 7
1.1.5 Exercícios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2 Continuidade . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.1 Limites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.2 Continuidade . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.2.3 Exercícios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.3 Derivabilidade . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3.1 Derivadas parciais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3.2 Derivada direccional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.3.3 Derivabilidade . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.3.4 Derivadas de ordem superior . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.3.5 Derivação de funções compostas . . . . . . . . . . . . . . . . . . . . . . . 24
1.3.6 Funções implícitas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.3.7 Aplicações geométricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.3.8 Fórmula de Taylor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
1.3.9 Exercícios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
1.4 Extremos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
1.4.1 Exercícios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2 Cálculo integral vectorial 42

2.1 Integral de Linha . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2.1.1 Primeiras noções . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2.1.2 Motivação . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.1.3 Propriedades gerais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.1.4 Integral de linha independente do caminho . . . . . . . . . . . . . . . . . 48
2.1.5 Exercícios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
2.2 Integral duplo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
2.2.1 Integral de Riemann . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
2.2.2 Integral repetido . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
2.2.3 Mudança de variáveis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
ii
iii
2.2.4 Aplicações . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
2.2.5 Teorema de Green . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
2.2.6 Exercícios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
2.3 Integral de superfície . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.3.1 Superfícies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.3.2 Plano tangente e normal a uma superfície . . . . . . . . . . . . . . . . . 63
2.3.3 Integrais de superfície . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
2.3.4 Outras noções de integral de superfície . . . . . . . . . . . . . . . . . . . 67
2.3.5 Aplicações . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
2.3.6 Exercícios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
2.4 Integral triplo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
2.4.1 Integral de Riemann . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
2.4.2 Integral repetido . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
2.4.3 Mudança de variáveis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
2.4.4 Aplicações . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
2.4.5 Teorema de Green . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
2.4.6 Exercícios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
Bibliografia 82

c Hermenegildo Borges de Oliveira, 2008/2009
Capítulo 1
Funções reais de várias variáveis
Neste primeiro capítulo, iremos dar início ao estudo das funções funções reias de várias variáveis
reias, isto é, funções reais com mais do que uma variável real. A análise que foi feita para
funções reais de apenas uma variável real será naturalmente extendida a dimensões superiores
para podermos trabalhar com funções de várias variáveis. Em função dos programas leccinados
nas disciplinas de Análise Matemática II para os diferentes cursos, parte deste capítulo pode já
ter sido estudada nessa disciplina. Nesse caso, este capítulo é uma introdução para o resto do
texto.
1.1 Introdução
1.1.1 Noções algébricas
Seja N ∈ N arbitrário. Denota-se por RN o conjunto de todas as sucessões finitas de N números
reais que se podem representar pelo N-uplo
x = (x1 , . . . , xN ) .
Em particular, R2 é o conjunto de todos os pares (x, y) de números reais e R3 é o conjunto de

todos os triplos (x, y, z) de números reais. Num sentido mais alargado, designamos os elementos
genéricos de RN por pontos. O real xk , para k = 1, . . . , N, designa-se por k-ésima componente,
ou k-ésima coordenada, do ponto x = (x1 , . . . , xN ). Dados dois elementos x = (x1 , . . . , xN ),
y = (y1 , . . . , yN ) ∈ RN , tem-se x = y se e só se cada componente de x é igual à respectiva
componente de y, isto é,
x = y ⇐⇒ xk = yk ∀ k = 1, . . . , N.
Em RN define-se a adição por
x + y = (x1 , . . . , xN ) + (y1 , . . . , yN ) = (x1 + y1 , . . . , xN + yN )
e a multiplicação por escalar por
αx = α(x1 , . . . , xN ) = (αx1 , . . . , αxN ),
1
2 1. FUNÇÕES REAIS DE VÁRIAS VARIÁVEIS
para quaisquer x, y ∈ RN e α ∈ R. O conjunto RN munido com estas duas operações tem uma
estrutura algébrica de espaço vectorial, ou espaço linear, isto é, RN é um grupo comutativo
para a adição e
RN 6= ∅ ⇐= (0, . . . , 0) ∈ RN ;
(αx + βy) ∈ RN ⇐= x, y ∈ RN , α, β ∈ R.
A noção de que RN é um grupo comutativo para a adição, quer dizer que, para quaisquer x,
y, z ∈ RN :
(i) x + y ∈ RN ;
(ii) x + (y + z) = (x + y) + z;
(iii) x + y = y + x;
(iv) existe o elemento neutro 0 = (0, . . . , 0) e x + 0 = x;
(v) todo x = (x1 , . . . , xN ) ∈ RN tem simétrico − x = (−x1 , . . . , −xN ) ∈ RN .
Neste sentido, os elementos de RN podem ser denominados vectores. Sejam, agora, x1 , x2 , . . . ,

xk vectores de RN . Diz-se que o vector x é uma combinação linear de x1 , x2 , . . . , xk , se
existirem k escalares α1 , α2 , . . . , αk tais que
x = α1 x1 + α2 x2 + · · · + αk xk .
Por outro lado, diz-se que x1 , x2 , . . . , xk ∈ RN são vectores linearmente independentes,

se
α1 x1 + α2 x2 + · · · + αk xk = 0 =⇒ α1 = α2 = · · · = αk = 0.
Denotamos por ek , com k ∈ {1, . . . , N}, o elemento de RN com as componentes todas nulas
excepto a k-ésima componente que é 1. Assim,
e1 = (1, 0, . . . , 0), e2 = (0, 1, 0, . . . , 0), . . . , eN = (0, . . . , 0, 1).
Usando esta notação, podemos escrever, para qualquer vector x = (x1 , . . . , xN ) ∈ RN ,
x = x1 e1 + · · · + xN eN .
Verifica-se que os vectores e1 , e2 , . . . , eN são linearmente independentes. As duas afirmações

anteriores exprimem o facto de e1 , e2 , . . . , eN ser uma base vectorial de RN . Isto significa
que qualquer elemento de RN se escreve de modo único como combinação linear de e1 , e2 , . . . ,
eN .
1.1.2 Noções topológicas

Definição 1.1.1 Sejam x = (x1 , . . . , xN ), y = (y1 , . . . , yN ) ∈ RN quaisquer. Definimos o
produto interno, ou produto escalar, de x por y por
x · y = (x1 , . . . , xN ) · (y1 , . . . , yN ) = x1 y1 + · · · + xN yN .

Podemos definir vários produtos internos em RN . Mas, qualquer definição de produto interno,
terá de satisfazer às propriedades seguintes:
(i) x · y = y · x;
(ii) (x + y) · z = x · z + y) · z e x · (y + z) = x · y + x) · z;
(iii) (αx) · y = x · (αy);
(iv) 0 · 0 = 0 e, para qualquer x =6 0, x · x > 0;
para quaisquer x, y, z ∈ RN e α ∈ R. O espaço vectorial RN munido do produto interno
anterior designa-se por espaço euclidiano. Este produto interno induz a norma seguinte,
denominada norma euclidiana.
Definição 1.1.2 Seja x = (x1 , . . . , xN ) ∈ RN arbitrário. Definimos a norma de x por
p q
kxk = k(x1 , . . . , xN )k = (x1 , . . . , xN ) · (x1 , . . . , xN ) = x21 + · · · + x2N .
Do mesmo que modo para o produto interno, podemos definir várias normas em RN . No
entanto, qualquer definição de norma, terá de satisfazer às propriedades seguintes:
(i) kxk > 0 se x 6= 0 e kxk = 0 se e só se x = 0;
(ii) kαxk = |α|kxk;
(iii) kx + yk ≤ kxk + kyk (desigualdade triangular);
para quaisquer x, y ∈ RN e qualquer α ∈ R. A propriedade (iii) acima é, ainda, válida para
um número finito de parcelas:
kx1 + · · · + xk k ≤ kx1 k + · · · + kxk k.
Deste modo, RN chama-se um espaço vectorial normado e pode-se provar a propriedade seguinte,
denominada desigualdade de Cauchy-Schwarz,
|x · y| ≤ kxk kyk,
válida para quaisquer x, y ∈ RN . Recordando que os elementos de RN são vectores, definimos
a projecção do vector x sobre o vector y como sendo o vector
x·y
p = p y, com p = .
kyk2
No caso particular de y = ek , a projecção de x sobre ek é p = xk ek . Define-se o ângulo entre
dois vectores não nulos x e y por θ, onde
x·y
cos(θ) = , θ ∈ [0, π].
kxk kyk
Sai desta definição que x e y são vectores perpendiculares, ou ortogonais, se
x · y = 0.
Para os elementos da base {e1 , e2 , . . . , eN } tem-se
1 , se i = j

ei · ej = , kei k = 1 ,
0 , se i 6= j,
e, assim, dizemos que {e1 , e2 , . . . , eN } é uma base ortonormada de RN .

Definição 1.1.3 Sejam x = (x1 , . . . , xN ) e y = (y1 , . . . , yN ) dois quaisquer elementos de RN .

Chama-se distância de x a y ao real (não negativo) seguinte:
p
d(x, y) = kx − yk = (x1 − y1 )2 + · · · + (xN − yN )2 .
Qualquer distância d que se possa definir em RN , tem de satisfazer às propriedades seguintes:
(i) d(x, y) ≥ 0 e d(x, y) = 0 se e só se x = y;
(ii) d(x, y) = d(y, x) (propriedade comutativa);
(iii) d(x, z) ≤ d(x, y) + d(y, z) (desigualdade triangular);
para quaisquer x, y, z ∈ RN .
Prenda-mo-nos agora com algumas noções topológicas propriamente ditas. Sejam a ∈ RN
e r > 0. Definem-se a bola aberta de centro a e raio r por
B(a, r) = {x ∈ RN : d(x, a) < r} = {x ∈ RN : kx − ak < r}
e a bola fechada de centro a e raio r por
B(a, r) = {x ∈ RN : d(x, a) ≤ r} = {x ∈ RN : kx − ak ≤ r} .
Na definição seguinte enumeramos as noções topológicas mais importantes.
Definição 1.1.4 Sejam A um subconjunto de RN , com N ∈ N, e a ∈ A.
1. Diz-se que a é um ponto interior do conjunto A, se existir uma bola aberta B(a, r) to-
talmente contida em A. O conjunto de todos os pontos interiores a A designa-se interior
do conjunto A e denota-se por intA.
2. Diz-se que a é um ponto exterior ao conjunto A, se existir uma bola aberta B(a, r) que
não intersecte o conjunto A. O conjunto de todos os pontos exteriores a A designa-se
exterior do conjunto A e denota-se por extA.
3. Diz-se que a é um ponto fronteiro ao conjunto A, se toda a bola aberta B(a, r) intersecta
o conjunto A bem como o seu conjunto complementar, RN \ A. O conjunto de todos os
pontos fronteiros a A designa-se fronteira do conjunto A e denota-se por frA.
4. Diz-se que a é um ponto aderente ao conjunto A, se for um ponto interior ou fronteiro
ao conjunto A. O conjunto de todos os pontos aderentes de A designa-se aderência, ou
fecho, do conjunto A, denota-se por adA e tem-se adA = intA ∪ frA.
5. Diz-se que a é um ponto de acumulação do conjunto A, se toda a bola aberta B(a, r)
contém, pelo menos, um ponto de A distinto de a. O conjunto de todos os pontos de
acumulação de A designa-se derivado do conjunto A e denota-se por A′ .
6. Diz-se que a é um ponto isolado do conjunto A, se existir uma bola aberta B(a, r) que
não intersecte nenhum outro ponto de A distinto de a.
Um subconjunto A de RN diz-se um conjunto aberto, se
A = intA ,
e diz-se um conjunto fechado, se
A = adA .

1.1.3 Sucessões de vectores

Nesta secção introduzimos as sucessões numéricas de termos vectoriais. Usamos aqui a notação
RM em vez de RN para, como vamos ver a seguir, não confundir o índice n habitual na notação
das sucessões com a dimensão de espaço euclidiano.
Definição 1.1.5 Uma sucessão numérica infinita de termos vectoriais é uma função de
variável natural e com valores vectoriais, digamos em RM . Usando a escrita habitual para as
funções, uma sucessão, digamos f, escreve-se da forma seguinte:
f : N −→ RM
n 7→ (f1 (n), . . . , fM (n)).
Por simplicidade de escrita, designamos apenas por sucessão vectorial qualquer sucessão numérica
infinita de termos vectoriais. Para distinguir das sucessões vectoriais, as sucessões numéri-
cas de termos reais serão referidas como sucessões escalares. As considerações feitas para as
sucessões escalares são também válidas para as sucessões vectoriais. Assim, os vectores f(1) =
(f1 (1), . . . , fM (1)), f(2) = (f1 (2), . . . , fM (2)), . . . , f(n) = (f1 (n), . . . , fM (n)), . . . denominam-se
termos da sucessão: primeiro termo, segundo termo, . . . , n-ésimo termo, . . . . O contra-
domínio da função f denomina-se conjunto dos termos da sucessão. Usando a notação
habitual de letras indexadas nos números naturais, denotamos os termos da sucessão acima por
u1 , u2 , . . . , un , . . . , cuja escrita na forma vectorial é a seguinte:
u1 = (u11 , u12 , . . . , u1M ) , u2 = (u21 , u22 , . . . , u2M ) , . . . , un = (un1 , un2 , . . . , unM ), . . . ;
onde unk designa a n-ésima componente do vector uk . Deste modo, verifica-se que cada sucessão
em RM de termo geral un determina M sucessões escalares de termos gerais
u nk , k = 1, . . . , M.
As sucessões unk designam-se sucessões componentes ou sucessões coordenadas da sucessão
un . Por exemplo, un3 é a sucessão terceira componente, ou a sucessão terceira coordenada, da
sucessão unM .
Exercício exemplo 1.1.1 Indique as sucessões componentes da sucessão vectorial seguinte:

1 n
un = , 2 , ln(n + 1) .
n
As operações algébricas de soma e subtracção de sucessões vectoriais un e vn definem-se,
como habitualmente, por:
(u + v)n = un + vn e (u − v)n = un − vn .
A multiplicação de uma sucessão vectorial un por um escalar α define-se por:
(α u)n = α un .
Muitas das definições e resultados sobre sucessões escalares podem ser estendidos às sucessões
vectoriais. No entanto, aqui, iremos apenas prender-nos com as noções e alguns resultados
sobre sucessão limitada e sucessão convergente.

Definição 1.1.6 Uma sucessão vectorial un é limitada, se

∃ C > 0 : kun k ≤ C.
Exercício exemplo 1.1.2 Verifique que a sucessão seguinte é limitada:

1 −n n
un = , e , (−1) .
n
Proposição 1.1.1 Uma sucessão vectorial un em RM , un = (un1 , . . . , unM ), é limitada se e só
se cada uma das suas sucessões componentes unk , com k = 1, . . . , M, é uma sucessão limitada.
SEM DEMONSTRAÇÃO: Ver, por exemplo, Campos Ferreira p. 30.
Exercício exemplo 1.1.3 Usando a proposição anterior, verifique que a sucessão seguinte
não é limitada:
√ ln(n + 1)

un = n, .
n
Definição 1.1.7 Uma sucessão vectorial un de termos em RM diz-se convergente para a ∈
RM , se
∀ ε > 0 ∃ p = p(ε) ∈ N : ( ∀ n ∈ N e n > p ) ⇒ kun − ak < ε.
O vector a da definição anterior chama-se limite da sucessão e, tal como para as sucessões
escalares, escrevemos
lim un = a.
n→+∞
Exercício exemplo 1.1.4 Usando a definição, mostre que

n−1

1
lim , 0, n = (1, 0, 0).
n→+∞ n 3
Proposição 1.1.2 Uma sucessão vectorial un em RM , un = (un1 , . . . , unM ), é convergente para
a ∈ RM , a = (a1 , . . . , aM ), se e só se cada uma das suas sucessões componentes unk converge
para ak , com k = 1, . . . , M.
Desta proposição decorre que

lim un = a ⇐⇒ lim un1 = a1 , lim un2 = a2 , . . . , lim unM = aM .
n→+∞ n→+∞ n→+∞ n→+∞
Exercício exemplo 1.1.5 Usando a proposição anterior, verifique que a sucessão seguinte é
divergente: n
−n 1 n
un = e , 1 + , (−1) .
n
Proposição 1.1.3 Sejam un e vn sucessões vectoriais convergentes em RM e αn uma sucessão
escalar convergente em R. Então as sucessões vectoriais un + vn , un − vn e αn un são conver-
gentes em RM e as sucessões escalares un · vn e kun k são convergentes em R.

1.1.4 Funções reais de várias variáveis

Este capítulo centra-se no estudo de funções do tipo
f : RN −→ R ,
com N ∈ N, isto é, funções cujos objectos têm várias variáveis e as imagens estão em R. Por
isto, estas funções recebem, habitualmente, o nome de campos escalares para as distinguir
dos campos vectoriais, cujas imagens estão em RM , com M ∈ N não necessariamente igual a
N. O domínio da função f , que iremos denotar por Df , é o subconjunto de RN onde a função
está definida. O contra-domínio é o subconjunto de R onde a função toma valores. O gráfico
da função f acima é o subconjunto de RN +1 seguinte
Graf (f ) = {(x, x) ∈ RN +1 : x = f (x), x = (x1 , . . . , xN ) ∈ Ω ⊆ RN } .
Como se depreende, é apenas possível esboçar de forma realística os domínios de funções reais de
duas ou três variáveis. O domínio de uma função de duas variáveis é esboçado no plano gerado
pelo sistema de eixos cartesianos das abcissas e das ordenadas, habitualmente denominados eixo
dos xx e eixo dos yy, respectivamente. O domínio de uma função de três variáveis é esboçado
no espaço gerado pelo sistema de eixos cartesianos das abcissas, das ordenadas e das cotas,
habitualmente denominados eixo dos xx, eixo dos yy e eixo dos zz, respectivamente.
Exercício exemplo 1.1.6 Determine e esboce os domínios das funções seguintes:
p
f (x, y) = x2 − y 2 e g(x, y, z) = ln[1 − (x2 + y 2 + z 2 )] .
Do mesmo modo, é apenas possível esboçar, com acuidade geométrica realística, os gráficos de
funções reais de duas variáveis apenas. O gráfico de uma função de duas variáveis esboça-se
num sistemas de três eixos cartesianos, fazendo corresponder cada ponto do domínio, desenhado
no plano gerado pelos eixos dos xx e dos yy, a um único ponto no eixo dos zz.
Exercício exemplo 1.1.7 Esboce os gráficos das funções seguintes:
p
2 2 1 − (x2 + y 2) x2 + y 2 ≤ 1
f (x, y) = x + y e g(x, y) =
0 x2 + y 2 > 1 .
Assim, iremos estudar, apenas, funções reais de duas ou, quanto muito, três variáveis, uma
vez que para mais do que três variáveis as interpretações geométricas já são mais complicadas.
Diga-se, em abono da verdade, que esta opção será, na maioria das situações, para simplificar
os cálculos.
Definição 1.1.8 Sejam f uma função real de N ∈ N variáveis reais e c uma constante real.
Chama-se conjunto de nível c da função f ao subconjunto de RN onde f atinge o valor c,
isto é,
Nc (f ) = {x ∈ RN : f (x) = c} ,
sendo c denominado nível da função.
Em particular, no caso de funções de duas variáveis, o conjunto de nível designa-se por linha
de nível e, no caso de funções de três variáveis, por superfície de nível.
Exercício exemplo 1.1.8 Determine e esboce os conjuntos de nível c = −9, c = −4, c = −1,
c = 0, c = 1, c = 4 e c = 9 da função f (x, y) = 4y 2 − x2 . A partir das linhas de nível, esboce o
gráfico da função.

1.1. INTRODUÇÃO 8
1.1.5 Exercícios
1. Considere a função seguinte:
x2 − y 3
f (x, y) = .
2xy
Determine:

1 1 1
a) f (y, x) ; b) f (−x, −y) ; c) f , ; d) .
x y f (x, y)
2. Determine o domínio de cada uma das funções seguintes e represente-o graficamente:

1 1 y
a) f (x, y) = √ +√ ; b) f (x, y) = arcsen ;
x+y x−y x

√ x−y
c) f (x, y) = xsen y ; d) f (x, y) = tg ;
x2 + y 2
√ √ √
e) f (x, y, z) = ln(x2 + y 2 − z) ; f) f (x, y, z) = x + y + z .
3. Para cada um dos conjuntos representados na questão anterior, indique:

a) o conjunto dos pontos interiores;
b) o conjunto dos pontos exteriores;
c) o conjunto dos pontos fronteiros;
d) o conjunto dos pontos de acumulação;
e) a aderência.
4. Represente graficamente as linhas ou superfícies de nível correspondentes aos valores c

indicados das funções seguintes:
a) f (x, y) = x2 + y 2 , para c = 0, 1, 4, 9;
b) f (x, y) = exy , para c = e−2 , e−1 , 1, e, e2 ;
c) f (x, y, z) = x + y + z, para c = −1, c = 0, c = 1; √
d) f (x, y, z) = sen(x + y + z ),
2 2 2
para c = −1, − 12 , 0, 22 , 1.
5. Determine os conjuntos de nível das funções seguintes:
a) f (x, y) = x2 + 2y 2 ; b) f (x, y) = x2 − y 2 ;
c) f (x, y, z) = 2x2 + y 2 + 3z 2 ; d) f (x, y, z) = x2 + 2y 2 .
6. A partir das linhas de nível, esboce os gráficos das funções seguintes:
a) f (x, y) = x2 + 2y 2 ; b) f (x, y) = x2 − y 2 ;
p p
c) f (x, y) = 1 − x2 + y 2 ; d) f (x, y) = 3 x2 + y 2 .
1.2. CONTINUIDADE 9
1.2 Continuidade
1.2.1 Limites
A noção de limite para funções de várias variáveis é em tudo idêntica à mesma para funções de
apenas uma variável. A principal diferença é que o módulo dos objectos é agora substituído pela
norma. Por outro lado, ainda relativamente aos objectos, os intervalos abertos são substituídos
por discos abertos em dimensão N = 2 e por esferas abertas em dimensão N = 3. A definição
seguinte é habitualmente deniminada por definição de limite segundo Cauchy.
Definição 1.2.1 Sejam f : RN −→ R uma função escalar e a um ponto de acumulação do

domínio Df . Diz-se que um ponto b ∈ R é o limite da função f no ponto a, ou quando x
tende para a, se
∀ ε > 0 ∃ δ > 0 : ( ∀ x ∈ Df e kx − ak < δ) =⇒ |f (x) − b| < ε .
No caso de existir tal b ∈ R, escrevemos
lim = b .
x−→a
Observemos que a tem de ser um ponto de acumulação do domínio Df , pois, caso contrário,
existiria sempre um δ > 0 para o qual
(B(a, δ) \ {a}) = ∅
e a implicação
( ∀ x ∈ Df e kx − ak < δ) =⇒ |f (x) − b| < ε
seria verdadeira para qualquer b ∈ R.
Para o cálculo de limites de funções escalares, procede-se como no caso de funções de uma
variável apenas. Começamos por tentar determinar o valor da função no ponto onde se pretende
calcular o limite. Se a função estiver definida nesse ponto, então o limite será o valor que a
função aí toma. Mesmo que a função não esteja definida nesse ponto, podemos determinar o
limite em R = [−∞, +∞], desde que não se obtenha uma indeterminação.
Exercício exemplo 1.2.1 Calcule o limite seguinte:

sen(x + y)
lim .
(x,y)→(0,π) x+y
Os casos em que se obtêm indeterminações são mais delicados, pois, tal como para funções de
uma só variável, temos de levantar essas indeterminações. Recordemos que as indeterminações
principais existentes em R = [−∞, +∞] são:
∞−∞; 0×∞; 1∞ ; ∞0 .
A análise da existência ou não destes limites, consiste em reduzir esse estudo a limites de funções
de uma só variável. Nos casos mais simples, podemos usar propriedades da função em estudo
para tornar o cálculo do limite mais simples.
2.2 CONTINUIDADE
c Hermenegildo Borges de Oliveira 10
Exercício exemplo 1.2.2 Calcule o limite seguinte:

x2
x+y
1
lim 1+ .
(x,y)→(∞,0) x
Nas situações mais complicadas, não é possível levantar a indeterminação e proceder como
anteriormente. Nestes casos, temos de analisar a forma como nos podemos aproximar do ponto
em que se pretende calcular o limite para obtermos alguma informação sobre a existência ou
não de limite. Se a função em estudo tem limite nesse ponto, para provar que o limite existe
e é determinada quantidade, temos de usar a definição. A informação de que o limite vai ser
determinada quantidade, é nos dada por um dos dois processos seguintes:
Limites direccionais - este processo consiste em calcular os denominados limites direccionais

de uma função escalar, isto é, em calcular os limites segundo as várias direcções possíveis de
aproximação ao ponto onde se pretende calcular o limite.
Exercício exemplo 1.2.3 Calcule os limites direccionais da função seguinte
2x3 − y 3
f (x, y) =
x2 + y 2
no ponto (x, y) = (0, 0).
Limites repetidos - este processo consiste em reduzir o cálculo do limite de uma função real
de N variáveis ao cálculo de N limites sucessivos (ou repetidos). Na verdade, este processo é
um caso particular do anterior, pois calculamos, apenas, os limites segundo as direcções dos
eixos das coordenadas.
Exercício exemplo 1.2.4 Calcule os limites repetidos da função seguinte
2x3 − y 3
f (x, y) =
x2 + y 2
no ponto (x, y) = (0, 0).
No entanto, convém frisar, que mesmo no caso dos limites anteriores existirem e serem iguais
a determinada quantidade, não podemos concluir que o limite existe e que é igual a essa
quantidade. Isto decorre do facto de num ponto em RN , com N ≥ 2, passarem infinitas
direcções, incluindo as dos eixos de coordenadas. E, por conseguinte, para uma quantidade
finita de direcções o limite pode existir e ser o mesmo, mas ficar sempre a possibilidade de
haver uma direcção onde o limite venha a não existir ou ser diferente. Nestes casos, o único
processo que nos permite concluir que tal limite existe, é usar a definição.
2x3 − y 3
lim = 0.
(x,y)→(0,0) x2 + y 2
2.2 CONTINUIDADE
De ora em diante e apenas por simplicidade de apresentação e de cálculo, todos os exemplos que
vamos considerar são de limites na origem. Para estudar a existência de limites noutros pontos,
procede-se de modo inteiramente análogo. Por outro lado, sempre que possível, podemos fazer
uma mudança de variável para recuperar o caso na origem. Convém, no entanto, realçar o caso
em que o ponto onde se pretende calcular o limite tem todas as componentes infinitas. Neste
caso, a definição de limite (Definição 1.2.1) sofre uma ligeira adaptação:

1
∀ ε > 0 ∃ δ > 0 : ∀ x ∈ Df e kxk > =⇒ |f (x) − b| < ε .
δ
Esta definição ainda é válida quando apenas uma das componentes do ponto onde se pretende
calcular o limite é infinito.
x+y
lim = 0.
(x,y)→(+∞,+∞) x2 + y 2
Muitas das propriedades dos limites de funções reais de variável real permanecem válidas para
os limites de funções escalares.
Proposição 1.2.1 O limite de uma função f : RN −→ R num ponto de acumulação a do
domínio Df quando existe, é único.
DEMONSTRAÇÃO - EXERCÍCIO: Análoga à demonstração para funções reais de variável real.
Este resultado de unicidade do limite é muito útil na prática. De facto, a proposição anterior
obriga a que, quando o limite exista, têm de ser iguais todos os limites direccionais e, por
maioria de razão, os limites repetidos. Assim, se verificarmos que para determinada função não
coincidem dois dos seus limites direccionais, podemos automaticamente concluir que a função
não tem limite no ponto em consideração.
Exercício exemplo 1.2.7 Mostre que a função seguinte não tem limite no ponto (0, 0):
x2 − y 2
f (x, y) = .
x2 + y 2
São, ainda, válidas as propriedades de adição, subtracção, multiplicação, divisão e multiplicação
por escalar de limites, tal como no caso de funções reais de variáveis reais. A proposição seguinte
dá-nos uma definição equivalente de limite, conhecida por definição de limite segundo Heine.
Proposição 1.2.2 Sejam f uma função escalar e a um ponto de acumulação de Df . Tem-se:
lim f (x) = b ⇐⇒ (∀ xn ∈ Df e xn 6= a : xn −→ a =⇒ f (xn ) −→ b) .
x→a
DEMONSTRAÇÃO: AULA TEÓRICA. 1
A proposição anterior tem, também, interesse prático no cálculo de limites. Tal como os limites
direccionais, este resultado é muito útil para mostrar que o limite não existe.
1
Observe-se que o natural n, índice da sucessão vectorial xn , não é necessariamente igual à dimensão do
conjunto de partida da função f .
2.2 CONTINUIDADE
Exercício exemplo 1.2.8 Usando a proposição anterior, considere as sucessões xn e yn a

seguir indicadas para mostrar que o limite seguinte não existe:
x2 − y 2

1 1 1 1
lim , xn = , , yn = √ , √ .
(x,y)→(0,0) x2 + y 2 n 2n n n+1
Outra propriedade dos limites que convém destacar, é a de limite de funções compostas.
Proposição 1.2.3 Sejam f : RN −→ R, g : R −→ R e a um ponto de acumulação de Df

suponhamos que existe o limite limx−→a f (x) = b. Então, se b é um ponto de acumulação de
Dg e se existir o limite limx−→b g(x), também existe o limite da função composta g ◦ f no ponto
ae
lim (g ◦ f )(x) = lim g(x) .
x→a x→b
DEMONSTRAÇÃO: AULA TEÓRICA.
Exercício exemplo 1.2.9 Usando a proposição anterior, mostre que

!
x2 − y 2
lim cos p = 1.
(x,y)→(0,0) x2 + y 2
1.2.2 Continuidade
A noção de continuidade está intimamente ligada à noção de limite. Tal como para os limites,
a continuidade de funções escalares é uma generalização da correspondente noção para funções
reais de uma só variável.
Definição 1.2.2 Sejam f : RN −→ R e a um ponto de acumulação de domínio Df . Diz-se

que a função f é contínua no ponto a, se
lim f (x) = f (a)

x→a
ou, equivalentemente, se
∀ ε > 0 ∃ δ > 0 : ( ∀ x ∈ Df e kx − ak < δ) =⇒ |f (x) − f (a)| < ε .
A Proposição 1.2.2 dá-nos a seguinte condição necessária e suficiente de continuidade.
Proposição 1.2.4 Sejam f uma função escalar e a um ponto de acumulação de Df . Tem-se:
f é contínua em a ⇐⇒ (∀ xn ∈ Df e xn 6= a : xn −→ a =⇒ f (xn ) −→ f (a)) .
DEMONSTRAÇÃO - EXERCÍCIO: Análoga à demonstração da Proposição 1.2.2 2 .
A função f diz-se contínua num conjunto A ⊆ Df , se for contínua em todos os pontos de

A. Diremos que é contínua, sem especificar onde, se for contínua em Df .
2
Ver nota de roda-pé da página 11.
2.2 CONTINUIDADE
As propriedades de continuidade das funções reais de variável real estendem-se naturalmente

às funções escalares. Assim, a soma, a subtracção, o produto e a multiplicação por escalar
de funções contínuas ainda é uma função contínua. A função identidade, as transformações
lineares e os polinómios de N variáveis são funções contínuas. Todas as funções elementares
são funções contínuas nos seus domínios de definição.
De um modo geral, todas as funções que encontramos são contínuas, excepto nos possíveis
pontos onde o denominador se anula ou onde a função não está definida. Na prática, justifica-
se que determinada função é contínua no seu domínio de definição, dizendo que as funções
elementares que a compõem são contínuas.
Por vezes, podemos prolongar por continuidade uma função até um ponto fora do seu domínio.
Mas, como nos diz a definição de limite, esse ponto terá de ser um ponto de acumulação do
domínio.
Exercício exemplo 1.2.10 Estude a função seguinte quanto à continuidade:

( x2 −y2
√ (x, y) 6= (0, 0)
f (x, y) = x2 +y 2
0 (x, y) = (0, 0) .
O exemplo seguinte mostra que nem sempre é possível prolongar, por continuidade, uma função
até um ponto fora do seu domínio.
Exercício exemplo 1.2.11 Verifique se é possível prolongar por continuidade a função seguinte
ao ponto (x, y) = (0, 0):
xy
f (x, y) = 2 .
x + y2
Tal como fizemos para os limites na Proposição 1.2.3, e pelo seu interesse prático, convém
destacar a propriedade seguinte sobre a continuidade de funções compostas.
Proposição 1.2.5 Sejam f : RN −→ R uma função contínua em a ∈ Df , g : R −→ R uma

função contínua em b = f (a). Então a função composta g ◦ f é contínua no ponto a.
DEMONSTRAÇÃO - EXERCÍCIO: Análoga à demonstração da Proposição 1.2.3.
Exercício exemplo 1.2.12 Usando a proposição anterior, estude a função seguinte quanto à
continuidade: 
x 2 −y 2
cos √ 2 2 (x, y) 6= (0, 0)

f (x, y) = x +y
1 (x, y) = (0, 0) .

2.2 CONTINUIDADE
1.2.3 Exercícios
1. Calcule os limites seguintes:
x
sen2 y

− 2 1 2
y x
a) lim lim e x (y−1) ; b) lim lim 1 + ; c) lim lim 1 + ;
x→0 y→2 y→2 x→+∞ x y→0 x→+∞ 3xy 2
 
" #4y  ( x )
 cos x 1
1 + xy − 1
d) lim lim 1 + ; e) lim lim ;
x→+∞ y→+∞ y  y→0 x→+∞ y
 x 
ln(1+y) 2
1+ − 1
( tgy )
 x 1
f) lim lim ; g) lim lim 1 + ln(1 + xy) ;
y→0 x→+∞ y  x→0 y→0 x
 
y2 sen( y1 )
tgx 1
1+ − 1

  ( )
 y
 1 y sen( x1 )
h) lim lim ; i) lim lim 1 + ln 1 + .
x→0 y→+∞
 senx 
 y→0 x→+∞ y x
 
2. Considere a função seguinte:

x−y
f (x, y) = .
x+y
Mostre que
h i
lim lim f (x, y) 6= lim lim f (x, y)
x→0 y→0 y→0 x→0
e conclua que f (x, y) não tem limite no ponto (0, 0).
3. Seja
x2 y 2
f (x, y) = .
x2 y 2 + (x − y)2
Mostre que h i
lim lim f (x, y) = lim lim f (x, y) = 0 ,
x→0 y→0 y→0 x→0
mas f (x, y) não tem limite no ponto (0, 0).
4. Discuta a existência dos limites seguintes, usando a definição para provar os que existem:
xy x+y 1 − cos(x2 + y 2 )
a) lim ; b) lim ; c) lim ;
(x,y)→(0,0) (x2 + y 2 )x2 y 2
p p
(x,y)→(0,0) x2 + y 2 (x,y)→(0,0) x2 + y 2
x2 − y 2 x2 x2 + y 2
d) lim ; e) lim ; f) lim ;
(x,y)→(0,0) x2 + y 2 (x,y)→(0,0) x2 + y 2 − 1 + (x − 1)2 (x,y)→(+∞,+∞) x4 + y 4

1 xy x+y
g) lim x + ysen ; h) lim 2 2
; i) lim .
(x,y)→(0,0) x (x,y)→(0,0) x + y (x,y)→(+∞,+∞) x − xy + y 2
2
2.3 DERIVABILIDADE
5. Calcule os limites seguintes:
− 1 ln(x + ey ) x2 + y 2
a) lim e x2 (y−1)2 ; b) lim ; c) lim .
(x,y)→(+∞,+∞) ex+y
p
(x,y)→(0,1) (x,y)→(1,2) x2 + y 2
6. Usando a composição de limites, calcule os limites seguintes:
sen(x + y) sen(xy) (x2 +y2 )2

a) lim ; b) lim ; c) lim x2 + y 2 .
(x,y)→(+∞,+∞) x+y (x,y)→(0,2) x (x,y)→(0,0)
7. Estude as funções seguintes quanto à continuidade:

1 1
a) f (x, y) = arctg 2
; b) f (x, y) = cos ; c) f (x, y) = ln |1 − xy| ;
x −y xy
( 2xy
xsen y1 se y 6= 0 x2 +y 2
se x2 + y 2 6= 0
d) f (x, y) = e) f (x, y) =
0 se y = 0 ; 0 se x2 + y 2 = 0 .
(
sen xy x2 + y 2 se x2 + y 2 ≤ 1

se (x, y) 6= (0, 0)
f) f (x, y) = g) f (x, y) =
1 se (x, y) = (0, 0) ; 2 − (x2 + y 2) se x2 + y 2 > 1 .
8. Determine os valores de α e β de modo que as funções seguintes sejam contínuas em R2 :

( ( x|y|
sen(x2 +y 2 )
se (x, y) 6= (0, 0) √ se (x, y) 6= (0, 0)
a) f (x, y) = x2 +y 2 b) f (x, y) = x2 +y 2 .
α se (x, y) = (0, 0) ; β se (x, y) = (0, 0)
9. Verifique se é possível prolongar por continuidade as funções seguintes:
y3 x3 − y 3
a) f (x, y) = , (x, y) 6= (0, 0) ; b) f (x, y) = , y 6= x .
x2 + y 2 x−y
1.3 Derivabilidade
1.3.1 Derivadas parciais
Consideremos uma função escalar f : RN −→ R, N natural superior a 1, e seja a ∈ int(Df ),
com a = (a1 , . . . , N). Fixando x2 = a2 , . . . , xN = aN , a função
f (x1 , a2 , . . . , aN )
passa a ser uma função de uma só variável independente x1 . Suponhamos que esta função é
derivável, isto é, existe e é finito o limite seguinte
f (x1 + h, a2 , . . . , aN ) − f (x1 , a2 , . . . , aN )
lim ,
h→0 h
2.3 DERIVABILIDADE
para qualquer x1 , tal que (x1 , . . . , xN ) ∈ Df . Então a derivada (ordinária) de f (x1 , a2 , . . . , aN )

em a1 recebe o nome de derivada parcial de f em relação à variável x1 no ponto
(a1 , . . . , aN ) e denotá-mo-la por um dos símbolos seguintes:
∂f
(a) ; fx′ 1 (a) ; (Dx1 f ) (a) .
∂ x1
Definição 1.3.1 Sejam f : RN −→ R uma função escalar e a ∈ int(Df ), com a = (a1 , . . . , N).
Suponhamos que a função f (a1 , . . . , ak−1 , xk , ak+1, . . . , aN ) é uma função derivável em xk = ak
e tal que (x1 , . . . , xN ) ∈ Df . Então, chama-se derivada parcial de primeira ordem em
relação à variável xk ao seguinte limite
∂f f (a1 , . . . , ak−1 , ak + h, ak+1, . . . , aN ) − f (a1 , . . . , ak−1 , ak , ak+1 , . . . , aN )

(a) = lim .
∂ xk h→0 h
Exercício exemplo 1.3.1 Usando a definição, calcule as duas derivadas parciais de primeira
ordem da função f (x, y) = xy + sen(x + y) no ponto (x, y) = (0, π).
Quando, para determinada função escalar f , existe a derivada parcial fx′ k , k = 1, . . . , N, num
ponto a = (a1 , . . . , ak , . . . , N) ∈ Df , esta coincide com a derivada ordinária da função real de
uma variável real apenas f (a1 , . . . , ak−1 , xk , ak+1 , . . . , aN ) no ponto xk = ak . Assim, podemos
utilizar os formulários de derivação das funções reais, de uma variável real apenas, para o cálculo
das derivadas parciais de funções escalares. Na prática, para calcular a derivada parcial de uma
função escalar f em ordem a determinada variável, consideramos todas as outras variáveis como
constantes e aplicamos as fórmulas de derivação para funções reais, de uma variável real apenas,
em relação a essa variável.
Exercício exemplo 1.3.2 Usando os formulários de derivação de funções reais de uma var-
iável real apenas, calcule as duas derivadas parciais de primeira ordem da função f (x, y) =
xy + sen(x + y) no ponto (x, y) = (0, π).
Interpretação geométrica
Consideremos o caso particular N = 2 de uma função real f de duas variáveis reais (x, y). No
caso de existirem, as duas derivadas parciais de primeira ordem, num ponto a = (a, b) ∈ int(Df ),
são dadas por:
∂f f (a + h, b) − f (a, b)
(a, b) = lim
∂x h→0 h
e
∂f f (a, b + h) − f (a, b)
(a, b) = lim .
∂y h→0 h
O gráfico de f (x, y) é uma superfície que se esboça no sistema de eixos cartesianos xyz. Seccio-
nando a superfície z = f (x, y) pelo plano y = b, obtemos uma função de uma variável apenas
(linha) de equação
z = f (x, b) = ϕ(x) .
2.3 DERIVABILIDADE
Esta linha esboça-se num plano paralelo ao plano definido pelos eixos dos xx e dos yy e tem-se
∂f
(a, b) = ϕ′ (a).
∂x
Do que se conhece para funções de uma variável apenas, fx′ (a, b) pode ser interpretada como o
declive da recta tangente ao gráfico da linha z = ϕ(x) no ponto x = a, isto é,
∂f
(a, b) = tg α ,
∂x
onde α é o ângulo formado pela tangente ao gráfico da linha z = ϕ(x) no ponto x = a e o plano
definido pelos eixos dos xx e dos yy. Analogamente,
∂f
(a, b) = tg β ,
∂y
sendo β o ângulo formado pela tangente ao gráfico da linha z = φ(y) no ponto y = b e o plano
definido pelos eixos dos xx e dos yy. E a linha z = φ(y) obtém-se, seccionando, no mesmo
ponto (a, b), a superfície z = f (x, y) pelo plano x = a,
1.3.2 Derivada direccional

Definição 1.3.2 Sejam f : RN −→ R uma função escalar, a ∈ int(Df ), a = (a1 , . . . , aN ), e u
um vector arbitrário de RN , u = (u1 , . . . , uN ). Chama-se derivada da função f no ponto a
dirigida segundo o vector u ao limite seguinte, quando existe,
f (a + hu) − f (a) f (a1 + h u1 , . . . , aN + h uN ) − f (a1 , . . . , aN )

lim ≡ lim .
h→0 h h→0 h
O limite anterior denota-se por um dos símbolos seguintes:
fu′ (a) ; (Du f ) (a) .
No caso de kuk = 1, a derivada da função f no ponto a dirigida segundo o vector u, chama-se

derivada direccional. No caso particular de u ser um vector unitário, digamos u = ek , a
derivada direccional coincide com a derivada parcial:
f (a + h ek ) − f (a) ∂f
fe′ k (a) = lim = (a) .
h→0 h ∂ xk
Portanto, o conceito de derivada parcial é um caso particular do conceito de derivada dirigida.
O significado geométrico da derivada direccional é análogo ao de derivada parcial. Por exemplo,
para uma função escalar f de duas variáveis reais (x, y), tem-se, para qualquer vector unitário
u ∈ R2 ,
fu′ (a, b) = tg γ ,
onde γ é o ângulo que a tangente ao gráfico de f , no ponto (a, b, f (a, b)) e na direcção de u,
forma com o plano definido pelos eixos dos xx e dos yy.
2.3 DERIVABILIDADE
Exercício exemplo 1.3.3 Calcule as derivadas dirigidas da função f (x, y) = xy + sen(x + y),
segundo os vectores u = (1, 1) e v = (0, 1), no ponto (x, y) = (0, π).
Definição 1.3.3 Sejam f : RN −→ R uma função escalar, a ∈ int(Df ). Suponhamos que
existem as derivadas parciais fx′ 1 (a),. . . ,fx′ N (a). Ao vector de RN

∂f ∂f
(a), . . . , (a)
∂ x1 ∂ xN
chama-se gradiente de f no ponto a e denota-se por um dos símbolos seguintes:
∇f (a) ; gradf (a) .
Se a função f tiver todas as derivadas parciais em todos os pontos interiores ao seu domínio
Df , chama-se função gradiente de f , à função vectorial seguinte:
∇f : int(Df ) −→ RN
(x1 , . . . , xN ) 7→ ∇f (x1 , . . . , xN ) .
Exercício exemplo 1.3.4 Calcule o gradiente da função f (x, y) = xy + sen(x + y) no ponto
(x, y) = (0, π) e indique a direcção de maior crescimento da função no ponto (x, y) = (0, π).
Proposição 1.3.1 Sejam f : RN −→ R uma função escalar, a ∈ int(Df ) e u um vector
arbitrário de RN . No caso de existirem, a derivada dirigida de f no ponto a segundo u e o
gradiente de f no mesmo ponto estão relacionados por:
fu′ (a) = ∇f (a) · u .
Exercício exemplo 1.3.5 Usando a proposição anterior, √ calcule
√ a derivada direccional da
função f (x, y) = xy+sen(x+y), segundo a direcção u = ( 2/2, 2/2), no ponto (x, y) = (0, π).
Interpretação geométrica do gradiente

Usando a proposição anterior e conhecimentos de Geometria Analítica, podemos escrever
fu′ (a) = ∇f (a) · u = k∇f (a)k kuk cos θ , θ = ∠ (u, ∇f (a)) ,
desde que ∇f (a) 6= 0. No caso de kuk = 1,
fu′ (a) = k∇f (a)k cos θ .
Temos cos θ = 1 quando θ = 0, isto é, o maior valor de cos θ é obtido quando seleccionamos u
com a mesma direcção e sentido de ∇f (a). Então, a derivada direccional tem o valor máximo,
fu′ (a) = k∇f (a)k ,
e u é a direcção de maior crescimento da função. Quando θ = π, cos θ = −1 e corresponde
ao caso de u ter a mesma direcção de ∇f (a), mas sentido oposto. Neste caso, a derivada
direccional tem o valor mínimo,
fu′ (a) = −k∇f (a)k ,
e u é a direcção de maior decrescimento da função. Assim, temos as interpretações ge-
ométricas seguintes para o gradiente:
2.3 DERIVABILIDADE
1. O gradiente de uma função escalar, calculado num ponto interior ao seu domínio, indica-
nos a direcção de maior crescimento (gradiente) da função, observada no seu domínio.
2. A norma do gradiente dá-nos a taxa de crescimento da função na sua direcção de maior
crescimento.
Como primeira consequência do exposto, verifica-se que o gradiente de uma função f num ponto
a ∈ Df é perpendicular ao conjunto de nível correspondente a esse ponto:
∇f (a) · Nc f = 0 , onde c = f (a) , ∇f (a) 6= 0.
Outra consequência, é a possibilidade de ver o comportamento gráfico de uma função de duas
variáveis a partir do conhecimento dos gradientes calculados em alguns pontos do domínio da
função.
Definição 1.3.4 Consideremos uma função f : R2 → R de domínio Df . Designa-se por

campo de direcções da função f à correspondência que a cada ponto (a, b) ∈ Df , com
∇f (a, b) 6= (0, 0), associa a direcção das rectas que passam (a, b) e de declive igual a
1
∇f (a, b) .
k∇f (a, b)k
O ponto (a, b) ∈ Df é representado no plano xy e a associação que define o campo de direcções
pode ser representada por uma seta com a direcção e sentido de ∇f (a, b)/k∇f (a, b)k. Como o
declive de qualquer recta no plano xy é dado por
dy
,
dx
em cada ponto (x, y) ∈ Df , o declive da seta que representa o campo de direcções corresponde
ao declive da curva y solução da equação diferencial seguinte:
∂f
dy ∂y
= ∂f
.
dx ∂x
Nos casos mais simples, esta equação diferencial pode ser resolvida separando as variáveis e
depois integrando a equação. Os casos mais complicados remetemos para o Curso de Equações
Diferenciais (ordinárias).
Exercício exemplo 1.3.6 Esboce os campos de direcções da função f (x, y) = 4 − 2x2 − y 2 .
1.3.3 Derivabilidade
Definição 1.3.5 Sejam f : RN → R uma função escalar e a ∈ int(Df ), a = (a1 , . . . , aN ).
Diz-se que f é uma função derivável no ponto a, se existir α = (α1 ,. . . ,αN ) constante tal
que, para qualquer h ∈ RN , h = (h1 , . . . , hN ), se tem
f (a + h) = f (a) + α · h + o(khk) ,
com h independente α e a + h ∈ Df .
2.3 DERIVABILIDADE
A notação o(khk) quer dizer que o(khk) é um infinitésimo quando comparado com khk, isto é,
o(khk)
lim = 0.
khk→0 khk
Na prática, torna-se útil escrever a expressão f (a + h) = f (a) + α · h + o(khk) na forma seguinte

q
f (a1 + h1 , . . . , aN + hN ) = f (a1 , . . . , aN ) + α1 h1 + · · · + αn hN + o h21 + · · · + h2N .
Exercício exemplo 1.3.7 Usando a definição, mostre que a função f (x, y, z) = x2 + 2x − y −

z 2 + 3 é derivável no ponto (0, 0, 0).
De maneira equivalente, diremos que f é derivável no ponto a, se existir α = (α1 ,. . . ,αN )

constante tal que
f (x) = f (a) + α · (x − a) + o(kx − ak) ,
em todo o ponto x ∈ Df . Dizer que f é derivável no ponto a equivale a dizer que o diferencial
f (a + h) − f (a)
pode ser aproximado por uma função linear de h. Deste modo, podemos dizer que f é derivável
no ponto a, se e só se existir uma única aplicação linear La : RN → R tal que
f (a + h) = f (a) + La (h) + o(khk) .
Chama-se derivada de f no ponto a, e denota-se por f ′ (a), à única aplicação linear nas
condições anteriores e, neste caso, podemos escrever
f (a + h) = f (a) + f ′ (a)(h) + o(khk) .
Proposição 1.3.2 Seja f : RN → R uma função derivável num ponto a ∈ int(Df ). Então:
1. A função f é contínua no ponto a;
2. Para todo o vector u ∈ RN , existe a derivada dirigida fu′ (a). Em particular, a função f
admite todas as derivadas parciais de primeira ordem e, tendo em conta as notações da
definição anterior,
∂f ∂f
(a) = α1 , . . . , (a) = αN ,
∂ x1 ∂ xN
e
f (a + h) = f (a) + ∇f (a) · h + o(khk) .
A primeira afirmação da proposição anterior permite-nos concluir que, se f não for contínua no
ponto a, então a função também não é derivável nesse ponto. A segunda afirmação permite-nos
concluir que, se não existir uma derivada dirigida de f no ponto a, então a função também não
é derivável no ponto a.
2.3 DERIVABILIDADE
Exercício exemplo 1.3.8 Usando a resolução do Exercício 1.3.7, indique, sem efectuar mais
cálculos, as derivadas parciais fx′ (0, 0, 0), fy′ (0, 0, 0) e fz′ (0, 0, 0).
A proposição seguinte dá-nos uma condição suficiente de derivabilidade, o que na prática
e para justificar a derivabilidade de uma função escalar, nos permite evitar a complicada
Definição 1.3.5.
Proposição 1.3.3 Sejam f : RN → R uma função escalar e a ∈ RN . Suponhamos que existem
e são contínuas, no ponto a, todas as derivadas parciais de primeira ordem, fx′ 1 , . . . fx′ N . Então
a função f é derivável no ponto a.
Exercício exemplo 1.3.9 Mostre que para a função seguinte existem todas as derivadas par-
ciais no ponto (x, y) = (0, 0), mas a função não é derivável nesse ponto:
(
xy 2
x2 +y 4
, x 6= 0
f (x, y) =
0, x = 0.
As hipóteses da proposição anterior podem ser enfraquecidas. Basta que N − 1 das derivadas
parciais fx′ 1 , . . . fx′ N sejam contínuas no ponto a. Uma função escalar f para a qual existam
e sejam contínuas, num ponto de acumulação do seu domínio, todas as derivadas parciais de
primeira ordem, diz-se função de classe C 1 , ou função continuamente derivável, nesse
ponto. De um modo geral, dado k ∈ N, diz-se que f é uma função de classe C k , ou que é
uma função continuamente derivável até à ordem k, se existirem e forem contínuas todas
as derivadas parciais de ordem k. Designam-se por funções de classe C 0 as funções contínuas
e as função de classe C ∞ , também designadas funções indefinidamente deriváveis, são
aquelas paras as quais existem e são contínuas todas as derivadas parciais de qualquer ordem.
Definição 1.3.6 Sejam f : RN → R uma função escalar, a um ponto de RN e u = (u1, . . . , uN )
um vector de RN . Designa-se por diferencial da função f no ponto a relativo ao vector u ao
escalar seguinte
∂f ∂f
(d f )u (a) = ∇f (a) · u = (a) u1 + · · · + (a) uN .
∂ x1 ∂ xN
Resulta da Proposição 1.3.1 que o diferencial de uma função escalar f relativo a um vector,
não é mais do que a derivada de f segundo esse mesmo vector, isto é,
(d f )u(a) = fu′ (a) .
Se, na definição anterior, escolhermos para vector u o vector h = (h1 , . . . , hN ) da Definição 1.3.5,
o qual podemos escrever como h = (d x1 , . . . , d xN ), e não fizermos referência ao ponto a,
obtemos o diferencial seguinte
∂f ∂f
df = d x1 + · · · + d xN .
∂ x1 ∂ xN
Esta última expressão do diferencial é habitualmente designada por diferencial total da
função f .
2.3 DERIVABILIDADE
Exercício exemplo 1.3.10 Para a função f (x, y) = ln(x2 + y 2), determine o seu diferencial
no ponto (x, y) = (1, 1) relativo ao vector u = (−1, 1). Determine, também, o diferencial total
da função f .
1.3.4 Derivadas de ordem superior

Se, para cada k = 1, . . . , N, a função f da Definição 1.3.1 tiver derivada parcial em ordem a
xk em todos os pontos interiores ao seu domínio Df , chama-se função derivada parcial de
f em ordem a xk à função seguinte:
fx′ k : int(Df ) −→ R
(x1 , . . . , xN ) 7→ fx′ k (x1 , . . . , xN ) .
Assim, as novas funções fx′ k podem admitir, por sua vez, derivadas parciais. Suponhamos que
fx′ k admite derivada parcial em ordem à variável xj , com j = 1, . . . , N, em todos os pontos do
domínio Dfx′ k . Chama-se (função) derivada parcial de segunda ordem de f em relação a
xk , primeiro, e a xj , depois, à derivada parcial de primeira ordem de fx′ k em relação a xj . Esta
derivada denota-se por um dos símbolos seguintes:
∂2 f

∂ ∂f ′
≡ ; fx′′k xj ≡ fx′ k x ; Dxk xj f ≡ Dxj (Dxk f ) .
∂ xj ∂xk ∂ xj ∂ xk j
No caso em que j = k, isto é, xj = xk , designamos a derivada parcial de segunda ordem de f

em relação a xk por derivada quadrada de f em ordem a xk e denota-se por um dos símbolos
seguintes:
∂2 f

∂ ∂f ′
2
≡ ; fx′′2 ≡ fx′ k x ; Dx2k f ≡ Dxk (Dxk f ) .
∂ xk ∂ xk ∂ xk k k
Se j 6= k, isto é, xj 6= xk , designamos a derivada parcial de segunda ordem de f em relação a

xk , primeiro, e a xj , depois, por derivada rectangular ou derivada cruzada de f em relação
a xk e a xj . No caso de existirem as derivadas parciais fx′ k xj em todos os pontos de Dfx′ k xj ,
podemos definir as derivadas parciais de terceira ordem
∂3 f ∂2 f

∂ ′
; fx′′′i xj xk ≡ fx′′k xj

≡ ; D xk xj xi f ≡ D xi D xk xj f .
∂ xi ∂xj ∂xk ∂ xi ∂ xj ∂ xk xi
No caso de algum dos xk , xj , xi ser igual a outro, podemos simplificar a escrita como no caso
das derivadas parciais de segunda ordem. O mesmo raciocínio pode ser aplicado às derivadas
parciais de quarta ordem e assim sucessivamente.
Exercício exemplo 1.3.11 Determine as (funções) derivadas parciais de primeira e segunda

ordem da função f (x, y) = arccos(x/y).
Para cada função escalar f de N variáveis vão, possivelmente, existir N p derivadas parciais de
ordem N. No caso particular de uma função escalar de duas variáveis, verifica-se que, em muitas
situações, as derivadas parciais cruzadas de segunda ordem são iguais, independentemente da
ordem por que variável se deriva.
2.3 DERIVABILIDADE
Proposição 1.3.4 (Schwarz) Sejam f uma função escalar de duas variáveis (x, y) e D um
subconjunto aberto de Df . Suponhamos que existem as funções derivadas parciais fx′ e fy′ em D
e que uma das duas derivadas parciais de segunda ordem, fx′′y ou fy′′x , existe e é contínua em
D. Então a outra derivada parcial de segunda ordem também existe e tem-se:
∂2 f ∂2 f
(a, b) = (a, b) ∀ (a, b) ∈ D .
∂ x ∂y ∂ y ∂x
O resultado anterior pode ser generalizado a qualquer função escalar f definida num domínio
de RN , com N ≥ 2. Por exemplo, prova-se que, para uma função f definida num domínio de
R3 e assumindo que as condições correspondentes da proposição anterior são satisfeitas,
∂4 f ∂4 f
= .
∂ x2 ∂y ∂z ∂ z ∂y ∂x2
Exercício exemplo 1.3.12 Mostre que para a função f (x, y, z) = sen(xyz) se tem sempre
∂3 f ∂3 f
= .
∂ x ∂y ∂z ∂ z ∂x ∂y
Uma outra noção que convém ser estendida a uma ordem superior, é a noção de diferencial.
Nas condições da Definição 1.3.6, sabemos que o diferencial da função f no ponto a, relativo
ao vector u, é o escalar definido por:
∂f ∂f
(d f )u(a) = (a) u1 + · · · + (a) uN .
∂ x1 ∂ xN
Conhecida esta expressão, podemos escrever o diferencial de primeira ordem do modo seguinte:

∂ ∂
(d f )(u) (a) = u1 + · · · + uN f (a) .
∂ x1 ∂ xN
Se existirem todas as derivadas parciais de ordem k > 1 da função f no ponto a, então podemos
definir o diferencial de ordem k de f no ponto a, relativo ao vector u, pela forma simbólica
seguinte:
k
k ∂ ∂
(d f )(u) (a) = u1 + · · · + uN f (a) .
∂ x1 ∂ xN
Por exemplo, no caso particular de k = 2 e de uma função de duas variáveis (x, y), a fórmula
anterior tem o significado seguinte:
∂2f ∂2f 2
2∂ f
(d 2 f )(u,v) (a, b) = u2 (a, b) + 2uv (a, b) + v (a, b) .
∂ x2 ∂ y∂ x ∂ y2
Procedendo de forma análoga e usando as notações da Definição 1.3.5, definimos o diferencial
total de ordem k da função f por:
k
k ∂ ∂
d f= d x1 + · · · + d xN f .
∂ x1 ∂ xN
Exercício exemplo 1.3.13 Determine o diferencial total de segunda ordem da função
2 +y 2 )
f (x, y) = e−(x .
2.3 DERIVABILIDADE
1.3.5 Derivação de funções compostas

Nesta subsecção vamos considerar a composição de um campo escalar f : RN → R com um
campo vectorial g : RM → RN , podendo se ter M = 1. A função composta f ◦ g é um campo
escalar e fica, assim, definida
f ◦ g : RM −→ R .
Caso de apenas uma variável independente - este é o caso mais simples, em que M = 1
e N é um natural qualquer maior do que 1.
Proposição 1.3.5 Sejam ϕ1 : R → R, . . . , ϕN : R → R N funções reais, de uma variável real

apenas, deriváveis em t e seja f : RN → R um campo escalar derivável em x ≡ (x1 , . . . , xN ) =
(ϕ1 (t), . . . , ϕN (t)). Então a função composta
u : R −→ R
t 7→ u(t) = f (ϕ1 (t), . . . , ϕN (t))
é derivável em t e a sua derivada é calculada pela fórmula

du ∂ f d ϕ1 ∂ f d ϕ2 ∂ f d ϕN
u′ (t) ≡ = + +···+ .
dt ∂ x1 d t ∂ x2 d t ∂ xN d t
SEM DEMONSTRAÇÃO: Ver, por exemplo, Dias Agudo p. 129.
A expressão da derivada de u pode ser simplificada, se considerarmos u como a função composta
z = f (x1 , . . . , xN ) , x1 = x1 (t) , . . . , xN = xN (t) .
Por isto e embora sendo um abuso de notação, podemos escrever a derivada de u na forma
dz ∂ z d x1 ∂ z d x2 ∂ z d xN
z ′ (t) ≡ = + +···+ .
dt ∂ x1 d t ∂ x2 d t ∂ xN d t
Usando a notação matricial, podemos ainda escrever a derivada anterior na forma
 d x1 
dt
 d x2 
dz ∂z ∂z ∂z
 dt

z ′ (t) ≡ = ··· .
 
dt ∂ x1 ∂ x2 ∂ xN  ..

 . 

d xN
dt
Exercício exemplo 1.3.14 Determine a derivada z ′ (t) da função seguinte, resultante da com-
posição das funções indicadas:
x
z= , x = e2t , y = ln t .
y
2.3 DERIVABILIDADE
As derivadas de ordem superior são calculadas aplicando a fórmula da Proposição 1.3.5 a cada
função que se obtém por derivação. Por exemplo, para a segunda derivada z ′′ (t) e usando a
notação mais simples z = f (x1 , . . . , xN ), obtemos a fórmula
d2 z
2
∂ z 2 d x2 ∂ z 2 d xN d x1 ∂ z d2 x1

d dz ∂ z d x1
≡ = + +···+ + +
d t2 dt dt ∂ x21 d t ∂ x2 x1 d t ∂ xn x1 d t dt ∂ x1 d t2
2
∂ z d x1 ∂ z 2 d x2 ∂ z 2 d xN d x2 ∂ z d2 x2

+ + · · · + + +
∂ x1 x2 d t ∂ x22 d t ∂ xn x2 d t dt ∂ x2 d t2
·· · · · · · · · · · · +
∂ 2 z d x1 ∂ z 2 d x2 ∂ z 2 d xN d xN ∂ z d2 xN

+ +···+ + .
∂ x1 xN d t ∂ x2 xN d t ∂ x2N d t dt ∂ xN d t2
Exercício exemplo 1.3.15 Determine a derivada z ′′ (t) da função resultante da composição

das funções indicadas no Exercício 1.3.14.
Caso de várias variáveis independentes - este é o caso geral para N e M naturais, su-
postamente maiores do que 1.
Proposição 1.3.6 Sejam ϕ1 : RM → R, . . . , ϕN : RM → R N funções escalares deriváveis

em t = (t1 , . . . , tM ) e seja f : RN → R um campo escalar derivável em
x = (ϕ1 (t), . . . , ϕN (t)) = (ϕ1 (t1 , . . . , tM ), . . . , ϕN (t1 , . . . , tM )) .
Então a função composta
u : RM −→ R
t 7→ u(t) = f (ϕ1 (t), . . . , ϕN (t))
onde, não é demais recordar,
t = (t1 , . . . , tM ) e u(t) = f (ϕ1 (t), . . . , ϕN (t)) = f (ϕ1 (t1 , . . . , tM ), . . . , ϕN (t1 , . . . , tM )) ,
é derivável em t = (t1 , . . . , tM ) e as suas derivadas parciais são calculadas pelas fórmulas
∂u ∂ f ∂ ϕ1 ∂ f ∂ ϕ2 ∂ f ∂ ϕN
= + +···+ ,
∂ t1 ∂ x1 ∂ t1 ∂ x2 ∂ t1 ∂ xN ∂ t1
∂u ∂ f ∂ ϕ1 ∂ f ∂ ϕ2 ∂ f ∂ ϕN
= + +···+ ,
∂ t2 ∂ x1 ∂ t2 ∂ x2 ∂ t2 ∂ xN ∂ t2
..
.
∂u ∂ f ∂ ϕ1 ∂ f ∂ ϕ2 ∂ f ∂ ϕN
= + +···+ .
∂ tM ∂ x1 ∂ tM ∂ x2 ∂ tM ∂ xN ∂ tM
2.3 DERIVABILIDADE
As expressões das derivadas parciais de u também podem ser simplificadas, neste caso, se
considerarmos u como a função composta
z = f (x1 , . . . , xN ) , x1 = x1 (t1 , . . . , tM ) , . . . , xN = xN (t1 , . . . , tM ) .
Novamente e tendo presente que se trata de um abuso de notação, podemos escrever as derivadas
parciais de u na forma
∂z ∂ z ∂ x1 ∂ z ∂ x2 ∂ z ∂ xN
= + +··· ,
∂ t1 ∂ x1 ∂ t1 ∂ x2 ∂ t1 ∂ xN ∂ t1
∂z ∂ z ∂ x1 ∂ z ∂ x2 ∂ z ∂ xN
= + +··· ,
∂ t2 ∂ x1 ∂ t2 ∂ x2 ∂ t2 ∂ xN ∂ t2
..
.
∂z ∂ z ∂ x1 ∂ z ∂ x2 ∂ z ∂ xN
= + +··· .
∂ tM ∂ x1 ∂ tM ∂ x2 ∂ tM ∂ xN ∂ tM
Usando a notação matricial, podemos também escrever a derivadas parciais anteriores na forma
 ∂ x1 ∂ x1
· · · ∂∂ txM1

∂ t1 ∂ t2
 ∂ x2 ∂ x2 ∂ x2 
 ∂t
1 ∂ t 2
· · · ∂ tM 

∂z ∂z ∂z
∂z ∂z ∂z

· · · ∂ tM = ∂ x1 ∂ x2 · · · ∂ xN  . .. .. .
.

∂ t1 ∂ t2
 . . . 
 
∂ xN ∂ xN ∂ xN
∂ t1 ∂ t2
· · · ∂ tM
Exercício exemplo 1.3.16 Determine as derivadas parciais de primeira ordem zx e zy da

função seguinte, resultante da composição das funções indicadas:
u
z = arctg , u = sen(xy) , v = cos(xy) .
v
As derivadas parciais de ordem superior são calculadas aplicando a fórmula da Proposição 1.3.6
a cada função que se obtém por derivação parcial. Por exemplo, para a segunda derivada parcial
cruzada em ordem a t1 e a t2 e usando a notação mais simples z = f (x1 , . . . , xN ), obtemos a
fórmula
∂2 z
2
∂ z 2 ∂ x2 ∂ z 2 ∂ xN ∂ x1 ∂ z ∂ 2 x1

∂ ∂z ∂ z ∂ x1
≡ = + + · · · + + +
∂ t2 t1 ∂ t2 ∂ t1 ∂ x21 ∂ t2 ∂ x2 x1 ∂ t2 ∂ xn x1 d t2 d t1 ∂ x1 ∂ t2 t1
2
∂ z ∂ x1 ∂ z 2 ∂ x2 ∂ z 2 ∂ xN ∂ x2 ∂ z ∂ 2 x2

+ + · · · + + +
∂ x1 x2 ∂ t2 ∂ x22 ∂ t2 ∂ xn x2 d t2 d t2 ∂ x2 ∂ t22
·· · · · · · · · · · · +
∂ 2 z ∂ x1 ∂ z 2 ∂ x2 ∂ z 2 ∂ xN ∂ xN ∂ z ∂ 2 xN

+ +···+ + .
∂ x1 xN ∂ t2 ∂ x2 xN ∂ t2 ∂ x2N d t2 d t1 ∂ xN ∂ t2 t1
Exercício exemplo 1.3.17 Determine as derivadas parciais de segunda ordem zxx , zyy e zxy
da função resultante da composição das funções indicadas no Exercício 1.3.16.
2.3 DERIVABILIDADE
1.3.6 Funções implícitas

Seja y uma função real de N variáveis reais (x1 , . . . , xN ) e suponhamos que a equação
F (x1 , . . . , xN , y) = 0
não é resolúvel em ordem a y. Então dizemos que a equação F (x1 , . . . , xN , y) = 0 define, de
forma implícita, y como função de (x1 , . . . , xN ). Um caso particular, é o da equação
F (x, y, z) = 0
definir implicitamente z como função de (x, y). Poder-se-á dar o caso da equação F (x, y, z) = 0
definir x ou y como funções implícitas de (y, z) ou (x, z), respectivamente.
Proposição 1.3.7 Seja

F : RN × R −→ R
(x1 , . . . , xN , y) 7→ F (x1 , . . . , xN , y)
uma função de classe C 1 num subconjunto aberto D ⊂ RN × R, com N ≥ 1. Suponhamos que

as seguintes condições são satisfeitas:
(i) (a1 , . . . , aN , b) é uma solução da equação F (x1 , . . . , xN , y) = 0 ;
(ii) Fy (a1 , . . . , aN , b) 6= 0 .
Então existem uma bola aberta B(a0 , ε) ⊂ RN , com a0 = (a1 , . . . , aN ), e uma função
ϕ: B(a0 , ε) −→ R
(x1 , . . . , xN ) 7→ ϕ(x1 , . . . , xN ) ,
tal que ϕ(a1 , . . . , aN ) = b, que define y como função implícita de (x1 , . . . , xN ):
y = ϕ(x1 , . . . , xN ) .
SEM DEMONSTRAÇÃO: Ver, por exemplo, Campos Ferreira, p. 130.
No caso particular da equação f (x, y) = 0 definir implicitamente y como uma função derivável
de x, y = ϕ(x), tem-se
dy f ′ (x, y)
y′ ≡ = − x′ .
dx fy (x, y)
Exercício exemplo 1.3.18 Verifique que a equação
π
5xey − sen (x + y) + x − 5 = 0
2
define implicitamente y como uma função derivável de x no ponto (x, y) = (1, 0) e calcule y ′(1).
Se a equação F (x1 , . . . , xN , y) = 0 definir implicitamente y como uma função de classe C 1 de
(x1 , . . . , xN ), y = ϕ(x1 , . . . , xN ), tem-se
∂y Fx′ (x1 , . . . , xN , y)
= − ′k
∂ xk Fy (x1 , . . . , xN , y)
para qualquer variável xk , com 1 ≤ k ≤ N.
2.3 DERIVABILIDADE
Exercício exemplo 1.3.19 Considere a equação
x2 y 2 + xz 3 − 2x4 + zy 3 = 6 .
1. Mostre que esta equação define implicitamente, tanto z como uma função derivável de
(x, y), y como uma função derivável de (x, z) ou x como uma função derivável de (y, z),
todas no ponto (x, y, z) = (1, −1, 2).
2. Calcule zx , yz e xy no ponto (x, y, z) = (1, −1, 2).
3. Verifique que, neste caso, se tem

∂x ∂y ∂z
= −1 .
∂y ∂z ∂x
1.3.7 Aplicações geométricas

As aplicações geométricas a que nos iremos referir dizem respeito, apenas, aos casos das dimen-
sões N = 2 ou N = 3. Comecemos por considerar um campo escalar
F : R3 −→ R
(x, y, z) 7→ F (x, y, z) .
A equação
F (x, y, z) = 0
representa, num sistema de eixos cartesianos (x, y, z) e de forma implícita, uma superfície em
R3 , que denotamos por S. Dizemos que um ponto P = (a, b, c) pertencente à superfície S
é ponto regular, se existem as derivadas parciais Fx , Fy e Fz , e são contínuas em P , e se,
além do mais, pelo menos, uma das derivadas parciais, Fx (P ), Fy (P ) ou Fz (P ) é não nula. Se
todas as derivadas parciais Fx (P ), Fy (P ) e Fz (P ) forem nulas, ou não existir, no ponto P , pelo
menos, uma das derivadas parciais, dizemos que P = (a, b, c) é um ponto crítico, ou ponto
singular, da superfície S.
Definição 1.3.7 Seja S uma superfície em R3 e P = (a, b, c) um ponto regular de S. Chama-

se plano tangente à superfície S no ponto P ao plano formado por todas as rectas tangentes
à superfície S no ponto P .
O ponto P = (a, b, c) da definição anterior é denominado ponto de contacto do plano tangente

com a superfície S.
Proposição 1.3.8 Sejam S uma superfície em R3 , representada de forma implícita pela equação
F (x, y, z) = 0, e P = (a, b, c) um ponto regular de S. Então:
(i) o vector ∇F (P ) é perpendicular ao plano tangente à superfície S no ponto P .
(ii) a equação cartesiana do plano tangente à superfície S no ponto P é dada por
∂F ∂F ∂F
(P )(x − a) + (P )(y − b) + (P )(z − c) = 0 .
∂x ∂y ∂z
2.3 DERIVABILIDADE
No caso da equação da superfície F (x, y, z) = 0 puder ser escrita na forma explícita

z = f (x, y) ,
a equação do plano tangente simplifica-se consideravelmente. Neste caso, a equação carte-
siana do plano tangente à superfície S no ponto P = (a, b, c), onde c = f (a, b), é escrita
na forma seguinte
∂f ∂f
z−c = (a, b)(x − a) + (a, b)(y − b) .
∂x ∂y
Exercício exemplo 1.3.20√ Determine
√ a equação cartesiana do plano tangente à superfície
√
x + y + z = 1 no ponto ( 3/3, 3/3, 3/3).
2 2 2
Um vector perpendicular ao plano tangente a uma superfície S, representada de forma implícita

pela equação F (x, y, z) = 0, num ponto regular P , designa-se por vector normal, se tiver
norma 1. Pela Proposição 1.3.8, verifica-se que o vector normal à superfície S no ponto P =
(a, b, c) é dado por

1 1 ∂F ∂F ∂F
N= ∇F (P ) =
∂ x (P ), ∂ y (P ), ∂ z (P ) .

k∇F (P )k ∂F ∂F ∂F
∂ x (P ), ∂ y (P ), ∂ z (P )
Se a superfície S puder ser escrita, de forma explícita, por z = f (x, y), o vector normal à
superfície no ponto P = (a, b, c), com c = f (a, b), escreve-se do modo seguinte

1 ∂f ∂f
N= − (a, b), − (a, b), 1 .
∂f
−
∂x (a, b), − ∂f
(a, b), 1
∂ x ∂ y
∂y
Definição 1.3.8 Seja S uma superfície em R3 e P = (a, b, c) um ponto regular de S. Chama-se

recta normal à superfície S à recta perpendicular ao plano tangente à superfície S no ponto
de contacto P .
Proposição 1.3.9 Sejam S uma superfície em R3 , representada de forma implícita pela equação
F (x, y, z) = 0, e P = (a, b, c) um ponto regular de S. Então, as equações cartesianas da
recta normal à superfície S no ponto P = (a, b, c) são dadas por
x−a y−b z−c
∂F
= ∂F
= ∂F
.
∂x
(P ) ∂y
(P ) ∂z
(P )
Resulta da demonstração da proposição anterior que as equações paramétricas da recta

normal, digamos r(t) = (x(t), y(t), z(t)), à superfície S no ponto P = (a, b, c), são dadas por:
x(t) = a + ∂∂ Fx (a, b, c) t




r(t) ≡ (x(t), y(t), z(t)) = P + t ∇ F (P ) ⇔ y(t) = b + ∂∂ Fy (a, b, c) t , com t ∈ R .


z(t) = c + ∂∂ Fz (a, b, c) t

2.3 DERIVABILIDADE
No caso da equação da superfície S puder ser escrita na forma explícita z = f (x, y), as equações
cartesianas da recta normal à superfície S no ponto P = (a, b, c), com c = f (a, b), simplificam-se
a
x−a y−b
−∂f = −∂f = z−c
∂x
(a, b) ∂y
(a, b)
e as correspondentes equações paramétricas reduzem-se a


 x(t) = a − ∂∂ fx (a, b) t

∂f ∂f 
r(t) ≡ (x(t), y(t), z(t)) = (a, b, c) + t − (a, b), − (a, b), 1 ⇔ y(t) = b − ∂∂ fy (a, b) t .
∂x ∂y 


 z(t) = c + t
Exercício exemplo 1.3.21 Determine as equações cartesianas

√ √ e√ as equações paramétricas da
recta normal à superfície x2 + y 2 + z 2 = 1 no ponto ( 3/3, 3/3, 3/3).
1.3.8 Fórmula de Taylor

Na Definição 1.3.5 e Proposição 1.3.2, vimos que para uma função escalar f : RN −→ R de
classe C 1 numa bola aberta B(a, ε) ⊂ Df , com a = (a1 , . . . , aN ), se tem
f (x) = f (a) + ∇f (a) · (x − a) + o(kx − ak) ,
para todo x ∈ B(a, ε). Notemos que o segundo termo do segundo membro da fórmula acima -
o diferencial de f , pode ser escrito na forma matricial seguinte
 
x1 − a1
 x2 − a2 
∇f (a) · (x − a) = fx1 (a) fx2 (a) · · · fxN (a)  .
 ··· 
xN − aN
A fórmula anterior dá-nos uma aproximação, em termos do diferencial de primeira ordem de

uma função escalar f e designá-mo-la por aproximação de primeira ordem da função f
em torno do ponto a.
Exercício exemplo 1.3.22 Determine a aproximação de primeira ordem da função f (x, y) =

ex sen y em torno do ponto (x, y) = (0, 0).
A proposição seguinte dá-nos uma aproximação mais precisa da função f , a qual já vai fazer
intervir o diferencial de segunda ordem. Designamos esta fórmula por aproximação de se-
gunda ordem da função f em torno do ponto a.
Proposição 1.3.10 Seja f : RN −→ R um campo escalar de classe C 2 numa bola aberta

B(a, ε) ⊂ Df , com a = (a1 , . . . , aN ). Então, para todo x ∈ B(a, ε),
1
f (x) = f (a) + ∇f (a) · (x − a) + [(x − a)H(a)] · (x − a) + o(kx − ak2 ) ,
2!
2.3 DERIVABILIDADE
onde  
fx1 x1 fx1 x2 · · · fx1 xN
 fx2 x1 fx2 x2 · · · fx2 xN 
H = .. .. .. ..
 
. . . .

 
fxn x1 fxn x2 · · · fxn xN
é a denominada matriz Hessiana da função f e a notação o(kx − ak2 ) significa
o(kx − ak2 )
lim = 0.
x→a kx − ak2
SEM DEMONSTRAÇÃO: Ver, por exemplo, Apostol p. 308.
Notemos, também aqui, que o terceiro termo do segundo membro da fórmula da proposição
anterior - o diferencial de segunda ordem de f , pode ser escrito com a notação matricial seguinte
1
[(x − a)H(a)] · (x − a) =
2!
  
fx1 x1 (a) fx1 x2 (a) ··· fx1 xN (a) x1 − a1
1  fx2 x1 (a) fx2 x2 (a)
 ··· fx2 xN (a)  x2 − a2 
x1 − a1 x2 − a2 · · · xN − aN  .. .. .. .. .. .
 
2! . . . . .

  
fxn x1 (a) fxn x2 (a) · · · fxn xN (a) xN − aN
Exercício exemplo 1.3.23 Determine a aproximação de segunda ordem da função f (x, y) =

1/[(1 − x)(1 − y)] em torno do ponto (x, y) = (0, 0).
Podemos continuar este processo de aproximação de uma função escalar f até ao diferencial
de uma ordem N qualquer, desde que a função tenha derivabilidade suficiente. No entanto,
para as aproximações terem mais rigor, convém estabelecer formas precisas de determinar o
infinitésimo o(kx − akN ).
Proposição 1.3.11 (Fórmula de Taylor) Seja f : RN −→ R um campo escalar de classe

C M +1 , com M ∈ N, numa bola aberta B(a, ε) ⊂ Df , com a = (a1 , . . . , aN ). Então, para todo
x ∈ B(a, ε),
′ 1 ′′ 1 (M )
f (x) = f (a) + f(x−a) (a) + f(x−a) (a) + · · · + f (a) + RM (x − a)
2! M! (x−a)
onde
1 (M +1)
RM (x − a) = f(x−a) (a + θ(x − a)) , 0 < θ < 1.
(M + 1)!
A fórmula escrita na proposição anterior, designa-se por Fórmula de Taylor de ordem N

em torno do ponto a e RM (x − a) é o resto de ordem M desta fórmula, denominado resto
de Taylor. No caso de a = 0, esta fórmula designa-se por Fórmula de Mac-Laurin de
2.3 DERIVABILIDADE
(k)
ordem N. A notação f(x−a) (a) indica a derivada de ordem k, com 1 ≤ k ≤ M + 1, da função f
no ponto a segundo o vector x − a. Usando o resultado expresso na Proposição 1.3.1, podemos
escrever a fórmula de Taylor, até aos termos de segunda ordem, usando a notação
′
f(x−a) (a) = ∇f (a) · (x − a) ,
′′
f(x−a) (a) = [(x − a)∇(∇f )(a)] · (x − a) = [(x − a)H(a)] · (x − a) ,
onde H(a) é a matriz Hessiana de f , calculada no ponto a, referida na Proposição 1.3.10.
No entanto, para as derivadas de ordem superior a dois, esta notação torna-se de tal forma
(k)
complicada que é preferível não a utilizar. Por outro lado, atendendo à Definição 1.3.6, f(x−a) (a)
também pode ser interpretada como o diferencial de ordem k da função f no ponto a relativo
ao vector x − a, (d f )(x−a) (a). Usando a notação dos diferenciais, a fórmula de Taylor pode ser
escrita na forma seguinte:
1 2 1
f (x) = f (a) + (d f )(x−a) (a) + (d f )(x−a) (a) + · · · + (d M f )(x−a) (a) + RM (x − a)
2! M!
onde
1
RM (x − a) = (d M +1 f )(x−a) (a + θ(x − a)) , 0 < θ < 1 .
(M + 1)!
Conhecida a expressão do diferencial de primeira ordem da função escalar f relativo ao vector
x − a, podemos escrevê-lo na forma simbólica seguinte

∂ ∂
(d f )(x−a) (a) = (x1 − a1 ) + · · · + (xN − aN ) f (a) .
∂ x1 ∂ xN
De forma análoga, podemos escrever, simbolicamente, a expressão do diferencial de qualquer
ordem k, com 1 ≤ k ≤ M + 1, de f relativo ao vector x − a por
k
k ∂ ∂
(d f )(x−a) (a) = (x1 − a1 ) + · · · + (xN − aN ) f (a) .
∂ x1 ∂ xN
Portanto, a escrita da fórmula de Taylor em termos dos diferenciais é mais simples, o que, na
resolução de exercícios práticos, se torna numa grande vantagem.
Exercício exemplo 1.3.24 Determine a Fórmula de Mac-Laurin de terceira ordem da função
f (x, y) = ex y .
1.3.9 Exercícios
Derivadas parciais
1. Determine as (funções) primeira e segundas derivadas parciais das funções seguintes:
senh(x + y)
a) f (x, y) = arctg(x + ey ) ; b) f (x, y) = ;
ln(x − y)
c) f (x, y, z) = zsen(xy) + x cos(yz) ; d) f (x, y) = [cotg(x + y)]ln(x+y) ;
√ tgh(xyz)
e) f (x, y) = xy arccos(x + y) + 3
xy ; f) f (x, y, z) = .
xyz
2.3 DERIVABILIDADE
2. Verifique que a função

z = ln(ex + ey )
satisfaz a equação
2
∂2z ∂2z ∂2z

− = 0.
∂x2 ∂y 2 ∂x∂y
3. Considere as funções seguintes:
( 3 3 (
x −y 1
x2 +y 2
se (x, y) 6= (0, 0) xy 2 sen y
se y 6= 0
f (x, y) = ; g(x, y) = .
0 se (x, y) = (0, 0) 0 se y = 0
a) Estude a continuidade de f e g.
b) Determine as primeiras derivadas parciais de f e g em todos os pontos de R2 .
∂2f ∂2g
c) Determine ∂x∂y e ∂y∂x em todos os pontos de R2 .
4. Calculepos gradientes das funções seguintes nos pontos indicados:

a) z = x2 + y 2 no ponto (0, 3);
b) u = z 3 cos(x − no ponto 3 π
;

y)
x + sen(x + y − z) 4
π, 4
+ 1, 2
y
c) v = arctg x no ponto (2, 2);
d) u = xyz no ponto (1, −1, 1);
e) u = 2x2 + 3y 2 − cos z no ponto 1, 0, π2 .

5. Determine o produto interno do gradiente da função u = x2 y 2 z 2 no ponto (1, 1, −1) com

o gradiente da função u = x + y + z no ponto (2, −1, −1).
Derivadas direccionais
1. Determine as derivadas
p dasfunções seguintes nos pontos e segundo os vectores ~v indicados:
a) f (x, y) = ln x2 + y 2 , P= (1, 1), ~v = (1, 2);
b) f (x, y, z) = 2x2 − 3xy + 6z 2 , P= (1, 1, 0), ~v = ~e1 + ~e2 ;
c) f (x, y, z) = z − e sen y,
x
P= (ln 3, 32 π, −3), ~v = (x, y, z);
→
d) f (x, y) = xyex , P= (1, 1), ~v =PQ, sendo Q= (4, −3);
2 2
e) f (x, y) = ex y , P= (0, 1), ~v tem a direcção da recta y = 2x.
2. Determine as derivadas direccionais das funções seguintes nos pontos e direcções indica-
dos:
a) f (x, y) = x3 − 2x2 y + xy 2 + 1 no ponto (1, 2) na direcção ~v = (3, 4);
2 3
b) f (x, y) = 1 + x2 + y3 no ponto (2, 3) na direcção ~v = (3, 4);
c) f (x, y, z) = π1 (x − y + z)tg(x + y − z) no ponto (π, π, π) na direcção definida pelo vector
~v = (1, 2, 2);
d) f (x, y, z) = xsen(y + z) + ysen(x + z) + zsen(x + y) no ponto ( π4 , − π4 , π2 ) na direcção
~v = (1, 2, 2); √
e) f (x, y, z) = ln cos(x2 ) + cos(y 2) + cos(z 2 ) + e − 3 2 2 no ponto
√ √ √
π
2
, 2π , 2π na direcção ~v = (1, 2, 2).
2.3 DERIVABILIDADE
Derivabilidade
1. Estude as funções seguintes quanto à derivabilidade:
a) f (x, y) = xy + 2x2 ; b) f (x, y, z) = z − ex sen y ;
√ ( x |y|
xy se xy ≥ 0 √ 2 2 se (x, y) 6= (0, 0)
c) f (x, y) = ; d) f (x, y) = x +y .
0 se xy < 0 0 se (x, y) = (0, 0)
2. Determine o diferencial total das funções seguintes:
x2 − y 2
a) f (x, y) = tg sen(x + y) + y 2 ; b) f (x, y) = cos(x2 + y 2 ) +

;
ex
xy
c) f (x, y, z) = (x + y + z)senh(xyz) ; d) f (x, y) = (ln(x + y))e .
Derivada da função composta

1. Determine f ′ (t), onde f (x, y) = x2 − y2, x = cos(t) e y = sen(t), de dois modos distintos:
a) Substituindo x = cos(t) e y = sen(t) em f (x, y) = x2 − y 2 e determinando a derivada;
b) Usando o Teorema de Derivação da Função Composta.
2. Resolva o exercício anterior, agora para f (x, y) = cos(xy), x = et e y = e−t .
3. Determine f ′ (t) usando o Teorema de Derivação da Função Composta:
a) f (x, y) = x2 − y 2 , x = senh t , y = cosh t ;
2
ln(1 + ex )
a) f (x, y) = , x = t , y = e−t ;
1 + ey2
c) f (x, y, z) = exyz cos(x + y + z) , x = ln t , y = cos t , z = sen t ;
4. Para cada uma das alíneas seguintes:
a) f (u, v) = u2 v 3 + 2v 2 u , u = x + y , v = x2 − y 2 ;
x+y
b) f (u, v) = ln(u + v) , u = ex , v = ;
xy
d) f (u, v) = eu+v , u = t3 s2 , v = 5s , s = 2xy 3 , t = sen(xy) ;
determine as derivadas parciais fx e fy de dois modos distintos:
(i) substituindo as expressões e determinado a derivada;
(ii) usando o Teorema de Derivação da Função Composta.
5. Considere a função

x
f (x, y) = xy + xg , com y 6= 0 ,
y
em que g : R → R é de classe C 2 e tal que g(1) = −1 e g ′(1) = 2.
a) Calcule ∇f (1, 1).
b) Determine a derivada direccional de f no ponto (1, 1) na direcção do ponto (0, 2).
2.3 DERIVABILIDADE
Funções implícitas
1. Mostre que a equação
1
x3 + 4x2 y − y =
3
define y como função de x numa vizinhança do ponto (0, −1/3) e calcule y ′ (1).
2. Seja y uma função de x determinada pela equação
p y
ln x2 + y 2 = arctan .
x
Calcule yx e yxx .
3. Mostre que a equação
xey + yex + z + cos(xz) = 0
define z como função de x e y numa vizinhança do ponto (0, 0, −1) e calcule zx (0, 0) e
zy (0, 0).
4. Considere a equação
x ln(yz) − y ln(xz) = 0 .
a) Mostre que, numa vizinhança do ponto (1, 1, 1), esta equação define y como função de
x e z.
b) Calcule yx (1, 1), yxx (1, 1) e yxz (1, 1).
5. Considere a equação
2 +y 2
z 3 − z 2 + 2ex = 2.
a) Determine z0 de modo que esta equação defina z como função de x e y numa vizinhaça
do ponto (0, 0, z0).
b) Calcule zx (0, 0), zy (0, 0), zxx (0, 0), zyy (0, 0) e zxy (0, 0).
Aplicações geométricas
1. Considere a função
f (x, y) = x2 + 2xy + 2y 4 .
a) Determine os pontos onde o gradiente de f é vertical;
b) Determine os pontos da curva C= {(x, y) ∈ R2 : f (x, y) = 1} onde a tangente é
horizontal.
2. Determine os planos tangentes às superfícies seguintes nos pontos indicados:
2 2 2
a) x25 − y16 + z9 = 0 no ponto (5, 4, 3);
b) x2 y 2 − 2x3 z − yz 2 = 1 no ponto (−1, 1, 2);
c) z 2 + arccos(z − y) − x2 = π2 no ponto (−1, 1, 1);
d) x + arccos(z − y) + 5z = 1 + arccos 5 no ponto 0, 45 , 15 ;
2 3

e) y 3 + 3x2 z + yz = 7x no ponto (2, 1, 1).

3. Verifique se 2x + y − z = 1 é o plano tangente à superfície y 2 + 2y 3 z 2 + yz = 11x2 no
ponto (1, 1, 2).
2.4 EXTREMOS
4. Considere a função:
x
r
f (x, y) = 1− .
x+y
Determine:
a) a equação da linha de nível que passa no ponto (0, 1);
b) a direcção de maior crescimento da função no ponto (0, 1);
c) a equação do plano tangente à superfície z = f (x, y) no ponto (0, 1);
d) as equações paramétricas da recta normal à superfície z = f (x, y) no ponto (0, 1).
5. Considere a função:
x2

f (x, y) = ln 1 − .
y
Determine:
a) a equação da linha de nível que passa no ponto (1, 2);
b) a direcção de maior crescimento da função no ponto (1, 2);
c) a equação do plano tangente à superfície z = f (x, y) no ponto (1, 2);
d) as equações paramétricas da recta normal à superfície z = f (x, y) no ponto (1, 2).
Fórmula de Taylor
1. Determine a Fórmula de Taylor, até aos termos de terceira ordem, das funções seguintes
em torno dos pontos (a, b) indicados.
a) f (x, y) = xy 2 , com (a, b) = (1, 2) ;
b) f (x, y) = xseny + ysenx com (a, b) = (0, 0).
2. Determine a Fórmula de Mac-Laurin, até aos termos de terceira ordem, das funções
seguintes:
a) f (x, y) = cos x cos y;
b) f (x, y) = y x+1 .
1.4 Extremos
Consideremos uma função escalar de classe C 1
f: RN −→ R
(x1 , . . . , xN ) 7→ f (x1 , . . . , xN ) .
e seja a = (a1 , . . . , aN ) um ponto interior ao seu domínio Df ⊂ RN . Dizemos que a função f

tem um máximo relativo, ou máximo local, no ponto a, se
f (a) ≥ f (x) ∀ x = (x1 , . . . , xN ) ∈ B(a, ε) ,
para algum ε > 0. A função f tem um mínimo relativo, ou mínimo local, no ponto a, se
f (a) ≤ f (x) ∀ x = (x1 , . . . , xN ) ∈ B(a, ǫ) ,

2.4 EXTREMOS
para algum ǫ > 0. Designamos por máximo ou mínimo relativo, ou local, da função f ao
valor da imagem f (a). Quando não houver necessidade de especificar se este valor é máximo
ou mínimo relativo, diremos que f (a) é um extremo relativo, ou extremo local, da função
f . Em muitas aplicações estamos interessados no maior ou menor valor possível de uma função
escalar. A função f tem um máximo absoluto, ou máximo global, no ponto a, se
f (a) ≥ f (x) ∀ x = (x1 , . . . , xN ) ∈ Df .
A função f tem um mínimo absoluto, ou mínimo global, no ponto a, se
f (a) ≤ f (x) ∀ x = (x1 , . . . , xN ) ∈ Df .
Ao valor f (a), nas condições anteriores, designamos por máximo ou mínimo absoluto, ou
global. No caso de não querermos especificar se se trata de máximo ou mínimo, designaremos,
apenas, por extremo absoluto ou global. Em geral, uma função escalar não tem necessari-
amente um máximo ou mínimo absoluto. No entanto, por um resultado da Análise que não
demonstraremos aqui 3 , sabemos que toda a função contínua num domínio fechado e limitado
tem máximo e mínimo absolutos.
Proposição 1.4.1 Sejam f : RN → R um campo escalar e a ∈ int(Df ) onde f é de classe C 1 .

Se f tem um extremo no ponto a = (a1 , . . . , aN ), então
∇ f (a1 , . . . , aN ) = (0, . . . , 0) .
Aos pontos a nas condições da Proposição 1.4.1, isto é, aos pontos a ∈ int(Df ) onde f é de
classe C 1 e tais ∇ f (a) = 0, vamos designar por pontos de estacionaridade da função f .
Recordamos que anteriormente já havíamos introduzido a noção de pontos críticos, ou pontos
singulares, pontos onde as derivadas parciais eram todas nulas ou onde, pelo menos, uma das
derivadas parciais não existia. Neste sentido, os pontos de estacionaridade de uma função
escalar são pontos singulares. Os pontos de estacionaridade de uma função escalar poderão
ser de dois tipos diferentes: extremos ou pontos de sela. Os extremos podem ser máximos ou
mínimos e os pontos de sela são pontos de estacionaridade que não são extremos.
Exercício exemplo 1.4.1 (AULA TEÓRICA) Determine o ponto de estacionaridade da função

2 2
f (x, y) = e−(x +y ) e concluia, apenas por análise da função f , qual o tipo de extremo.
Proposição 1.4.2 Sejam f : RN → R um campo escalar de classe C m , com m > 1, e a =

(a1 , . . . , aN ) ∈ int(Df ). Suponhamos que a é um ponto de estacionaridade de f e m é a mais
baixa ordem das derivadas direccionais de f que não se anulam em a. Nestas condições:
(m)
1. se m é par e fu (a) > 0 para qualquer direcção u emergente de a, então f (a) é um
mínimo local;
(m)
2. se m é par e fu (a) < 0 para qualquer direcção u emergente de a, então f (a) é um
máximo local;
3
Ver, por exemplo, Campos Ferreira, p. 55 (Teorema de Weierstrass).
2.4 EXTREMOS
(m)
3. se m é ímpar, ou m é par e fu (a) não tem sinal definido para alguma direcção u
emergente de a, então f (a) é um ponto de sela;
(m)
4. se m é par e fu (a) tem sinal constante, com excepção de um número finito de direcções
(m)
u emergentes de a para as quais fu (a) = 0, nada se pode concluir sobre f (a).
Em termos práticos, o ponto 4 desta proposição é o de mais difícil resolução, tendo-se, na maioria
dos casos, de efectuar um estudo local, ou determinar as derivadas direccionais de f de ordem
superior a m, para averiguar a natureza do ponto de estacionaridade a. Os resultados expressos
na proposição anterior são particularmente úteis quando m = 2. Por isso, é importante ter
presente a caracterização da derivada direccional fu2 (a), com u = (u1 , . . . , uN ), através da
respectiva matriz Hessiana H(a):
  
fx1 x1 (a) · · · fx1 xN (a) u1
fu2 (a) = [u H(a)] · u = u1 · · · uN 
 .. .. ..   .. 
. . .  .  .
fxn x1 (a) · · · fxn xN (a) uN
Observemos que a matriz Hessiana H(a) é uma matriz simétrica de ordem N × N com as
entradas todas reais. Portanto, em termos da Álgebra Linear, fu2 (a) é uma forma bilinear
quadrática associada à matriz Hessiana H(a), que pode ser caracterizada como definida pos-
itiva ou negativa, semi-definida positiva ou negativa e, ainda, indefinida.
Definição 1.4.1 Consideremos a forma bilinear quadrática
fu2 (a) = [u H(a)] · u = uT H(a) u , uT = u1 · · · uN .

1. Diz-se que fu2 (a) é definida positiva, se uT H(a) u > 0 para todo u ∈ RN \ {0}.
2. Diz-se que fu2 (a) é definida negativa, se uT H(a) u < 0 para todo u ∈ RN \ {0}.
3. Diz-se que fu2 (a) é semi-definida positiva, se uT H(a) u ≥ 0 para todo u ∈ RN .
4. Diz-se que fu2 (a) é semi-definida negativa, se uT H(a) u ≤ 0 para todo u ∈ RN .
Se fu2 (a) não satisfaz nenhum dos casos da definição anterior, dizemos que fu2 (a) é uma forma
quadrática indefinida. Deste modo, a Proposição 1.4.2 pode ser escrita na forma algébrica
seguinte, recorrendo às caracterizações da forma quadrática fu2 (a).
Proposição 1.4.3 Sejam f : RN → R um campo escalar de classe C m , com m > 1, e a =

(a1 , . . . , aN ) ∈ int(Df ). Suponhamos que a é um ponto de estacionaridade de f e m = 2 é a
mais baixa ordem das derivadas direccionais de f que não se anulam em a. Nestas condições:
1. se fu2 (a) é definida positiva, então f (a) é um mínimo local;
2. se fu2 (a) é definida negativa, então f (a) é um máximo local;
3. se fu2 (a) é indefinida, então f (a) corresponde a um ponto de sela.;

2.4 EXTREMOS
4. se fu2 (a) é semi-definida positiva ou negativa, nada se pode concluir sobre f (a).
Tal como na Proposição 1.4.2, o ponto 4 da Proposição 1.4.3 é o de mais difícil resolução,
tendo-se, também, de efectuar um estudo local, ou determinar as derivadas direccionais de f
de ordem superior a 2, para averiguar a natureza do ponto de estacionaridade a. Na resolução
de exercícios práticos sobre extremos, é mais simples caracterizar a forma bilinear quadrática
fu2 (a) através da matriz Hessiana H(a) que lhe está associada. O resultado seguinte da Álgebra
Linear simplifica-nos bastante essa caracterização.
Proposição 1.4.4 Consideremos a forma bilinear quadrática
fu2 (a) = uT H(a) u , uT = u1 · · · uN .

Seja
 
fx1 x1 (a) · · · fx1 xN (a)
fx1 x1 (a) fx1 x2 (a) .. .. ..
△1 = fx1 x1 (a), △2 = det , ··· , △N = det  . . .
 
fx2 x1 (a) fx2 x2 (a)

fxn x1 (a) · · · fxn xN (a)
a cadeia de menores principais da matriz Hessiana H(a).

1. Se △k > 0 para todo k = 1, . . . , N, então fu2 (a) é definida positiva, pelo que f (a) é um
mínimo local;
2. Se △k < 0 para todo k ímpar e △k > 0 para todo k par, com k = 1, . . . , N, então fu2 (a)
é definida negativa, pelo que f (a) é um máximo local;
3. Se △k > 0 para algum k ímpar ou △k < 0 para algum k par, com k = 1, . . . , N, então
fu2 (a) é indefinida, pelo que f (a) corresponde a um ponto de sela;
4. Se △j = 0 para algum j = 1, . . . , N − 1 e △k =
6 0 para k > j, com k = 1, . . . , N, nada se
2
pode concluir sobre a forma quadrática fu (a) e, por consequência, sobre f (a).
SEM DEMONSTRAÇÃO: Ver, por exemplo, Dias Agudo, p. 186.
Se conjugarmos os casos 3 e 4 da Proposição 1.4.4, ainda podemos tirar as conclusões seguintes:
3-4 (i) se △j = 0 para algum j ímpar e △k < 0 para algum k > j par, então fu2 (a) é
indefinida, pelo que f (a) corresponde a um ponto de sela;
3-4 (ii) se △j = 0 para algum j par e △k > 0 para algum k > j ímpar, então fu2 (a) é
indefinida, pelo que f (a) corresponde a um ponto de sela.
Se △N = 0, então a característica4 da matriz Hessiana H(a) é menor do que N, digamos

c(H(a)) = r e r < N. Neste caso, o ponto 4 da Proposição 1.4.4, ainda, nos permite as
4
A característica de uma matriz A denota-se por c(A) e, por definição, é o maior número de linhas, ou
colunas, de A linearmente independentes.
2.4 EXTREMOS
afirmações seguintes, apesar de serem inconclusivas quanto à natureza dos pontos de esta-
cionaridade:
4 (i) se △k > 0 para todo k = 1, . . . , r e △k = 0 para todo k = r + 1, . . . , N, então fu2 (a) é

semi-definida positiva e nada se pode concluir f (a);
4 (ii) se, para todo k = 1, . . . , r, △k < 0 para k ímpar e △k > 0 para k par, e se, para todo
k = r + 1, . . . , N, △k = 0, então fu2 (a) é semi-definida negativa e nada se pode concluir sobre
f (a).
Exercício exemplo 1.4.2 Determine os pontos de estacionaridade da função a seguir indi-

cada e classifique-os como extremos:
f (x, y, z) = x2 + y 2 + 3z 2 + yz + 2xz − xy .
A Proposição 1.4.4 conjugada com a Proposição 1.4.3 simplifica bastante o estudo dos extremos
de funções escalares de duas variáveis.
Proposição 1.4.5 (O caso particular N = 2) Sejam f : R2 → R um campo escalar de

classe C 2 e
fxx (a, b) fxy (a, b)
H(a, b) = .
fyx (a, b) fyy (a, b)
a matriz Hessiana de f calculada no ponto (x, y) = (a, b) . Suponhamos que (a, b) é um ponto
de estacionaridade de f . Nestas condições:
1. se det[H(a, b)] > 0 e fxx (a, b) > 0, então f (a, b) é um mínimo local;
2. se det[H(a, b)] > 0 e fxx (a, b) < 0, então f (a, b) é um máximo local;
3. se det[H(a, b)] < 0, então f (a, b) corresponde a um ponto de sela;
4. se det[H(a, b)] = 0, nada se pode concluir.
Exercício exemplo 1.4.3 Determine os pontos de estacionaridade da função a seguir indi-

cada e classifique-os como extremos:
f (x, y) = x3 − 3x2 + y 2 .
1.4.1 Exercícios
1. Determine e classifique os pontos de estacionaridade da funções seguintes:
1
a) f (x, y) = x2 + y 4 ; b) f (x, y) = x2 − 2xy + y 3 − 3y ;
3
c) f (x, y) = x2 + y 2 + 4x − 2y + 5 ; d) f (x, y) = xy(x − 1) ;
e) f (x, y) = x4 + y 4 − 2x2 ; f) f (x, y) = xy ln(x2 + y 2) ;
2 +y 2 +z 2
g) f (x, y, z) = ex ; h) f (x, y, z) = x2 + y 2 + 3z 2 + yz + 2xz − xy .
2.4 EXTREMOS
2. Determine os extremos da função y definida implicitamente por
y 3 − 3x2 y + x3 = 3 .
3. Determine os extremos da função z definida implicitamente por
x + 2xy + 3z 2 + x2 z = 1 .
Capítulo 2
Cálculo integral vectorial
2.1 Integral de Linha

2.1.1 Primeiras noções
Nesta secção vamos desenvolver a noção de integral ao longo de uma curva que é descrita
num espaço euclidiano. Por simplicidade de escrita e das aplicações práticas, vamos consid-
erar apenas o plano R2 e o espaço R3 , os quais iremos designar apenas por plano e espaço,
respectivamente.
Definição 2.1.1 Seja I um intervalo fechado de R. Uma linha, ou curva, é um conjunto de

pontos cujas coordenadas cartesianas são dadas por:

x = x(t)
, no caso do plano;
y = y(t)

 x = x(t)
y = y(t) , no caso do espaço;
z = z(t)

onde x e y, ou x, y e z são funções contínuas da variável t ∈ I.
Por vezes, denominamos as linhas por curvas e, nesse sentido, é habitual denotar as linhas por
C, notação que iremos seguir neste curso. Outra notação frequente para as linhas é Γ. As linhas
são habitualmente escritas na forma paramétrica:
r(t) = (x(t), y(t)), t ∈ I, no caso do plano;
ou
r(t) = (x(t), y(t), z(t)), t ∈ I, no caso do espaço.
No entanto podem, também, aparecer escritas na forma cartesiana:
y = f (x), x ∈ I, no caso do plano;
42
2.4 EXTREMOS
onde f é uma função contínua dada, ou

y = f (x, z)
, x ∈ I, no caso do espaço,
z = g(x, y)
onde f e g são funções contínuas dadas. Podem, ainda, ser escritas por outras formas. Por ex-
emplo, fazendo uso das coordenadas polares no caso do plano R2 , ou das coordenadas cilíndricas
ou esféricas no caso do espaço R3 . Mas neste curso vamos estudar, apenas, curvas definidas na
forma paramétrica ou na forma cartesiana.
Definição 2.1.2 Designa-se por caminho qualquer linha a que lhe está adstrita um sentido
em que é descrita, uma origem, uma extremidade e a cada ponto uma multiplicidade.
Os caminhos são, também, denominados por trajectórias e são representados por aplicações
vectoriais contínuas. Sendo [a, b] ⊂ R, com a < b, temos as representações seguintes:
r : [a, b] −→ R2
⇔ r(t) = (x(t), y(t)), t ∈ [a, b];
t 7→ (x(t), y(t))
r : [a, b] −→ R3
⇔ r(t) = (x(t), y(t), z(t)), t ∈ [a, b];
t 7→ (x(t), y(t), z(t))
nos casos do plano e espaço, respectivamente. Designam-se por origem do caminho ao ponto
A = r(a) e por extremidade do caminho ao ponto B = r(b). Neste caso, dizemos que
o caminho está orientado de A para B. No caso de A = B, i.e. r(a) = r(b), diz-se que
o caminho é fechado. Caso contrário, diz-se que o caminho é aberto. Dizemos que um
ponto C tem multiplicidade k, se corresponder a k valores distintos t1 , . . . , tk ∈ [a, b], i.e., se
r(t1 ) = · · · = r(tk )=C.
Por simplicidade de escrita, mas tendo em conta que se trata de um abuso de linguagem, iremos
designar os caminhos por linhas e continuaremos a utilizar a letra C para denotar os caminhos.
Exemplo 2.1.1 Esboce os caminhos seguintes em sistemas de eixos cartesianos apropriados:

a) r(t) = (cos(t), sen(t)), com t ∈ [0, π/2].
Resolução: Temos x(t) = cos(t) e y(t) = sen(t), pelo que r(t) = (cos(t), sen(t)) corresponde ao
arco da circunferência x2 + y 2 = 1 entre os pontos de origem A = (cos(0), sen(0)) = (1, 0) e
extremidade B = (cos(π/2), sen(π/2)) = (0, 1).
b) r(t) = (cos(t), sen(t), 3t), com t ∈ [0, 2π].
Resolução: Temos x(t) = cos(t), y(t) = sen(t) e z(t) = 3t. Raciocinando de modo análogo
ao da alínea anterior, r(t) = (cos(t), sen(t), 3t) vai corresponder ao arco da hélice x2 + y 2 = 1
entre os planos z = 0 e z = 6π. Tem por origem o ponto A = (cos(0), sen(0), 3 × 0) = (1, 0, 0)
e a sua extremidade é B = (cos(2π), sen(2π), 3 × 2π) = (0, 1, 6π).
Definição 2.1.3 Sejam F uma função vectorial e C um caminho parametrizado por r(t), com
t ∈ [a, b] e a < b. Define-se o integral de linha do campo vectorial F sobre o caminho C
por Z b
dr
Z Z
F ≡ F(r) · dr = F(r(t)) · dt.
C C a dt
2.4 EXTREMOS
No caso do plano, se o caminho C está parametrizado por r(t) = (x(t), y(t)) e

F: R2 −→ R2
,
(x, y) 7→ (F1 (x, y), F2(x, y))
i.e. F(x, y) = (F1 (x, y), F2(x, y)), a expressão do integral de linha simplifica-se a
Z Z b
F= (F1 x′ + F2 y ′ ) dt.
C a
Se, no espaço, o caminho C está parametrizado por r(t) = (x(t), y(t), z(t)) e
F: R3 −→ R3
,
(x, y, z) 7→ (F1 (x, y, z), F2 (x, y, z), F3 (x, y, z))
i.e. F(x, y, z) = (F1 (x, y, z), F2 (x, y, z), F3 (x, y, z)) a expressão simplifica-se a
Z Z b
F= (F1 x′ + F2 y ′ + F3 z ′ ) dt.
C a
Se C for um caminho fechado, usamos a notação

I Z
F em vez de F.
C C
Observemos que as expressões dos segundos membros das equações anteriores estão escritas de
uma forma abreviada. Convém ter em mente que todas as funções aí envolvidas vêm escritas
em função de t. Por exemplo, F1 = F1 (x(t), y(t)) no caso do plano, ou F1 = F1 (x(t), y(t), z(t))
no caso do espaço e x′ = x′ (t).
Exemplo 2.1.2 Calcule o integral de linha C F nas condições seguintes:

R
a) F(x, y) = (−y, −xy) e C é o caminho parametrizado por r(t) = (cos t, sent), com t ∈ [0, π/2].
Resolução: Sabendo que x(t) = cos(t) e y(t) = sen(t), F1 (x(t), y(t)) = −y(t) = −sen(t) e
F2 (x(t), y(t)) = −x(t)y(t) = − cos(t)sen(t), temos
Z Z π Z π
2 2
′ ′
F= (F1 x + F2 y ) dt = (sen2 (t) − sen(t) cos2 (t)) dt =
C 0 0
t= π2
1 1 1 π 1
t − sen(2t) + cos3 (t) = − .
2 4 3 t=0 4 3
b) F(x, y, z) = (z, x, y) e C é o caminho parametrizado por r(t) = (cos t, sent, 3t), com t ∈
[0, 2π].
Resolução: Neste caso, x(t) = cos(t), y(t) = sen(t) e z(t) = 3t, F1 (x(t), y(t), z(t)) = z(t) = 3t,
F2 (x(t), y(t), z(t)) = x(t) = cos(t), F3 (x(t), y(t), z(t)) = y(t) = sen(t). Então
Z Z 2π Z 2π
′ ′ ′
F= (F1 x + F2 y + F3 z ) dt = (−3tsen(t) + cos2 (t) + 3sen(t)) dt =
C 0 0
t=2π
1 1
3t cos(t) − sen(t) + t + sen(2t) = 7π .
2 4 t=0
2.4 EXTREMOS
No caso particular da função vectorial F ter apenas uma componente não nula, obtemos uma
nova noção de integral de linha. Designaremos este novo integral por integral de linha relati-
vamente a uma variável.
Definição 2.1.4 Sejam f uma função escalar e C um caminho parametrizado por r(t) =
(x(t), y(t)), no caso do plano, ou r(t) = (x(t), y(t), z(t)), no caso do espaço, com t ∈ [a, b]
e a < b. O integral de linha do campo escalar f ao longo do caminho C relativo à variável
x, y ou z é dado, respectivamente, por
Z Z b Z Z b Z Z b
′ ′
f dx = f x dt, f dy = f y dt, f dz = f z ′ dt.
C a C a C a
Os casos mencionados na definição anterior, correspondem a considerar na Definição 2.1.3,

respectivamente, F = (f, 0) e F = (0, f ) no caso do plano, ou F = (0, 0, f ). F = (f, 0, 0),
F = (0, f, 0) e F = (0, 0, f ) no caso do espaço.
Exemplo 2.1.3 Calcule os integrais de linha

Z Z Z
f dx, f dy e f dz,
C C C
onde f (x, y, z) = xy + y 2 − xyz e C é o arco da parábola dado por y = x2 e z = 0 entre os

pontos A = (−1, 1, 0) e B = (2, 4, 0).
Resolução: Observemos que o caminho em questão é um ramo da parábola y = x2 inscrita no
plano z = 0. Uma parametrização possível para este caminho é r(t) = (t, t2 , 0), com t ∈ [−1, 2] e
onde fizemos x(t) = t, y(t) = x(t)2 = t2 e z(t) = 0. Por outro lado, tem-se f (x(t), y(t), z(t)) =
x(t)y(t) + y(t)2 − x(t)y(t)z(t) = t3 + t4 . Logo:
2 2 t=2
1 1 207
Z Z Z
′
f dx = f (x(t), y(t), z(t)) x (t) dt = (t + t ) dt = t4 + t5
3 4
= ;
C −1 −1 4 5 t=−1 20
2 2 t=2
1 5 1 6 213
Z Z Z
′ 4 5
f dy = f (x(t), y(t), z(t)) y (t) dt = 2 (t + t ) dt = 2 t + t = ;
C −1 −1 5 6 t=−1 15
Z Z 2 Z 2
′
f dy = f (x(t), y(t), z(t)) z (t) dt = (t3 + t4 ) × 0 dt = 0 .
C −1 −1
Definição 2.1.5 Sejam f uma função escalar e C um caminho parametrizado por r(t) =
(x(t), y(t)), no caso do plano, ou r(t) = (x(t), y(t), z(t)), no caso do espaço, com t ∈ [a, b]
e a < b. O integral de linha do campo escalar f ao longo do caminho C relativamente ao
comprimento de arco, denotado por s, é definido por
Z b
dr
Z
f ds = f (r(t))
dt dt,

C a
sendo ds = kdr/dtk dt designado o elemento de arco.

2.4 EXTREMOS
No caso do plano e do espaço, a expressão anterior simplifica-se, respectivamente a

Z Z b p Z Z b p
f ds = f (x′ )2 + (y ′ )2 dt e f ds = f (x′ )2 + (y ′ )2 + (z ′ )2 dt
C a C a
Exemplo 2.1.4 Calcule o integral de linha C f ds, onde:

R
a) f (x, y) = x2 y e C é o caminho parametrizado por r(t) = (sen cos t), com t ∈ [0, π/2];
Resolução: Observe-se que, neste caso, o caminho em questão é o ramo da circunferência x2 +
y 2 = 1 entre os pontos A = (sen(0), cos(0)) = (0, 1), de origem, e B = (sen(π/2), cos(π/2)) =
(1, 0), extremidade (percorrido no sentido contrário ao dos ponteiros do relógio). Aqui, x(t) =
sen(t), y(t) = cos(t) e f (x(t), y(t)) = x(t)2 y(t) = sen2 (t) cos(t). Então
Z Z π
2 p
f ds = f (x(t), y(t)) (x′ (t))2 + (y ′(t))2 dt =
C 0
π
1 3 t= π2 1
Z
2 p
sen2 (t) cos(t) cos2 (t) + (−sen(t))2 dt = sen (t) t=0 = .
0 3 3
b) f (x, y, z) = 1/(x2 + y 2 + z 2 ) e C é o caminho parametrizado por r(t) = (cos(t), sen(t), t),
com t ∈ [0, 1].
Resolução: Agora, o caminho é o ramo de hélice x2 + y 2 = 1 entre os planos z = 0 e z = 1.
A parametrização é dada por x(t) = cos(t), y(t) = sen(t) e z(t) = t. Logo, f (x(t), y(t), z(t)) =
(cos2 (t) + sen2 (t)2 + t2 )−1 = (1 + t2 )−1 e
Z Z 1 p
f ds = f (x(t), y(t), z(t)) (x′ (t))2 + (y ′(t))2 + (z ′ (t))2 dt =
C 0
√
√ Z 1
dt √ t=1 2
2 2
= 2 [arctgt]t=0 = π.
0 1+t 4
2.1.2 Motivação
O conceito de integral de linha tem bastantes aplicações em diversos campos das ciências. Va-
mos aqui considerar, apenas como exemplos de aplicação, um exemplo físico e outro geométrico.
Definição 2.1.6 Sejam C um caminho parametrizado por r(t), onde t ∈ [a, b] e F um campo
de forças. Definimos o trabalho realizado pela força F quando o seu ponto de aplicação
percorre o caminho C por Z Z
W= F≡ F(r(t)) · dr.
C C
Deste modo, verificamos que o trabalho realizado por um campo de forças ao longo de uma
curva é dado pela Definição 2.1.3.
Exemplo 2.1.5 Calcule o trabalho realizado pelo campo de forças F(x, y, z) = (x, −xy, z 2 ),
sobre uma partícula que se move sobre a hélice r(t) = (cos(t), sen(t), t), desde o ponto A =
(1, 0, 0) até ao ponto B = (−1, 0, 3π).
2.4 EXTREMOS
Resolução: A parametrização do caminho é x(t) = cos(t), y(t) = sen(t) e z(t) = t e, pela última
componente, vemos que t ∈ [0, 3π]. Então, F(x(t), y(t), z(t)) = (cos(t), − cos(t)sen(t), t2 ) e
Z Z 3π Z 3π
′ ′ ′
− cos(t)sen(t) − cos2 (t)sen(t) + t2 dt

W= F= (F1 x + F2 y + F3 z ) dt =
C 0 0
t=3π
1 1 1 2
= cos2 (t) + cos3 (t) + t3 = 9π 3 − .
2 3 3 t=0 3
Definição 2.1.7 Seja C um caminho parametrizado por r(t), onde t ∈ [a, b]. Definimos o
comprimento do caminho C por
Z b
dr
Z
s= ds = dt.
dt
C a
Por abuso de linguagem, poderemos falar em comprimento da curva C em vez de compri-

mento do caminho C. Tal como vimos acima, também aqui, a expressão para o comprimento
da curva C pode ser escrita nas formas mais simples
Z bp Z bp
s= ′ 2 ′ 2
(x ) + (y ) dt e s= (x′ )2 + (y ′ )2 + (z ′ )2 dt,
a a
nos casos do plano e do espaço, respectivamente.
Exemplo 2.1.6 Calcule o comprimento das curvas seguintes:

a) circunferência (x + 2)2 + (y − 1)2 = 9;
Resolução: Uma parametrização possível é x(t) = 3 cos(t)−2 e y(t) = 3sen(t)+1, com t ∈ [0, π].
O comprimento desta circunferência é, então, dado por:
Z Z 2π p Z 2π p
ds = ′ 2 ′ 2
(x (t)) + (y (t)) dt = (−3sen(t))2 + (3 cos(t))2 dt = 6π .
C 0 0
b) hélice, parametrizada por r(t) = (t, sen(t), cos(t)), desde o ponto A = (0, 0, 1) até ao ponto
B = (3π, 0, −1).
Resolução: A parametrização indicada é x(t) = t, y(t) = sen(t) e z(t) = cos(t), e, pela primeira
componente, vemos que t ∈ [0, 3π]. Assim, o comprimento da hélice é dado por:
Z Z 3π p Z 3π p √
ds = ′ 2 ′ 2 ′ 2
(x (t)) + (y (t)) + (z (t)) dt = 1 + (cos(t))2 + (−sen(t))2 dt = 3 2π .
C 0 0
2.1.3 Propriedades gerais

Proposição 2.1.1 Sejam F e G campos vectoriais, α uma constante real e C um caminho.
Então:
1. C (F + G) = C F + C G;
R R R
2. C αF = α C F;
R R
2.4 EXTREMOS
3. Se C é a união de dois caminhos, digamos C1 e C2 , então

R R R
C
F= C1
F+ C2
F.
DEMONSTRAÇÃO - EXERCÍCIO:
Exemplo 2.1.7 Calcule o integral de F(x, y) = (y, x) ao longo do caminho rectilíneo C que
une a origem de referencial ao ponto (3, 1), passando pelo ponto (3, 0).
Antes de passarmos ao resultado seguinte, convém referir que se entende por caminho oposto
de um dado caminho C como sendo o caminho percorrido ao longo da curva C, mas em sentido
contrário. Denotamos o caminho oposto por C − .
Proposição 2.1.2 Seja C um caminho parametrizado por r(t), com t ∈ [a, b] e C − o seu
caminho oposto. O caminho C − é parametrizado por r− (t) = r(a + b − t) e temos
Z Z
F=− F
C− C
Exemplo 2.1.8 a) Calcule o integral de F(x, y) = (x2 − 2xy, y 2 − 2xy) ao longo da parábola
y = x2 desde o ponto A = (−1, 1) até ao ponto B = (1, 1).
b) Calcule o integral anterior, mas no caminho oposto.
2.1.4 Integral de linha independente do caminho

Como temos vindo a ver nas secções precedentes, um dos principais problemas para o cálculo
dos integrais de linha prende-se com a determinação de uma parametrização para o caminho.
Contudo, existem integrais de linha para os quais não é necessário conhecer como é descrito o
caminho, isto é, a parametrização. Basta conhecermos a origem e a extremidade do caminho.
Recordemos que se define e denota o gradiente de uma função escalar f por

∂f ∂f
∇f = , , no caso do plano,
∂x ∂y

∂f ∂f ∂f
∇f = , , , no caso do espaço.
∂x ∂y ∂z
Por vezes, o gradiente de uma função f é denotado por grad f .
Proposição 2.1.3 Sejam C um caminho que une o ponto A ao ponto B e F um campo vectorial
tal que F = ∇f para alguma função escalar f . Então
Z
F = f (B) − f (A).
C
2.4 EXTREMOS
A proposição anterior diz-nos que, a existir tal campo escalar, o integral de linha vai ser inde-
pendente do caminho que é percorrido entre o ponto A e o ponto B. Isto é, o integral de linha
sobre caminhos distintos que unam o ponto A ao B tem sempre o mesmo valor. Em particular,
se C for um caminho fechado, A = B e, por consequência C F = 0. Deste modo, vamos denotar
R
o integral independente do caminho por

Z B
F.
A
Em termos físicos, podemos interpretar o resultado anterior do modo seguinte. Quando um

potencial (a função escalar f ) existe, o trabalho realizado por uma força (a função vectorial F)
para mover uma partícula ao longo de uma curva que una o ponto A a B, é a diferença entre
os potenciais de B e A.
Na resolução de exercícios práticos, há que determinar a função escalar f , para podermos

calcular o integral de linha.
Exemplo 2.1.9 Calcule o integral de linha de F(x, y, z) = (2xy 3 z, 3x2 y 2 z, x2 y 3) ao longo do

caminho que une o ponto A = (1, −1, 2) a B = (−3, 2, 5).
Existem casos em que o enunciado é omisso quanto ao integral ser independente do caminho
ou não. Nestes casos, necessitamos de condições que nos garantam que determinado campo
vectorial corresponde ao gradiente de alguma função escalar. O próximo resultado fornece-nos
uma condição necessária, mas não suficiente, para verificarmos isso.
Proposição 2.1.4 Seja F um campo vectorial. Suponhamos que existe uma função f tal que
F = ∇f . Então:
∂F2 ∂F1
= , no caso do plano; (2.1.1)
∂x ∂y
ou
∂F3 ∂F2 ∂F1 ∂F3 ∂F2 ∂F1
= , = , = , no caso do espaço. (2.1.2)
∂y ∂z ∂z ∂x ∂x ∂y
Este resultado é muito útil pelo que diz a sua afirmação contra-recíproca. Isto é, se a condição
(2.1.1), ou a condição (2.1.2), não é satisfeita, então não existe um campo escalar f tal que
F = ∇f .
Exemplo 2.1.10 Considere o integral de linha de F = (5z, xy, x2 z) ao longo de uma curva
que une o ponto A = (0, 0, 0) a B = (1, 1, 1).
a) Verifique que este integral não é independente do caminho.
b) Para corroborar a), calcule este integral de linha quando o caminho é parametrizado por:
(i) r1 (t) = (t, t, t), com t ∈ [0, 1];
(ii) r2 (t) = (t, t, t2 ), com t ∈ [0, 1].
2.4 EXTREMOS
A questão que agora se põe é a de saber em que condições a afirmação recíproca da expressa
na Proposição 2.1.4 é suficiente para o integral de linha ser independente do caminho.
Proposição 2.1.5 Seja C um caminho contido num domínio simplesmente conexo1 D de R2
ou R3 que contém o caminho C. Suponhamos que F(x, y) = (F1 (x, y), F2(x, y)), ou F(x, y, z) =
(F1 (x, y, z), F2 (x, y, z), F3 (x, y, z)) são tais que, respectivamente,
∂F2 ∂F1
= , no caso do plano;
∂x ∂y
ou
∂F3 ∂F2 ∂F1 ∂F3 ∂F2 ∂F1
= , = , = , no caso do espaço.
∂y ∂z ∂z ∂x ∂x ∂y
Então o integral C F é independente do caminho C.
R
2.1.5 Exercícios
1. Esboce as curvas seguintes no plano:
a) r(t) = (t, 2t√
− 1), com t ∈ [0, 1];
b) r(t) = (t, − 1 − t2 ), com t ∈ [−1, 1];
c) r(t) = (t, 1/t), com t ∈ (−∞, 0);
3
d) r(t) = 2 cos(t) , 3sen(t) , com t ∈ [0, 2π].

2. Esboce o caminho rectilíneo fechado C parametrizado por:

se 1 ≤ t < 2

 (t, t)
r(t) = (4 − t, t) se 2 ≤ t < 3 .
(1, 6 − t) se 3 ≤ t ≤ 5

3. Esboce o caminho fechado C parametrizado por:

(−t, 1 − t2 ) −1 ≤ t < 1

r(t) ≡ (x(t), y(t)) =
(t − 2, 0) 1≤t≤3
4. Esboce as curvas seguintes no espaço:

a) r(t) = (cos t, sen (2t), 0), com 0 ≤ t ≤ 2π;
b) r(t) = (t, t2 , t3 ), com 0 ≤ t ≤ 2.
5. Calcule o integral de linha C F nos casos seguintes:
R

a) F(x, y) = 4x−y 2 +y 2 , 4x2 +y 2 , onde C é o caminho dado no Exercício 1-d).
x
b) F(x, y) = (−y, x), onde C é o caminho dado no Exercício 3.

c) F(x, y, z) = (y 2 − z 2 , 2yz, −x2 ) e C é o caminho parametrizado por r(t) = (t, t2 , t3 ),
com t ∈ [0, 1];
d) F(x, y, z) = (cosh x, senh y, ez ) e C é o caminho parametrizado por r(t) = (t, t2 , t3 ),
com t ∈ [0, 2];
1
Um conjunto D diz-se um domínio simplesmente conexo, se qualquer caminho fechado C contido em D puder
ser encolhido até um ponto no interior de D sem intersectar pontos exteriores ao conjunto D.
2.4 EXTREMOS
6. Sendo C o caminho parametrizado por r(t) = (2 cos t, sen t), calcule os integrais de linha
seguintes:
Z Z Z
2 2
a) y dx; b) (x + y ) dy; c) y 2 dx + x2 dy.
C C C
7. Calcule o integral seguinte, onde C é o caminho do Exercício 2:

Z
2(x2 + y 2 ) dx + 2(x + y)2 dy.
C
8. Calcule o integral C y dx + z dy + x dz, onde C é o caminho percorrido no sentido positivo

R
e que resulta da intersecção das superfícies x + y = 2 e x2 + y 2 + z 2 = 2(x + y).

9. Calcule os integrais de linha C f ds, onde:
R
a) f (x, y) = x + y e C é o caminho rectilíneo fechado, percorrido no sentido positivo, que

passa pelos pontos A = (0, 0), B = (1, 0) e C = (0, 1);
b) f (x, y, z) = z e C é o caminho parametrizado por r(t) = (t cos t, t sen t, t), com 0 ≤ t ≤
2π.
10. Calcule os comprimentos das curvas seguintes:
a) r(t) = (cos t, sen t, t), com t ∈ [0, 2π];
b) asteróide x2/3 + y 2/3 = 1;
c) curva parametrizada pelo caminho dado no Exercício 2.
11. Mostre que a função vectorial envolvida nos integrais seguintes é o gradiente de um
potencial e calcule o valor dos respectivos integrais:
Z (3,π/2) Z (2,4,0)
x 2
a) e (cos y dx − sen y dy); b) ex−y+z (dx − dy + 2z dz);
(0,π) (0,−1,1)
Z (2,1,1)
c) yz dx + xz dy + xy dz.
(1,1,1)
12. Calcule os integrais de linha seguintes:

Z (4,3) Z (1,1,1)
2
a) (3z dx+ 6xz dz); b) [yz senh (xz) dx + cosh(xz) dy + xy senh (xz) dz] ;
(−1,5) (0,2,3)
Z (π,0,2)
c) (z 2 − ysen(x)) dx + (cos(x) − 2z) dy + (2xz − 2y + z) dz .
( π
4 )
,1,1
2.2 Integral duplo

Nesta secção vamos estender a noção de integral já conhecida para funções escalares de duas
variáveis. As regiões de integração vão ser agora subconjuntos de R2 . Primeiro, consideramos
regiões de integração rectangulares e depois consideramos regiões mais gerais com fronteiras
curvilíneas.
2.4 EXTREMOS
2.2.1 Integral de Riemann

Comecemos por estender a noção de partição de um intervalo limitado de R a um rectângulo
(limitado) contido em R2 . Para isso, recordemos que uma partição de um intervalo [a, b] ⊂
R é um conjunto finito de pontos, digamos x0 , x1 , . . . , xn , que divide o intervalo [a, b] em
subintervalos tais que a = x0 < x1 < · · · < xn−1 < xn = b, onde n é um número natural
arbitrário. O comprimento do i-ésimo intervalo é △xi = xi − xi−1 . A localização dos pontos
x0 , x1 , . . . , xn , e a consequente divisão do intervalo [a, b], é arbitrária. Em particular, os
subintervalos [xi−1 , xi ] não têm necessariamente o mesmo comprimento. A partição, desta
forma definida, vai ser denotada por P e escreve-se do modo seguinte:
P : a = x0 < x1 < . . . xn−1 < xn = b.
Definição 2.2.1 Sejam a, b, c e d números reais tais que a < b e c < d. Consideremos o
rectângulo
R = (x, y) ∈ R2 : a ≤ x ≤ b, c ≤ y ≤ d

e as partições dos intervalos [a, b] e [c, d], definidas, respectivamente, por
P1 : a = x0 < x1 < . . . xn−1 < xn = b e P2 : c = y0 < x1 < . . . ym−1 < xm = b,
onde n e m são números naturais arbitrários. Designa-se por partição do rectângulo R ao

conjunto seguinte
P = {(xi , yj ) ∈ R : 0 ≤ i ≤ n, 0 ≤ j ≤ m} .
Uma partição P do rectângulo R tal como acima definimos, determina mn subrectângulos de

R:
Rij = (xi , yj ) ∈ R2 : xi−1 ≤ x ≤ xi , yi−1 ≤ y ≤ yi ,

cujas áreas são dadas por

△xi △yj = (xi − xi−1 )(yj − yj−1).
A colecção destes subrectângulos Rij , é o que verdadeiramente constitui a partição P .
Definição 2.2.2 Seja f uma função definida num rectângulo R ⊂ R2 . Designamos por soma
de Riemann da função f no rectângulo R à quantidade seguinte:
n,
X m n X
X n
f (x∗ij )△xi △yj ≡ f (x∗ij )△xi △yj ≡ f (x∗11 )△x1 △y1 + · · · f (x∗nm )△xn △ym ;
i=0, j=0 j=0 i=0
onde x∗ij são pontos seleccionados aleatoriamente nos subrectângulos Rij respectivos.
Para a noção de integral duplo, interessa-nos que as partições sejam muito finas. Definimos a
quantidade que define a finura de dada partição P de um rectângulo R ⊂ R2 por
q
|P | = max (△xi )2 + (△yj )2 .
i, j
Observemos que |P | é o comprimento da maior diagonal de todos os subrectângulos Rij con-

siderados na partição P .
2.4 EXTREMOS
Definição 2.2.3 Sejam f uma função definida num rectângulo R ⊂ R2 e
P = {(xi , yj ) ∈ R : 0 ≤ i ≤ n, 0 ≤ j ≤ m}
uma partição arbitrária de R. Diz-se que a função f é integrável (à Riemann) no rectângulo

R, se existir (e for finito)o limite seguinte:
n X
X n
lim f (x∗ij )△xi △yj ,
|P |→0
j=0 i=0
independentemente de como a partição P do rectângulo R é formada, ou de como os pontos x∗ij

pertencentes aos subrectângulos Rij são escolhidos.
No caso de existir, o limite da definição anterior designa-se por integral da função f e denota-
se por Z Z Z Z
f (x, y) dx dy ou f (x, y) dy dx,
R R
onde dx dy, ou dy dx, indica um elemento de área ao qual ainda não está subjacente nenhuma
ordem de integração.
A noção de função integrável que acabamos de introduzir, extende-se a qualquer função definida
num conjunto limitado D ⊂ R2 que não seja propriamente um rectângulo. Apenas temos de
considerar um rectângulo R que contenha D e aí fazer a análise anterior. O único cuidado
a tomar para a definição fazer sentido, é fixar o valor de f (x∗ij ) igual a zero quando x∗ij não
pertencer a D.
Proposição 2.2.1 Seja D um subconjunto de R2 limitado e cuja fronteira é constituída pela

união de curvas (contínuas). Consideremos uma função f definida em D. Se f é contínua em
D excepto, quanto muito, num conjunto de medida nula2 , então f é integrável em D.
SEM DEMONSTRAÇÃO: Ver, por exemplo, Dias Agudo, Volume I, p. 236.
A proposição anterior diz-nos basicamente que todas as funções elementares que conhecemos
serão integráveis em domínios limitados.
As propriedades do análogo integral definido em R, serão válidas mutatis mutandis para o

integral duplo. Em particular, o integral duplo é um operador linear e satisfaz a propriedade
aditiva dos integrais.
2.2.2 Integral repetido

Em integração de funções reais de uma variável real apenas, o denominado Teorema Funda-
mental do Cálculo Integral dá-nos um modo muito prático de calcular integrais. A proposição
seguinte vai-nos permitir calcular alguns integrais duplos usando integrações repetidas de
funções reais de uma variável real.
2
Um subconjunto E de R2 diz-se de medida nula, se existir uma quantidade finita de rectângulos cuja união
contém E e de tal modo que a soma das áreas é tão pequena quanto se queira.
2.4 EXTREMOS
Proposição 2.2.2 (Fubini) Seja f uma função integrável num rectângulo

R = (x, y) ∈ R2 : a ≤ x ≤ b, c ≤ y ≤ d , onde a < b e c < d.

Rb Rd
Suponhamos que a f (x, y) dx existe para qualquer y ∈ [c, d], e que c f (x, y) dy existe para
qualquer x ∈ [a, b]. Então
Z Z Z d Z b Z b Z d
f (x, y) dx dy = f (x, y) dx dy = f (x, y) dy dx.
R c a a c
Os integrais expressos no segundo membro da equação dada na proposição anterior designam-se

por integrais repetidos ou iterados. Observe-se que na aplicação do Teorema Fundamental
do Cálculo Integral ao cálculo do integral entre parêntesis rectos, se primitiva a função f em
relação à variável aí referida, fixando a outra variável como se fosse constante.
Exemplo 2.2.1 Esboce a região de integração e calcule o integral respectivo:

Z πZ π
2
sen(2x − 3y) dx dy.
0 0
A Proposição 2.2.2 pode ser generalizada a qualquer domínio limitado D ⊂ R2 . Isto é, o cálculo
de um integral duplo num domínio limitado, resume-se ao cálculo de integrais repetidos.
Proposição 2.2.3 1. Sejam g e h duas funções reais de uma variável real, contínuas num
intervalo [a, b] ⊂ R, com a < b, e tais que, para cada a ≤ x ≤ b, g(x) ≤ h(x). Consideremos
uma função f contínua no domínio
D = (x, y) ∈ R2 : a ≤ x ≤ b, g(x) ≤ y ≤ h(x) .

Então Z b "Z #
Z Z h(x)
f (x, y) dx dy = f (x, y) dy dx.
D a g(x)
2. Sejam i e j duas funções reais de uma variável real, contínuas num intervalo [c, d] ⊂ R, com
c < d, e tais que, para cada c ≤ y ≤ d, i(y) ≤ j(y). Consideremos uma função f contínua no
domínio
D = (x, y) ∈ R2 : c ≤ y ≤ d, i(y) ≤ x ≤ j(y) .

Então "Z #
Z Z Z d j(y)
f (x, y) dx dy = f (x, y) dx dy.
D c i(y)
SEM DEMONSTRAÇÃO: Consequência da Proposição 2.2.2.
Exemplo 2.2.2 Esboce as regiões de integração e calcule os integrais respectivos:

Z π Z sen(y) Z 1 Z 1−x2
2
x
a) e cos(y) dx dy e b) x2 dy dx.
0 0 0 1−x
2.4 EXTREMOS
Na maioria das situações, podemos integrar tanto primeiro em ralação a uma variável, digamos
x, e em seguida relativamente à outra, y, como pela ordem inversa. No entanto, em qualquer
exercício prático, uma das alternativas anterior é mais fácil de calcular do que a outra. Mas,
existem situações em que, por diversas razões, é manifestamente impossível calcular o integral
primeiro relativamente a uma das variáveis, sendo a mais comum a impossibilidade de deter-
minar a primitiva da função dada em relação a essa variável. Nestes casos, somos obrigados a
calcular o integral como um integral repetido, mas com uma única ordem de integração possível.
No caso dos limites de integração serem constantes, então pela Proposição 2.2.2, podemos in-
verter a ordem de integração da forma que nos for mais conveniente. Contudo, se os limites de
integração não são constantes, a inversão da ordem de integração vai afectar também as funções
que limitam a região de integração. Na prática, ao inverter a ordem de integração destes inte-
grais, vamos considerar as inversas das funções que limitam o domínio, fazendo depois a devida
mudança de variável.
Exemplo 2.2.3 Inverta a ordem de integração do integral duplo seguinte, onde f (x, y) é uma
função arbitrária integrável no domínio indicado:
Z 1 Z − 12 ln y
f (x, y) dx dy.
e−2 ln y
Exemplo 2.2.4 Comece por verificar que, usando os métodos de primitivação conhecidos, é
impossível calcular o integral seguinte pela ordem de integração que está escrito. Inverta a
ordem de integração e calcule o seu valor:
Z 3Z 9
3
x3 ey dy dx .
0 x2
2.2.3 Mudança de variáveis

Em muitos exercícios práticos, e quando é possível, a simples inversão de ordem de integração
não simplifica muito o cálculo do integral duplo. Tal como acontece em integração de funções
reais de uma só variável real, também aqui podemos fazer uma transformação do integral duplo,
de modo ao seu cálculo ser mais simples. Esta transformação corresponde a uma mudança de
variáveis. Neste caso, não só a expressão designatória da função a integrar vem alterada, assim
como a própria região de integração passa a ser diferente.
Proposição 2.2.4 (Mudança de variáveis) Sejam D e Ω dois subconjuntos abertos de R2

limitados por funções diferenciáveis e seja
ϕ: Ω −→ D
(u, v) 7→ (x, y) = ϕ(u, v) ≡ (ϕ1 (u, v), ϕ2(u, v))
uma função bijectiva com derivadas parciais contínuas. Então, temos,

∂x ∂x

∂(x, y)
Z Z Z
f (x, y) dx dy = f (ϕ(u, v)) |J| du dv, J = det ≡ det ∂u ∂v .
∂y ∂y
D ϕ−1 (D) ∂(u, v) ∂u ∂v
2.4 EXTREMOS
SEM DEMONSTRAÇÃO: Ver, por exemplo, Dias Agudo, Volume I, p. 258
Observe-se que, por ϕ ser bijectiva, ϕ−1 (D) = Ω. A quantidade J expressa na proposição ante-
rior designa-se por jacobiano da mudança de variáveis (x, y) para (u, v) e a matriz [∂(x, y)/∂(u, v)]
é a denominada matriz jacobiana dessa transformação. Esta matriz poderá, ainda, ser escrita
da forma seguinte: ∂ϕ1 ∂ϕ1
∂(x, y)
≡ ∂ϕ∂u
2
∂v
∂ϕ2 .
∂(u, v) ∂u ∂v
Exemplo 2.2.5 Fazendo a mudança de variáveis u = x − y e v = x + y, calcule o integral

seguinte:
cos(x − y)
Z Z
onde D = (x, y) ∈ R2 : 1 ≤ x + y ≤ 2, x ≥ 0, y ≥ 0 .

dx dy,
D sen(x + y)
No integral duplo, o exemplo de mudança de variáveis com mais interesse prático é a trans-
formação para coordenadas polares. Recordemos o que são as coordenadas polares de um
ponto no plano. Seja P um ponto em R2 (plano) cujas coordenadas rectangulares num sistema
de eixos cartesiano são dadas por (x, y). Fixemos o triângulo rectângulo de vértices (0, 0), (x, 0)
e P = (x, y). Designemos a hipotenusa deste rectângulo por r e seja θ o ângulo formado entre
o semi-eixo positivo dos xx e a semi-recta com origem no ponto de coordenadas (0, 0) e que
passa por P = (x, y). Da trigonometria elementar, temos o seguinte:
cateto adjacente

 cos(θ) = = xr
hipotenusa


x = r cos(θ)
⇔ (2.2.3)
 cateto oposto y
y = rsen(θ) .
 sen(θ) =

=r
hipotenusa
Deste modo, podemos definir o ponto P num sistema de eixos cartesiano à custa das variáveis
(r, θ), as quais se designam por coordenadas polares do ponto P .
Figura 2.1: Coordenadas polares.
A mudança de variáveis para coordenadas polares é particularmente importante quando a região

de integração, vista em termos das coordenadas polares, tem fronteiras ao longo das quais r ou
θ é constante.
2.4 EXTREMOS
Proposição 2.2.5 O módulo do jacobiano da mudança de variáveis para coordenadas polares,

conforme (2.2.3), é |J| = r.
Exemplo 2.2.6 Fazendo a mudança de variáveis para coordenadas polares, calcule o integral
seguinte: Z Z
2 2
ex +y dx dy, onde D = (x, y) ∈ R2 : x2 + y 2 ≤ 1 .

D
2.2.4 Aplicações
Nesta secção vamos restringir-nos às aplicações geométricas, alertando apenas o leitor de exis-
tirem muitas outras aplicações, como o cálculo do centro de massa de um objecto, bem como
o cálculo dos seus momentos de inércia.
Seja D um domínio limitado de R2 . A área do domínio D é determinada através do integral

duplo, calculado sobre o domínio D, da função identicamente igual a 1. Ou seja,
Z Z
Área(D) = dx dy.
D
Exemplo 2.2.7 Usando integrais duplos, calcule a área do domínio D seguinte:
D = (x, y) ∈ R2 : −1 ≤ x ≤ 1, x ≤ y ≤ −x2 + x + 1 .

No caso do domínio D ser a reunião de dois domínios, digamos D1 e D2 , então:

Z Z Z Z
Área(D) = dx dy + dx dy.
D1 D2
Por outro lado, se o domínio D resultar da diferença de dois domínios, digamos D = D1 \ D2 ,

temos: Z Z Z Z
Área(D) = dx dy − dx dy.
D1 D2
Exemplo 2.2.8 Usando integrais duplos, calcule a área do domínio D seguinte:
D = (x, y) ∈ R2 : 1 ≤ x2 + y 2 ≤ 4 .

O outro exemplo de aplicação geométrica, é o cálculo de volumes de sólidos limitados por

superfícies z = f (x, y), onde f é uma função definida num domínio limitado D ⊂ R2 . Por
simplicidade de exposição, admitamos que f é uma função não negativa em D. O volume do
corpo
Ω = (x, y, z) ∈ R3 : (x, y) ∈ D, 0 ≤ z ≤ f (x, y)

é determinado pelo integral duplo de f , calculado sobre o domínio D, isto é,

Z Z
Volume(Ω) = f (x, y) dx dy .
D
2.4 EXTREMOS
No caso do corpo ser dado por
Ω = (x, y, z) ∈ R3 : (x, y) ∈ D,

g(x, y) ≤ z ≤ f (x, y) ,
então: Z Z
Volume(Ω) = [f (x, y) − g(x, y)] dx dy .
D
Exemplo 2.2.9 Usando integrais duplos, calcule o volume do corpo seguinte:
Ω = (x, y, z) ∈ R3 : 1 ≤ x2 + y 2 + z 2 ≤ 4, z ≥ 0 .

2.2.5 Teorema de Green

Nesta secção, vamos abordar a questão importante de como se relaciona o integral duplo com
o integral de linha. A resposta a esta questão, vai ser dada pelo Teorema de Green no plano.
Proposição 2.2.6 (Teorema de Green) Sejam f e g funções reais de duas variáveis reais,
com derivadas parciais contínuas num domínio aberto Ω de R2 . Sejam, ainda, C um caminho
simples, fechado e totalmente contido em D, onde D = C ∪ int(C) ⊂ Ω. Então:
Z Z
∂g ∂f
Z
f dx + g dy = − dx dy.
C D ∂x ∂y
SEM DEMONSTRAÇÃO: Ver, por exemplo, Dias Agudo, Volume 2, p. 103.
A equação integral do Teorema de Green, é muitas vezes designada por Fórmula de Riemann.
Exemplo 2.2.10 Verifique a Fórmula de Riemann do Teorema de Green para o campo vecto-
rial F(x, y) = (y 2 − 7y, 2xy + 2x) na circunferência x2 + y 2 = 1.
O Teorema de Green é uma ferramenta muito poderosa, não só em termos práticos de reso-
lução de alguns exercícios, assim como na demonstração de alguns resultados teóricos, alguns
deles já estabelecidos. Por exemplo, fica trivial mostrar que o integral de linha de um campo
conservativo calculado sobre um caminho simples fechado é zero.
Existem outras formas de apresentar o Teorema de Green no plano. Os dois resultados que
apresentamos a seguir num único enunciado, são casos particulares, no plano, de resultados
mais gerais.
Proposição 2.2.7 Sejam Ω, D e C como no enunciado da Proposição 2.2.6 e seja F um campo

vectorial cujas suas componentes, F1 e F2 , têm derivadas parciais contínuas em Ω. Admitamos
que o caminho C resulta de uma parametrização derivável e sejam n e t os vectores unitários
normal e tangente, respectivamente, à curva C. Então:
1. Teorema da Divergência
Z Z Z
div F dx dy = F · n ds;
D C
2.4 EXTREMOS
2. Teorema de Stokes Z Z Z
rot F · z dxdy = F · t ds,
D C
onde z é o vector normal unitário ao plano xy, isto é, z = (0, 0, 1).
SEM DEMONSTRAÇÃO: : É uma consequência imediata do Teorema de Green no plano dado

acima.
Como vimos, aquando do estudo dos integrais de linha, os integrais dos segundos membros das
equações integrais da proposição anterior, são integrais de linha relativamente ao comprimento
de arco, sendo ds o denominado elemento de arco. Designa-se por divergência de um campo
vectorial F(x, y) = (F1 (x, y), F2(x, y)), ao campo escalar definido por:
∂F1 ∂F2
div F = + .
∂x ∂y
Observe-se que qualquer vector, definido num espaço vectorial qualquer, se pode decompor em
duas componentes - uma normal e outra tangente. Em Análise Matemática, normal tem o
significado de perpendicular. Assim, para um campo vectorial qualquer F, tem-se:
F = Fn n + Ft t , onde Fn = F · n , Ft = F · t e n · t = 0 .
O vector unitário tangente à curva C é dado por

1 1
t= r′ (t) = p (x′ (t), y ′(t)) .
kr′ (t)k x′ (t)2 + y ′ (t)2
Por consequência, o vector unitário normal à mesma curva pode ser dado por
1
n= p (y ′ (t), −x′ (t)) .
x′ (t)2 + y ′ (t)2
Deste modo, o Teorema da Divergência escrito acima resume-se a

Z Z
∂F1 ∂F2
Z
+ dx dy = (F1 dy − F2 dx) ,
D ∂x ∂y C
onde usamos a igualdade formal n ds = (dy, −dx). Esta igualdade corresponde à Fórmula de
Riemann com f = −F2 e g = F1 . Por outro lado, sabemos que rot F denota o rotacional do
campo vectorial F = (F1 , F2 ) e, neste caso, reduz-se a
∂ F2 ∂ F1
rot F = − .
∂x ∂y
Do mesmo modo, o Teorema de Stokes escrito acima resume-se a
Z Z
∂ F2 ∂ F1
Z
− dxdy = (F1 dx + F2 dy) ,
D ∂x ∂y C
2.4 EXTREMOS
onde usamos, agora, a igualdade formal t ds = (dx, dy) e, na verdade, fizemos F = (F1 , F2 , 0) e

∂ F2 ∂ F1
rot F = 0, 0, − .
∂x ∂y
Verificamos pois que o que obtemos, é a Fórmula de Riemann com f = F1 e g = F2 .
A principal aplicação prática do Teorema de Green, é a possibilidade de podermos calcular a

área de determinado domínio plano, interior a um caminho, através do conhecimento apenas
do integral de linha.
Proposição 2.2.8 Sejam C e D, respectivamente, um caminho e um domínio nas condições

da Proposição 2.2.6. Então:
1
Z
Área(D) = x dy − y dx.
2 C
Exemplo 2.2.11 Usando o Teorema de Green no plano, calcule a área do círculo x2 + y 2 ≤ 1.
Se, na proposição anterior, fizermos uma mudança de variáveis para coordenadas polares, x =
r cos(θ) e y = rsen(θ), obtemos a fórmula seguinte:
1
Z
Área(D) = r 2 dθ.
2 C
Esta fórmula é útil para o cálculo da área interior a uma curva escrita em coordenadas polares.
2.2.6 Exercícios
1. Descreva as regiões de integração e calcule os integrais respectivos:
Z 2 Z 4 Z 2 Z x3
2 2
a) (x + y ) dx dy ; b) x, dy dx ;
0 0 1 x2
π
Z
2
Z y Z 1 Z |x|
c) sen x dx dy ; d) dy dx.
0 −y −1 0
2. Calcule os integrais das funções seguintes sobre as regiões indicadas:

a) f (x, y) = x cos(x + y) sobre o triângulo de vértices A = (0, 0), B = (π, 0) e C = (π, π).
b) f (x, y) = x sobre a região limitada por y = x2 e y = x3 .
c) f (x, y) = x − y sobre a região limitada por y = sen x e pelo eixo dos xx entre os pontos
x = 0 e x = π.
3. Inverta a ordem de integração dos integrais duplos seguintes:

Z 1Z 1 Z 2 Z 9−x
a) f (x, y) dx dy ; b) f (x, y) dy dx .
0 y 0 x2
2.4 EXTREMOS
4. Comece por verificar que, usando os métodos de primitivação conhecidos, é impossível

calcular os integrais seguintes pela ordem de integração indicadas. Inverta a ordem de
integração e calcule os integrais:
Z 4Z 2 Z 1Z 1p Z 8Z 2 √
1+x3 2
a) √
e dx dy ; b) y + 1 dy dx ; c) √
x4 + 1 dx dy .
0 y 0 x 0 3 y
5. Calcule os integrais seguintes usando as mudanças de variáveis indicadas:

Z Z
y−x
a) e y+x dx dy u = y − x, v = y + x,
D
onde D é a região limitada pela recta x + y = 2 e pelos dois eixos de coordenadas;

Z Z
b) (x2 + y 2) dx dy u = x + y, v = x − y,
D
onde D é a região limitada pelas rectas y = x, y = −x, y = x − 2 e y = 2 − x.
6. Calcule os integrais seguintes usando a mudança de variáveis para coordenadas polares:

Z Z
2 2 2
(x2 + y 2 )e(x +y ) dx dy , D = (x, y) ∈ R2 : x2 + y 2 ≤ 1 ;

a)
R
Z Z p
D = (x, y) ∈ R2 : 2 ≤ x2 + y 2 ≤ 4 .

b) x2 + y 2 dx dy ,
R
7. Calcule os integrais duplos seguintes usando a mudança de variáveis para coordenadas

polares:
Z 2 Z √4−y2 Z 3 Z √9−x2
3 2 2 2
a) (4 − x2 − y 2 ) 2 dx dy ; b) √
e−(x +y ) dy dx .
0 0 0 3x
8. Usando integrais duplos, calcule a área dos domínios planos seguintes:

a) D = {(x, y) ∈ R2 : 0 ≤ x ≤ π/4, senx ≤ y ≤ cos x};
b) D = {(x, y) ∈ R2 : −2 ≤ x ≤ x, −|x| ≤ y ≤ |x|};
c) D = {(x, y) ∈ R2 : 0 ≤ x ≤ 1, x2 − y 2 ≥ 0}.
9. Usando integrais duplos, calcule volume dos domínios espaciais seguintes:

a) D é o domínio do primeiro octante limitado pelas superfícies y = 1 − x2 e z = 1 − x2 ;
b) D é o domínio do primeiro octante interior ao cilindro x2 + z 2 = 4 e limitado pelo
plano x = y.
10. Usando o Teorema de Green, calcule os integrais de linha C F nas situações seguintes:
R
a) F(x, y) = (x2 ey , y 2 ex ) e C é o rectângulo de vértices A = (0, 0), B = (2, 0), C = (2, 3)

e D = (0, 3).
b) F(x, y) = (ey /x, ey ln x + 2x) e C é a curva dada por 1 + x2 ≤ y ≤ 2.
2.4 EXTREMOS
11. Usando o Teorema de Green, calcule os integrais de linha seguintes:

Z Z
2
a) (2x − y ) dx + (xy − 1) dy , b) (x4 − y 3) dx + (x3 + y 5) dy ,
C1 C2
onde C1 é o triângulo de vértices A = (0, 0), B = (2, 0) e C = (0, 1), e C2 é a circunferência

x2 + y 2 = 1.
12. Usando o Teorema de Green no plano, calcule as áreas interiores às curvas fechadas
parametrizadas por:
(−t, 1 − t2 ) −1 ≤ t < 1

sen t cos t
a) r(t) = , , t ∈ [0, 2π] ; b) r(t) =
3 2 (t − 2, 0) 1 ≤ t ≤ 3.
13. Usando o Teorema de Green, calcule as áreas das figuras geométricas limitadas pelas
curvas seguintes:
a) elipse x2 + 4y 2 = 4;
b) asteróide x2/3 + y 2/3 = 1.
2.3 Integral de superfície

2.3.1 Superfícies
Por superfície entende-se, de um dodo geral, como sendo o gráfico de uma função real de duas
variáveis reais e que é contínua em ambas as variáveis. Este modo de definir uma superfície
corresponde à denominada representação cartesiana da superfície. E, desta forma, dizemos
habitualmente que uma superfície está representada, de forma implícita, pela equação cartesiana
F (x, y, z) = 0.
No entanto, para definirmos o integral de superfície, esta noção de superfície não é satisfatória.
Tal como para as linhas que definidos aquando da secção sobre o integral de linha, vai ser
necessário definir o que se entende por representação paramétrica de uma superfície.
Definição 2.3.1 Seja D um subconjunto conexo3 de R2 . Chama-se superfície ao conjunto de

pontos de R3 cujas coordenadas cartesianas são dadas por:

 x = x(u, v)
y = y(u, v) ,
z = z(u, v)

onde x, y e z são funções contínuas das variáveis (u, v) ∈ D.
As superfícies são habitualmente denotadas pela letra S e podem ser escritas na forma vectorial
seguinte:
r(u, v) = (x(u, v), y(u, v), z(u, v));
onde r é uma aplicação vectorial definida de D ⊂ R2 em R3 .
3
2.4 EXTREMOS
Exemplo 2.3.1 Um cilindro de raio 2 definido pela equação cartesiana x2 + y 2 = 4 e com

−3 ≤ z ≤ 4 é representado parametricamente por:

 x = 2 cos(u)
y = 2sen(u) , onde (u, v) ∈ [0, 2π] × [−3, 4].
z=v

Exemplo 2.3.2 A superfície de uma esfera de raio 3 definida pela equação cartesiana x2 +
y 2 + z 2 = 9 é representada parametricamente por:

 x = 3 cos(v) cos(u) h π πi
y = 3 cos(v)sen(u) , onde (u, v) ∈ [0, 2π] × − , .
2 2
z = 3sen(v)

Exemplo 2.3.3 1. Determine uma equação cartesiana que represente a superfície definida
parametricamente por 
 x = u + 2v
y = 2u + v + 1 .
z =u+v+2

2. Determine uma representação paramétrica que represente a superfície definida pela equação
cartesiana 3x + 4y + 6z = 24.
2.3.2 Plano tangente e normal a uma superfície

Comecemos por recordar as noções de plano tangente e normal a uma superfície S que é
representada, de forma implícita, pela equação cartesiana:
F : R3 −→ R
F (x, y, z) = 0, onde
(x, y, z) 7→ F (x, y, z) .
Sendo P = (a, b, c) um ponto regular pertencente à superfície S, sabemos que:
(i) o vector ∇F (P ) é perpendicular ao plano tangente à superfície S no ponto P .
(ii) a equação cartesiana do plano tangente à superfície S no ponto P é dada por
∂F ∂F ∂F
(P )(x − a) + (P )(y − b) + (P )(z − c) = 0 .
∂x ∂y ∂z
Observemos que o ponto P é o denominado ponto de contacto da superfície com o plano
tangente. O ponto P diz-se regular, se existem as derivadas parciais Fx , Fy e Fz , e são
contínuas em P , e se, além do mais, pelo menos, uma das derivadas parciais, Fx (P ), Fy (P ) ou
Fz (P ) é não nula. Se todas as derivadas parciais Fx (P ), Fy (P ) e Fz (P ) forem nulas, ou não
existir, no ponto P , pelo menos, uma das derivadas parciais, dizemos que P = (a, b, c) é um
ponto singular da superfície S. O vector normal à superfície S no ponto P = (a, b, c) é
dado por

1 1 ∂F ∂F ∂F
n= ∇F (P ) = (P ), (P ), (P ) .
k∇F (P )k ∂ x (P ), ∂∂ Fy (P ), ∂∂ Fz (P ) ∂ x
∂F ∂y ∂z
A ideia agora, é a de definir, de forma equivalente, o plano tangente e a normal à superfície S

no ponto de contacto P = (a, b, c), mas usando a representação paramétrica da superfície.
2.4 EXTREMOS
Proposição 2.3.1 Seja S uma superfície parametrizada por:


 x = x(u, v)
y = y(u, v) ,
z = z(u, v)

tal, que em cada ponto P = (a, b, c) de S, os vectores derivadas parciais

∂r ∂x ∂y ∂z ∂r ∂x ∂y ∂z
= , , e = , ,
∂u ∂u ∂u ∂u ∂v ∂v ∂v ∂v
são contínuos e são linearmente independentes. Então, em cada ponto de S, tem-se:
(i) existe uma única normal cuja direcção depende de forma contínua dos pontos de S e que é
dada por:
1 ∂r ∂r
n = ∂r ∂r × ; (2.3.4)
× ∂u ∂v
∂u ∂v
(ii) existe um único plano tangente à superfície S no ponto P = (a, b, c) que é gerado pelos
∂r ∂r
vectores ∂u e ∂v e cuja equação é dada por:

∂r ∂r
(x − a, y − b, z − c) · (u0 , v0 ) × (u0 , v0 ) = 0;
∂u ∂v
onde (u0 , v0 ) é um ponto em D tal que r(u0, v0 ) = (a, b, c) ≡ P .
Figura 2.2: Plano tangente e recta normal.

∂r ∂r
Como sabemos da Álgebra Linear, os vectores ∂u
e
são linearmente independentes se e só se
∂v
 
e1 e2 e3
 
∂r ∂r ∂r ∂r  ∂x ∂y ∂z 
× 6= 0, onde × ≡ det 
 ∂u ∂u ∂u  .

∂u ∂v ∂u ∂v  
∂x ∂y ∂z
∂v ∂v ∂v
Exemplo 2.3.4 Usando a proposição anterior, determine a normal √ e√
a equação do plano tan-
gente à superfície S, definida por y 2 + z 2 = 4, no ponto P = (3, − 2, 2).
2.4 EXTREMOS
2.3.3 Integrais de superfície

Neste momento, estamos em condições de poder definir o integral de superfície. No que se segue,
diremos que S é uma superfície de classe C 1 se for parametrizada por uma função vectorial
r(u, v) = (x(u, v), y(u, v), z(u, v)) tal que os vectores derivadas parciais ∂u
∂r
e ∂v
∂r
existirem e
forem contínuos no subdomínio de R onde se define a parametrização.
2
Definição 2.3.2 Sejam f um campo escalar definido e contínuo num subconjunto Ω de R3 , S,

de classe C 1 , uma superfície parametrizada por r(u, v) = (x(u, v), y(u, v), z(u, v)), com (u, v) ∈
D e tal que r (D) ⊆ Ω. Define-se o integral de superfície do campo escalar f sobre a
superfície S, relativamente à área S, por:

∂r ∂r
Z Z Z Z
f (x, y, z) dS = f (r(u, v))
× du dv.
S D ∂u ∂v
Observe-se que esta noção de integral de superfície é a versão análoga para superfícies do
integral de linha relativamente ao comprimento de arco. Por outro lado, note-se que ∂u
∂r ∂r
× ∂v
dá-nos a recta normal à superfície S. O elemento

∂r ∂r
∂u × ∂v du dv

designa-se por elemento de área e, é comum denotá-lo por dS.
Exemplo 2.3.5 Sendo f (x, y, z) = x + y + z, calcule o integral, relativamente à área, sobre a

superfície definida pelo rectângulo 0 ≤ x ≤ 1, 0 ≤ y ≤ 1 e z = 12 .
A proposição seguinte permite-nos calcular o integral de superfície de uma função escalar no

caso de a superfície puder ser escrita, de modo explícito, na forma cartesiana.
Proposição 2.3.2 Sejam f um campo escalar definido e contínuo num subconjunto Ω de R3 ,

S uma superfície definida pela equação cartesiana z = ϕ(x, y), com (x, y) ∈ D e ϕ é uma função
de classe C 1 em D. Então, o integral de f , relativamente à área, sobre a superfície S é dado
por: s 2 2
∂ϕ ∂ϕ
Z Z Z Z
f (x, y, z) dS = f (x, y, ϕ(x, y)) 1 + + dx dy.
S D ∂x ∂y
No caso da superfície S ser definida por y = ϕ(x, z) ou x = ϕ(y, z), tem-se, respectivamente,
s 2 2
∂ϕ ∂ϕ
Z Z Z Z
f (x, y, z) dS = f (x, ϕ(x, z), z) 1 + + dx dz
S D ∂x ∂z
ou s 2 2
∂ϕ ∂ϕ
Z Z Z Z
f (x, y, z) dS = f (ϕ(y, z), y, z) 1 + + dy dz.
S D ∂y ∂z
2.4 EXTREMOS
Exemplo 2.3.6 Usando a proposição anterior, calcule o integral de f (x, y, z) = x + y + z sobre

a superfície do cubo 0 ≤ x ≤ 1, 0 ≤ y ≤ 1 e 0 ≤ z ≤ 1.
Usando a noção de integral de superfície, relativamente à área S de uma função escalar, podemos
definir o que se entende por integral de superfície, relativamente à área S de uma função
vectorial.
Definição 2.3.3 Sejam F um campo vectorial definido e contínuo num subconjunto Ω de R3 , S
uma superfície, de classe C 1 , parametrizada por r(u, v) = (x(u, v), y(u, v), z(u, v)), com (u, v) ∈
D e tal que r (D) ⊆ Ω. Define-se o integral de superfície do campo vectorial F sobre a
superfície S, relativamente à área S, por:

∂r ∂r
Z Z Z Z Z Z

F≡ F · n dS = F(r(u, v)) · n(u, v) ×
∂u ∂v du dv.

S S D
A definição anterior é uma caso particular da Definição 2.3.5, onde se considera

f (r(u, v)) = F(r(u, v)) · n(u, v).
Exemplo 2.3.7 Calcule o integral do campo vectorial F(x, y, z) = (y, −x, z 2 ) sobre o parabolóide
z = x2 + y 2, com 0 ≤ z ≤ 1.
Tal como para o integral de linha, no caso particular da função vectorial F ter apenas uma
componente não nula, obtemos uma nova noção de integral de superfície. Designaremos este
novo integral por integral de superfície relativamente a duas das três variáveis x, y e z.
Definição 2.3.4 Sejam f um campo escalar definido e contínuo num subconjunto Ω de R3 , S
uma superfície de classe C 1 parametrizada por r(u, v) = (x(u, v), y(u, v), z(u, v)), com (u, v) ∈
D e tal que r (D) ⊆ Ω. Define-se o integral de superfície do campo escalar F relativamente
às variáveis (x, y), (x, z) e (y, z), respectivamente, por:
∂x ∂x
∂(x, y) ∂(x, y)
Z Z Z Z
f dx dy = f (r(u, v)) du dv, onde = det ∂u ∂y
∂v
∂y ;
S D ∂(u, v) ∂(u, v) ∂u ∂v
∂x ∂x
∂(x, z) ∂(x, z)
Z Z Z Z
f dx dz = f (r(u, v)) du dv, onde = det ∂u ∂z
∂v
∂z ;
S D ∂(u, v) ∂(u, v) ∂u ∂v
∂y ∂y
∂(y, z) ∂(y, z)
Z Z Z Z
f dy dz = f (r(u, v)) du dv, onde = det ∂u ∂z
∂v
∂z .
S D ∂(u, v) ∂(u, v) ∂u ∂v
Como referimos acima, estas noções de integral de superfície correspondem a considerar na

Definição 2.3.3 F = (f, 0, 0), F = (0, f, 0) e F = (0, 0, f ), respectivamente. Observe-se que, tal
como no integral duplo, dx dy, dx dz e dy dz não indicam nenhuma ordem de integração. Resulta
das propriedades dos determinantes que, para os integrais de superfície da Definição 2.3.4, se
tenha:
Z Z Z Z Z Z Z Z Z Z Z Z
f dx dy = − f dy dx; f dx dz = − f dz dy; f dy dz = − f dz dy.
S S S S S S
Tendo em conta que estamos perante integrais de superfície, podemos escrever estas igualdades
na forma abreviada seguinte:
dx dy = −dy dx; dx dz = −dz dy; dy dz = −dz dy.
2.4 EXTREMOS
Exemplo 2.3.8 Calcule os integrais seguintes sobre a superfície esférica x2 + y 2 + z 2 = 1

correspondente ao hemisfério superior:
Z Z Z Z Z Z
a) z dx dy; b) y dx dz; c) x dy dz.
S S S
Usando estas noções de integral de superfície relativamente a duas das três variáveis x, y e z,
podemos então reescrever a igualdade integral da Definição 2.3.3.
Proposição 2.3.3 Nas condições da Definição 2.3.3, escrevendo o campo vectorial na forma
F(x, y, z) = (F1 (x, y, z), F2 (x, y, z), F3(x, y, z)), temos:
Z Z Z Z
F · n dS = F1 dy dz + F2 dz dx + F3 dx dy.
S S
Exemplo 2.3.9 Calcule o integral seguinte sobre a superfície esférica x2 + y 2 + z 2 = 1 corre-

spondente ao hemisfério superior:
Z Z
x dy dz + y dz dx + z dx dy.
S
2.3.4 Outras noções de integral de superfície

Nesta secção vamos generalizar algumas das definições e alguns resultados aposentados na
secção anterior. Comecemos por observar que em cada ponto de uma superfície S, supostamente
de classe C 1 , existem duas normais unitárias com sentidos opostos. Portanto, em determinadas
situações, convém ter presente em que face estamos a considerar o integral de superfície.
Pelo acima exposto, vejamos que a noção de integral de superfície de um campo vectorial
vem um pouco modificada. Pela definição da normal n dada em (2.3.4), podemos escrever a
identidade integral da Definição 2.3.3 na forma seguinte:

∂r ∂r
Z Z Z Z Z Z
F≡ F · n dS = ± F(r(u, v)) · × du dv;
S S D ∂u ∂v
sendo o sinal desta identidade integral escolhido de acordo com a face da superfície que se
considera.
No caso da superfície S puder ser escrita na forma cartesiana através da equação
F (x, y, z) = 0,
a proposição seguinte permite-nos escrever a identidade integral da Definição 2.3.3 numa forma
que faz intervir a função escalar F que define a superfície.
2.4 EXTREMOS
Proposição 2.3.4 Nas condições da Definição 2.3.3, supondo que a superfície S é definida
pela equação cartesiana F (x, y, z) = 0, temos:
Z Z Z Z Z Z
F≡ F · n dS = (F1 cos(α) + F2 cos(β) + F3 cos(γ)) dS;
S S S
onde
1 ∂F 1 ∂F 1 ∂F
cos(α) = , cos(β) = , cos(γ) =
D ∂x D ∂y D ∂z
e s 2 2 2
∂F ∂F ∂F
D=± + + .
∂x ∂y ∂z
Os cossenos cos(α), cos(β) e cos(γ), são designados cossenos directores da normal n, isto
é,
n = (cos(α), cos(β), cos(γ)) .
A escolha do sinal de D deve ser feita de acordo com a face da superfície que se considera.
Exemplo 2.3.10 Usando a proposição anterior, calcule o integral de F(x, y, z) = (x, y, z) sobre
a face exterior da superfície esférica x2 + y 2 + z 2 = 1 correspondente ao hemisfério superior.
Pelas Proposições 2.3.3 e 2.3.4, e sempre que faça sentido, podemos escrever a igualdade integral
seguinte:
Z Z Z Z
(F1 cos(α) + F2 cos(β) + F3 cos(γ)) dS = F1 dy dz + F2 dz dx + F3 dx dy.
S S
2.3.5 Aplicações
A principal aplicação dos integrais de superfície é de cariz geométrico e corresponde ao cálculo
da área de uma superfície.
Definição 2.3.5 Seja S uma superfície de classe C 1 parametrizada por r(u, v) = (x(u, v),
y(u, v), z(u, v)), com (u, v) ∈ D e tal que r (D) ⊆ Ω. Define-se a área da superfície de S por:
Z Z
∂r ∂r
Z Z
A(S) = dS ≡ ∂u × ∂v du dv.

S D
Exemplo 2.3.11 Calcule a área da superfície lateral do cilindro x2 + y 2 = 4, com 0 ≤ z ≤ 10.
No caso da superfície puder ser escrita, de modo explícito, na forma cartesiana z = ϕ(x, y),
com (x, y) ∈ D e ϕ uma função de classe C 1 em D, então, de acordo com a Proposição 2.3.2, a
área da superfície S pode ser calculada da forma seguinte:
s 2 2
∂ϕ ∂ϕ
Z Z
A(S) = 1+ + dx dy.
D ∂x ∂y
2.4 EXTREMOS
De modo análogo, se a superfície S puder ser escrita na forma y = ϕ(x, z) ou x = ϕ(y, z),
tem-se, respectivamente,
s 2 2 s 2 2
∂ϕ ∂ϕ ∂ϕ ∂ϕ
Z Z Z Z
A(S) = 1+ + dx dz, A(S) = 1+ + dy dz.
D ∂x ∂z D ∂y ∂z
Exemplo 2.3.12 Usando a equação cartesiana da esfera x2 + y 2 + z 2 = 1, calcule a sua área

de superfície.
Os integrais de superfície têm também muita aplicação em Mecânica dos Fluidos. Nesta área
da Mecânica é importante determinar a quantidade de massa de fluido que passa através de
uma superfície S por unidade de tempo. Esta quantidade é designada em Mecânica dos Fluidos
por fluxo ou vazão de um fluido.
Definição 2.3.6 Seja v o campo de velocidades (vector) associado a um fluido. O fluxo, ou

vazão, do fluido que atravessa uma superfície S é dado por:

∂r ∂r
Z Z Z Z
F (S) = v · n dS ≡ v(r(u, v)) · n(u, v)
× du dv.
S D ∂u ∂v
Exemplo 2.3.13 Determine o fluxo de água que atravessa o cilindro parabólico y = x2 , com
0 ≤ x ≤ 2, 0 ≤ z ≤ 3, sabendo que o seu campo de velocidades é dado por v = (3z 2 , 6, 6xz)
ms−1 .
2.3.6 Exercícios
1. Determine uma equação cartesiana que represente as superfícies definidas parametrica-
mente por
 
 x=1  x=√u
a) y = u , 0 ≤ u ≤ 1, 0 ≤ v ≤ 1; b) y = 1 − u2 − v 2 , u2 + v 2 ≤ 1.
z=v z=v
 
2. Determine uma representação paramétrica para as superfícies definidas pelas equações

cartesianas seguintes:
a) x2 +p(y − 1)2 + (z + 2)2 = 4 (esfera);
b) z = x2 + y 2 , com 0 ≤ z ≤ 3 (cone);
c) 9x2 + 4y 2 = 36 (cilindro elíptico).
3. Usando as representações paramétricas das superfícies em causa, determine as normais e

as equações dos planos tangentes às superfícies seguintes nos pontos P indicados:
a) 2x + 3y − 5z = 7, P √ = (3,
√2, 1);
b) x + y = 4, P = ( 2, − 2, 0);
2 2
c) x2 + y 2 + z 2 = 25, P = (3, 4, 0).
4. Sendo f (x, y, z) = x+y +z, calcule o integral de f , relativamente à área, sobre a superfície
do cubo 0 ≤ x ≤ 1, 0 ≤ y ≤ 1 e 0 ≤ z ≤ 1.
2.4 EXTREMOS
5. Calcule o integral de superfície, relativamente à área, seguinte:

Z Z
(x2 + y 2) dS, S : x2 + y 2 + z 2 = 1.
S
6. Usando as equações cartesianas das superfícies indicadas, calcule os integrais de superfície

seguintes:
1
Z Z
a) (x + y + z) dS, S : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, z = ;
S 2
Z Z
b) (x2 + y 2 ) dS, S : x2 + y 2 + z 2 = 1.
S
7. Calcule os integrais dos campos vectoriais seguintes sobre as superfícies indicadas:

a) F(x, y, z) = (x2 , 0, 3y 2), S é a porção do plano x + y + z = 1 que está no primeiro
octante.
b) F(x, y, z) = (0, y, −z), S é dada pelas superfícies do parabolóide y = x2 + z 2 , com
0 ≤ y ≤ 1, e do disco x2 + z 2 ≤ 1, com y = 1. (Resposta: 0)
8. Calcule os integrais seguintes:
Z Z
a) yz dy dz+xz dz dx+xy dx dy , S é a face externa de x + y + z = 1 no 1o octante ;
S
Z Z
b) xz dy dz + yz dz dx + x2 dx dy , S é a parte exterior da esfera x2 + y 2 + z 2 = 4 .
S
9. Calcule as área das superfícies S seguintes:

a) S é a superfície face externa do tetraedro x + y + z = 1 que está no 1o octante ;
b) S é a parte do parabolóide z = x2 + y 2 interior ao cilindro x2 + y 2 = 4 ;
c) S é a superfície total (incluindo a base) do cone z 2 = x2 + y 2 , com 0 ≤ z ≤ 4.
10. Esboce as superfícies descritas pelas equações paramétricas (ou equação vectorial) seguintes
e calcule as áreas respectivas:
a) r(u, v) = (cos(u), sen(u), v) , com u ∈ [0, 2π] e v ∈ [−1, 1];

√
b) r(u, v) = (u cos(v), u sen(v), u2), com 0 ≤ u ≤ 1, 0 ≤ v ≤ 2π , R: π(5 5 − 1)/6;
√
c) r(u, v) = (u, v, u2 + v 2 ), com 0 ≤ u2 + v 2 ≤ 1 , R: π(5 5 − 1)/6.
11. Considere a superfície S de R3 parametrizada por
r(u, v) = (cos(u), sen(u), v) , com u ∈ [0, 2π] e v ∈ (−1, 1).
a) Esboce a superfície S.
b) Usando integrais de superfície, calcule a área da superfície S.
c) Tomando S como a parte exterior da superfície dada, calcule o integral de superfície
Z Z
x2 dy dz + y 2 dz dx + z 2 dx dy .
S
2.4 EXTREMOS
12. Seja
q
E(x, y, z) = (x, y, z)
(x2
+ + z 2 )3/2
y2
o campo eléctrico criado por uma carga q localizada na origem. Calcule o fluxo de E
através da superfície esférica de raio 2 e centrada na origem, com a normal apontando
para fora da esfera. Observe que este fluxo não depende do raio da esfera.
(Resposta: 4πq)
2.4 Integral triplo

Depois de introduzido o integral duplo, torna-se mais fácil estender as noções de integração para
funções escalares de três variáveis. Portanto, as regiões de integração são agora subconjuntos
de R3 . Vamos, primeiro, considerar regiões de integração paralelipípedicas e, depois, iremos
considerar regiões mais gerais cujas fronteiras são limitadas por superfícies não propriamente
planares.
2.4.1 Integral de Riemann

De igual modo como fizemos para o integral duplo, comecemos por estender a noção de partição
para um paralelipípedo (limitado) contido em R3 .
Definição 2.4.1 Sejam a1 , a2 , b1 , b2 , c1 e c2 números reais tais que a1 < a2 , b1 < b2 , c1 < c2 .
Consideremos o paralelipípedo
P = (x, y, z) ∈ R3 : a1 ≤ x ≤ a2 , b1 ≤ y ≤ b2 , c1 ≤ z ≤ c2

e as partições dos intervalos [a1 , a2 ], [b1 , b2 ] e [c1 , c2 ], definidas, respectivamente, por
P1 : a1 = x0 < x1 < · · · < xn−1 < xn = a2 , P2 : b1 = y0 < y1 < · · · < ym−1 < ym = b2
P3 : c1 = z0 < z1 < · · · < zl−1 < zl = c2 , onde m, n e l são números naturais arbitrários.
Designa-se por partição do paralelipípedo P ao conjunto seguinte
P = {(xi , yj , zk ) ∈ P : 0 ≤ i ≤ n, 0 ≤ j ≤ m, 0 ≤ k ≤ l} .
Uma partição P do paralelipípedo P , tal como a definida acima, determina mnl paralelipípedos
contidos em P :
Pijl = (xi , yj , zk ) ∈ R3 : xi−1 ≤ x ≤ xi , yi−1 ≤ y ≤ yi, zk−1 ≤ x ≤ zk , ,

cujos volumes são dados por
△xi △yj △zk = (xi − xi−1 )(yj − yj−1)(zk − zk−1 ).
Em termos objectivos, a colecção formada por estes paralelipípedos Pijl , é o que verdadeira-
mente constitui a partição P.
2.4 EXTREMOS
Definição 2.4.2 Seja f uma função definida num paralelipípedo P ⊂ R3 . Designamos por
soma de Riemann da função f no paralelipípedo P à quantidade seguinte
n, m, l l X
m X
n
X X
f (x∗ijl )△xi △yj △zl ≡ f (x∗ijl )△xi △yj △zl
i=1, j=1, k=1 k=1 j=1 i=1
≡ f (x∗111 )△x1 △y1 △zl + · · · + f (x∗nml )△xn △ym △zl ,

onde x∗ijl são pontos seleccionados aleatoriamente nos paralelipípedos respectivos Pijl .
Para a noção de qualquer integral, interessa-nos que as partições sejam muito finas. Definimos
a quantidade que define a finura de dada partição P de um paralelipípedo P ⊂ R3 por
q
|P| = max (△xi )2 + (△yj )2 + (△zk )2 .
i, j, k
No caso do integral triplo, |P| corresponde ao comprimento da maior diagonal de todos os

paralelipípedos Pijk considerados na partição P.
Definição 2.4.3 Sejam f uma função definida num paralelipípedo P ⊂ R3 e
P = {(xi , yj , zk ) ∈ P : 0 ≤ i ≤ n, 0 ≤ j ≤ m, 0 ≤ k ≤ l}
uma partição arbitrária de P . Diz-se que a função f é integrável (à Riemann) no par-

alelipípedo P , se existir (e for finito) o limite seguinte:
l X
X m X
n
lim f (x∗ij )△xi △yj △zk ,
|P|→0
k=1 j=1 i=1
independentemente de como a partição P do paralelipípedo P é formada, ou de como os pontos

x∗ijl pertencentes aos subparalelipípedo Pijl são escolhidos.
No caso de existir, o limite da definição anterior designa-se por integral da função f e denota-
se por uma das formas seguintes:
Z Z Z Z Z Z Z Z Z
f (x, y, z) dx dy dz, f (x, y, z) dx dz dy, f (x, y, z) dy dx dz,
P P P
Z Z Z Z Z Z Z Z Z
f (x, y, z) dy dz dx, f (x, y, z) dz dx dy, f (x, y, z) dz dy dx;
P P P
onde dx dy dz, dx dz dy, dy dx dz, dy dz dx, dz dx dy ou dz dy dx indica um elemento de volume
ao qual ainda não está subjacente nenhuma ordem de integração.
Esta noção de função integrável que acabamos de introduzir, estende-se a qualquer função
definida num conjunto limitado Ω ⊂ R3 que não seja propriamente um paralelipípedo. Apenas
temos de considerar um paralelipípedo P que contenha Ω e aí fazer a análise anterior. Neste
caso, para a definição precedente fazer sentido, fixamos o valor de f (x∗ijl ) igual a zero quando
x∗ijl não pertencer a Ω.
2.4 EXTREMOS
Proposição 2.4.1 Seja Ω um subconjunto de R3 limitado e cuja fronteira é constituída pela

união de superfícies (contínuas). Consideremos uma função f definida em Ω. Se f é contínua
em Ω excepto, quanto muito, num conjunto de medida nula4 , então f é integrável em Ω.
A proposição anterior diz-nos, então, que todas as funções elementares de três variáveis que
conhecemos serão integráveis em domínios limitados de R3 .
As propriedades do integral de Riemann definido em R, serão válidas mutatis mutandis para o
integral triplo. Em particular, o integral triplo é um operador linear e satisfaz a propriedade
aditiva dos integrais.
2.4.2 Integral repetido

O processo de integração repetida do integral duplo generaliza-se a qualquer integral em dimen-
são superior, em particular também para o integral triplo. Podemos, assim, reduzir o cálculo
de um integral triplo a um integral duplo. Em particular, o Teorema de Fubbini (ver Secção 2)
generaliza-se de modo imediato ao integral triplo.
Proposição 2.4.2 Sejam f uma função contínua num conjunto limitado Ω ⊂ R3 e D ⊂ R2
um conjunto fechado e limitado.
1. Sejam g(x, y) e h(x, y) duas funções contínuas em D tais que, para quaisquer (x, y) ∈ D,
g(x, y) ≤ h(x, y), e seja Ω = {(x, y, z) ∈ R3 : g(x, y) ≤ z ≤ h(x, y), (x, y) ∈ D}. Então:
Z Z Z Z Z "Z h(x,y) #
f (x, y, z) dx dy dz = f (x, y, z) dz dx dy.
Ω D g(x,y)
2. Sejam g(x, z) e h(x, z) duas funções contínuas em D tais que, para quaisquer (x, z) ∈ D,
g(x, z) ≤ h(x, z), e seja Ω = {(x, y, z) ∈ R3 : g(x, z) ≤ y ≤ h(x, z), (x, z) ∈ D}. Então:
Z Z Z Z Z "Z h(x,z) #
f (x, y, z) dx dy dz = f (x, y, z) dy dx dz.
Ω D g(x,z)
3. Sejam g(y, z) e h(y, z) duas funções contínuas em D tais que, para quaisquer (y, z) ∈ D,
g(y, z) ≤ h(y, z), e seja Ω = {(x, y, z) ∈ R3 : g(y, z) ≤ x ≤ h(y, z), (y, z) ∈ D}. Então:
Z Z Z Z Z "Z h(y,z) #
f (x, y, z) dx dy dz = f (x, y, z) dx dy dz.
Ω D g(y,z)

Exemplo 2.4.1 Descreva a região de integração e calcule o respectivo integral:
Z 1Z 2Z 3
xy 2 z dx dy dz.
0 0 1
As considerações sobre o cálculo do integral repetido e inversão de ordem de integração feitas
para o integral duplo, podem ser adaptadas para o integral triplo.
4
2.4 EXTREMOS
2.4.3 Mudança de variáveis

Tal como acontece no cálculo de outros integrais, por vezes, também no cálculo do integral
triplo se torna necessário fazer mudança de variáveis para se poder calcular os integrais.
Proposição 2.4.3 (Mudança de variáveis) Sejam Ω e ∆ dois subconjuntos abertos de R3

limitados por funções diferenciáveis e seja
ϕ: ∆ −→ Ω
(u, v, w) 7→ (x, y, z) = ϕ(u, v, w)
uma função bijectiva com derivadas parciais contínuas. Então, se f (x, y, z) é uma função
integrável em ∆, temos:
Z Z Z Z Z Z
f (x, y, z) dx dy dz = f (ϕ(u, v, w)) |J| du dv dw ,
Ω ϕ−1 (Ω)
onde  
∂x ∂x ∂x
∂u ∂v ∂w
∂(x, y, z)  ∂y ∂y ∂y

J = det ≡ det  ∂u ∂v ∂w
.
∂(u, v, w)  
∂z ∂z ∂z
∂u ∂v ∂w
SEM DEMONSTRAÇÃO: Ver, por exemplo, Dias Agudo, Volume I, p. 258
Tendo em conta que ϕ é uma aplicação bijectiva, ϕ−1 (Ω) = ∆. Tal como no integral duplo,
a quantidade J continua a designar-se por jacobiano da mudança de variáveis (x, y, z) para
(u, v, w) e a matriz [∂(x, y, z)/∂(u, v, w)] é a denominada matriz jacobiana dessa transfor-
mação. Mais, tendo em conta que podemos escrever
(x, y, z) = ϕ(u, v, w) ≡ (ϕ1 (u, v, w), ϕ2(u, v, w), ϕ3(u, v, w)),
a matriz jacobiana pode ser escrita na forma seguinte:

 ∂ϕ ∂ϕ ∂ϕ1

1 1
∂u ∂v ∂w
∂(x, y, z)  ∂ϕ ∂ϕ ∂ϕ2

≡ 2
 ∂u ∂v
2
∂w
.
∂(u, v, w) 
∂ϕ3 ∂ϕ3 ∂ϕ3
∂u ∂v ∂w
Exemplo 2.4.2 Fazendo a mudança de variáveis u = x + y − z e v = x + 2y + z e w = z,

calcule o integral seguinte:
sen(x + y − z)
Z Z Z
dx dy dz ,
Ω x + 2y + z
onde n π o
Ω = (x, y, z) ∈ R3 : 1 ≤ x + 2y + z ≤ 2, 0 ≤ x + y − z ≤ , 0 ≤ z ≤ 1 .
4
2.4 EXTREMOS
No integral triplo, usamos com frequência uma mudança de variáveis correspondente à transfor-
mação para coordenadas polares do integral duplo. No integral triplo esta mudança de variáveis,
recebe o nome de transformação para coordenadas cilíndricas. Usando a trigonometria no es-
paço, qualquer ponto P de coordenadas paralelipípedicas (x, y, z) pode ser escrito em termos
das coordenadas cilíndricas: 
 x = r cos(θ)
y = r sen(θ) ;
z=z

onde r é a distância radial do ponto P à origem do referencial (x, y, z) = (0, 0, 0), θ ∈ [0, 2π] e z
é a cota do ponto P . θ é o ângulo formado pelo semi-eixo positivo dos xx e pela semi-recta com
origem em (x, y, z) = (0, 0, 0) e que passa pela projecção do ponto P no plano xy, e é medido
a partir do semi-eixo positivo dos xx até 360o .
Figura 2.3: Coordenadas cilíndricas.
Exemplo 2.4.3 Calcule o integral seguinte fazendo mudança de variáveis para coordenadas
cilíndricas:
Z Z Z
(x2 + y 2) dx dy dz , Ω = (x, y, z) ∈ R3 : x2 + y 2 ≤ 2z, z ≤ 2 .

Ω
Tal como as coordenadas polares no caso do integral duplo, a mudança de variáveis para coorde-
nadas cilíndricas no integral triplo é particularmente importante quando a região de integração
tem fronteiras ao longo das quais r ou θ é constante.
Outra mudança de variáveis muito comum em integrais triplos, é a transformação para coor-
denadas esféricas: 
 x = r cos(φ) cos(θ)
y = r cos(φ) sen(θ) ; (2.4.5)
z = r sen(φ)

onde r é a distância radial do ponto P à origem do referencial (x, y, z) = (0, 0, 0), θ ∈ [0, 2π],
φ ∈ [−π/2, π/2] e z é a cota do ponto P . θ é o ângulo formado pelo semi-eixo positivo dos xx
e pela semi-recta com origem em (x, y, z) = (0, 0, 0) e que passa pela projecção do ponto P no
plano xy, e é medido a partir do semi-eixo positivo dos xx até 360o . φ é o ângulo formado pelo
2.4 EXTREMOS
plano z = 0 e pela semi-recta com origem em (x, y, z) = (0, 0, 0) e que passa pelo ponto P , e é
medido a partir do plano z = 0 até ±90o .
Podemos introduzir as coordenadas esféricas equivalentes da forma seguinte:


 x = r sen(φ) cos(θ)
y = r sen(φ) sen(θ) ; (2.4.6)
z = r cos(φ)

onde, agora, φ é o ângulo formado pelo semi-eixo positivo dos zz e pela semi-recta com origem
em (x, y, z) = (0, 0, 0) e que passa pelo ponto P , e é medido a partir do semi-eixo positivo dos
zz até ±180o .
Figura 2.4: Coordenadas esféricas 2.4.5 (esquerda) e 2.4.6 (direita).
Proposição 2.4.4 Os jacobianos das transformações para coordenadas esféricas seguintes, são
os indicados:
 
 x = r cos(φ) cos(θ)  x = r sen(φ) cos(θ)
a) y = r cos(φ) sen(θ) ⇒ J = ±r cos(φ); b) y = r sen(φ) sen(θ) ⇒ J = ±r 2 sen(φ).
2
z = r sen(φ) z = r cos(φ)
 
A escolha do sinal + ou − está relacionada com a ordem pela qual escrevemos as colunas da
matriz jacobiana. Se escolhermos a ordem (r, φ, θ), virá − em a) e + em b). Se a ordem for
(r, θ, φ), então os sinais vêm trocados. As duas transformações para coordenadas esféricas são
equivalentes, temos apenas de ter o cuidado que o domínio de variação do ângulo φ é diferente
e que, na substituição do integral a calcular, deverá aparecer |J|. Assim, se usarmos a forma
das coordenadas esféricas expressa em a) da Proposição 2.4.4, obtemos:
h π πi
2
|J| = r cos(φ), pois cos(φ) ≥ 0 para φ ∈ − , .
2 2
Se, porventura, usarmos as coordenadas esféricas de b), temos:
|J| = r 2 sen(φ), pois sen(φ) ≥ 0 para φ ∈ [0, π] .

2.4 EXTREMOS
Exemplo 2.4.4 Calcule o integral seguinte fazendo mudança de variáveis para coordenadas
esféricas:
Z Z Z p
Ω = (x, y, z) ∈ R3 : x2 + y 2 + z 2 ≤ 4 .

x2 + y 2 + z 2 dx dy dz ,
Ω
2.4.4 Aplicações
Nesta secção vamos apenas considerar a aplicação de integrais triplos ao cálculo de volumes de
corpos.
Seja Ω um domínio limitado de R3 . O volume do domínio Ω é determinado através do integral

triplo, calculado sobre o domínio Ω, da função identicamente igual a 1. Ou seja,
Z Z Z
Volume(Ω) = dx dy dz.
Ω
Exemplo 2.4.5 Usando integrais triplos, calcule o volume do paralelipípedo

√ com vértices
√ V1 =
(0, 0, √
0), V2 = (1, 0, 0),
√ V3 = (0, 1, 0), V4 = (1, 1, 0), V5 = (0, 0, 2), V6 = (1, 0, 2), V7 =
(0, 1, 2), V8 = (1, 1, 2).
Se o domínio Ω for a reunião de dois domínios, digamos Ω1 e Ω2 , então:

Z Z Z Z Z Z
Volume(Ω) = dx dy dz + dx dy dz.
Ω1 Ω2
Por outro lado, se o domínio Ω resultar da diferença de dois domínios, digamos Ω = Ω1 \ Ω2 ,

temos: Z Z Z Z Z Z
Volume(Ω) = dx dy dz − dx dy dz.
Ω1 Ω2
Exemplo 2.4.6 Usando integrais triplos, calcule o volume do domínio Ω compreendido entre
o rectângulo, do plano z = 0, com vértices A = (0, 0), B = (3, 0), C = (3, 2) e D = (0, 2), e a
superfície z = 4x2 + 9y 2.
2.4.5 Teorema de Green

Nesta secção apresentamos duas extensões para o espaço R3 , do Teorema de Green no plano R2 :
o Teorema da Divergência e o Teorema de Stokes. O Teorema da Divergência está associado
à divergência de um campo vectorial e o Teorema de Stokes está associado com o rotacional.
Comecemos por definir o que se entende por divergência de um campo vectorial e recordar a
noção de rotacional.
Definição 2.4.4 Seja F um campo vectorial de R3 em R3 de classe C 1 , tal que F(x, y, z) =

(F1 (x, y, z), F2 (x, y, z), F3 (x, y, z)). Chama-se divergência do campo vectorial F ao campo
escalar denotado por div F e definido por
∂ F1 ∂ F2 ∂ F3
div F = + + .
∂x ∂y ∂z
2.4 EXTREMOS
Por vezes, também se usa a notação ∇ · F para denotar a divergência de F, onde ∇ é um

operador diferencial e que é dado por:

∂ ∂ ∂
∇= , , .
∂x ∂y ∂z
Usando esta notação, podemos escrever

∂ ∂ ∂
div F = ∇ · F ≡ , , · (F1 , F2 , F3 ).
∂x ∂y ∂z
Uma interpretação física da divergência é a seguinte. Se F for o campo de velocidades de um
fluido que atravessa uma superfície na direcção da normal unitária, então div F dá-nos a taxa
de variação de massa de fluido (por unidade de massa e por unidade de tempo) em cada ponto
do domínio limitado pela superfície. Em particular, se div F = 0, esta taxa de variação é nula
e o fluido diz-se incompressível.
Definição 2.4.5 Seja F um campo vectorial de R3 em R3 de classe C 1 , tal que F(x, y, z) =

(F1 (x, y, z), F2 (x, y, z), F3 (x, y, z)). Chama-se rotacional do campo vectorial F ao campo vec-
torial denotado por rot F e definido por
 
e1 e2 e3
 ∂ ∂ ∂ 
 ∂x ∂y ∂z  ,
rot F = det  
F1 F2 F3
onde (e1 , e2 , e3 ) é a base canónica de R3 .
O rotacional pode, também, ser definido usando o operador diferencial ∇:

∂ ∂ ∂
rot F = ∇ × F ≡ , , × (F1 , F2 , F3 ).
∂x ∂y ∂z
À luz da Mecânica dos Fluidos, também podemos dar uma interpretação física do rotacional.
Se F é o campo de velocidades de um fluido, então rot F é a quantidade, ou rotação, pela qual
o fluido gira em torno de determinado ponto. Neste contexto, o rotacional é habitualmente
designado por vorticidade.
Proposição 2.4.5 Seja Ω uma região fechada e limitada de R3 cuja fronteira S é uma super-
fície orientada5 limitada por uma curva simples fechada C. Seja F um campo vectorial cujas
suas componentes, F1 , F2 , F3 têm derivadas parciais contínuas em Ω e suponhamos que a su-
perfície S resulta de uma parametrização derivável. Sejam, ainda, n o vector unitário normal
à superfície S e t o vector unitário tangente à curva C. Então:
1. Teorema da Divergência
Z Z Z Z Z
div F dx dy z = F · n dS;
Ω S
5
Grosso modo, uma superfície diz-se orientada, se tiver dois lados.
2.4 EXTREMOS
2. Teorema de Stokes Z Z Z
rot F · n dS = F · t ds.
S C
SEM DEMONSTRAÇÃO: Ver, por exemplo, Dias Agudo, Volume I, pp. 126 e 130.
Recordemos, que, se a superfície S for parametrizada por r(u, v), então o vector unitário normal
à superfície S é dado por:
∂r
∂u
× ∂∂ vr
n = ∂ r ∂ r
.
∂u
× ∂ v
Por outro lado, se a curva C for parametrizada por r(t), então o vector unitário tangente à
curva C é dado por:
r′ (t)
t= ′ .
kr (t)k
Mais, se usarmos a notação dos integrais de superfície e de linha introduzida anteriormente,
podemos escrever as identidades integrais da proposição anterior nas formas abreviadas seguintes
mais simples:
• Teorema da Divergência
Z Z Z Z Z
div F dx dy z = F;
Ω S
• Teorema de Stokes Z Z Z
rot F = F.
S C
Podemos, ainda, escrevê-las nas formas seguintes mais intuitivas:
• Teorema da Divergência
Z Z Z
∂ F1 ∂ F2 ∂ F3
Z Z
+ + dx dy dz = F1 dydz + F2 dzdx + F3 dxdy ;
Ω ∂x ∂y ∂z S
• Teorema de Stokes
Z Z
∂ F3 ∂ F2 ∂ F1 ∂ F3 ∂ F2 ∂ F1
− dydz + − dzdx + − dxdy
S ∂y ∂z ∂z ∂x ∂x ∂y
Z
= F1 dx + F2 dy + F3 dz .
C
O Teorema da Divergência é, por vezes, designado por Teorema de Gauss e, habitualmente,

é designado na literatura russa por Teorema de Ostrogradsky-Gauss. Aquando do estudo do
integral duplo, vimos como o então também Teorema da Divergência relacionava o integral
duplo com o integral de linha. O Teorema da Divergência desta secção relaciona o integral
triplo com o integral de superfície.
2.4 EXTREMOS
Exemplo 2.4.7 Usando o Teorema da Divergência, calcule o integral de superfície

Z Z
x3 dydz + x2 y dzdx + x2 z dxdy
S
onde S é a superfície lateral do cilindro x2 + y 2 = 4 e 0 ≤ z ≤ 1 (superfície lateral e bases).

O Teorema de Stokes relaciona o integral de superfície com o integral de linha e tem a inter-
pretação física seguinte. O fluxo do rotacional de F através de S coincide com circulação de F
ao longo da fronteira de S.
Exemplo 2.4.8 Usando o Teorema de Stokes, calcule o integral de superfície
Z Z
rot F , F(x, y, z) = (z − y, x + z, −x − y),
S
onde S é a superfície z = 4 − x2 − y 2 e 0 ≤ z ≤ 4.
O Teorema de Stokes é um resultado essencialmente do espaço. Aquele de que falamos na secção
do integral duplo, é um caso particular do desta secção, quando se considera uma superfície S
plana.
2.4.6 Exercícios
1. Descreva as regiões de integração e calcule os respectivos integrais:
Z 1 Z 1−x Z 1−x−y Z 1Z 1Z 1√
a) dz dy dx; b) 1 − z 2 dx dy dz;
0 0 0 0 0 0
Z 1Z xZ xy Z πZ sen(y) Z cos(y)
c) x3 y 2 z dz dy dx; d) (π − z) dx dz dy.
0 0 0 0 0 0
2. Calcule os integrais seguintes usando as mudanças de variáveis indicadas:

Z Z Z
y−x+z
a) e y+x−z dx dy dz , u = y − x + z, v = y + x − z, w = z ,
Ω
onde Ω = {(x, y, z) : R3 : 0 ≤ y − x + z ≤ y + x − z, 0 ≤ y + x − z ≤ 1, 1 ≤ z ≤ 2}.

cos(x + 2y + z)
Z Z Z
b) dx dy dz , u = x + 2y + z, v = x + y − z, w = z ,
Ω x+y−z
onde Ω = (x, y, z) ∈ R3 : 1 ≤ x + y − z ≤ 2, 0 ≤ x + 2y + z ≤ π4 , 0 ≤ z ≤ 1 .

3. Calcule os integrais seguintes usando a mudança de variáveis para coordenadas cilíndricas:

Z Z Z
a) xyz dx dy dz ,
Ω
n p p o
onde Ω = (x, y, z) ∈ R3 : −1 ≤ y ≤ 1, 0 ≤ x ≤ 1 − y 2 , x2 + y 2 ≤ z ≤ x2 + y 2 ;
Z Z Z
2 2
ex +y dx dy dz , onde Ω = (x, y, z) ∈ R3 : x2 + y 2 ≤ 1, 0 ≤ z ≤ 1 .

b)
Ω
2.4 EXTREMOS
4. Calcule os integrais seguintes usando a mudança de variáveis para coordenadas esféricas:

Z Z Z
a) (x2 + y 2 + z 2 ) dx dy dz , onde
Ω
n p p p o
Ω = (x, y, z) ∈ R3 : 0 ≤ y ≤ 3, 0 ≤ x ≤ 9 − y 2, x2 + y 2 ≤ z ≤ 16 − (x2 + y 2 ) ;
Z Z Z
onde Ω = (x, y, z) ∈ R3 : 1 ≤ x2 + y 2 + z 2 ≤ 4, z ≥ 0 .

b) z dx dy dz ,
Ω
5. Usando integrais triplos, calcule o volume dos domínios Ω seguintes,:

a) Ω = {(x, y, z) ∈ R3 : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1 − x, 0 ≤ z ≤ 1 − x − y};
b) Ω = {(x, y, z) ∈ R3 : x + y + z ≤ 2, x ≥ 0, y ≥ 0, z ≥ 0};
c) Ω = {(x, y, z) ∈ R3 : x2 + y 2 ≤ 2z, z ≤ 2};
d) Ω = {(x, y, z) ∈ R3 : 1 ≤ x2 + y 2 + z 2 ≤ 4, z ≥ 0}.
6. Considere o elipsóide
2
(y − 1)2

3 (x − 1)

2
Ω = (x, y, z) ∈ R : + + (z − 1) = 1 .
4 9
a) Mostre que o módulo do jacobiano da transformação indicada a seguir é |J| = 6r 2 sen(φ):

 x = 1 + 2rsen(φ) cos(θ)
y = 1 + 3rsen(φ)sen(θ)
z = 1 + r cos(φ) .

b) Usando integrais triplos, e tendo em conta a transformação anterior com θ ∈ [0, 2π] e
φ ∈ [0, π], calcule o volume do elipsóide Ω.
7. Verifique a igualdade do Teorema da Divergência no caso do cone x2 +y 2 ≤ z 2 e 0 ≤ z ≤ 2.
8. Usando o Teorema da Divergência, calcule o integral de superfície F nos casos
RR
S
seguintes:
a) F(x, y, z) = (ex , ey , ez ), S : −1 ≤ x ≤ 1, −1 ≤ y ≤ 1, −1 ≤ z ≤ 1;
b) F(x, y, z) = (x3 , y 3 , z 3 ), S : x2 + y 2 + z 2 = 9.
9. Usando o Teorema da Divergência, calcule os integrais de superfície seguintes:
Z Z Z Z
3 2 2
a) x dydz + x y dzdx + x z dxdy ; b) (z − y) dy dz + y 3 dz dx + 2z 3 dx dy ;
S1 S2
onde S1 é a superfície do cubo 0 ≤ x ≤ 1, 0 ≤ y ≤ 1 e 0 ≤ z ≤ 1 e S2 é a superfície

(total) do cilindro x2 + y 2 = 4 com −3 ≤ z ≤ 3.
10. Verifique a igualdade do Teorema de Stokes no caso do parabolóide S : z = 1 − (x2 + y 2 )
e z ≥ 0; C o bordo deste parabolóide, i.e., x2 + y 2 = 1 e F(x, y, z) = (y, z, x) .
11. Usando o Teorema de Stokes, calcule o integral de superfície rot F nos casos seguintes:
RR
S
a) F(x, y, z) = (z , 5x, 0), S : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, z = 1;
2
b) F(x, y, z) = (ez , ez sen(y), ez cos(y)), S : z = y 2 , 0 ≤ x ≤ 4, 0 ≤ y ≤ 2.

Bibliografia
[1] T.M. Apostol. Calculus. John Wiley & Sons, New York, 1969.
[2] W.M. Boyce, R.C. DiPrima. Calculus. John Wiley & Sons, New York, 1988.
[3] B. Demidovitch (sob a redacção). Problemas e Exercícios de Análise Matemática. Editora Mir,
Moscovo, 1987.
[4] F.R. Dias Agudo. Análise Real. Escolar Editora, Lisboa, 1989.
[5] H.L. Guiadoras. Um Curso de Cálculo. Livros Técnicos e Científicos, Rio de Janeiro, 2002.
[6] M. Krasnov, A. Kiselev, G. Makarenko e E. Shikin. Mathematical Analysis for Engineers. Mir
Publishers, Moscow, 1989.
[7] V. Kravchenko. Apontamentos das aulas teóricas de Análise Matemática III. Universidade do
Algarve, 2009.
[8] E. Kreyszig. Advanced Engineering Mathematics. John Wiley & Sons, New York, 1999.
[9] S. Lang. Calculus of Several Variables. Springer-Verlag, New York, 1987.
[10] J.E. Marsden e M.J. Hoffman. Elementary Classical Analysis. Second Edition. W.E. Freeman and
Company, New York, 1995.
[11] W. Rudin. Principles of Mathematical Analysis. McGraw-Hill International Editions, 1976.
[12] S. Samko. Apontamentos das aulas teóricas de Análise Matemática III. Universidade do Algarve,
2008.
[13] V. Zorich. Mathematical Analysis. Springer-Verlag, Berlin Heidelberg, 2004.
82

Apontamentos Análise Matemática-III 2008-2009

Enviado por

Direitos autorais:

Formatos disponíveis

Apontamentos Análise Matemática-III 2008-2009

Enviado por

Dados do documento

Descrição 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

Apontamentos Análise Matemática-III 2008-2009

Enviado por

Direitos autorais:

Formatos disponíveis

Apontamentos de Análise Matemática III

Hermenegildo Borges de Oliveira

1 Funções reais de várias variáveis 1

2 Cálculo integral vectorial 42

Funções reais de várias variáveis

Em particular, R2 é o conjunto de todos os pares (x, y) de números reais e R3 é o conjunto de

Em RN define-se a adição por

x + y = (x1 , . . . , xN ) + (y1 , . . . , yN ) = (x1 + y1 , . . . , xN + yN )

e a multiplicação por escalar por

αx = α(x1 , . . . , xN ) = (αx1 , . . . , αxN ),

Neste sentido, os elementos de RN podem ser denominados vectores. Sejam, agora, x1 , x2 , . . . ,

Por outro lado, diz-se que x1 , x2 , . . . , xk ∈ RN são vectores linearmente independentes,

e1 = (1, 0, . . . , 0), e2 = (0, 1, 0, . . . , 0), . . . , eN = (0, . . . , 0, 1).

Usando esta notação, podemos escrever, para qualquer vector x = (x1 , . . . , xN ) ∈ RN ,

Verifica-se que os vectores e1 , e2 , . . . , eN são linearmente independentes. As duas afirmações

1.1.2 Noções topológicas

Definição 1.1.3 Sejam x = (x1 , . . . , xN ) e y = (y1 , . . . , yN ) dois quaisquer elementos de RN .

1.1.3 Sucessões de vectores

Exercício exemplo 1.1.1 Indique as sucessões componentes da sucessão vectorial seguinte:

Definição 1.1.6 Uma sucessão vectorial un é limitada, se

Exercício exemplo 1.1.4 Usando a definição, mostre que

Desta proposição decorre que

1.1.4 Funções reais de várias variáveis

2. Determine o domínio de cada uma das funções seguintes e represente-o graficamente:

3. Para cada um dos conjuntos representados na questão anterior, indique:

4. Represente graficamente as linhas ou superfícies de nível correspondentes aos valores c

5. Determine os conjuntos de nível das funções seguintes:

c) f (x, y, z) = 2x2 + y 2 + 3z 2 ; d) f (x, y, z) = x2 + 2y 2 .

6. A partir das linhas de nível, esboce os gráficos das funções seguintes:

Definição 1.2.1 Sejam f : RN −→ R uma função escalar e a um ponto de acumulação do

∀ ε > 0 ∃ δ > 0 : ( ∀ x ∈ Df e kx − ak < δ) =⇒ |f (x) − b| < ε .

No caso de existir tal b ∈ R, escrevemos

Exercício exemplo 1.2.1 Calcule o limite seguinte:

Exercício exemplo 1.2.2 Calcule o limite seguinte:

Limites direccionais - este processo consiste em calcular os denominados limites direccionais

Exercício exemplo 1.2.3 Calcule os limites direccionais da função seguinte

no ponto (x, y) = (0, 0).

Exercício exemplo 1.2.4 Calcule os limites repetidos da função seguinte

no ponto (x, y) = (0, 0).

Exercício exemplo 1.2.5 Usando a definição, mostre que

DEMONSTRAÇÃO: AULA TEÓRICA. 1

Exercício exemplo 1.2.8 Usando a proposição anterior, considere as sucessões xn e yn a

Proposição 1.2.3 Sejam f : RN −→ R, g : R −→ R e a um ponto de acumulação de Df

DEMONSTRAÇÃO: AULA TEÓRICA.

Exercício exemplo 1.2.9 Usando a proposição anterior, mostre que

Definição 1.2.2 Sejam f : RN −→ R e a um ponto de acumulação de domínio Df . Diz-se

lim f (x) = f (a)

∀ ε > 0 ∃ δ > 0 : ( ∀ x ∈ Df e kx − ak < δ) =⇒ |f (x) − f (a)| < ε .

A Proposição 1.2.2 dá-nos a seguinte condição necessária e suficiente de continuidade.

Proposição 1.2.4 Sejam f uma função escalar e a um ponto de acumulação de Df . Tem-se:

f é contínua em a ⇐⇒ (∀ xn ∈ Df e xn 6= a : xn −→ a =⇒ f (xn ) −→ f (a)) .

DEMONSTRAÇÃO - EXERCÍCIO: Análoga à demonstração da Proposição 1.2.2 2 .

A função f diz-se contínua num conjunto A ⊆ Df , se for contínua em todos os pontos de

As propriedades de continuidade das funções reais de variável real estendem-se naturalmente

Exercício exemplo 1.2.10 Estude a função seguinte quanto à continuidade:

Proposição 1.2.5 Sejam f : RN −→ R uma função contínua em a ∈ Df , g : R −→ R uma

DEMONSTRAÇÃO - EXERCÍCIO: Análoga à demonstração da Proposição 1.2.3.

2. Considere a função seguinte:

e conclua que f (x, y) não tem limite no ponto (0, 0).