Introducao Elementar As Tecnicas Do Calculo

TERCEIRA EDIÇÃO
INTRODUÇÃO ELEMENTAR ÀS TÉCNICAS DO
CÁLCULO DIFERENCIAL E
INTEGRAL
f (x)
b
Z
A = f (x)dx
a
a b x
Carlos E. I. Carneiro, Carmen P. C. Prado e Silvio R. A. Salinas
INSTITUTO DE FÍSICA DA UNIVERSIDADE DE SÃO PAULO

i
DOI:10.11606/9788529200095
Carlos E. I. Carneiro, Carmen P. C. Prado e Silvio R. A. Salinas
INTRODUÇÃO ELEMENTAR ÀS TÉCNICAS DO

CÁLCULO DIFERENCIAL E INTEGRAL
TERCEIRA EDIÇÃO
São Paulo
Instituto de Fı́sica USP
2018
ii
Este trabalho está licenciado com uma Licença

Creative Commons - Atribuição 4.0 Internacional. Esta licença
permite que outros distribuam, remixem, adaptem e criem a partir deste trabalho, mesmo
para fins comerciais, desde que atribuam aos autores o devido crédito pela criação original.
Capa: Carlos E. I. Carneiro
FICHA CATALOGRÁFICA
Preparada pelo Serviço de Biblioteca e Informação
do Instituto de Fı́sica da Universidade de São Paulo
Carneiro, Carlos E. I..Prado, Carmen P. C..Salinas, Silvio R. A.
Introdução elementar às técnicas do calculo diferencial e integral.

São Paulo, Instituto de Fı́sica 2018. 3.ed. 91p.
ISBN:978-85-292-0009-5
DOI:10.11606/9788529200095
1. CÁLCULO DIFERENCIAL E INTEGRAL (TEXTOS

ELEMENTARES); 2. Prado, Carmen P. C.; 3. Salinas, Silvio R. A.;
USP/IF/SBI-053/2018
Sobre os autores: Carlos E. I. Carneiro; Carmen P. C. Prado e Silvio R. A. Salinas são

professores do Instituto de Fı́sica da Universidade de São Paulo.
Sumário
1 Introdução 1
2 Limites 3
2.1 Limite de uma função . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2 Definição mais precisa de limite . . . . . . . . . . . . . . . . . . . . . . . . 5
3 Derivadas 7
3.1 Definição de derivada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.2 Propriedades mais comuns das derivadas . . . . . . . . . . . . . . . . . . . 10
3.3 Interpretação geométrica da derivada . . . . . . . . . . . . . . . . . . . . . 14
4 Integrais 18
4.1 O conceito de integral . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
4.2 Propriedades das integrais definidas . . . . . . . . . . . . . . . . . . . . . . 21
4.3 Teorema Fundamental do Cálculo (TFC) . . . . . . . . . . . . . . . . . . . 23
4.3.1 Demonstração pouco rigorosa do TFC . . . . . . . . . . . . . . . . 23
4.4 Integrais indefinidas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
4.5 Cálculo de integrais definidas . . . . . . . . . . . . . . . . . . . . . . . . . 28
4.6 As funções logaritmo e exponencial . . . . . . . . . . . . . . . . . . . . . . 30
4.7 Algumas técnicas de integração . . . . . . . . . . . . . . . . . . . . . . . . 32
4.7.1 Integral de uma derivada . . . . . . . . . . . . . . . . . . . . . . . . 32
4.7.2 Integração por partes . . . . . . . . . . . . . . . . . . . . . . . . . . 33
4.7.3 Mudança de variável de integração . . . . . . . . . . . . . . . . . . 35
4.8 O que fazer quando nada funciona? . . . . . . . . . . . . . . . . . . . . . . 41
5 Vetores 43
5.1 Conceito de vetor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
5.2 Componentes e módulo de um vetor; versor . . . . . . . . . . . . . . . . . . . 43
5.3 Operações com vetores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
5.3.1 Soma ou subtração . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
5.3.2 Produto de vetores . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
5.4 Funções vetoriais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
5.5 Sistema de coordenadas polares . . . . . . . . . . . . . . . . . . . . . . . . 52
iii
SUMÁRIO iv
6 Expansões em séries de potências. 56

6.1 Definições, séries geométrica e de Taylor . . . . . . . . . . . . . . . . . . . 56
6.2 A exponencial complexa . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
7 Equações diferenciais simples 63

7.1 Equações diferenciais de primeira ordem . . . . . . . . . . . . . . . . . . . 63
7.2 Equações diferenciais de segunda ordem . . . . . . . . . . . . . . . . . . . . 71
7.2.1 O oscilador harmônico . . . . . . . . . . . . . . . . . . . . . . . . . 71
7.2.2 Solução de equações diferenciais com a exponencial complexa . . . 73
7.2.3 Oscilações amortecidas . . . . . . . . . . . . . . . . . . . . . . . . 75
7.2.4 Oscilações amortecidas e forçadas . . . . . . . . . . . . . . . . . . . 78
7.2.5 Operadores diferenciais com coeficientes constantes . . . . . . . . . 80
A Soluções dos exercı́cios 83

A.1 Limites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
A.2 Derivadas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
A.3 Integrais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
A.4 Vetores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
A.5 Expansões em séries de potências . . . . . . . . . . . . . . . . . . . . . . . 87
A.6 Equações Diferenciais . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
Prefácio da segunda edição
Na sua primeira edição este livro pretendia apresentar as idéias básicas do cálculo diferen-
cial e integral aos estudantes do primeiro ano dos cursos de fı́sica e engenharia. Com isto
em mente, tentamos nos restringir ao essencial para que o aluno adquirisse uma visão de
conjunto, embora superficial, do cálculo, além de uma habilidade operacional mı́nima no
manejo de derivadas e integrais. Definir o que é essencial não é uma tarefa fácil e alguns
tópicos importantes não foram incluı́dos na primeira edição.
Nesta segunda edição tentamos preencher algumas lacunas. Incluimos uma seção
sobre as funções logaritmo e exponencial, que aparecem de forma essencial e intensiva na
fı́sica, matemática e engenharia. Incluimos também algumas técnicas para a resolução de
integrais, entre elas a integração por partes e a mudança da variável de integração. Sem
o conhecimento destas duas técnicas o número de integrais que podem ser calculadas é
muito limitado.
A existência de belos textos introdutórios de fı́sica, como as aulas de Richard Feynman,
que utilizam variáveis complexas para tratar oscilações e circuitos com corrente alternada,
levou-nos a adicionar uma seção sobre a exponencial complexa e outra sobre sua utilização
na resolução de equações diferenciais. Incluı́mos também um capı́tulo sobre a resolução
de equações diferenciais de primeira ordem com a propriedade de separação de variáveis.
Evidentemente estes tópicos não são tão elementares como o resto do material coberto
neste livro. Porém, fizemos um grande esforço para torná-los bem acessı́veis. Quem quiser
apenas uma visão geral do cálculo, mais próxima em espı́rito da primeira edição, pode ler
o livro até o item 5.1 sobre séries de Taylor e omitir o restante.
Finalmente, nesta edição foram acrescentadas listas de exercı́cios no final de cada
capı́tulo. As soluções de todos os exercı́cios, algumas em detalhe, são apresentadas nos
apêndices. As listas não são extensivas e visam apenas verificar se os conceitos básicos
foram bem assimilados.
São Paulo, agosto de 2011
Prefácio da terceira edição

Nesta terceira edição fizemos algumas mudanças na seção 3.7.3 (Mudança de variável de
integração). Também acrescentamos um capı́tulo sobre equações diferenciais ordinárias
lineares de segunda ordem onde analisamos sistemas que executam oscilações amortecidas
e forçadas. Em particular, estudamos as equações diferenciais do sistema massa mola
e do circuito RLC enfatizando a notável analogia entre elas. É na solução deste tipo
de equação que a técnica da exponencial complexa mostra toda a sua versatilidade. O
capı́tulo 6 sobre equações diferenciais é mais avançado do que os anteriores. Porém, como
as equações diferenciais aparecem em todos os domı́nios da fı́sica é importante fornecer
uma introdução ao assunto. Com este acréscimo, esperamos cobrir uma parte razoável
dos requisitos de cálculo para os cursos de fı́sica básica.
São Paulo, maio de 2018
v
Capı́tulo 1
Introdução
Essas notas introdutórias sobre cálculo diferencial e integral, elaboradas há vários anos,
eram distribuı́das em forma manuscrita para complementar a primeira disciplina de fı́sica,
sobre fenômenos mecânicos, oferecida aos alunos dos cursos de ciências básicas e enge-
nharia da USP em São Paulo. Elas devem ser lidas e estudadas em paralelo com um
bom texto introdutório de fı́sica, como o primeiro volume da coleção de H. Moysés Nus-
senzveig [1]. Nessas notas são introduzidos os conceitos de limite, derivada e integral, de
maneira intuitiva e sem qualquer preocupação com o rigor matemático. Esses conceitos
– e os teoremas matemáticos pertinentes – serão apresentados com mais rigor e detalhe
nas disciplinas de cálculo ou análise matemática.
Não temos a pretensão de substituir as aulas de cálculo. No entanto, num curso univer-
sitário é importante desde o inı́cio trabalhar com as leis do movimento utilizando recursos
do cálculo diferencial e integral, que foram inventados pelo próprio Newton para formular
a mecânica há cerca de trezentos anos. Nos textos de fı́sica os conceitos de derivada e de
velocidade instantânea são inseparáveis. Como encontrar a equação horária, mesmo dos
movimentos mais simples, ou calcular o trabalho de uma força sem introduzir o conceito
de primitiva, ou integral, de uma função? Nessas notas também vamos apresentar um
apanhado de noções básicas sobre vetores (velocidade, aceleração e força são grandezas ve-
toriais), incluindo as definições de produto escalar e produto vetorial, bem como a análise
do movimento circular em coordenadas polares. Na seção final discutimos a expansão de
Taylor, ferramenta importante nas aplicações do cálculo diferencial.
Muitas vezes certos conceitos ou técnicas básicas são ensinados em diversas disciplinas,
em épocas distintas, em nı́veis diferentes. Na nossa opinião isso apenas reforça o processo
de aprendizagem. É possı́vel que alguns alunos já tenham estudado as operações básicas
do cálculo diferencial ou já tenham sido expostos às noções de limite e derivada. Nesse
caso o nosso texto vai ser um mero reforço operacional. Pretendemos ensinar apenas
as idéias mais intuitivas e algumas técnicas simples. Todo o cálculo necessário para
o acompanhamento dos problemas da fı́sica deve ser visto no decorrer das disciplinas
introdutórias dos cursos de ciência ou tecnologia. Temos a esperança de que as aulas de
mecânica possam contribuir para um primeiro contato com esses métodos matemáticos.
O historiador da ciência Alexandre Koyré, enfatizando a conexão entre a fı́sica moderna
e a matemática, escreve que “um experimento é uma pergunta que fazemos à natureza
e que deve ser formulada numa linguagem apropriada. A revolução galileana pode ser
resumida na descoberta dessa linguagem, na descoberta de que as matemáticas são a
gramática da ciência fı́sica. Foi esta descoberta da estrutura racional da natureza que
1
CAPÍTULO 1. INTRODUÇÃO 2
formou a base a priori da ciência experimental moderna e tornou possı́vel a sua cons-
tituição.” A ciência moderna representa antes de tudo um profundo rompimento com
as idéias do cosmos aristotélico. O universo moderno é aberto, indefinido e até infinito,
unificado e governado pelas mesmas leis naturais. Desaparecem da perspectiva cientı́fica
todas as considerações baseadas no valor, na perfeição, na harmonia, na significação ou
no desı́gnio. É nesse contexto que a matemática se faz realidade e que as leis da fı́sica
clássica encontram valor e aplicação.
O estabelecimento de modelos matemáticos - no estilo das antigas leis de Kepler ou de
Galileu - ganhou uma dimensão extraordinária na segunda metade do século XX. Antes se
observava, classificava e especulava. Agora se dá ênfase ao valor de teorias ou de modelos,
em geral formulados com o auxı́lio da matemática, em vários ramos das ciências e das
tecnologias, incluindo até as ciências humanas como a sociologia e a psicologia. Esperamos
que o nosso estudo das equações do movimento constitua um exemplo de alcance mais
geral, e que a conexão entre a matemática, as ciências e a tecnologia esteja presente em
todo o ensino moderno das engenharias.
São Paulo, janeiro de 2007
Capı́tulo 2
Limites
2.1 Limite de uma função

O conceito de limite de uma função vai ser tratado com rigor nas disciplinas de análise
matemática. Vamos apresentar aqui um resumo com algumas idéias que serão úteis no
cálculo de derivadas e de integrais em problemas de interesse fı́sico.
O limite de uma função f (x) quando seu argumento x tende a x0 é o valor L para o
qual a função se aproxima quando x se aproxima de x0 (Note que a função não precisa
estar definida em x0 .).
Se f (x) está definida em x0 e seu gráfico não apresenta descontinuidades nem oscilações
muito fortes (como ocorre com a função sen (1/x) próximo de x = 0) é natural escrever
lim f (x) = L ≡ f (x0 ),

x→x0
ou seja, o limite é igual ao valor da função em x0 .
Exemplos
(a) Com f (x) = 3x3 + 2x + 4, temos
lim f (x) = f (2) = 3 × 23 + 2 × 2 + 4 = 32,

x→2
lim f (x) = f (0) = 4.
x→0
(b) Com f (x) = 2 sen x + 3 cos 3x, temos
lim f (x) = 2 sen 0 + 3 cos 0 = 3,

x→0
lim f (x) = 2 sen π + 3 cos 3π = −3.
x→π
Em alguns casos, no entanto, a função não é bem definida e pode haver problemas
sérios. Por exemplo, para que valor tende a função
x3 − 27
f (x) = , x 6= 3,
x−3
3
CAPÍTULO 2. LIMITES 4
quando x → 3? Esse tipo de limite pode parecer um tanto artificial, mas é exatamente
o tipo de problema que temos que resolver para calcular a velocidade ou a aceleração
instantâneas, que são expressas por uma fração em que tanto o numerador quanto o
denominador vão para zero. Se calcularmos f (3) obteremos
33 − 27 0
f (3) = =
3−3 0
e temos problemas sérios pela frente. No entanto, com uma calculadora de bolso é possı́vel
traçar um gráfico de f (x) contra x nas vizinhanças de x = 3. A partir dos valores da
tabela abaixo, traçamos o gráfico da figura 2.1.
f (x)
x f (x)
2, 997 26, 973 27,027
2, 998 26, 982
2, 999 26, 991 27,000
3, 000 −
3, 001 27, 009 26,991
3, 002 27, 018
3, 003 27, 027 x
2,999 3,000 3,001
Figura 2.1
Observando os valores numéricos dessa tabela, apesar da função não estar definida
para x = 3, dá para desconfiar que à medida que nos aproximamos de x = 3,
x3 − 27
f (x) = −→ 27,
x−3
não sendo necessário que a função esteja definida em x = 3. Observe que, para x 6= 3,
sempre podemos escrever
x3 − 27 (x − 3)(x2 + 3x + 9)
f (x) = = = x2 + 3x + 9.
x−3 (x − 3)
Como x2 + 3x + 9 = 27 para x = 3, há uma motivação muito forte para escrever L = 27.
Para lidar com situações como essa criou-se uma definição de limite onde o que acontece
exatamente no ponto em que se deseja calcular o limite não é importante. Importa apenas
o que ocorre nas vizinhanças desse ponto. Isso permite o cancelamento dos fatores comuns
no numerador e no denominador como acabamos de fazer. Assim,
x3 − 27
lim f (x) = lim = lim (x2 + 3x + 9) = 27.
x→3 x→3 x − 3 x→3
As funções da mecânica clássica são em geral muito bem comportadas. Porém, há
situações em que as funções são tão mal comportadas que o limite não existe mesmo. Por
exemplo, vamos considerar a função

+1 para x > 1,
f (x) =
−1 para x < 1,
que pode ser representada pelo gráfico da figura 2.2.
f (x)
+1
0 1 x
−1
Figura 2.2
Quando x → 1 o limite é claramente indefinido. Se fizermos x → 1 por valores maiores

do que 1 obtemos +1; se fizermos x → 1 por valores menores do que 1 obtemos −1. Note
que f (x) pode ser escrita na forma
x−1
f (x) = = sgn(x − 1),
|x − 1|
em que sgn é a “função sinal”, sgn (x) = 1 se x > 0 e sgn (x) = −1 se x < 0. Ainda é
possı́vel trabalhar com casos desse tipo, pois não há muitos dificuldades na presença de
uma descontinuidade isolada.
2.2 Definição mais precisa de limite

Depois desses exemplos intuitivos e meio óbvios, vale a pena apresentar uma definição mais
formal de limite, com todos os épsilons e deltas. Considere uma função f (x) definida no
domı́nio x1 < x < x0 e x0 < x < x2 (não precisando, portanto, estar definida no ponto
x0 ).
A função f se aproxima do limite L próximo de x0 (limx→x0 f (x) = L) se

para todo ǫ > 0 existe um δ > 0 tal que, para todo x, se 0 < |x − x0 | < δ
então |f (x) − L| < ǫ.
Na prática isso significa que, quando o limite existe, se for dado um “limite de to-
lerância” ε > 0 arbitrariamente pequeno (tão pequeno quanto se queira), podemos sempre
encontrar um outro número δ > 0 tal que, para qualquer valor de x entre x0 − δ e x0 + δ ,
o valor da função estará dentro do “limite de tolerância”, ou seja, f (x) estará entre L − ε
e L + ε. Isto funciona quer f (x) seja definida ou não no ponto x0 .
Por exemplo, vamos considerar o limite da função f (x) = 2x2 para x → x0 = 3. É
claro que
lim (2x2 ) = 18 ≡ L.
x→3
Dando o “limite de tolerância” ε = 0, 1 temos
2x21 = 18 − 0, 1 −→ x1 = 2, 991655... = 3 − 0, 008344...,

2x22 = 18 + 0, 1 −→ x1 = 3, 008321... = 3 + 0, 008321....
É claro que, para δ = 0, 008, temos
|f (x) − 18| < ε = 0, 1 quando 0 < |x − 3| < δ = 0, 008.

Adotando um “limite de tolerância” menor obviamente vamos ter que encontrar um valor
de δ menor também (para ε = 0, 01 é fácil perceber que tudo funciona com δ = 0, 0008, e
assim por diante).
Finalmente, a partir da idéia de limite podemos definir continuidade de uma função.
A função f é contı́nua em x0 se
lim f (x) = f (x0 ).
x→x0
Exercı́cios
(1) Mostre que
x2 − 5x + 6
(a) lim = −1, (b) lim x (x − 2) = 0,
x→2 x−2 x→2
y 2 − 25
(c) lim x3 − 4x + 1 = 1, (d) lim = 10.
x→0 y→5 y − 5
(2) Calcule o valor dos seguintes limites:
x−5 x2 − x − 2 x+5
(a) lim , (b) lim , (c) lim ,
x→5 x2 − 25 x→1 x2 − 1 x→5 x2 − 25
5x3 + 8x2 1
(d) lim , (e) lim .
x→0 3x4 − 16x2 x→0 x2
Note que há limites que não existem ou que vão para o infinito 1 .
1
Dizemos que limx→a f (x) = ∞ se para qualquer N existe um δ > 0 tal que para todo x, se 0 <
|x − a| < δ então f (x) > N . Analogamente, limx→a f (x) = −∞ se para qualquer N existe um δ > 0 tal
que para todo x, se 0 < |x − a| < δ então f (x) < N .
Capı́tulo 3
Derivadas
3.1 Definição de derivada

O cálculo diferencial foi inventado por Leibnitz e Newton, que sempre disputaram a
primazia das suas propostas! Embora usasse uma notação um tanto complicada, Newton
desenvolveu o conceito de derivada e percebeu a sua utilidade na formulação matemática
da mecânica. Tecnicamente a derivada de uma função não passa de um caso especial
de limite. A velocidade instantânea (que é a derivada da posição em relação ao tempo)
corresponde ao limite da velocidade média para um intervalo de tempo muito pequeno
(que tende a se anular).
Para calcular a derivada de uma função f (x) num certo ponto x0 , nós inicialmente
damos um acréscimo ∆x em x0 e calculamos a diferença
∆f = f (x0 + ∆x) − f (x0 )
e a razão
∆f f (x0 + ∆x) − f (x0 )
= .
∆x ∆x
A derivada no ponto x0 , designada por (df /dx)x=x0 , é dada pelo limite

df ∆f f (x0 + ∆x) − f (x0 )
= lim = lim .
dx x=x0
∆x→0 ∆x ∆x→0 ∆x
Exemplos
(a) Calcular a derivada de f (x) = 3x2 num certo ponto x0 .

Temos
f (x0 ) = 3x20
e
f (x0 + ∆x) = 3(x0 + ∆x)2 ,
de onde vem que
∆f ≡ f (x0 + ∆x) − f (x0 ) = 3(x0 + ∆x)2 − 3x20 = 6x0 ∆x + 3∆x2
7
CAPÍTULO 3. DERIVADAS 8
e
∆f 6x0 ∆x + 3∆x2
= = 6x0 + 3∆x.
∆x ∆x
Portanto,
df ∆f
= lim = 6x0 .
dx x=x0 ∆x→0 ∆x
Como x0 é um ponto genérico, também podemos escrever de forma genérica, para
qualquer ponto x,
df
f (x) = 3x2 =⇒ = 6x.
dx
(b) Calcular a derivada de f (x) = Axn , onde A é uma constante e n é um número

inteiro.
No ponto x0 temos
f (x0 ) = Axn0
e
f (x0 + ∆x) = A(x0 + ∆x)n .
Usando a expansão do binômio de Newton, podemos escrever

n n n n n−1 n n−2 2 n
(x0 + ∆x) = x0 + x0 ∆x + x0 (∆x) + ... + (∆x)n ,
0 1 2 n
com a notação combinatorial

n n!
= .
p (n − p)!p!
Então
n(n − 1) n−2 n(n − 1)(n − 2) n−3
∆f = A n x0n−1 ∆x + A x0 (∆x)2 + A x0 (∆x)3 + ...,
2 6
de onde obtemos
∆f n(n − 1) n−2
= A n x0n−1 + A x0 (∆x) + ....
∆x 2
No limite ∆x → 0 , somente o primeiro termo sobrevive, pois todos os outros têm
pelo menos um fator ∆x. Então temos a derivada

df ∆f
= lim = A n x0n−1 .
Esses exemplos nos conduzem a uma regra de derivação importantı́ssima, que é

preciso saber de cor:
df
f (x) = Axn =⇒ = Anxn−1 .
dx
Na realidade essa regra é válida para qualquer valor de n (negativo, fracionário, etc).
Note que, para n = 0, f (x) = A, em que A é uma constante, e df /dx = 0. Vamos
ver que essa regra vai permitir o cálculo da derivada de qualquer função polinomial.
(c) Calcular a derivada de f (x) = cos x no ponto x0 .

No ponto x0 temos
f (x0 ) = cos x0
e
f (x0 + ∆x) = cos(x0 + ∆x) = cos x0 cos ∆x − sen x0 sen ∆x,
de onde vem que
∆f cos x0 cos ∆x − sen x0 sen ∆x − cos x0
= .
∆x ∆x
Tomando o limite de forma ingênua, isto é, colocando diretamente ∆x = 0 nessa
expressão, tem-se uma indeterminação (zero sobre zero!). No entanto, para ∆x
muito pequeno temos
cos ∆x ≈ 1 e sen ∆x ≈ ∆x,
que conduz ao resultado
∆f
≈ − sen x0 .
∆x
Temos então a regra de derivação
df
f (x) = cos x =⇒ = − sen x.
dx
Também é fácil obter uma regra para a função seno,
df
f (x) = sen x =⇒ = cos x.
dx
Uma justificativa mais adequada para essas fórmulas será apresentada durante o
primeiro curso de cálculo (veja também a seção final, sobre séries de Taylor).
Ao invés de continuar com mais exemplos, vamos dar uma relação das derivadas mais
comuns no curso de fı́sica (veja a tabela 3.1). Certamente você vai aprender a justificativa
de todas essas fórmulas nas disciplinas de cálculo. Note que as funções exponencial,
f (x) = exp (x), e logaritmo natural, f (x) = ln x, são muito usadas em fı́sica.
Tabela 3.1: Algumas derivadas fundamentais
f (x) df /dx
constante zero
xn nxn−1
sen x cos x
cos x − sen x
x
e ex
ln x 1/x
cosh x senh x
senh x cosh x
O cosseno hiperbólico e o seno hiperbólico são dados por

1 x
cosh x = e + e−x
2
e
1 x
senh x = e − e−x .
2
É muito importante conhecer os gráficos de todas essas funções. Em particular, na figura
3.1 esboçamos os gráficos da função exponencial, ex ≡ exp x, e da função logaritmo, ln x.
Vamos estudá-las em mais detalhe no capı́tulo 3 (Integrais).
y y
1 ex
ln x
x x
1
Figura 3.1
Sugerimos agora que você trace gráficos de mais algumas funções: (i) f (x) = exp(−x);
(ii) f (x) = tan x = sen x/ cos x, que é a função tangente trigonométrica; e (iii) f (x) =
tanh x = senh x/ cosh x, que é a função tangente hiperbólica. Todas essas funções são
muito úteis em fı́sica.
3.2 Propriedades mais comuns das derivadas

Vamos relacionar algumas propriedades, facilmente demonstráveis, que simplificam enor-
memente o cálculo das derivadas.
(1) Dada a função f (x) = c g(x), onde c é uma constante e g(x) é outra função (bem
comportada), temos
df dg
=c .
dx dx
(2) Dada a funcão f (x) = af1 (x) + bf2 (x) , onde f1 (x) e f2 (x) são funções e a e b são
constantes, temos
df df1 df2
=a +b ,
dx dx dx
ou seja, a derivada é uma “operação linear”.
(3) Dado o produto de funções, f (x) = f1 (x)f2 (x), temos

df df1 df2
= f2 + f1 .
dx dx dx
Essa regra do produto é tão importante que vale a pena ser demonstrada. De fato,
temos
f (x0 ) = f1 (x0 )f2 (x0 )
e
f (x0 + ∆x) = f1 (x0 + ∆x)f2 (x0 + ∆x).
Portanto,
∆f = f1 (x0 + ∆x)f2 (x0 + ∆x) − f1 (x0 )f2 (x0 ).
Somando e subtraindo f1 (x0 )f2 (x0 + ∆x), temos
∆f = f1 (x0 + ∆x)f2 (x0 + ∆x) − f1 (x0 )f2 (x0 + ∆x) +

+ f1 (x0 )f2 (x0 + ∆x) − f1 (x0 )f2 (x0 ).
Então
∆f = [f1 (x0 + ∆x) − f1 (x0 )] f2 (x0 + ∆x) +
+ f1 (x0 ) [f2 (x0 + ∆x) − f2 (x0 )] .
Portanto,
∆f f1 (x0 + ∆x) − f1 (x0 )
= f2 (x0 + ∆x) +
∆x ∆x
f2 (x0 + ∆x) − f2 (x0 )
+ f1 (x0 ) .
∆x
Tomando o limite ∆x → 0, finalmente temos

df df1 df2
= f2 (x0 ) + f1 (x0 ) .
dx x=x0 dx x=x0 dx x=x0
Em termos gerais, podemos escrever

df d(f1 f2 ) df1 df2
= = f2 + f1 .
dx dx dx dx
(4) Dado o quociente de duas funções, f (x) = f1 (x)/f2 (x), temos

df 1 df1 df2
= 2 f2 − f1 .
dx f2 dx dx
Faça um esforço para demonstrar essa propriedade.
(5) Muitas vezes temos que calcular a derivada de uma “função de função”. Vamos
considerar a função f = f (y) onde y = y(x). A derivada é dada pela “regra da
cadeia”,
df df dy
= .
dx dy dx
Vale a pena demonstrar essa regra. Vamos então escrever
∆f = f [y(x0 + ∆x)] − f [y(x0 )] .
Somando e subtraindo y(x0 ) também temos
∆f = f [y(x0 + ∆x) − y(x0 ) + y(x0 )] − f [y(x0 )] .
Agora vamos usar a notação abreviada y(x0 ) = y0 e y(x0 +∆x)−y(x0 ) = ∆y. Então
∆f f [y0 + ∆y] − f [y0 ] f [y0 + ∆y] − f [y0 ] ∆y

= = .
∆x ∆x ∆y ∆x
Tomando o limite ∆x → 0, vem

df df dy
= ,
dx x=x0 dy y=y0 dx x=x0
ou seja,
df df dy
= ,
dx dy dx
como querı́amos demonstrar.
Na tabela 3.2, onde a, b e c são constantes, apresentamos um resumo destas proprie-

dades.
Tabela 3.2: Propriedades importantes das derivadas
d df
(1) {cf } = c
dx dx
d df1 df2
(2) {af1 + bf2 } = a +b
dx dx dx
d df1 df2
(3) {f1 f2 } = f2 + f1
dx dx dx

d f1 1 df1 df2
(4) = 2 f2 − f1
dx f2 f2 dx dx
d df dy
(5) {f (y (x))} =
dx dy dx
Exemplos
(a) f (x) = Ax4 + Bx2 + Cx + D.
O cálculo da derivada é imediato,

df
= 4Ax3 + 2Bx + C.
dx
(b) f (x) = (Ax4 + Bx2 + Cx + D)3 .
Basta fazer f (y) = y 3 , y = Ax4 + Bx2 + Cx + D e aplicar a “regra da cadeia”:

df df dy
= = 3y 2 4Ax3 + 2Bx + C =
dx dy dx
2
= 3 Ax4 + Bx2 + Cx + D 4Ax3 + 2Bx + C .
(c) f (t) = exp (at2 + b).
Basta fazer f (y) = exp (y) com y = at2 + b. Então

df df dy
= = ey (2at) = 2at exp at2 + b .
dt dy dt
−3
(d) f (t) = (t + 1)2 (t2 + 2t) .
Basta fazer f (t) = f1 (t)f2 (t), com

−3
f1 (t) = (t + 1)2 ; f2 (t) = t2 + 2t .
Então,
df1 df2 −4
= 2 (t + 1) ; = −3 t2 + 2t (2t + 2) ,
dt dt
de onde obtemos
df −3 h −4 i
2 2 2
= 2 (t + 1) t + 2t + (t + 1) −3 t + 2t (2t + 2) .
dt
3.3 Interpretação geométrica da derivada

A derivada de uma função y = y(x) num certo ponto x0 corresponde ao valor da tangente
da curva y contra x no ponto x0 . Isso pode ser facilmente visualizado através de um
argumento gráfico (veja a figura 3.2).
y
y1
y2 θ1
θ2
y3 θ3
y0
x0 x3 x2 x1 x
Figura 3.2
Considere três pontos obtidos através de acréscimos em x0 . Para o maior acréscimo,

x1 − x0 , temos
∆y y1 − y0
= = tan θ1 .
∆x x1 − x0
Para x2 , temos
∆y y2 − y0
= = tan θ2 .
∆x x2 − x0
Da mesma forma, para x3 escrevemos

∆y y3 − y0
= = tan θ3 .
∆x x3 − x0
A partir dessa construção, fica óbvio que, à medida que ∆x diminui, a razão ∆y/∆x vai
se aproximando de tan θ, onde θ é o ângulo formado entre a reta tangente à curva y(x),
passando pelo ponto de coordenadas (x0 , y0 ), e o eixo x. Então temos

dy ∆y
= lim = tan θ.
Essa interpretação da derivada tem inúmeras utilidades. Por exemplo, dada a equação
de uma trajetória unidimensional, x = x(t), a velocidade é dada por v = dx/dt. Gra-
ficamente isto significa que a velocidade é a tangente da curva num gráfico de x contra
t. Dada a velocidade em função do tempo, v = v(t), a aceleração é definida através da
derivada a = dv/dt, que pode ser interpretada como a tangente da curva no gráfico de v
contra o tempo t.
Exemplo: movimentos retilı́neos
Vamos utilizar o conceito de derivada para obter as equações dos movimentos retilı́neos
mais simples.
(a) No movimento retilı́neo uniforme (MRU) a velocidade é constante,
v = v0 = constante.
Temos então a aceleração,

dv
a= = 0.
dt
A equação horária do MRU é dada pela expressão
x = v0 t + constante.
Portanto, podemos verificar a expressão da velocidade,

dx d
= [v0 t + constante] = v0 .
dt dt
Fazendo t = 0 na equação horária, percebemos que a constante é a posição inicial.
Em resumo, o MRU com x = x0 no instante inicial t = 0 é caracterizado pelas
equações
x = v 0 t + x0 ; v = v 0 e a = 0 .
(b) O movimento retilı́neo uniformemente variado (MRUV) é definido por uma

aceleração constante,
a = a0 = constante.
No ensino médio os alunos devem ter aprendido que a equação horária do MRUV é
dada pela expressão
1
x = x 0 + v 0 t + a0 t 2 ,
2
em que a constante v0 é a velocidade inicial (no instante de tempo t = 0) e a
aceleração a0 é uma contante. Então é fácil obter a velocidade,

dx d 1 2
v= = x 0 + v 0 t + a0 t = v 0 + a0 t
dt dt 2
e verificar que a aceleração é dada pela constante a0 ,

dv d
a= = [v0 + a0 t] = a0 .
dt dt
(c) Vamos agora considerar agora um movimento descrito pela equação horária
1
x = c t 3 + v 0 t + x0 ,
6
a velocidade instantânea de um corpo que executa esse movimento é dada por
dx 1
v= = c t 2 + v0
dt 2
e a aceleração por
dv
a= = c t.
dt
Essa é, portanto, a equação horária de um corpo sujeito a uma aceleração
que varia linearmente com o tempo. Note que as constantes x0 e v0 são a
posição e a velocidade no instante inicial.
Exercı́cio
Calcule a derivada em relação a x ou a t das funções abaixo, onde a, b, c, ω e φ são

constantes.
y1 (x) = (x2 + 5)8 ;
y2 (t) = cos(ωt + φ);
y3 (t) = [cos(ωt)]2 ;
y4 (x) = sen (ax2 + bx);
y5 (t) = exp(ωt);
y6 (x) = exp (ax2 + bx) ;
y7 (x) = ln (ax2 + bx + c) ;
√
y8 (x) = ax2 + bx + c;
1
y9 (x) = √ .
ax2 + c
Capı́tulo 4
Integrais
4.1 O conceito de integral

Dada a equação horária x = x(t), já vimos que é possı́vel obter a velocidade instantânea
v(t) tomando a derivada de x em relação a t, isto é,
d
v(t) = x(t).
dt
Frequentemente temos que resolver o problema inverso: dada a velocidade v = v(t),
precisamos calcular o espaço percorrido entre um instante inicial ti e um instante final tf ,
isto é x(tf ) − x(ti ) ≡ xf − xi . Esse problema tem uma solução gráfica muito simples, que
conduz ao conceito de integral.
Vamos considerar o gráfico de v contra t indicado na figura 4.1.
v
∆t2 ∆t3
v̄3 11
00
0000
1111
0000
1111
00
11
0000
1111
0000
1111
0000
1111
∆t1 0000
1111
0000
1111
0000
1111
∆t4
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
0000
1111
ti = t0 t1 t2 t3 tf = t4 t
Figura 4.1
Se o gráfico de v contra t fornecesse a velocidade média v ti →tf a solução do problema

seria trivial; nesse caso o espaço percorrido seria dado por
xf − xi = v ti →tf (tf − ti ).
Mas a velocidade média não é conhecida. No entanto, se o intervalo tf − ti fosse bem
pequeno, a velocidade média v ti →tf seria muito aproximadamente igual a qualquer valor
18
CAPÍTULO 4. INTEGRAIS 19
da velocidade v(t) nesse intervalo. Usualmente o intervalo entre ti e tf não é pequeno, mas
sempre pode ser dividido num certo número de subintervalos (na figura 4.1 escolhemos
apenas quatro intervalos menores, de comprimentos ∆t1 = t1 − t0 , ∆t2 = t2 − t1 , ∆t3 =
t3 −t2 e ∆t4 = t4 −t3 ). Num subintervalo genérico (entre t2 e t3 , por exemplo), a velocidade
média é definida como
x(t3 ) − x(t2 ) x3 − x2
v t2 →t3 = = ,
t3 − t2 ∆t3
de onde vem que
x3 − x2 = v t2 →t3 ∆t3 .
Vamos simplificar um pouco a notação, escrevendo
v t2 →t3 = v 3 .
Então
x3 − x2 = v 3 ∆t3 .
A partir dessas considerações, é fácil perceber que a distância
x(tf ) − x(ti ) ≡ xf − xi ≡ x4 − x0
será dada por
xf − xi = (x4 − x3 ) + (x3 − x2 ) + (x2 − x1 ) + (x1 − x0 )

= v 4 ∆t4 + v 3 ∆t3 + v 2 ∆t2 + v 1 ∆t1 ,
onde xi = x0 e xf = x4 . Essa expressão pode ser escrita numa forma bem mais compacta,
4
X
xf − xi = v j ∆tj .
j=1
Ao invés de considerar apenas 4 subdivisões, poderı́amos ter subdividido o intervalo

tf − ti em N subintervalos bem menores. Nesse caso terı́amos
N
X
xf − xi = v j ∆tj .
j=1
Note que esta expressão é exata. O problema é que não conhecemos as velocidades médias
v j . Porém, no limite de N muito grande e quando o maior ∆tj vai a zero 1 , v j tende à
velocidade instantânea vj ≡ v(tj ), e a soma das inúmeras parcelas, que se denomina
integral definida, costuma ser escrita na forma
N
X N
X Z tf
xf − xi = lim v j ∆tj = lim vj ∆tj ≡ v(t)dt,
∆t→0 ∆t→0 ti
j=1 j=1
em que ∆t → 0 significa que todos os intervalos ∆tj vão a zero. Observe que a soma
P
virou um S estilizado. Ao invés de um ı́ndice j que assume valores discretos, há
uma variável de integração contı́nua em t (a velocidade instantânea vj é substituı́da pela
1
Consequentemente, todos os ∆tj vão a zero.
velocidade instantânea v(t), e ∆tj passa a ser um “intervalo infinitesimal” dt). À medida
PN
que N aumenta, é fácil perceber que a soma j=1 vj ∆tj corresponde cada vez mais
fielmente à área sob a curva do gráfico de v contra t. Nesse limite a soma, ou melhor,
a integral definida, corresponde exatamente à área sob a curva da função v = v(t) entre
t = ti e t = tf . Dessa forma, a integral é o caso particular de um limite — é um tipo de
limite em que as parcelas de uma soma tendem a zero, mas o número de parcelas tende
a infinito.
Os matemáticos são mais cuidadosos. Para funções contı́nuas, ao invés de usarem a
velocidade média na somatória, eles definem duas somas: a primeira soma usando o menor
valor da velocidade em cada subintervalo ∆tj ; a segunda soma usando o maior valor da
velocidade em cada subintervalo ∆tj . No limite em que N → ∞ e ∆tj → 0 para todos os
subintervalos, se essas duas somas convergirem para um mesmo valor, fica então definida
a integral de Riemann desta função.
Exemplos
(a) Dada a velocidade v(t) = v0 (constante), qual o espaço percorrido entre ti e tf ?

Graficamente, v(t) é dada pela figura 4.2.
v
v0 11111111111
00000000000
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
00000000000
11111111111
ti tf t
Figura 4.2
Então Z tf
xf − xi = v(t)dt = área hachurada = v0 (tf − ti ).
ti
Tomando xi = 0 para ti = 0 e um ponto genérico xf = x para tf = t, temos
x − x0 = v0 t,
ou seja,
x = x0 + v0 t,
que é a conhecidı́ssima equação do MRU.
(b) Dada a velocidade v(t) = v0 +at, onde v0 e a são constantes, qual o espaço percorrido
entre os instantes ti e tf ?
Vamos observar o gráfico da figura 4.3, em que
d
tan θ = v(t) = a.
dt
v 1111111111
0000000000
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
θ 0000000000
1111111111
0000000000
1111111111
v0 0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
ti tf t
Figura 4.3
Então Z tf
xf − xi = v(t)dt = área hachurada.
ti
Usando a fórmula da área de um trapézio, temos
v(tf ) + v(ti )
xf − xi = área hachurada = (tf − ti )
2
1 1
= [v0 + atf + v0 + ati ](tf − ti ) = v0 (tf − ti ) + a t2f − t2i .
2 2
Tomando de novo xi = 0 para ti = 0 e um ponto genérico xf = x para tf = t,
recuperamos a famosa equação horária do MRUV
1
x = x0 + v0 t + at2 .
2
Até agora vimos dois exemplos muito simples, sem nenhuma dificuldade para calcular
a “área sob a curva” (e encontrar o valor da integral definida). As situações práticas, no
entanto, podem ser bem mais complicados. Há poucas figuras geométricas cujas áreas
podem ser calculadas tão facilmente. Na grande maioria das vezes temos que utilizar
algumas propriedades gerais e um arsenal de truques para calcular diretamente as integrais
definidas (e obter, portanto, o valor das “áreas sob as curvas”). Com este objetivo,
vamos apresentar algumas propriedades muito simples das integrais e enunciar o “teorema
fundamental do cálculo”.
4.2 Propriedades das integrais definidas

É interessante apontar as seguintes propriedades das integrais definidas:
(1) Quebra dos limites de integração:

Z tf Z tm Z tf
f (t) dt = f (t) dt + f (t) dt.
ti ti tm
Para ti ≤ tm ≤ tf , essa propriedade é meio óbvia. Para se convencer disso, basta

observar a figura 4.4; a área total é a soma das áreas.
111111
000000
000001
111110
00000
11111
00000
111110
1
0
1
0
1
000001
11111
000000000
111111111
0
1 0
1 00000
11111
000000000
111111111
0 0
1
00000
111110
1
0
1
000000000
111111111
0
1 0
1
00000
111110
1
000000000
111111111
0
1 0
1
00000
11111
000000000
111111111 0
1
0
1 0
1
00000
11111
000000000
111111111 0
1
0
1 0
1
00000
11111
000000000
111111111 0
1
0
1 0
1
00000
11111
000000000
111111111
0
1 0
1
0
1
00000
11111
000000000
111111111
0
1 0
1 0
1
00000
11111
000000000
111111111
0
1 0
1 0
1
00000
11111
000000000
111111111
0
1 0
1 0
1
1 00000
11111
000000000
111111111
0 0
1
00000
111110
1
0
1
000000000
111111111
0
1 0
1
000001
111110
1
000000000
111111111
0
1 0
1
00000
11111
000000000
111111111 0
0
1 0
1
00000
11111
000000000
111111111 0
1
0
1
ti tm tf t
Figura 4.4
(2) Inversão dos limites:

Z b Z a
f (t) dt = − f (t) dt.
a b
Trocando os limites de integração a integral muda de sinal.

Como a integral provém de uma soma, o lado esquerdo representa o limite da soma,
desde t = a até t = b , com b > a , sendo todos os subintervalos ∆tj positivos
(∆tj > 0). No entanto, indo de b para a (com b > a ), todos os subintervalos ∆tj
serão negativos. Não há dúvidas, portanto, que a troca dos limites de integração
acarreta apenas a multiplicação por −1.
É fácil verificar que essa propriedade acaba garantindo a validade da propriedade (1)
para qualquer tm (isto é, mesmo para tm > tf ). Portanto, no cálculo das integrais
definidas é preciso levar em conta o sinal algébrico das “áreas sob as curvas”.
(3) Mesmos limites de integração:

Z a
f (t)dt = 0.
a
Essa propriedade é obvia, pois não há área sob “um único ponto da curva”.
(4) Multiplicação por uma constante:

Z t2 Z t2
Af (t)dt = A f (t)dt,
t1 t1
onde A é uma constante (ou uma função independente de t). Essa propriedade
também é obvia, bastando considerar uma soma em que todas as parcelas estejam
multiplicadas pela constante A.
(5) Linearidade:
Z t2 Z t2 Z t2
[Af (t) + Bg(t)] dt = A f (t)dt + B g(t)dt,
t1 t1 t1
onde A e B são constantes. Essa propriedade também é óbvia, pois a soma é

associativa, isto é, as parcelas sempre podem ser agregadas. Essa propriedade indica
que a integração é uma “operação linear”.
4.3 Teorema Fundamental do Cálculo (TFC)

Como o espaço percorrido é dado pela integral da velocidade (que, por sua vez, é a
derivada do espaço), a integração deve corresponder a uma operação inversa da derivação.
O “teorema fundamental do cálculo” torna esta idéia mais precisa.
Teorema: Se a função F (x) for dada por

Z x
F (x) = f (t)dt,
a
onde a é uma constante arbitrária, então

dF (x)
= f (x).
dx
A função F (x) se chama primitiva de f (x). A sua derivada coincide com o integrando
f (t) no ponto t = x. Estamos tomando bastante cuidado com a notação—como o intervalo
de integração vai de a até x, estamos usando o sı́mbolo t como variável de integração (no
extremo inferior, t = a; no extremo superior, t = x ). é claro que poderı́amos ter escolhido
qualquer outra letra (y, z, w, etc) como variável de integração.
4.3.1 Demonstração pouco rigorosa do TFC

Dada a expressão Z x
F (x) = f (t)dt,
a
temos Z Z Z
x+∆x x x+∆x
F (x + ∆x) = f (t)dt = f (t)dt + f (t)dt.
a a x
Então Z x+∆x
F (x + ∆x) − F (x) = f (t)dt.
x
Esta última integral é a área sob a curva do gráfico de f (t) contra t entre t = x e
t = x + ∆x. Para ∆x muito pequeno, temos
Z x+∆x
f (t)dt = área ≈ f (x)∆x.
x
Então
dF (x) F (x + ∆x) − F (x)
= lim = f (x),
dx ∆x→0 ∆x
como querı́amos demonstrar. ❚
Vamos verificar como este teorema funciona em dois casos conhecidos.
Exemplos
(a) Seja f (t) = A, com A constante.
f (t)
A 1111111111
0000000000
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
a x t
Figura 4.5
Nesse caso basta calcular a área sob o gráfico da função constante f (t) = A, entre
t = a e t = x, Z x Z x
F (x) = f (t) dt = A dt = A(x − a).
a a
Então é fácil verificar que
dF (x)
= A = f (x).
dx
(b) Seja f (t) = A + Bt, com A e B constantes.
f (t)
A + Bx 1111111111
0000000000
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
A + Ba 0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
0000000000
1111111111
a x t
Figura 4.6
Nesse caso, considerando a figura 4.6, temos

Z x Z x
F (x) = f (t) dt = (A + Bt)dt = área do trapézio =
a a
(A + Bx) + (A + Ba) 1
= (x − a) = A(x − a) + B x2 − a2 .
2 2
Então é claro que
dF
= A + Bx = f (x).
dx
Agora vai ser fácil descobrir o que acontece num caso mais complicado, em que não
seja trivial obter a área. Vamos tomar, por exemplo,
f (t) = A + Bt + Ct2 ,
com A, B e C constantes. Temos então que calcular

Z x Z x
F (x) = f (t)dt = (A + Bt + Ct2 )dt.
a a
Obviamente fica complicado apelar para uma fórmula que dê a área sob a curva de
f (t) entre t = a e t = x. Mas, a partir do “teorema fundamental do cálculo” temos
dF (x)
f (x) = = A + Bx + Cx2 .
dx
Então é possı́vel usar as “regras de derivação ao contrário” para “garimpar” a função
F (x). De fato, é simples verificar que
1 1
F (x) = Ax + Bx2 + Cx3 + k,
2 3
onde k é uma constante arbitrária (pois a derivada de uma constante é sempre nula).
Para encontrar a constante k é muito fácil. Basta notar que
Z a
F (x = a) = f (t)dt = 0.
a
Então
1 1
Aa + Ba2 + Ca3 + k = 0,
2 3
de onde finalmente obtemos
1 1 1 1
F (x) = Ax + Bx2 + Cx3 − Aa − Ba2 − Ca3 .
2 3 2 3
4.4 Integrais indefinidas

Para cada valor da constante k na expressão
1 1
G(x) = Ax + Bx2 + Cx3 + k
2 3
temos uma função G(x) diferente. Cada uma dessas funções é chamada primitiva da
função
f (x) = A + Bx + Cx2 ,
pois dG/dx = f (x) para qualquer valor de k.
Essas primitivas formam uma famı́lia de funções que são normalmente simbolizadas
como Z
G(x) = f (x) dx ,
sem a preocupação de especificar os limites de integração. Isso é o que se chama integral

indefinida. Como dG/dx = f (x), é claro que
Z x
F (x) = f (x) dx = G(x) − G(a),
a
com dF (x)/dx = dG(x)/dx = f (x).
O “teorema fundamental do cálculo” pode então ser reescrito na forma

Z x
G(x) − G(a) = f (t) dt ,
a
em que G(x) é uma primitiva genérica de f (x). Note que a constante aditiva k, distin-
guindo as diferentes primitivas, desaparece quando se faz a diferença G(x) − G(a).
Exemplos
(a) Dada a função f (x) = sen x, calcular a sua primitiva G(x). Temos
Z
G(x) = sen x dx.
Portanto, olhando a tabela de derivação “ao contrário”, obtemos

G(x) = − cos x + k.
É claro que
G(x) − G(a) = − cos x + cos a.
(b) Dada a função f (x) = cos x, calcular G(x). É simples perceber que
Z
G(x) = cos x dx = sen x + k.
(c) Dada a função f (x) = x4 − 10x2 , calcular G(x). Temos

Z Z
1 10
G(x) = f (x)dx = x4 − 10x2 dx = x5 − x3 + k.
5 3
(d) No MRUV, é dada a aceleração a(t) = a (constante). Obter v(t) e x(t), com as
condições iniciais v(t0 ) = v0 e x(t0 ) = x0 .
A velocidade será dada por
Z
v(t) = a dt = at + k1 .
Como v(t0 ) = v0 , temos

at0 + k1 = v0 ,
de onde obtemos a constante k1 . Podemos então escrever a expressão da velocidade
em termos mais usuais,
v(t) = v0 + a(t − t0 ).
A equação horária é obtida a partir de uma integração da velocidade,
Z
1
x(t) = v(t) dt = v0 t + at2 − at0 t + k2 .
2
A constante k2 é definida pela condição inicial,
1
x0 = v0 t0 + a t20 − a t20 + k2 .
2
Portanto, podemos escrever a equação horária do MRUV na forma bem conhecida,
1
x(t) = x0 + v0 (t − t0 ) + a(t − t0 )2 .
2
Na tabela abaixo registramos algumas integrais indefinidas razoavelmente simples que

vão aparecer em problemas de fı́sica. Note que a e k são constantes arbitrárias. Não deixe
de verificar que está tudo correto, conferindo com a tabela das derivadas.
Tabela 4.1: Algumas integrais indefinidas
R
f (x) G(x) = f (x)dx
1
xn , n 6= −1 xn+1 + k
n+1
1
sen (ax) − cos(ax) + k
a
1
cos(ax) sen (ax) + k
a
1
exp(ax) exp(ax) + k
a
1
ln x + k
x
1
senh (ax) cosh(ax) + k
a
1
cosh(ax) senh (ax) + k
a
4.5 Cálculo de integrais definidas

A partir do “teorema fundamental do cálculo” temos
Z x
F (x) = G(x) − G(a) = f (t) dt,
a
onde G(x) é uma primitiva genérica de f (x). Então, para x = b, vem

Z b
G(b) − G(a) = f (t) dt.
a
Para calcular uma integral definida basta achar uma primitiva G(x) e encontrar os
seus valores nos extremos do intervalo de integração.
É comum utilizarmos a notação
Z b
b
G(b) − G(a) = G(x) | =| dG,
a a
Com essa nova notação, temos

Z b
b
f (x)dx = G(x) || = G(b) − G(a).
a a
Exemplos
(a) Calcular a área hachurada na figura 4.7, definida pela função

2
f (x) = ,
x2
entre x = 1 e x = 2.
f (x)
2 11111111
00000000
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
1 00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
0.5 00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
00000000
11111111
1 2 x
Figura 4.7
Essa área hachurada será dada por

Z 2 Z 2
2
A= f (x)dx = dx.
1 1 x2
Mas a primitiva de f (x) = 2/x2 é G(x) = −2/x + k. Note que no cálculo das
integrais definidas podemos omitir a constante k, que vai ser sempre cancelada na
diferença G(b) − G(a), onde a e b são os limites de integração. Então temos
2
2 2
A= − + k = − + k − (−2 + k) = 1.
x 1 2
(b) Calcular o trabalho executado pela força F (x) = −4x+x2 no percurso entre x1 = −1
e x2 = +1 (por simplididade, não estamos nos preocupando com as unidades corretas
de distância, trabalho, etc).
Nessas condições, o trabalho é dado por
Z x2 Z +1
W1→2 = F (x)dx = (−4x + x2 )dx
x1 −1
+1
2 1 3 1 1 2
= −2x + x = −2 + − −2 − = .
3 −1 3 3 3
(c) Calcular a área entre a parábola y(x) = 4 − x2 e o eixo x no intervalo −2 ≤ x ≤ 2.

A área vai ser dada pela expressão
Z +2 +2

2 1 3

A= 4 − x dx = 4x − x =
−2 3 −2

8 8 32
= 8− − −8 + = .
3 3 3
(d) A equação de estado de um mol de um fluido é P V = RT , onde P é a pressão, V

é o volume, R = 8, 3 J/(K·mol) é a constante universal dos gases. Num processo
termodinâmico isotérmico, à temperatura T = 300 K, o gás se expande de um
volume inicial VA = 2 ℓ até um volume final VB = 4 ℓ. Calcule o trabalho WA→B
realizado para ir de A até B.
ZVB ZVB
RT VB VB
WA→B = P dV = dV = RT ln(V ) || = RT ln
V VA VA
VA VA
3
= 2, 49 × 10 ln(2) J .
4.6 As funções logaritmo e exponencial

A partir do que apresentamos até agora não é difı́cil concluir que integrar é bem mais
difı́cil do que derivar. De fato, a derivada de qualquer combinação de funções simples
(funções trigonométricas, potências, etc) é sempre uma combinação de funções simples
e pode ser calculada facilmente usando as regras de derivação que apresentamos nas
tabelas 3.1 e 3.2 nas páginas 10 e 13, respectivamente. Derivada é o nome de um tipo
de limite, e calculamos seu valor calculando o valor desse limite. Em contraste, não
existem regras como as apresentadas nas tabelas 3.1 e 3.2 para calcular integrais. Integral
também é o nome dado a um tipo de limite, mas em geral é impossı́vel calcular o valor
desse limite diretamente. Calculamos uma integral quando somos capazes de encontrar
a sua primitiva. A maioria das técnicas de integração apenas transformam o integrando
original de modo inteligente, conveniente, para que se torne mais simples vislumbrar uma
primitiva. Além disso, a maioria das integrais com integrandos que são combinações de
funções simples não possuem primitivas que são combinações de funções simples. Isto não
quer dizer que a integral não existe (ou seja, que o limite não existe), mas apenas que ela
não pode ser expressa como uma combinação de funções simples. Se uma integral deste
tipo aparece com muita frequência, então ela recebe um nome especial. Este é o caso da
função logaritmo.
A função logaritmo na base e, ln(x), também conhecida como logaritmo neperiano, ou
simplesmente logaritmo, é definida através da integral
Zx
1
ln(x) = dt, com x > 0. (4.1)
t
1
Da expressão acima decorre imediatamente que
ln(1) = 0.
Do “teorema fundamental do cálculo” vem
d ln(x) 1
= (4.2)
dx x
Olhando o gráfico da função ln(x) na página 10 não é difı́cil ver que esta função
tem uma inversa (geometricamente, a inversa pode ser obtida rebatendo-se o gráfico da
função através da bissetriz do primeiro e terceiro quadrandes). A função exponencial
exp(x), também escrita como ex , é por definição a inversa da função logaritmo:
exp ≡ ln−1 . (4.3)

Vamos calcular a derivada da exponencial. Pela definição de inversa, temos
d ln(ex )
ln(ex ) = x =⇒ = 1.
dx
Usando a “regra da cadeia” (propriedade (5) na pág. 13), podemos calcular a derivada,
d ln(ex ) d ln(y) | dex 1 dex dex

= = = 1 =⇒ = ex . (4.4)
dx dy |y=ex dx ex dx dx
A função exponencial e sua inversa, a função logaritmo, são extremamente importan-
tes. Elas aparecem com muita frequência em todas as áreas da fı́sica e da matemática. É
importante que você se familiarize com seus gráficos (veja a página 10) e as propriedades
listadas na tabela 4.2.
Vamos demonstrar que o logaritmo do produto é a soma dos logaritmos.
Logaritmo do produto.
ln(xz) = ln(x) + ln(z). (4.5)
Demonstração:
Defina a função f (x) ≡ ln(xz), onde a variável z é mantida fixa e só x pode variar.
Temos
df (x) d ln(xz) d ln(y) | d(xz) 1 1

= = = z = .
dx dx dy |y=xz dx xz x
Lembre que d ln(x)/dx = 1/x, ou seja ln(x) e f (x) têm a mesma derivada. Duas
funções que têm a mesma derivada ou são iguais ou diferem apenas por uma
constante aditiva. Portanto,
Tabela 4.2: Algumas propriedadades das funções ln(x) e ex
Para x, y, z > 0 e a, b reais quaisquer
ln(xz) = ln(x) + ln(z),

ln(xa ) = a ln(x),
ln(1/x) = − ln(x),
ln(1) = 0
ea+b = ea · eb ,
1/ea = e−a ,
(ea )b = eab
e0 = 1
ln(ea ) = a,
eln(x) = x.
f (x) ≡ ln(xz) = ln(x) + const. (4.6)

Para descobrir o valor da constante basta lembrar que ln(1) = 0. Colocando x = 1
na equação acima vem
ln(z) = ln(1) + const. = const. =⇒ const. = ln(z). (4.7)

Substituindo a equação (4.7) na equação (4.6) completamos a demonstração.
4.7 Algumas técnicas de integração

4.7.1 Integral de uma derivada
Do “teorema fundamental do cálculo” decorre imediatamente que
Zb
dF (x) b
dx = F (x) || = F (b) − F (a),
dx a
a
uma vez que segundo este teorema o integrando é a derivada da primitiva.

Esta é uma maneira ligeiramente diferente de escrever o TFC. É claro que o TFC
garante que o integrando é sempre a derivada de uma função. Porém, como já vimos,
nem sempre esta função pode ser escrita em termos de combinações de funções simples. Se
conseguirmos escrever o integrando como uma derivada de uma função simples, a integral
se torna trivial.
Se a integral for indefinida temos simplesmente
Z
dF (x)
dx = F (x),
dx
onde omitimos, como de costume, a constante de integração.
Exemplo
No cálculo do potencial eletrostático sobre o eixo z devido a um disco de raio R,

uniformemente carregado com densidade superficial de carga σ, colocado sobre o
plano xy e com centro na origem do sistema de coordenadas, a integral relevante é
ZR
σ s
V (z) = √ ds.
2ǫ0 z 2 + s2
0
Neste caso o integrando é a derivada de uma função simples:
ZR √
σ d ( z 2 + s2 ) σ √ 2 R σ √ 2
V (z) = ds = z + s2 || = ( z + R2 − |z|).
2ǫ0 ds 2ǫ0 0 2ǫ0
0
4.7.2 Integração por partes

Vamos usar agora a fórmula da derivada de um produto de funções para obter mais uma
técnica útil para o cálculo de integrais. Já mostramos que
d(f1 (x) f2 (x)) df1 (x) df2 (x)

= f2 (x) + f1 (x) ,
dx dx dx
portanto
Zb Zb Zb
df1 (x) d(f1 (x) f2 (x)) df2 (x)
f2 (x) dx = dx − f1 (x) dx.
dx dx dx
a a a
Note que o integrando da primeira integral do lado direito da equação é uma derivada.
Assim
Zb Zb
df1 (x) | b df2 (x)
f2 (x) dx = f1 (x) f2 (x) | − f1 (x) dx.
dx a dx
a a
Para que esta propriedade seja útil precisamos escrever o integrando como um produto
(df1 /dx) f2 de tal forma que (1) a derivada de df2 /dx seja mais simples do que f2 e que
consigamos integrar facilmente df1 /dx afim de obter f1 , e (2) a integral do lado direito da
equação acima seja mais simples do que a integral original.
A propriedade acima também vale para integrais indefinidas. Neste caso temos (omi-
tindo constantes de integração)
Z Z
df1 (x) df2 (x)
f2 (x) dx = f1 (x) f2 (x) − f1 (x) dx.
dx dx
Exemplos
(a) Calcule a integral Z

x3 ln(x) dx.
Colocamos f2 (x) = ln(x) e df1 (x)/dx = x3 =⇒ f1 (x) = x4 /4. Assim, lembrando

que d ln(x)/dx = 1/x, obtemos
Z Z 4
3 x4 x 1 x4 x4
x ln(x) dx = ln(x) − dx = ln(x) −
4 4 x 4 16
(b) Calcule a integral Z

x sen (x) dx.
Escolhemos f2 (x) = x e df1 (x)/dx = sen (x) =⇒ f1 (x) = − cos(x) . Portanto,

Z Z
x sen (x) dx = −x cos(x) + cos(x) dx = −x cos(x) + sen (x).
(c) Calcule a integral Z

x3
dx.
(x2 + 4)3/2
x 1
Escolhemos f2 (x) = x2 e df1 (x)/dx = =⇒ f1 (x) = − 2 . Por-
(x2 + 4) 3/2 (x + 4)1/2
tanto,
Z Z
x3 x2 2x
2 3/2
dx = − 2 1/2
+ dx
(x + 4) (x + 4) (x + 4)1/2
2
x2
= − 2 1/2
+ 2(x2 + 4)1/2 ,
(x + 4)
onde calculamos a integral no lado direito da equação acima percebendo que o seu
integrando (sem o fator 2) é a derivada de (x2 + 4)1/2 .
(d) Considere a distribuição de probabilidades
P (x) = Ax exp(−2x),
definida para valores não negativos da variável aleatória x (isto é, para 0 ≤ x < ∞).
Qual é o valor da constante A para que essa distribuição seja normalizada?
Para a distribuição ser normalizada,
Z∞ Z∞
P (x)dx = 1 ⇐⇒ Ax exp(−2x) dx = 1.
0 0
1
Escolhemos f2 (x) = Ax e df1 (x)/dx = exp(−2x) =⇒ f1 (x) = − exp(−2x). Por-
2
tanto,
Z∞ Z∞
Ax ∞ A
Ax exp(−2x) dx = 1 ⇐⇒ − exp(−2x) || + exp(−2x) dx = 1
2 0 2
0 0
A ∞ A
⇐⇒ − exp(−2x) || = 1 ⇐⇒ = 1 =⇒ A = 4.
4 0 4
4.7.3 Mudança de variável de integração

Integrais definidas
Uma das técnicas mais versáteis para calcular integrais é a mudança de variável de inte-
gração. Ela é baseada na igualdade
Zy(b) Zb
dy(x)
f (y) dy = f (y(x)) dx. (4.8)
dx
y(a) a
Demonstração:
Seja F (y) a primitiva de f (y). Isto significa que
dF (y)
= f (y). (4.9)
dy
Se conhecemos a primitiva de f (y), o TFC nos fornece imediatamente o valor da
integral do lado esquerdo da equação (4.8):
Zy(b)
y(b)
f (y) dy = F (y) || = F (y(b)) − F (y(a)). (4.10)
y(a)
y(a)
Para calcular a integral do lado direito da equação (4.8) precisamos achar a primitiva
da função f (y(x)) dy(x)/dx. Vamos mostrar que esta primitiva é F (y(x)), onde F (y) é a
primitiva de f (y). De fato, usando a fórmula da derivada de função de função (regra da

cadeia) obtemos
dF (y(x)) dF (y) dy(x) dy(x)

= = f (y) ,
dx dy dx dx
onde usamos a equação (4.9). Portanto,
Zb
dy(x) b
f (y(x)) dx = F (y(x)) || = F (y(b)) − F (y(a)). (4.11)
dx a
a
Comparando a equação (4.10) e a equação (4.11) vemos que as integrais na equação

(4.8) são iguais, como querı́amos demonstrar. ❚
Observações:
Obviamente a equação (4.8),
Zy(b) Zb
dy(x)
f (y) dy = f (y(x)) dx,
dx
y(a) a
pode ser usada tanto da direita para a esquerda como da esquerda para a direita.
É conveniente usar a equação (4.8) da direita para a esquerda quando for fácil ver que
o integrando tem a forma f (y) dy/dx e identificar y(x). Neste caso, após mudarmos os
limites de integração: a → y(a) e b → y(b), fazemos as substituições
dy(x)
y(x) → y e dx → dy , (4.12)
dx
veja os exemplos (a) e (b) abaixo.
Geralmente a equação (4.8) é usada da esquerda para a direita. Neste caso, precisamos
escolher uma função y(x) que leve a uma integral mais simples. Além disto, precisamos
ser capazes de achar a inversa y −1 da função y(x) escolhida e expressar x em função de
y para determinar os limites de integração da integral em x que aparece no lado direito
da equação (4.8). Após mudarmos os limites de integração: y(a) → a = y −1 (y(a)) e
y(b) → b = y −1 (y(b)), fazemos as substituições
dy(x)
y → y(x) e dy → dx , (4.13)
dx
veja os exemplos (c) e (d) abaixo.
Na prática trabalhamos de uma maneira mais informal e substituimos as equações
(4.12) e (4.13) respectivamente por
dy(x)
y(x) = y, dx = dy, limites de integração: y(a) e y(b) , (4.14)
dx
dy(x)
y = y(x), dy = dx, limites de integração: a = y −1 (y(a)) e b = y −1 (y(b)) , (4.15)
dx
como mostramos nos exemplos abaixo.
Exemplos
(a) Calcule a integral

Zb
sen 3 (x) cos(x)dx.
a
Observe que cos(x) = d sen (x)/dx. Isto sugere colocar y(x) = sen (x), usar a
equação (4.8) da direita para a esquerda fazendo sen (x) → y e cos(x)dx → dy.
Zb sen(b)
Z
3 y 4 |sen(b) sen 4 (b) sen 4 (a)
sen (x) cos(x)dx = y 3 dy = = − .
4 |sen(a) 4 4
a sen(a)
Vamos refazer o cálculo usando a equação (4.14):
d sen (x)
sen (x) = y, dx = cos(x)dx = dy, limites: y(a) = sen (a), y(b) = sen (b),
dx
Zb sen(b)
Z
y 4 |sen(b) sen 4 (b) sen 4 (a)
=⇒ sen 3 (x) cos(x)dx = y 3 dy = = −
4 |sen(a) 4 4
a sen(a)
e reobtemos o resultado anterior.
(b) Calcule a integral

Z3
ln(x)
dx.
x
2
Note que d ln(x)/dx = 1/x. Assim, uma boa escolha consiste em colocar y(x) =
ln(x), usar a equação (4.8) da direita para a esquerda e a equação (4.14).
d ln(x) 1
ln(x) = y, dx = dx = dy, limites: ln(2), ln(3),
dx x
Z3 ln(3)
Z
ln(x) 1 2 |ln(3) 1
=⇒ dx = y dy = y | = [ln2 (3) − ln2 (2)].
x 2 ln(2) 2
2 ln(2)
(c) Calcule a integral

Z1
1
p dy.
1 − y2
1/2
Uma boa escolha é y(x) = sen (x) porque a igualdade 1− sen 2 (x) = cos2 (x) permite
simplificar o denominador. Usamos a equação (4.8) da esquerda para a direita e a
equação (4.15).
d sen (x)
y = sen (x), dy = dx = cos(x)dx,
dx
x = arcsen (y), limites: arcsen (1/2) = π/6 e arcsen (1) = π/2
Z1 Zπ/2 Zπ/2
1 1 1
=⇒ p dy = p cos(x) dx = cos(x) dx
1 − y2 1 − sen 2 (x) cos(x)
1/2 π/6 π/6
Zπ/2
π/2 π
= dx = x || = ,
π/6 3
π/6
(d) Calcule a integral

Z5
1
dy.
y ln(y)
3
Escolhemos y(x) = ex , usamos a equação (4.8) da esquerda para a direita e a equação

(4.15).
dex
y = ex , dy = dx = ex dx, x = ln(y), limites: ln(3) e ln(5)
dx
Z5 ln(5)
Z ln(5)
Z
1 1 1 ln(5)
=⇒ dy = ex dx = dx = ln(x) ||
y ln(y) ex x x ln(3)
3 ln(3) ln(3)
= ln(ln(5)) − ln(ln(3)).
Note que neste exemplo também é fácil usar a equação (4.8) da direita da esquerda,
como fizemos no exemplo (a). Deixaremos esta resolução como um exercı́cio.
Por razões didáticas, quando usamos a equação (4.8) da esquerda para a direita
colocamos a variável de integração da integral inicial igual a y. Quando trabalhamos
da direita para a esquerda colocamos a variável de integração igual x. Obviamente, o
nome dado à variável de integração é irrelevante. Com a prática você vai usar a equação
(4.8) sem se preocupar com estes detalhes.
Integrais indefinidas
Para terminar, apresentamos a fórmula de mudança de variável de integração para inte-
grais indefinidas.
Se você está usando a equação (4.8) da esquerda para a direita, você sai de uma
integral em y e deve chegar no final numa integral em y. Assim, é conveniente escrever
Zy Zx(y)
′ ′ dy(x′ ) ′
f (y ) dy = f (y(x′ )) dx , (4.16)
dx′
onde x(y) é a inversa de y(x), isto é y(x(y)) = y.

Se você está usando a equação (4.8) da direita para a esquerda, a integral de saı́da é
em x e a integral de chegada também deve ser em x. É conveniente colocar
Zx Zy(x)
dy(x′ ) ′
f (y(x′ )) dx = f (y ′ ) dy ′ . (4.17)
dx′
Os limites superiores são os valores onde calculamos as primitivas. Note que nas
equações (4.16) e (4.17) o limite superior na integral do lado direito é igual à variável nova
expressa em termos da variável de saı́da. As linhas em x e y foram colocadas para deixar
clara a diferença entre a variável de integração e o valor onde calculamos a primitiva. Na
prática estas linhas podem ser omitidas, como fazemos nos exemplos abaixo.
É possı́vel usar as equações (4.16) e (4.17) sem os limites superiores. Neste caso, fica
implı́cito que no final dos cálculos você deve voltar à variável de saı́da (veja o exemplo
(e) abaixo). Preferimos escrever o limite superior nas integrais porque isto ajuda a seguir
as mudanças de variável de integração, especialmente no caso em que várias mudanças de
variável são feitas em seguida.
Demonstrações:
A demonstração da equação (4.16) é análoga à da equação (4.8). Já mostramos que
se F (y) é a primitiva de f (y) então F (y(x)) é a primitiva de f (y(x)) dy(x)
dx
(veja a página
35). Portanto,
Zy
y
f (y ′ ) dy ′ = F (y ′ ) || = F (y),
Zx(y)
dy(x′ ) ′ |x(y) = F (y(x(y))) = F (y).
f (y(x′ )) dx = F (y(x ′
)) |
dx′
Assim, os dois lados da equação (4.16) são iguais. ❚

Para demonstrar a equação (4.17) basta verificar que a derivada em x da integral no
lado direito é igual ao integrando da integral no lado esquerdo. Isto significa que a integral
do lado direito é a primitiva de f (y(x))dy(x)/dx. Usando o TFC e a “regra da cadeia”
obtemos
 y(x)   z 
Z Z
d  ′ ′ d  ′

′  dy(x) dy(x)
f (y ) dy = f (y ) dy = f (y(x)) .
dx dz dx dx
z=y(x)
Esta derivada é igual ao integrando da integral do lado esquerdo da equação (4.17), como
querı́amos mostrar. ❚
Exemplos
(e) Calcule a integral Z

1
p dy.
1 − y2
dy(x)
Colocamos y(x) = sen (x) =⇒ x(y) = sen −1 (y) = arcsen (y); dy → dx
dx =
d sen (x)
dx
dx = cos(x) dx. Portanto,
Zy arcsen
Z (y)
1 1
p dy = p cos(x) dx
1 − y2 1 − sen 2 (x)
arcsen
Z (y) arcsen
Z (y)
1 arcsen (y)
= cos(x) dx = dx = x || = arcsen (y).
cos(x)
Se você tivesse omitido o limite de integração você teria obtido

Z Z
1 1
p dy = p cos(x) dx
1 − y2 1 − sen 2 (x)
Z Z
1
= cos(x) dx = dx = x,
cos(x)
que só faz sentido se você substituir x por arcsen (y) no final.
(f ) Calcule a integral Z
sen 3 (x) cos(x)dx.
dy(x)
Definimos y(x) = sen (x), portanto cos(x) dx = dx → dy e
dx
Zx sen
Z (x)
y 4 | sen (x) sen 4 (x)
sen 3 (x) cos(x)dx = y 3 dy = = .
4 | 4
(g) Mostre que

Z
1 1
dy = ln(ay + b) .
ay + b a
x b dy(x) 1
Definimos x(y) = ay + b =⇒ y = − e dy → dx = dx
a a dx a
Zy Z
ay+b
1 1 1 1 ay+b 1
dy = dx = ln(x) || = ln(ay + b).
ay + b x a a a
4.8 O que fazer quando nada funciona?

Apresentamos neste livro as principais técnicas de integração. Você irá aprender muitas
outras nos cursos de Cálculo. Porém, muitas vezes você não conseguirá calcular a primitiva
de uma integral. Como já dissemos, a maioria das integrais não pode ser calculada em
termos de funções simples. Como ter certeza de que não existe um truque que permite o
cálculo da integral? A rigor, é impossı́vel ter esta certeza. Porém, quando tudo o mais
falha vale a pena consultar uma tabela de integrais, onde estão listadas as primitivas de
várias funções e os valores de várias integrais definidas que aparecem com frequência. As
tabelas variam muito em tamanho. Uma tabela que contém uma boa escolha de integrais
e que vai ser útil durante toda a graduação é a da Coleção Schaum [6]. Uma das tabelas
mais completas é a de Gradshteyn e Ryzhik [7] cuja sexta edição tem mais de mil páginas!
Se uma integral não puder ser encontrada no Gradshteyn e Ryzhik muito possivelmente
não existe uma expressão analı́tica para ela. Na internet, o excelente site da Wolfram
Alpha [8] permite calcular interativamente integrais definidas e indefinidas.
Mesmo não achando sua integral nas tabelas, ainda existe um último recurso—o cálculo
numérico. Quase sempre queremos calcular uma integral definida. Neste caso podemos
explorar a idéia de que a integral é essencialmente a área entre o eixo das abscissas e o
gráfico da função. É o que fazem as técnicas numéricas. Elas são extremamente eficientes
e fornecem rapidamente o resultado da integral com precisão desejada.
Apesar das soluções analı́ticas serem muito mais elegantes do que as soluções numéricas,
você não deve ter preconceitos contra o cálculo numérico. À medida em que você for
avançando em seus estudos você perceberá que os métodos computacionais desempenham
papel essencial na ciência e na tecnologia.
Exercı́cios
(1) Calcule as integrais

Z Z Z
1
(2x + x4 + e−x )dx; cos(ωt + φ)dt; dx;
Z p Z Z a + bx
z 1
y 1 + y 2 dy; dz; dt;
Z Z z2
+1 2
Z t +a
2
1 x
dt; dx; ln(x)dx.
t 2 − a2 (x2 + a2 )3/2
(2) Calcule o trabalho da força F (x) = −kx para deslocar um corpo que se move apenas
ao longo do eixo x da posição x = x0 até a posição x = x1 .
(3) Calcule a área delimitada pelas funções y1 (x) = x2 e y2 (x) = −x entre x = 0 e

x = 1.
Sugestão: primeiramente esboce os gráficos de y1 (x) e y2 (x).
(4) Um corpo que se desloca apenas sobre o eixo dos x tem aceleração que varia linear-
mente com o tempo, a(t) = ct, onde c é uma constante. Sabendo-se que em t = 0 o
corpo está em x = x0 com velocidade v = v0 , determine a função x(t) que fornece a
posição do corpo no instante t.
Sugestão: primeiramente integre a aceleração para determinar a velocidade, impo-
nha a condição inicial e integre novamente para determinar x(t).
(5) Calcule o trabalho necessário para levar uma partı́cula de massa m da superfı́cie da
Terra até o infinito (com velocidade zero no infinito).
Dado: a força que a Terra exerce sobre a partı́cula é F (r) = GM m/r2 , onde M
é a massa da Terra e G é a constante gravitacional; considere o raio da Terra RT
conhecido.
(6) Considere uma partı́cula de massa m que se move ao longo do eixo x sujeita à ação
de uma força F (t). Através de uma mudança de variável apropriada mostre que o
impulso da força F (t) é igual à variação do momento da partı́cula, isto é,
Zt2
F (t)dt = mv(t2 ) − mv(t1 ) ,
t1
onde v(ti ) é a velocidade no instante t = ti .

Sugestão: lembre que a aceleração é a = dv/dt.
Capı́tulo 5
Vetores
5.1 Conceito de vetor

Existem muitas grandezas fı́sicas que não podem ser completamente descritas por um
simples número. Para descrever essas grandezas (como forças, deslocamentos, velocida-
des, etc) precisamos especificar um número (módulo), uma direção e um sentido. Essas
grandezas são denominadas vetoriais, em contraposição às grandezas escalares, que po-
dem ser caracterizadas por um número (como temperatura, energia, etc). O vetor é uma
entidade matemática associada a um módulo, uma direção e um sentido. Na natureza
existem também outras grandezas que são muito mais complexas, necessitando para sua
caracterização de entidades matemáticas mais complicadas que os vetores (por exemplo,
a tensão ou as deformações de um sólido anisotrópico são descritas por entidades deno-
minadas tensores).
Vetores já devem ter sido vistos no ensino médio. Não vamos rever as idéias mais
intuitivas. Também não vamos nos preocupar com um tratamento rigoroso ou detalhado
(que deve ser apresentado numa disciplina especı́fica). Vamos apenas rever os conceitos
básicos e introduzir uma notação mais prática a fim de facilitar as operações (soma,
multiplicações, derivação) com vetores.
5.2 Componentes e módulo de um vetor; versor

Apesar dos vetores serem independentes dos sistemas de coordenadas, é muito interessante
escrevê-los em termos de suas componentes num determinado sistema. A figura 5.1 indica
um vetor F~ que tem uma componente F~x ao longo da direção x.
F~
α
F~x x
Figura 5.1
43
CAPÍTULO 5. VETORES 44

A notação F~ designa o módulo de um vetor F~ . Então, é claro que

~ ~
Fx = F cos α.
Adotando um sistema de eixos cartesianos (isto é, de eixos ortogonais, como mostrado
na figura 5.2), todo vetor planar pode ser decomposto em componentes ao longo das
direções x e y.
y
F~y F~
0 F~x x
Figura 5.2
Usando a “regra da soma do paralelogramo” temos
F~ = F~x + F~y ,
ou seja, um vetor sempre pode ser escrito como a soma de suas componentes num sistema
de eixos ortogonais. Pelo teorema de Pitágoras temos
2 2 2
~
F = F~x + F~y .
Embora tenha sido apresentado um exemplo bidimensional, é claro que isto tudo
também funciona em três dimensões (para um sistema de eixos cartesianos x − y − z).
Para tornar mais simples a representação de um vetor é interessante introduzir a noção
de versor, que é um vetor de módulo unitário, funcionando como uma espécie de “unidade
de direção”. Na figura 5.3 representamos os versores ~ı e ~, que são vetores unitários nas
direções x e y respectivamente (na direção z costuma-se usar o sı́mbolo ~k). Note que
|~ı| = |~| = 1.
O vetor F~x pode então
ser
~
escrito como Fx = Fx~ı, onde Fx é um escalar, cujo módulo
~
corresponde ao módulo Fx . Da mesma forma temos
F~y = Fy~.
Portanto,
F~ = Fx~ı + Fy~.
Em três dimensões terı́amos
F~ = Fx~ı + Fy~ + Fz~k.
~
~ı x
Figura 5.3
Exemplo
Escrever em termos dos versores cartesianos os dois vetores ~v1 e ~v2 , indicados na
figura 5.4.
y y
2
1 ~v1 1 2 3 4
x x
1 2 3
-1
~v2
-2
Figura 5.4
é fácil perceber que

~v1 = 3~ı + 2~; ~v2 = 4~ı − 2~
e que
|~v1 |2 = 32 + 22 = 13; |~v2 |2 = 42 + 22 = 20.
√ √ √
Então, |~v1 | = 13 e |~v2 | = 20 = 2 5.
5.3 Operações com vetores

5.3.1 Soma ou subtração
Dados os vetores
~a = ax~ı + ay~ + az~k
e
~b = bx~ı + by~ + bz~k,
o vetor soma ou subtração é definido como
~a ± ~b = (ax ± bx )~ı + (ay ± by ) ~ + (az ± bz ) ~k.
Exemplo
Vamos considerar os vetores
~v1 = 3~ı + 2~; ~v2 = 4~ı − 2~.
A soma é dada por

~v1 + ~v2 = 7~ı.
é fácil verificar que os mesmos resultados poderiam ter sido obtidos através da “regra
do paralelogramo” (ou de qualquer outra regra geométrica desse tipo).
5.3.2 Produto de vetores

Há pelo menos três tipos de produtos envolvendo vetores:
(i) produto de um número escalar por um vetor, dando como resultado um vetor;
(ii) produto de um vetor por outro vetor, dando como resultado um escalar (é o
chamado produto escalar);
(iii) produto de um vetor por outro vetor, dando como resultado um terceiro vetor (é
o chamado produto vetorial).
Vamos considerar cada um desses casos.
(i) Produto de um vetor por um escalar

Dado o vetor
~a = ax~ı + ay~ + az~k
e o escalar A, temos
A~a = Aax~ı + Aay~ + Aaz~k.
Como exemplo, vamos considerar o vetor ~v1 = 3~ı + 2~. Multiplicando por 5, temos
15~ı + 10~; multiplicando por −6, obtemos −18~ı − 12~.
(ii) Produto escalar entre dois vetores

Dados os vetores
~a = ax~ı + ay~ + az~k
e
~b = bx~ı + by~ + bz~k,
definimos o produto escalar,
~a · ~b = ax bx + ay by + az bz .
Por exemplo, o produto escalar dos vetores ~v1 = 3~ı + 2~ e ~v2 = 4~ı − 2~ é dado por
~v1 · ~v2 = 3 × 4 + 2 × (−2) = 12 − 4 = 8.
Há uma forma alternativa, muito conveniente, de escrever o produto escalar entre
dois vetores. Dados ~a e ~b, formando um ângulo θ, é fácil mostrar que

~
~a · b = |~a| ~b cos θ.
Como ~a e ~b (para θ 6= 0 ou π) definem um plano, que pode ser chamado plano xy,
basta demonstrar esta relação no espaço cartesiano bidimensional. Vamos considerar
a figura 5.5.
y
~b
θ2
~a
θ1
x
Figura 5.5
O produto escalar entre ~a e ~b é dado por
~a · ~b = ax bx + ay by .
Mas
ax = a cos θ1 ; ay = a sen θ1 ;
bx = b cos θ2 ; by = b sen θ2 ,

com a = |~a| e b = ~b. Então
~a · ~b = ab cos θ1 cos θ2 + ab sen θ1 sen θ2 == ab [cos θ1 cos θ2 + sen θ1 sen θ2 ] =

= ab cos (θ1 − θ2 ) = ab cos (θ2 − θ1 ) = ab cos θ.
Por exemplo, vamos calcular o ângulo entre os vetores ~v1 = 3~ı + 2~ e ~v2 = 4~ı − 2~.
O produto escalar é dado por
~v1 · ~v2 = 3 × 4 + 2 × (−2) = 12 − 4 = 8.

Os módulos desse dois vetores são dados por

√ √
|~v1 | = (9 + 4)1/2 = 13 e |~v2 | = (16 + 4)1/2 = 2 5.
Então
~v1 · ~v2 8
cos θ = = √ = 0.496...,
|~v1 | |~v2 | 2 13 × 5
de onde vem θ ≈ 60, 255◦ .
O produto escalar entre dois vetores tem uma série de propriedades facilmente de-
monstráveis:

(a) ~a · b + ~c = ~a · ~b + ~a · ~c,
~
conhecida como associatividade;
(b) ~a · ~b = ~b · ~a,
que é conhecida como comutatividade;
(c) ~a · ~b = 0 =⇒ ~a ⊥ b,
~
indicando que o produto escalar é nulo quando os vetores forem perpendicula-
res;

(d) ~a · b = |~a| ~b =⇒ ~a k b,
~ ~
mostrando que se o produto escalar for dado pelo produto dos módulos (cos θ =
1) então os dois vetores são paralelos.
(iii) Produto vetorial entre dois vetores

Dados os vetores
~a = ax~ı + ay~ + az~k
e
~b = bx~ı + by~ + bz~k,
o produto vetorial é um terceiro vetor definido como
~a × ~b = (ay bz − az by )~ı + (az bx − ax bz ) ~ + (ax by − ay bx ) ~k.
Há várias formas de se lembrar dessa definição. É fácil verificar que essa expressão
do produto vetorial pode ser obtida através de um determinante simbólico,

~ı ~ ~k

~a × ~b = ax ay az .
bx by bz
Exemplo
Vamos considerar de novo os vetores ~v1 = 3~ı + 2~ e ~v2 = 4~ı − 2~. O produto vetorial
é dado por

~ı ~ ~k

~v1 × ~v2 = 3 2 0 = 0~ı + 0~ + (−6 − 8) ~k = −14~k.
4 −2 0
No caso do produto vetorial também há uma definição alternativa, muito conveniente
em aplicações fı́sicas, mas um pouco mais complicada. O ~
vetor produto ~a × b deve

ser especificado por: (i) um módulo, que é dado por |~a| ~b sen θ, onde θ é o ângulo
entre ~a e~b; (ii) uma direção, que é perpendicular ao plano
definido
por ~a e ~b; ou seja,
~a e ~b são ambos perpendiculares ao produto vetorial ~a × ~b ; e (iii) um sentido, que
é dado pela “regra do saca-rolhas” ou “regra da mão direita” (com a mão direita
acompanhamos o vetor ~a e tentamos atingir a ”ponta”do vetor ~b; o dedo polegar da
mão direita vai apontar no sentido do produto vetorial).
Considerando o plano xy formado pelos vetores ~a e ~b, é fácil verificar essa definição
alternativa. De fato, observando a figura 5.6, temos
y
~b
θ2
~a
θ1
x
Figura 5.6

~ı ~ ~k

~a × ~b = ax ay 0 = (ax by − ay bx ) ~k

bx by 0
= (ab cos θ1 sen θ2 − ab sen θ1 cos θ2 ) ~k
= ab (cos θ1 senθ2 − sen θ1 cos θ2 ) ~k = ab sen (θ2 − θ1 ) ~k,
de onde vem que

~a × ~b = |~a| ~b sen θ ~k.
Como o eixo z é perpendicular aos eixos x e y, e deve estar orientado para “fora do
papel”, não há dúvidas de que os itens (i), (ii) e (iii) da definição alternativa vão
ser devidamente satisfeitos.
Há uma série de propriedades imediatas do produto vetorial:

(a) ~a × ~b = −~b × ~a.

Cuidado: o produto vetorial não é uma operação comutativa;
(b) ~a × ~b = 0 =⇒ ~a k ~b,
ou seja, dois vetores paralelos têm produto vetorial nulo;

(c) ~a × ~b = ab =⇒ ~a ⊥ ~b,
ou seja, nesse caso ~a e ~b são dois vetores perpendiculares;

(d) ~a × ~b + ~c = ~a × ~b + ~a × ~c,
continuando válida a propriedade associativa em relação à soma;

(e) ~a × ~b × ~c = ~b (~a · ~c) − ~c ~a · ~b ,
que é uma propriedade mais complicada, que está sendo dada apenas para
registro, mas que não é difı́cil verificar (podemos usar um desses artifı́cios
mnemônicos para lembrar do resultado final: “bac” menos “cab”;

(e) ~a · ~b × ~c = ~b · (~c × ~a) = ~c · ~a × ~b ,
que se chama produto misto, e que também está sendo apenas registrado. Note
a propriedade cı́clica dessas espressões (em relação à ordem da letras a, b e c).
Note também que aparece um sinal − associado a cada inversão de ordem do
produto vetorial.
5.4 Funções vetoriais

Uma função vetorial é um vetor que depende de uma ou mais variáveis. Por exemplo, o
vetor posição ~r ou o vetor velocidade ~v são em geral funções do tempo t. Então temos as
funções vetoriais ~r = ~r (t) e ~v = ~v (t).
Movimento circular uniforme (MCU) em coordenadas cartesianas

Nesse movimento a trajetória é um cı́rculo, de raio constante R, e “ângulos iguais são
percorridos em tempos iguais”, ou seja, a “taxa de variação do ângulo com o tempo”, ou
derivada, dθ/dt, é constante.
Utilizando a notação da figura 5.7, o MCU (movimento circular uniforme) é definido
por
|~r| = R(constante)
e
dθ
θ = θ0 + ωt, com = ω,
dt
onde a constante ω é a “velocidade angular”. O vetor posição em coordenadas cartesianas
é dado por
~r = R cos (θ0 + ωt)~ı + R sen (θ0 + ωt) ~.
P
~r
θ
0 x
Figura 5.7
Note que apenas as componentes rx = rx (t) = R cos(θ0 + ωt) e ry = ry (t) = R sen (θ0 + ωt)
dependem do tempo (os versores ~ı e ~ são obviamente fixos, constantes).
A derivada de uma função vetorial é definida como um vetor formado pela derivada
de cada uma das suas componentes cartesianas. Assim temos
d~r drx (t) dry (t)
~v = = ~ı + ~
dt dt dt
= −R ω sen (θ0 + ωt)~i + R ω cos (θ0 + ωt)~,
que é o vetor velocidade, e
d~v dvx (t) dvy (t)
~a = = ~ı + ~
dt dt dt
= −R ω 2 cos (θ0 + ωt)~i − R ω 2 sen (θ0 + ωt) ~,
que é o vetor aceleração.
A partir das expressões de ~a e de ~r temos
~a = −ω 2~r,
que é a famosa aceleração centrı́peta, cujo módulo é dado por ω 2 |~r|, mas que está orientada
na direção radial, contrariamente a ~r. Também é fácil calcular o módulo da velocidade.
De fato, temos
|~v |2 = [−R ω sen (θ0 + ωt)]2 + [R ω cos(θ0 + ωt)]2 = R2 ω 2 .
Então
|~v | = v = ω R,
que é outra expressão famosa do MCU (v é a velocidade tangencial e ω é a velocidade
angular).
Vamos agora mostrar que ~v e ~r são perpendiculares, ou seja, que a velocidade ~v é
tangente à trajetória. Para isso, basta calcular o produto escalar,
~r · ~v = rx vx + ry vy = [R cos (θ0 + ωt)] [−R ω sen (θ0 + ωt)]
+ [R sen (θ0 + ωt)] [R ω cos(θ0 + ωt)] = 0.
Como ~r · ~v = 0, o vetor velocidade é tangente à curva da trajetória.
5.5 Sistema de coordenadas polares

O sistema cartesiano é muito prático e simples, mas em determinadas situações pode ser
mais conveniente usar outros sistemas de referência. Vale a pena enfatizar que no sistema
cartesiano os versores básicos ~ı, ~ e ~k são constantes, em módulo, direção e sentido,
facilitando as operações de derivação. Em outros sistemas a situação pode ser mais
complicada.
No caso do movimento circular é muito conveniente considerar um sistema de coordenadas
polares. De acordo com a figura 5.7, um ponto P pode ser representado pelo par (x, y)
de coordenadas cartesianas ou, alternativamente, pelo par (r, θ) de coordenadas polares.
y
~eθ
y ~er
P
θ
r
~ θ
0 ~ı x x
Figura 5.8
Como está indicado nesse gráfico, r é a distância entre o ponto P e a origem O.

Portanto, r é um número positivo, que varia de 0 até ∞. O ângulo θ, que o vetor posição
~r, dado por OP ~ , faz com o eixo x, varia entre 0 e 2π.
Nessa figura também estão representados os versores~ı e ~ (em coordenadas cartesianas)
e ~er e ~eθ (em coordenadas polares). Observe que ~er tem a direção e o sentido do vetor
posição ~r; ~eθ é normal a ~er , orientado no sentido de θ crescente. Aqui há uma questão de
notação. Nós estamos preferindo a notação ~er e ~eθ para designar os vetores unitários em
coordenada polares, mas há quem prefira simplesmente rb e θ. b Convidamos o(a) leitor(a)
a fazer a sua escolha.
Em coordenadas cartesianas é muito simples obter os deslocamentos elementares: fi-
xando x, variamos y e obtemos dy; fixando y, variamos x e obtemos dx. A área elementar
é dx dy. Em coordenadas polares é um pouco mais complicado: fixando θ, variamos r
e obtemos dr; fixando r, variamos θ e obtemos o deslocamento elementar rdθ. A área
elementar é dada por rdrdθ.
Considerando a figura 18, podemos projetar ~er e ~eθ nos eixos x e y para escrever
~er = cos θ~ı + sen θ~,

~eθ = − sen θ~ı + cos θ~.
Invertendo essas equações, também temos
~ı = cos θ~er − sen θ~eθ ,

~ = sen θ~er + cos θ~eθ .
Movimento circular uniforme (MCU) em coordenadas polares

A posição ~r é dada por
~r = R~er .
Só que o vetor unitário ~er é uma função do tempo,
~er = cos (θ0 + ωt)~ı + sen (θ0 + ωt) ~.
Podemos agora calcular a velocidade,

d~r d~er
~v = =R ,
dt dt
com
d~er
= −ω sen (θ0 + ωt)~ı + ω cos (θ0 + ωt) ~ = ω~eθ .
dt
Portanto, temos
~v = ωR~eθ ,
que é uma expressão bem mais compacta do que em coordenadas cartesianas. Além disso,
essa expressão mostra que a velocidade ~v é tangente à trajetória.
A aceleração é dada por
d~v d~eθ
~a = = ωR .
dt dt
Mas
d~eθ
= −ω cos (θ0 + ωt)~ı − ω sen (θ0 + ωt) ~ = −ω~er .
dt
Então
~a = −ω 2 R~er ,
que é a expressão da aceleração centrı́peta, como nós já tı́nhamos obtido antes.
Descrição de um movimento circular arbitrário

Um movimento circular arbitrário é descrito pelo vetor posição
~r = R~er ,
onde |~r| = R é uma constante (o raio do cı́rculo) e
~er = [cos θ (t)]~ı + [ sen θ (t)] ~,
em que θ (t) é uma função qualquer do tempo.

No movimento circular uniforme (MCU), θ(t) = θ0 + ωt, onde ω é constante. No
movimento circular uniformemente acelerado (MCUA),
1
θ(t) = θ0 + ω0 t + αt2 ,
2
onde ω0 e α são constantes (α é a aceleração angular constante). No caso geral, dθ/dt =
ω (t) e dω/dt = α (t) são funções arbitrárias do tempo.
A velocidade é dada por

d~r d~er
~v = =R ,
dt dt
onde
d~er d d dθ dθ
= [cos θ(t)]~ı + [ sen θ(t)] ~ = − sen θ ~ı + cos θ ~
dt dt dt dt dt
= − [ω(t) sen θ(t)]~ı + [ω(t) cos θ(t)] ~ = ω (t) ~eθ .
Então
~v = R ω(t) ~eθ ,
mostrando que a velocidade permanece tangente à trajetória.
A aceleração é dada pela derivada dessa última expressão,
d~v dω d~eθ
~a = = R ~eθ + R ω .
dt dt dt
Mas
d~eθ d d dθ dθ
= [− sen θ(t)]~ı + [cos θ(t)] ~ = − cos θ ~ı − sen θ ~
dt dt dt dt dt
dθ dθ
= − [cos θ~ı + sen θ~] = − ~er .
dt dt
Então
dω
~a = R ~eθ − R ω 2~er .
dt
Portanto, quando dω/dt = α 6= 0, além da componente centrı́peta há também uma
componente tangencial da aceleração.
Exercı́cios
1. Desenhe, em um sistema cartesiano, os vetores ~u = 3~ı + 2 ~ e ~v = −3~ı − ~ .
2. Dados os vetores ~a = 4~ı + 8 ~ e ~b = 3~ı − 3 ~ − ~k, calcule:

1
(a) ~a;
8
(b) a componente y do vetor ~b;
(c) ~a + ~b e ~a − ~b;
(d) o módulo do vetor ~b;
(e) o produto escalar ~a . ~b ;
(f) o ângulo entre ~a e ~b.
3. Calcule a primeira e a segunda derivada dos vetores abaixo:
(a) ~y (t) = (−3t3 )~ı + (2t) ~;

(b) ~r(t) = −2 ~k;
(c) ~u(t) = A cos(ωt)~ı.
4. Calcule a integral Z t
~v (t′ ) dt′
0
quando (a) ~v (t) = A t2 ~ı e (b) ~v (t) = −ω sen (ωt)~ı + ω cos(ωt) ~.
5. Mostre que se o módulo da velocidade é constante, ou o vetor velocidade é constante,

ou o vetor aceleração é perpendicular ao vetor velocidade.
Capı́tulo 6
Expansões em séries de potências.
6.1 Definições, séries geométrica e de Taylor

Desde a Antiguidade, muitos matemáticos, fı́sicos e astrônomos se preocuparam com
problemas relacionados com sequências e séries. Talvez a primeira referência a esses
problemas seja o “paradoxo do movimento”, ou paradoxo de Aquiles que não consegue
alcançar a tartaruga, proposto por Zenão de Eléia, que viveu na Grécia há cerca de
450 anos AC. Essas questões já envolviam as idéias de limite, soma de séries infinitas e
quantidades infinitesimais. O desenvolvimento dos conceitos relacionados com sequências
e séries, em particular a soma de séries infinitas, ocupou inúmeros matemáticos desde
então, confundindo-se com o desenvolvimento do conceito de limite e do próprio cálculo
diferencial e integral.
Uma sequência é um conjunto de números, em geral relacionados através de uma certa
regra. Um exemplo famoso é a “sequência de Fibonacci”, que é uma sequência de números
inteiros na qual cada número, a partir do terceiro, é igual à soma dos dois números ante-
riores (1, 1, 2, 3, 5, 8, ...). Esta sequência, que tem muitas propriedades curiosas, continua
sendo utilizada em várias áreas da matemática e da ciência.
Uma série é um conjunto ordenado de infinitos termos, relacionados entre si por algum
tipo de operação. Se a diferença entre as sucessivas parcelas de uma série for constante,
temos uma série aritmética, ou PA (progressão aritmética) ilimitada, que é um exemplo fa-
miliar, certamente estudado no ensino médio. Outro exemplo muito conhecido é uma série
na qual a razão entre dois termos consecutivos é constante. Essa é uma série geométrica
ou PG (progressão geométrica) infinita. Existem vários outros tipos de séries, importan-
tes tanto para matemáticos como para fı́sicos ou engenheiros. Dizemos que uma série
a1 , a2 , a3 , ... é convergente quando a sequência de somas parciais, S1 = a1 , S2 = a1 + a2 ,
S3 = a1 + a2 + a3 , ... converge para um valor definido S. À medida que n cresce, caso
a sequência de somas parciais Sn seja oscilante, ou se aproxime de ± ∞, dizemos que a
série é divergente.
No ensino médio aprende-se que a soma de todos os (infinitos) termos de uma PG
infinita é dada por
X∞
a1
S= an = ,
n=1
1−q
desde que a razão q = an+1 /an seja em módulo menor do que a unidade (−1 < q < +1).
Quando |q| ≥ 1 a PG diverge.
56
CAPÍTULO 6. EXPANSÕES EM SÉRIES DE POTÊNCIAS. 57
Uma expansão em série de potências de uma função f (x) é uma representação desta
função como uma soma, em geral com infinitas parcelas, onde cada parcela é uma potência
da variável x. Em termos mais gerais, as parcelas de uma expansão em série podem ser
outras funções de x, em geral elementares.
No exemplo particular da soma da PG, lembrando que
an = a1 q n−1 ,
temos a série de potências

X ∞
1
f (x) = = xn−1 = 1 + x + x2 + x3 + ....
1 − x n=1
A possibilidade de representar algumas funções (e não apenas números!) através de

séries de potências ou de algumas outras séries infinitas foi reconhecida e trabalhada pela
tradicão matemática indiana. No inı́cio do século XVIII, o matemático inglês Taylor
desenvolveu um método geral para construir a expansão em série de potências de uma
função genérica bem comportada conhecida como série de Taylor. Em muitos casos uma
função f (x) infinitamente diferenciável e definida num intervalo aberto (a − ∆, a + ∆),
centrado no ponto a, pode ser escrita como uma série de potências infinita, dada pela
expressão
X∞
1 (n)
f (x) = f (a) (x − a)n , (6.1)
n=0
n!
onde
d n
f (x)
f (n) (a) = (6.2)
dxn x=a
é a n-ésima derivada de f (x) no ponto a, e
n! = n × (n − 1) × (n − 2) × · · · × 1
é o fatorial de n. Para a = 0, temos um caso especial conhecido como série de Maclaurin.

Além de perceber que muitas funções podem ser representadas como um polinômio de
grau infinito, A + B x + C x2 + D x3 + . . . , como os seus predecessores também sabiam,
Taylor descobriu a fórmula geral para encontrar os infinitos coeficientes A, B, C, . . . , desse
polinômio. As condições necessárias para que a série de Taylor convirja serão estudas nos
cursos de Cálculo.
Exemplo: a função cos(x)
Como exemplo concreto, vamos obter a expansão em série de Taylor da função f (x) =
cos x, para a = 0. Precisamos então calcular as deridadas f (n) (0) para vários valores de
n. De acordo com essa notação, f (0) (0) é a própria função calculada para x = 0, ou seja,
f (0) (0) = 1. Além disso, temos

(1) df
f (0) = = − sen x |x=0 = 0,
dx x=0

d2 f
f (2) (0) = = − cos x |x=0 = −1,
dx2 x=0

d3 f
f (3) (0) = = sen x |x=0 = 0,
dx3 x=0

d4 f
f (4) (0) = = cos x |x=0 = 1,
dx4 x=0
e assim por diante. Então temos

1 2 1 4
f (x) = cos x = 1 − x + x − ...
2! 4!
Portanto, na vizinhança do ponto x = 0, a função f (x) = cos x pode ser aproximada por
polinômios de grau n, com uma precisão que aumenta com o valor de n. A figura 6.1
mostra os gráficos da função f (x) = cos x e dos polinômios gerados pela série de potências
para diversos valores de n (n = 2, n = 4, n = 6, n = 10 e n = 20).
cos x
n=2
2 n=4
n=6
n = 10
n = 20
−3π −2π −π 0 π 2π 3π
−1
−2
Figura 6.1
Uma expansão em série convergente pode ser usada para que se obtenha uma forma
aproximada de uma função, eventualmente complicada, em termos de um polinômio de
grau relativamente baixo. Quando dizemos, por exemplo, que para θ pequeno, sen θ ≈ θ,
estamos aproximando a função pelo primeiro termo da série de Taylor nas vizinhanças da
origem! Levando em conta os termos seguintes dessa série,
1 3 1 5
sen θ = θ − θ + θ + ...,
6 120
podemos ir além dessa aproximação de primeira ordem. Quando θ for muito pequeno
(θ << 1), θ3 /6 vai ser muito menor do que θ, e o resultado aproximado em primeira
ordem,
sen θ ≈ θ
deve funcionar muito bem. No entanto, se θ não for tão pequeno, talvez seja necessário
adotar a aproximação seguinte,
1
sen θ ≈ θ − θ3 ,
6
ou até mesmo apelar para polinômios de grau superior a 3.
Como polinômios são fáceis de derivar, integrar ou mesmo interpretar, a substituição
de uma função por suas expansões em série tornou-se recurso muito importante. Esse
é um método interessante para calcular integrais complicadas, para resolver equações
diferenciais, para entender o comportamento assintótico de uma função, para visualizar
como ela cresce ou decresce, ou para encontrar a solução aproximada de algum problema
de interesse. No decorrer do seu curso devem surgir muitas aplicações das séries de Taylor.
Registramos agora algumas séries de Taylor bem conhecidas:
(1) exponencial de x,
X∞
1 1 xn
e x = 1 + x + x2 + x3 + . . . = , para qualquer x;
2 6 n=0
n !
(2) logaritmo de 1 + x,
X∞
1 2 1 3 (−1)n n+1
ln(1 + x) = x − x + x − . . . = x , para |x| < 1;
2 3 n=0
(n + 1)
(3) série geométrica,

∞
X
1
= xn , para |x| < 1;
1 − x n=m
(4) 1 + x elevado a um α real qualquer,

1
(1 + x)α = 1 + αx + α(α − 1)x2 + . . . , para |x| < 1;
2
(5) seno de x,
∞
X (−1)k 2k+1 1 1
sen x = x = x − x3 + x5 − . . . , para qualquer x;
k=0
(2k + 1) ! 3! 5!
(6) cosseno de x,
X∞
(−1)k 2k 1 1
cos x = x = 1 − x2 + x4 − . . . , para qualquer x.
k=0
(2k) ! 2! 4!
6.2 A exponencial complexa

Representação polar dos números complexos
Usando um sistema de coordenadas
√ cartesianas num plano, podemos associar ao número
complexo z = x + iy, onde i = −1, o ponto de coordenadas (x, y), conforme a figura
6.2.
z = x + iy
y
2
y
+
√ x2
|z |=
θ
0 x
Figura 6.2
p
Definindo a distância do número complexo z até a origem como |z| = x2 + y 2 e
o ângulo θ = arctan(y/x) (veja a figura 6.2) em analogia com as coordenadas polares,
podemos escrever

x = |z| cos θ
=⇒ z = x + iy = |z| cos θ + i|z| sen θ = |z|(cos θ + i sen θ).
y = |z| sen θ
É conveniente definir a exponencial complexa eiθ através da expansão de Taylor de ex ,

que vimos na seção anterior, colocando x = iθ:
X∞
iθ 1 2 1 3 1 4 1 5 (iθ)n
e ≡ 1 + iθ + (iθ) + (iθ) + (iθ) + (iθ) . . . = .
2! 3! 4! 5! n=0
n!
Note que i0 = 1, i1 = i, i2 = −1, i3 = −i, i4 = 1, i5 = i, i6 = −1, i7 = −i, i8 = 1, etc.
O padrão 1, i, -1, −i se repete sempre nesta ordem. As potências pares de i produzem 1 e
-1 alternadamente, enquanto que as potências ı́mpares produzem i e −i alternadamente.
Separando os termos sem i dos termos com i obtemos

iθ 1 2 1 4 1 3 1 5
e = 1 − (θ) + (θ) − . . . + i θ − (θ) + (θ) − . . . .
2! 4! 3! 5!
Os termos sem i representam a expansão do cosseno enquanto que os termos multipli-
cados por i representam a expansão do seno. Portanto,
eiθ = cos θ + i sen θ . (6.3)

Esta última equação permite demonstrar que eiθ tem as propriedades da exponencial
usual. Por exemplo,
eiθ1 eiθ2 = (cos θ1 + i sen θ1 )(cos θ2 + i sen θ2 )
= cos θ1 cosθ2 − sen θ1 sen θ2 + i(cos θ1 sen θ2 + sen θ1 cos θ2 )
= cos(θ1 + θ2 ) + i sen (θ1 + θ2 ) = ei(θ1 +θ2 ) .
Finalmente, podemos escrever um número complexo na forma polar:
z = x + iy = |z| cos θ + i|z| sen θ = |z|eiθ . (6.4)

Como veremos no capı́tulo 6, a exponencial complexa será extremamente útil na
solução de equações diferenciais ordinárias lineares de segunda ordem com coeficiente
constantes.
Exercı́cios
(1) Expanda até ordem x5 a função

1 1+x
f (x) = ln .
2 1−x
Sugestão: lembre que ln(a/b) = ln(a) − ln(b)
(2) Escreva as expansões de senh (x) e de cosh(x) em torno de x = 0.

Sugestão: lembre que d senh (x)/dx = cosh(x), d cosh(x)/dx = senh (x), senh (0) = 0,
e cosh(0) = 1. Compare seus resultados com as expansões para sen (x) e cos(x), dadas
neste capı́tulo.
(3) Um dipolo elétrico é um sistema com carga total zero constituı́do por uma carga
positiva +q e uma carga negativa −q, separadas por uma distância 2a, conforme a figura
6.3. O campo elétrico das duas cargas num ponto P do eixo x é dado por
y
−q +q P
a a x
Figura 6.3

~ = q 1 1
E − ~ı,
4πǫ0 (x − a)2 (x + a)2
onde x é a coordenada do ponto P . Mostre que para x >> a o campo elétrico é propor-
cional a 1/x3 .
Sugestão: escreva
1 1 1 1 1 a
2
= 2 2
= 2 (1 ± a/x)±2 = 2 (1 ± δ)−2 onde δ ≡
(x ± a) x (1 ± a/x) x x x
e em seguida expanda as funções (1 ± δ)−2 até primeira ordem em δ.
(4) Um quadrupolo elétrico é um sistema com carga total zero e momento de dipolo
elétrico também zero. Uma realização possı́vel é mostrada na figura 6.4. O campo elétrico
y
+q −2q +q P
a a x
Figura 6.4
das três cargas num ponto P do eixo x é dado por

~ q 1 2 1
E= − + ~ı,
4πǫ0 (x − a)2 x2 (x + a)2
onde x é a coordenada do ponto P . Mostre que para x >> a o campo elétrico é propor-
cional a 1/x4 .
Sugestão: faça como no exercı́cio anterior do dipolo, mas desta vez expanda as funções
(1 ± δ)−2 , onde δ ≡ a/x, até segunda ordem em δ.
(5) Um capacitor e um resistor ligados em série são descritos pela equação diferencial
dQ Q
R =− ,
dt C
onde Q é a carga na placa positiva do capacitor e −dQ/dt é a corrente no circuito.
Suponha que em t = 0 a carga no capacitor seja Q(0) = Q0 .
Sugestão: suponha uma solução do tipo q(t) = exp (pt), onde p é uma constante. Após
determinar p, use o fato da equação ser homogênea para escrever a solução geral como
Q(t) = Aq(t) = A exp (pt). Finalmente, determine a constante A usando Q(0) = Q0 .
Observe que como a equação diferencial é de primeira ordem no tempo podemos resolver o
problema diretamente no campo dos reais, não sendo necessário usar números complexos.
Capı́tulo 7
Equações diferenciais simples
7.1 Equações diferenciais de primeira ordem

Equações que envolvem derivadas de funções, denominadas equações diferenciais, estão
presentes em todos os ramos da fı́sica. Nós já vimos algumas equações diferenciais muito
simples. Por exemplo, a equação diferencial de um movimento retilı́neo uniforme (MRU),
com velocidade v constante, é dada por
dx
= v. (7.1)
dt
Se v é derivada de x o “teorema fundamental do cálculo” (TFC) garante que x é a integral
de v, Z
x= v dt =⇒ x = vt + c.
Obtemos assim a equação horária do MRU, x = x(t), a menos de uma constante c.

Essa equação diferencial é ordinária (não envolve derivadas parciais), linear (não envolve
potências, do tipo (dx/dt)2 ou x2 , ou de ordem superior) e de primeira ordem (porque
só envolve uma derivada primeira). Equações diferenciais de primeira ordem sempre
produzem uma solução que depende de uma “constante de integração” (mais adiante
vamos ver que as soluções das equações de segunda ordem dependem de duas constantes
de integração). Para determinar a constante c, podemos dar, por exemplo, uma condição
inicial, x(t = 0) = x0 . Com essa condição inicial temos a equação horária usual do MRU,
x = x0 + vt.
O movimento retilı́neo uniformememente variado (MRUV), com aceleração constante,

proporciona mais alguns exemplos de equações diferenciais muito simples. Por exemplo,
no MRUV temos
dv
= a,
dt
em que a aceleração a é constante. Usando o TFC obtemos
Z
v (t) = a dt = at + c1 = at + v0 ,
63
CAPÍTULO 7. EQUAÇÕES DIFERENCIAIS SIMPLES 64
onde a constante c1 é igual à velocidade inicial v0 . Para obter a equação horária x = x (t),
usamos novamente o TFC,
Z
1
x (t) = (at + v0 ) dt = at2 + v0 t + c2 ,
2
onde a constante c2 deve ser escolhida como a posição inicial, c2 = x0 , resultando na
equação horária usual x = x (t) do MRUV. É claro que poderı́amos ter escrito uma
equação diferencial de segunda ordem,
d2 x
= a, (7.2)
dt2
e integrado duas vezes em seguida1 para obter a forma geral da equação horária,
1
x (t) = at2 + c1 t + c2 ,
2
em que as constantes de integração c1 e c2 são identificadas como a velocidade e a posição
inicial, respectivamente, c1 = v0 e c2 = x0 .
Infelizmente, bem poucas equações diferenciais podem ser integradas diretamente como
as equações (7.1) e (7.2) usando o TFC. Além disto, não existe uma método capaz de
resolver todos os tipos de equações diferenciais. Neste capı́tulo vamos apresentar o método
de separação de variáveis que pode ser aplicado a uma classe de equações diferenciais que
aparece com frequência na fı́sica.
Vamos primeiramente mostrar como a separação de variáveis funciona com as equações
diferenciais do MRU e do MRUV que acabamos de resolver com o TFC. Depois enunciare-
mos mais precisamente quais equações diferenciais podem ser resolvidas com esta técnica
e finalmente mostraremos porque ela funciona.
Começamos colocando toda a dependência em x de um lado da equação diferencial e
toda a dependência em t do outro lado.
dx
MRU: = v −→ dx = v dt,
dt
dx
MRUV: = at + v0 −→ dx = (at + v0 ) dt.
dt
Depois de separarmos as variáveis, integramos os dois membros da equação
Z Z
MRU: dx = v dt =⇒ x = vt + c,
Z Z
at2
MRUV: dx = (at + v0 ) dt =⇒ x = + v0 t + c2 .
2
Os resultados coincidem com os obtidos com o TFC. Evidentemente, este procedimento
precisa ser justificado. Qual é o sentido de tratar o sı́mbolo da derivada como uma
fração, deixar o numerador desta fração num lado da equação, o denominador no outro
e em seguida integrar cada lado da equação em relação a variáveis diferentes? Antes de
1
Para calcular a primeira integral colocamos d2 x/dt2 = dv/dt.
mostrarmos que isto pode ser feito, vamos ver para que tipos de equações diferenciais este
método funciona.
O método de separação de variáveis pode ser aplicado a equações diferenciais do tipo
dy(x)
= F (x)G(y(x)). (7.3)
dx
Nesta equação y(x) é a função que queremos determinar, F (x) e G(y) são funções ar-
bitrárias de x e y, respectivamente. Neste tipo de equação a dependência em x pode ser
separada da dependência em y. Isto vai permitir integrar a equação diferencial direta-
mente.
Esta técnica é muito simples na prática. Começamos por separar as dependências
em x e em y na equação (7.3), jogando tudo que depende de y para o lado esquerdo da
equação e tudo que depende de x para o lado direito:
dy
= F (x)dx. (7.4)
G(y)
Em seguida, integramos o lado esquerdo em y e o lado direito em x, obtendo
Z y(x) Z x
dy ′
= F (x′ )dx′ , y0 ≡ y(x0 ). (7.5)
y0 G(y ′ ) x0
Através desta equação determinamos y em função de x. Note que a condição y(x0 ) = y0

está automaticamente satisfeita (basta colocar x = x0 na equação (7.5) para verificar que
ambos os membros da equação se anulam).
Os limites inferiores na equação (7.5) produzem uma constante aditiva. Isto permite
escrever a equação (7.5) de uma maneira mais informal, como fizemos na discussão do
MRU e do MRUV,
Z Z
dy
= F (x)dx + const. (7.6)
G(y)
A equação (7.6) é mais simples do que a equação (7.5) mas, ao contrário desta, não satisfaz
a condição y(x0 ) = y0 automaticamente. É necessário usar a constante de integração para
impor esta condição.
Demonstração:
Vamos começar reescrevendo a equação (4.8) de uma maneira ligeiramente diferente,
Z y(x) Z x
′ ′ dy(x′ ) ′
f (y ) dy = f (y(x′ )) dx , (7.7)
y(x0 ) x0 dx′
onde colocamos b = x, a = x0 e mudamos os nomes das variáveis de integração (x → x′ e

y → y′).
Usando a equação (7.7) com f = 1/G podemos escrever
Z y(x) Z x
1 1 dy(x′ ) ′
dy ′ = dx . (7.8)
y(x0 ) G(y ′ ) x0 G(y(x′ )) dx′
A equação diferencial (7.3) fornece
1 dy(x′ )
= F (x′ ). (7.9)
G(y(x′ )) dx′
Substituindo a equação (7.9) na equação (7.8) obtemos

Z y(x) Z x
1
dy ′ = F (x′ ) dx′ . (7.10)
y(x0 ) G(y ′ ) x0
Do lado esquerdo da equação (7.10) temos uma função de y, do lado direito uma
função de x. Se conseguirmos resolver y como função de x o problema está resolvido.
Se não conseguirmos, teremos y como uma função implı́cita de x e sempre podemos usar
técnicas numéricas para obter y(x). ❚
Note que é preciso conhecer a função y(x) em algum ponto x0 . Caso contrário, o
problema fica determinado a menos de uma constante.
Exemplos
(a) Decaimento radiativo. O processo de decaimento de um núcleo radiativo (carbono-

14, por exemplo) pode ser representrado pela equação diferencial
dN
= −λN,
dt
em que N = N (t) é o número de átomos de carbono-14 no instante t e a “constante
de decaimento” λ > 0 é intrepretada como o inverso da vida média do núcleo radi-
ativo (a vida média do carbono-14 é de 8033 anos). Para estabelecer essa equação
fizemos a hipótese (muito razoável) de que, em cada instante, a taxa de decaimento,
∆N/∆t, é proporcional ao número de átomos radiativos N (t).
De novo, fica muito simples separar as variáveis e calcular as integrais (indefinidas)
dN
= −λ dt; (7.11)
N
de onde vem
ln N = −λt + c,
em que c é uma constante de integração. Agora é conveniente redefinir a constante
e escrever a solução geral N (t) na forma
N (t) = C exp (−λt) ,
em que a nova constante C ≡ exp(c) deve ser interpretada como o número de átomos
radiativos no instante inicial,
N (t) = N0 exp (−λt) .
É claro que N (t) → 0 para t → ∞, quando não há mais núcleos de carbono-14.
Alternativamente, poderı́amos ter usado integrais definidas na integração da equação

(7.11),
Z N Z t
dN ′ ′ ′ |
N
′ |
t
= − λ dt =⇒ ln(N ) | N0 = −λt |0
N0 N
′
0
N
=⇒ ln = −λ t =⇒ N (t) = N0 e−λt ,
N0
que coincide com o resultado anterior. O uso de integrais definidas ou indefinidas é

uma questão de preferência.
(b) Atrito viscoso. Um corpo de massa m caindo em um meio fluido, com uma
velocidade não muito alta, sofre uma força de atrito viscoso, contrária à velocidade,
F~atrito = −γ ~v ,
em que γ > 0 é a constante de atrito viscoso e ~v é a velocidade. Utilizando como

sistema de referência um eixo vertical apontando para baixo, podemos escrever a
equação diferencial de movimento

dv mg
m = −γ v + mg = −γ v − , (7.12)
dt γ
em que tanto a velocidade v quanto a força peso mg são escritas com sinal positivo
devido à escolha do referencial. A constante mg/γ tem dimensão de velocidade2 .
Definindo
mg
vl ≡ ,
γ
a equação diferencial (7.12) pode ser escrita na forma
dv
m = −γ (v − vl ) ,
dt
Separando as variáveis obtemos
dv γ
= − dt.
v − vl m
Integrando a equação (veja o exemplo (g) na página 41) e colocando v(t = 0) ≡ v0

chegamos a
Zv(t) Zt
dv ′ γ v(t) γ
=− dt′ =⇒ ln (v ′ − vl ) || =− t
′
v − vl m v0 m
v0 0

v(t) − vl γ v(t) − vl
=⇒ ln = − t =⇒ = e−γt/m (7.13)
v0 − vl m v0 − vl
2
Veremos mais adiante que mg/γ é a velocidade limite do corpo. Ela é atingida quando a força viscosa
se torna igual à força peso.
Note que tanto para a velocidade inicial menor do que a velocidade limite (v0 ≤
v(t) < vl ) como para a velocidade inicial maior do que a velocidade limite (v0 ≥
v(t) > vl ) a razão (v(t) − vl )/(v0 − vl ) > 0 e o argumento do ln na equação (7.13) é
positivo, como deveria ser. Finalmente,
v(t) = vl + (v0 − vl ) e−γt/m , (7.14)
onde vemos que para t → ∞ a exponencial se anula e a velocidade do corpo tende
para a velocidade limite vl , aumentando quando v0 < vl , ou diminuindo quando
v0 > vl .
(c) Força de arrasto. Sabe-se que um corpo caindo no ar, com uma velocidade sufi-
ciente alta, sofre uma força de arrasto, dada por
1
F = CρAv 2 ,
2
em que A é a seção reta efetiva do corpo, ρ é a densidade do ar, v é a veloci-
dade de queda e C é uma constante adimensional que depende da forma do objeto
(tipicamente varia no intervalo 0,5—1,0).
Utilizando o mesmo referencial do problema anterior em que o eixo vertical aponta
para baixo, a equação de movimento é dada por

dv 1 2 1 2 2mg
m = − CρAv + mg = − CρA v − , (7.15)
dt 2 2 CρA
onde m é a massa do corpo e g é a aceleração da gravidade.
A grandeza 2mg/(CρA) tem dimensão de velocidade ao quadrado. Definindo a
constante r
2mg
vl ≡
CρA
que corresponde, como veremos mais adiante, à velocidade limite podemos reescrever
a equação (7.15) como
dv 1
m = − CρA v 2 − vl2 .
dt 2
Estamos agora diante de uma equação diferencial ordinária, de primeiro grau, mas
que não é mais linear (devido à presença do termo dependente de v 2 ). Não há
técnicas de solução simples para equações não lineares, que em geral demandam
enorme esforço. Nesse caso, devido ao fato das variáveis poderem ser separadas,
não é muito difı́cil resolvê-la.
Separando as variáveis, obtemos
dv CρA
2
=− dt. (7.16)
v2 − vl 2m
Vamos agora usar a propriedade

1 1 1 1
2 2
= −
v − vl 2vl v − vl v + vl
e introduzir uma constante de tempo

m
τ≡
CρAvl
para reescrever a equação (7.16) como

1 1 1
− dv = − dt
v − vl v + vl τ
Integrando os dois lados, temos

Z v Z t
1 1 ′ 1
− dv = − dt′
v0 v ′ − vl v ′ + vl τ 0
v v t
=⇒ ln(v ′ − vl ) || − ln(v ′ + vl ) || =−
v0 v0 τ

v − vl v + vl t
=⇒ ln − ln =− .
v0 − vl v0 + vl τ
Note que tanto para v, v0 < vl como para v, v0 > vl o argumento do primeiro
logaritmo na equação acima é positivo como deve ser.

vl − v vl + v0 t vl − v vl − v0 −t/τ
ln = − =⇒ = e .
vl + v vl − v0 τ vl + v vl + v0
Após um pouco de álgebra isolamos v,



vl − v0 −t/τ
 1 − vl + v0 e 
v(t) = vl 
 ,
vl − v0 −t/τ 
1+ e
vl + v0
onde é fácil ver que v(t = 0) = v0 e que para t → ∞, v(t) → vl exponencialmente
rápido.
(d) Oscilador harmônico. O método que apresentamos não se aplica diretamente ao
oscilador harmônico, cuja equação de movimento é dada por
dv
m = −kx, (7.17)
dt
onde k é a constante elástica da mola e x mede o quanto o corpo de massa m preso
na mola se afastou da posição de equilı́brio da mola (x = 0). O problema é que
aparecem três variáveis: v, x e t. Porém, usando a conservação da energia é possı́vel
integrar a equação (7.17) em duas etapas.
Primeiramente, vamos usar a equação (7.17) para mostrar que a energia se conserva.
Multiplicando os dois membros da equação (7.17) por v = dx/dt obtemos
dv dx
mv = −kxv = −kx . (7.18)
dt dt
Como vdv/dt = (1/2)d(v 2 )/dt e xdx/dt = (1/2)d(x2 )/dt, podemos reescrever a

equação (7.18) como

m dv 2 k dx2 d mv 2 kx2
=− =⇒ + = 0. (7.19)
2 dt 2 dt dt 2 2
A equação (7.19) mostra que a derivada em relação ao tempo da grandeza mv 2 /2 +

kx2 /2 é nula. Em outras palavras, ela independe do tempo. Esta grandeza é a
energia, e podemos escrever
mv 2 kx2
+ = E = const. (7.20)
2 2
É conveniente reescrever a energia E de outra forma. Quando a mola está distendida
ao máximo, x = xmáx ≡ A, a velocidade da massa é nula e toda a energia está sob
a forma de energia potencial da mola, E = kA2 /2.
A equação (7.20) mostra que existe uma relação entre x e v. Esta relação pode ser
usada para reduzir o número de variáveis na equação (7.17). Temos
r
mv 2 kx2 kA2 k
+ = =⇒ v = ± (A2 − x2 ). (7.21)
2 2 2 m
A equação (7.21) está bem definida uma vez que A2 ≥ x2 . Vamos substituir a
solução v ≥ 0 na equação (7.17) (a solução v ≤ 0 fornece o mesmo resultado e será
deixada como exercı́cio),
r
dv k x dx
m = −m √ = −kx. (7.22)
dt m A2 − x2 dt
r
x dx k 1 dx
=⇒ √ = x =⇒ x √ −ω =0 (7.23)
A2 − x2 dt m A2 − x2 dt
r
k
onde definimos ω0 ≡ . Desprezando a solução trivial x = 0, ficamos com a
m
equação
1 dx
√ − ω0 = 0. (7.24)
A2 − x2 dt
Separando as variáveis e integrando os dois lados da equação obtemos
Z x Z t
dx
√ = ω0 dt. (7.25)
x0 A2 − x2 0
Fazendo a mudança de variável de integração
x = A sen θ =⇒ θ = arcsen (x/A) e dx → A cos θdθ
na equação (7.25) obtem-se

Z arcsen (x/A) Z arcsen (x/A)

cos θ dθ
√ = ω0 t =⇒ dθ = ω0 t
arcsen (x0 /A) 1 − sen 2 θ arcsen (x0 /A)
=⇒ arcsen (x/A) = ω0 t + arcsen (x0 /A) =⇒ x = A sen (ω0 t + φ),
onde definimos arcsen (x0 /A) ≡ φ.
7.2 Equações diferenciais de segunda ordem

Estudaremos apenas as equações diferenciais ordinárias lineares de segunda ordem que
descrevem oscilações amortecidas e forçadas. Nestas equações que são extremamente im-
portantes na fı́sica aparecem no máximo derivadas de segunda ordem. No que vai ser feito
abaixo, só não vamos demonstrar o fato de que a solução geral destas equações depende de
duas constantes arbitrárias. Isto é razoável, basta lembrar que para determinar a solução
da equação de movimento (segunda lei de Newton), que é uma equação diferencial de
segunda ordem no tempo, precisamos fixar a posição e a velocidade iniciais.
Quando há dissipação de energia, a técnica que usamos na seção anterior para resolver
a equação do oscilador harmônico—baseada na conservação de energia— não funciona e
precisamos utilizar outros métodos. Por razões didáticas, vamos introduzir estes novos
métodos usando a equação diferencial do oscilador harmônico sem dissipação de energia
e aproveitaremos para explicar porque esta equação é tão importante na fı́sica.
7.2.1 O oscilador harmônico

Geralmente, associamos à expressão “oscilador harmônico” a imagem de uma massa presa
a uma mola oscilando em torno de uma posição de equilı́brio. Este sistema parece ser
muito particular para merecer maior interesse. Porém, a equação diferencial do oscilador
harmônico possui inúmeras aplicações. Por exemplo, pequenas oscilações de uma partı́cula
em torno de um ponto de equilı́brio estável de um potencial arbitrário podem ser bem
aproximadas por oscilações harmônicas. Há inúmeras situações na fı́sica e na engenharia
em que se torna importante considerar os efeitos de pequenas oscilações.
Considere uma partı́cula de massa m que executa pequenas oscilações em torno de um
ponto de equilı́brio estável, com coordenada y = a, de um potencial U (y). Para pequenas
oscilações, podemos expandir U (y) em torno de y = a. Termos proporcionais a (y − a)n
com n > 2 são muito pequenos. Assim, basta expandir U (y) até ordem (y − a)2 . Usando
as equações (6.1) e (6.2) obtemos
dU (y) | 1 d2 U (y) |
U (y) ≈ U (a) + (y − a) + (y − a)2 . (7.26)
dy |y=a 2 dy 2 |y=a
Como y = a é um ponto de equilı́brio estável a derivada primeira se anula e a derivada

segunda é positiva, d2 U (y)/dy 2 |x=a ≡ k > 0. Portanto,
1
U (y) ≈ U (a) + k(y − a)2 . (7.27)
2
De acordo com a segunda lei de Newton, md 2 y/dt2 = F (y) = −dU (y)/dy. Assim, temos
a equação de movimento
d2 y
m 2 = −k(y − a). (7.28)
dt
Definimos y(t) ≡ x(t) + a. Substituindo y por esta expressão na equação (7.28) e lem-
brando que a derivada da constante a é zero obtemos uma equação para x(t),
d2 x
m = −kx, (7.29)
dt2
que é a forma mais conhecida da equação do oscilador harmônico. Finalmente, é usual
definirmos ω02 = k/m. Podemos então reescrever a equação (7.29) na forma
d2 x
= −ω02 x. (7.30)
dt2
Note que se tivéssemos mantido termos de ordem mais alta na expansão (7.26) a
equação de movimento deixaria de ser linear. Por exemplo, se mantivéssemos o termo de
terceira ordem, proporcional a (y − a)3 , terı́amos na equação de movimento (7.28) um
termo adicional proporcional a (y − a)2 . Equações não lineares são muito difı́ceis de serem
resolvidas. Em várias situações é possı́vel considerar o efeito destes termos de ordem mais
alta como uma perturbação em torno do caso linear. Este tema é tratado na teoria de
perturbações que utiliza como ponto de partida a solução do problema linear.
Além das aplicações em problemas que envolvem pequenas oscilações, vale a pena
enfatizar a notável analogia entre o oscilador harmônico e circuitos elétricos. Do ponto
de vista matemático, um circuito elétrico constituı́do de um capacitor com capacitância
C ligado em série a um indutor com indutância L é equivalente a um sistema onde uma
massa m é presa a uma mola com constante elástica k, conforme a figura.
ch 0
1
+Q0
1
0
0
1
0
1
0
1
k m
0
1
C
−Q0
L < > 0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
R
Figura 7.1: Circuito LC e sistema massa mola
Ao fecharmos a chave ch no circuito, a carga positiva inicial Q0 ≡ Q(t = 0) começa a

passar para a placa negativa diminuindo a carga em ambas as placas (a soma das cargas
nas duas placas é sempre zero por conservação de carga). Se não fosse pelo indutor o
processo terminaria quando as cargas nas placas superior e inferior se anulassem. Devido
ao indutor, cargas positivas continuam sendo retiradas da placa superior do capacitor
mesmo depois dela ficar neutra. Assim, a placa superior começa a ficar negativa e a
inferior positiva. Quando a carga da placa superior for igual a −Q0 e a da placa inferior
Q0 , recaimos na situação inicial mas com os sinais das cargas nas placas trocados. Nesse
instante, a carga positiva, agora na placa inferior, começa a retornar para a placa superior
até voltarmos à situação inicial e o ciclo se repete. Como não há dissipação de energia a
carga positiva fica passando de uma placa para outra indefinidamente. No caso do sistema
massa mola, a situação inicial análoga à do capacitor consiste em deslocar a massa da
posição de equilı́brio e soltá-la. Tirar a massa da posição de equilı́brio equivale a carregar
o capacitor. Soltar a massa equivale a fechar a chave ch do circuito. Como não há atrito
a massa vai oscilar em torno de uma posição de equilı́brio indefinidamente. A analogia
entre o sistema elétrico e mecânico fica ainda mais clara quando escrevemos a equação
diferencial da carga Q(t) no circuito LC,
d2 Q(t) 1
L = − Q(t). (7.31)
dt2 C
Comparando as equações (7.29) e (7.31) fica claro que a indutância L desempenha um
papel análogo ao da massa m e representa a “inércia elétrica” do circuito. Ela mede a
resistência que o circuito oferece a mudanças na corrente, assim como a massa mede a
resistência que o corpo oferece a mudanças na velocidade. O inverso da capacitância C
é o análogo da constante elástica k da mola. Além disto, como as equações são análogas
podemos imediatamente escrever a solução da equação do circuito LC e desenvolver uma
boa intuição sobre seu funcionamento. De fato, dividindo ambos os membros da equação
(7.31) por L obtemos
d2 Q(t) 1
2
=− Q(t) ≡ −ω02 Q(t). (7.32)
dt LC
Fazendo Q(t) → x(t) reobtemos a equação (7.30). Portanto, podemos escrever ime-
diatamente a solução geral da equação do circuito LC:
1
Q(t) = A cos(ω0 t + φ), com ω0 = √ , (7.33)
LC
As constantes A e φ são fixadas através das condições iniciais Q(t = 0) ≡ Q0 e
I(t = 0) = −dQ(t)/dt|t=0 ≡ I0 , onde I(t) é a corrente através do circuito. Estas expressões
são análogas às condições iniciais do sistema massa mola: x(t = 0) ≡ x0 e v(t = 0) =
dx(t)/dt|t=0 ≡ v0 .
Este exemplo mostra um aspecto importante das equações diferenciais na fı́sica. A
comparação das equações diferenciais do circuito LC e do sistema massa mola nos permite
identificar o papel desempenhado por cada elemento do circuito e, em particular, entender
o significado fı́sico da indutância.
7.2.2 Solução de equações diferenciais com a exponencial com-

plexa
Vamos explicar esta técnica analisando primeiramente a equação do oscilador harmônico
sem dissipação. Como já vimos este sistema é descrito pela equação diferencial
d2 x
2
= −ω02 x, (7.34)
dt
onde ω02 = k/m.
Observe que esta equação é linear, pois todos os termos dependem de x e não de
x , x3 ou potências mais altas de x; também não aparece nenhuma função de x. Além
2
disto, a equação é homogênea. Isto significa que todos os termos são proporcionais a
uma derivada de x ou ao próprio x e nenhum termo contém apenas uma função de t. É
imediato verificar que, graças à homogeneidade, se x(t) é uma solução da (7.34), então
c x(t), onde c é uma constante, também é solução.
Vamos agora supor que a solução da equação (7.34) seja da forma
x(t) = ept . (7.35)

√
p2 ept = −ω02 ept =⇒ p2 = −ω02 =⇒ p = ± −1 ω0 = ±i ω0 . (7.36)
Supor que a solução da equação (7.34) é da forma x(t) = ept permitiu transformar a
equação diferencial do oscilador harmônico em uma equação algébrica (p2 = −ω02 ). Esta
técnica funciona com qualquer equação diferencial ordinária (sem derivadas parciais) linear
e homogênea.
Encontramos soluções do tipo x(t) = ept = e±iω0 t . Apesar de termos visto no capı́tulo
anterior o significado matemático de uma exponencial complexa, seu significado fı́sico
ainda precisa ser analisado. Como chegamos a uma solução complexa tendo partido da
equação diferencial (7.34) que só apresenta quantidades reais? É possı́vel extrair um
significado fı́sico desta solução complexa?
Ao transformarmos uma equação diferencial em uma equação algébrica podemos ob-
ter soluções complexas, bastando lembrar que as raı́zes da equação do segundo grau
ax2 + bx + c = 0 com b2 − 4ac < 0 são complexas. Para que a técnica de transfor-
mar equações diferenciais em equações álgébricas seja consistente precisamos resolver a
equação diferencial logo de saı́da no corpo dos números complexos. Desta forma, podemos
aceitar soluções complexas. Assim, substituimos a equação diferencial (7.34) por
d2 z
= −ω02 z onde z(t) = x(t) + iy(t). (7.37)
dt2
Devido à linearidade da equação, x — a parte real de z — nunca se mistura com y
— a parte imaginária de z — e ambas satisfazem a mesma equação diferencial. De fato,
substituindo z por x + iy na (7.37) obtemos
2
d2 x 2 dy 2
+ ω0 x + i + ω0 y = 0. (7.38)
dt2 dt2
Para que a equação (7.37) seja satisfeita é preciso que x(t) e y(t) satisfaçam simultane-
amente a equação diferencial que queremos resolver (para que um número complexo seja
zero é necessário que tanto a sua parte real como a sua parte imaginária sejam nulas).
Portanto, se encontrarmos um z(t) que satisfaz a equação diferencial que queremos
resolver então tanto sua parte real x(t), como sua parte imaginária y(t), satisfazem a
mesma equação diferencial e são soluções do problema. Em outras palavras, resolvemos
o problema no corpo dos complexos e no final dos cálculos tomamos a parte real, ou a
parte imaginária, do z encontrado.
Para o oscilador harmônico, partimos da equação (7.37) e supomos que a solução seja
do tipo z = ept . Isto leva, como já mostramos, à equação algébrica p2 = −ω02 que produz
a solução z = e±iω0 t . Tomando a parte real de z obtemos
x(t) = Re(z(t)) = Re(e±iω0 t ) = Re[cos(ω0 t) + i sen (ω0 t)] = cos(ω0 t),

onde usamos a equação (6.3). A solução obtida, cos(ω0 t), satisfaz a equação do oscilador
mas não é a solução mais geral. Para obter a solução geral lembramos que a equação (7.37)
é uma equação homogênea. Isto significa que se z(t) é uma solução então C z(t), onde C
é uma constante, também é solução. Como estamos resolvendo o problema no corpo dos
complexos a constante mais geral é complexa, C = Aeiφ , onde A e φ são constantes reais,
e usamos a forma polar do número complexo C. Finalmente,
x(t) = Re(Cz(t)) = Re(Aeiφ eiω0 t ) = Re(Aei(ω0 t+φ) ) = A cos(ω0 t + φ). (7.39)
Esta técnica também é muito conveniente para resolver as equações do circuito RLC e
do oscilador harmônico com amortecimento viscoso. Vamos resolver a equação do circuito
RLC.
7.2.3 Oscilações amortecidas

Para descrever sistemas elétricos um pouco mais realistas precisamos levar em conta a
dissipação de energia. Esta dissipação ocorre na forma de calor devido à resistência
elétrica R do circuito. Podemos levar este efeito em conta adicionando ao circuito LC
uma resistência R obtendo assim o circuito RLC em série, conforme a figura.
ch 00
11
11
00
00
11
00
11 k m
+Q0 00
11
00
11
00
11
C
−Q0
L < > 00
11
00
11
00
11
00
11
00
11
00
11
00
11
R 00
11 ρ
Figura 7.2: Circuito RLC e sistema massa mola num meio viscoso
A equação diferencial do circuito RLC expressa o fato de que a diferença de potencial

no indutor é igual à soma das diferenças de potencial no resistor e no capacitor (você vai
aprender no curso de Fı́sica 3 as convenções para os sinais destas diferenças de potenciais).
A equação é
dI(t) 1 d2 Q(t) dQ(t) 1

−L = RI − Q(t) ⇐⇒ L 2
+R + Q(t) = 0 , (7.40)
dt C dt dt C
onde usamos I = −dQ/dt.
O sistema mecânico equivalente ao circuito RLC é o composto de uma massa e uma
mola imersas num meio viscoso que oferece uma força resistiva Fvisc proporcional à velo-
cidade,
dx
Fvisc = −ρv = −ρ . (7.41)
dt
Levando este termo em conta ao escrever a equação de movimento para a massa m

obtemos
d2 x dx d2 x dx
m 2
= −kx − ρ ⇐⇒ m 2
+ ρ + kx = 0 . (7.42)
dt dt dt dt
Qualitativamente, sabemos o que vai acontecer com o sistema massa mola no meio
viscoso quando deslocamos a massa de sua posição de equilı́brio. Se o parâmetro ρ for
pequeno a massa vai oscilar algumas vezes com amplitude cada vez menor, devido à dis-
sipação de energia, até parar na posição de equilı́brio. Se o parâmetro ρ for grande, a
massa vai diretamente para a sua posição de equilı́brio sem oscilações. Como a equação
do circuito RLC é igual à equação deste sistema, a fı́sica é análoga. Esta análise quali-
tativa do circuito RLC é interessante mas não fornece muitas informações; por exemplo
não conseguimos determinar quando se dá a transição do regime oscilante para o regime
puramente amortecido. Para isto precisamos resolver a equação (7.40), ou equivalente-
mente a equação (7.42). Vamos mostrar como resolver este tipo de equação fazendo uso
da exponencial complexa que introduzimos na seção 7.2.2.
Primeiramente, dividimos todos os membros da equação (7.40) por L,
d2 Q(t) R dQ(t) 1 d2 Q(t) dQ(t)

+ + Q(t) = 0 ⇐⇒ + γ + ω02 Q(t) = 0, (7.43)
dt2 L dt LC dt2 dt
onde
R 1
γ≡ e ω02 ≡ . (7.44)
L LC
Analogamente ao que fizemos para o oscilador harmônico, supomos que a solução da
equação (7.43) tem a forma
Q(t) = ept . (7.45)

Substituindo e equação (7.45) na equação (7.43) obtemos
(p2 + γp + ω02 )ept = 0 =⇒ p2 + γp + ω02 = 0, (7.46)

uma vez que ept 6= 0. As soluções da equação de segundo grau em p são
p r
−γ ± γ 2 − 4ω02 γ γ2
p± = =− ± − ω02 . (7.47)
2 2 4
p
Observe que podemos ter soluções complexas se γ/2 < ω0 ⇔ R < 2 L/C (veja a
equação (7.44)). Há três tipos de solução.
p
(1) Amortecimento subcrı́tico para γ/2 < ω0 ⇔ R < 2 L/C.
Neste caso p± é√complexo e é conveniente fatorar −1 na raiz da equação (7.47).
Lembrando que −1 = i obtemos
r
γ γ2 γ
p± = − ± i ω02 − ≡ − ± iω1 , (7.48)
2 4 2
onde definimos
r
γ2
ω1 ≡ ω02 −, (7.49)
4
como γ = R/L, ω1 = ω0 quando a resistência R = 0 e reobtemos o circuito LC.
Podemos obter a solução geral como fizemos para o oscilador harmônico. Escolhemos
uma das soluções, por exemplo ep+ t , multiplicamos por uma constante complexa
C = Aeiφ e tomamos a parte real da expressão resultante. Então
Q(t) = Re[Ae−iφ e−γt/2+iω1 t ] = Ae−γt/2 cos(ω1 t + φ). (7.50)

Como a solução encontrada depende de duas constantes arbitrárias, A e φ, pode-
mos afirmar que ela é de fato a solução geral para o amortecimento subcrı́tico. Se
tivéssemos escolhido ep− t ao invés de ep+ t terı́amos obtido a mesma solução.
p
(2) Amortecimento supercrı́tico para γ/2 > ω0 ⇔ R > 2 L/C
Agora as duas raı́zes da equação (7.46) são reais:
r
γ γ2 γ
p± = − ± − ω02 ≡ − ± β, (7.51)
2 4 2
onde definimos
r
γ2
β≡ − ω02 . (7.52)
4
A solução geral é obtida tomando uma combinação linear das duas soluções ep+ t e
ep− t . Assim temos
Q(t) = Ae−(γ/2−β)t + Be−(γ/2+β)t . (7.53)

Note que se tivéssemos tomado uma das soluções, por exemplo ep+ t , multiplicado por
uma constante complexa C, como fizemos no amortecimento subcrı́tico, obterı́amos
A cos φ ep+ t que não é a solução geral uma vez que A cos φ ≡ A′ é na verdade uma
única constante arbitrária.
p
(3) Amortecimento crı́tico para γ/2 = ω0 ⇔ R = 2 L/C
Neste caso a equação (7.46) para p tem uma única raiz:
γ
p=− . (7.54)
2
e obtemos a solução Q(t) = Ae−γt/2 . Como esta solução só contém uma constante
arbitrária ela não é a solução geral. Uma maneira de obter outra solução indepen-
dente é usar a seguinte combinação linear das soluções indepentes encontradas para
o amortecimento supercrı́tico:
e−(γ/2−β)t − e−(γ/2+β)t
Q(t) = lim = te−γt/2 , (7.55)
β→0 2β
onde expandimos e±βt até primeira ordem em β para deduzir o resultado (Por que
não é necessário ir além da primeira ordem?). Assim, a solução geral é
Q(t) = (A + Bt)e−γt/2 . (7.56)
Este modo de encontrar a outra solução da equação do amortecimento crı́tico é

muito engenhoso mas pouquı́ssimo intuitivo. Uma maneira mais longa mas aplicável
a outros problemas é apresentada na seção (7.2.5).
A forma das soluções para os amortecimentos subcrı́tico, crı́tico e supercrı́tico é mos-

trado na figura 7.3. As três curvas foram obtidas mudando-se apenas a resistência R
do circuito (equivalente à viscosidade ρ no sistema massa mola). Nos três casos a carga
inicial no capacitor Q(t = 0) antes de fecharmos a chave ch é a mesma.
Q
Supercrı́tico
Supercritico
Crı́tico
Critico
Subcrı́tico
Subcritico
0
t
Figura 7.3: Soluções do circuito RLC com amortecimentos subcrı́tico, crı́tico e supercrı́tico
7.2.4 Oscilações amortecidas e forçadas

Finalmente, vamos analisar o oscilador harmônico amortecido forçado e o circuito elétrico
equivalente. As equações diferenciais destes dois sistemas são exemplos de equações
diferenciais ordinárias de segunda ordem lineares e não homogêneas.
Conforme a figura, no sistema massa mola com amortecimento viscoso aplicamos uma
força harmônica F0 cos ωt com frequência f = ω/2π. No circuito RLC inserimos uma
fonte de tensão que fornece uma voltagem alternada V0 cos ωt (poderia também ser uma
antena recebendo o sinal de uma estação de rádio com frequência f = ω/2π).
A equação diferencial do circuito RLC com uma fonte de voltagem alternada é dada
por
d2 Q(t) dQ(t) 1 d2 Q(t) dQ(t) V0

L 2
+ R + Q(t) = V 0 cos ωt ⇐⇒ 2
+γ + ω02 Q(t) = cos ωt,
dt dt C dt dt L
(7.57)
em que ω0 e γ são definidos da forma usual (veja a equação (7.44)).
V0 cos ωt
0
1
1
0
0
1
0
1
0
1
k m
0
1
0
1
C L < > 0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
F0 cos ωt
R ρ
Figura 7.4: Circuito RLC e sistema massa mola num meio viscoso forçados.
A solução geral da equação (7.57), exatamente como a da equação homogênea, depende

de duas constantes arbitrárias. Suponha que conseguimos achar uma solução particular
Qpart (t) da equação (7.57), isto é
d2 Qpart (t) dQpart (t) 2 V0

+ γ + ω0 Q part (t) = cos ωt. (7.58)
dt2 dt L
Geralmente, como o seu nome indica, a solução particular não depende de nenhuma
constante arbitrária. Porém, se tivermos a solução geral Qhom (t) da equação homogênea—
-como é o nosso caso—, a solução geral da equação não homogênea é simplesmente
Q(t) = Qhom (t) + Qpart (t). (7.59)

Vamos mostrar primeiramente que Q(t) = Qhom (t) + Qpart (t) é solução da equação não
homogênea (7.57). De fato, substituindo a equação (7.59) na equação (7.57) obtemos

d2 Qhom (t) dQhom (t) 2
+γ + ω0 Qhom (t)
dt2 dt
2
d Qpart (t) dQpart (t) 2 V0
+ +γ + ω0 Qpart (t) − cos ωt = 0, (7.60)
dt2 dt L
onde usamos a equação (7.43) com Q(t) = Qhom (t) e a equação (7.58) para mostrar que os
dois termos entre colchetes são nulos. Como Qhom (t) contém duas constantes arbitrárias,
Q(t) = Qhom (t) + Qpart (t) é a solução geral. Assim, tudo que precisamos agora é encontrar
uma solução particular.
Nos três casos de amortecimento: subcrı́tico, crı́tico e supercrı́tico as soluções Qhom (t)
vão exponencialmente a zero à medida que t aumenta. Por esta razão Qhom (t) é chamada
solução transiente. A solução dominante para tempos razoavelmente longos, Qpart (t), é a
solução estacionária.
A maneira mais conveniente de obter Qpart (t) consiste em trabalhar no corpo dos
complexos. Defino z(t) tal que Re(z(t)) = Qpart (t). Lembrando que cosωt = Re(eiωt ),
vemos que a equação (7.57) é a parte real da equação
d2 d V0
2
z(t) + γ z(t) + ω02 z(t) = eiωt . (7.61)
dt dt L
Vamos usar o sistema massa mola para guiar nossa intuição. Neste sistema, depois que
a solução transiente yhom (t) se torna desprezı́vel, a massa oscila com uma frequência igual
à da força aplicada. Portanto, a solução estacionária ypart (t) deve ser do tipo A cos(ωt+ϕ).
Como veremos abaixo, A e ϕ são constantes reais determinadas pela equação diferencial,
não são constantes arbitrárias. As duas constantes arbitrárias estão contidas em yhom (t).
Esta análise sugere fazer a substituição
z(t) → Ceiωt ≡ Aeiϕ eiωt , (7.62)

na equação (7.61). Note que Re(Ceiωt ) = A cos(ωt + ϕ) como queremos. Após a substi-
tuição, obtemos uma equação para a constante complexa C.
V0 V0
−ω 2 + iγω + ω02 Ceiωt = eiωt =⇒ C = 2
. (7.63)
L L(ω0 − ω 2 + iγω)
Para colocar a constante C na forma Aeiϕ colocamos ω02 − ω 2 + iγω na forma polar
(veja a equação (6.4)). Então
q
γω
ω02 2
− ω + iγω = (ω02 − ω 2 )2 + γ 2 ω 2 e , iα
α = arctan 2
. (7.64)
ω0 − ω 2

V0 V0
C= p
2
e−iα ≡ p 2 eiϕ , (7.65)
2 2 2
L (ω0 − ω ) + γ ω 2 L (ω0 − ω 2 )2 + γ 2 ω 2
onde definimos ϕ = −α. Assim,

V0 iωt+iϕ γω
z(t) = p e , ϕ = − arctan 2
. (7.66)
L (ω02 − ω 2 )2 + γ 2 ω 2 ω0 − ω 2
Finalmente, a solução estacionária é

V0 γω
Qpart (t) = Re(z(t)) = p cos(ωt + ϕ), ϕ = − arctan 2
.
L (ω02 − ω 2 )2 + γ 2 ω 2 ω0 − ω 2
(7.67)
O problema dos osciladores amortecidos e forçados é tratado em detalhe nas referências
[1] e [2]. Aqui nos preocupamos principalmente com os aspectos matemáticos da solução
das equações diferenciais destes sistemas.
7.2.5 Operadores diferenciais com coeficientes constantes

Começamos reescrevendo a equação diferencial (7.43) do circuito RLC como
2
d d 2
+ γ + ω0 Q(t) = 0. (7.68)
dt2 dt
b atuando sobre Q(t) como
Definimos o operador diferencial O
b d2 d
O ≡ 2 + γ + ω02 . (7.69)
dt dt
b e equação (7.68) se escreve como

Em termos de O
b Q(t) = 0.
O (7.70)
b pode
No caso do amortecimento crı́tico, ω02 = γ 2 /4 e é fácil verificar que o operador O
ser fatorado como o produto de dois operadores d/dt + γ/2 e a equação (7.68) é reescrita
como

d γ d γ
+ + Q(t) = 0, (7.71)
dt 2 dt 2
onde aplicamos primeiro o operador d/dt + γ/2 mais à direita sobre Q(t) e em seguida o
operador d/dt + γ/2 mais à esquerda sobre o resultado da primeira aplicação. Definimos

d γ
y(t) ≡ + Q(t). (7.72)
dt 2
Em termos de y(t) a equação (7.71) se escreve como

d γ
+ y(t) = 0. (7.73)
dt 2
A solução geral desta equação ordinária de primeira ordem é facilmente encontrada usando
os métodos desenvolvidos na seção 7.1:
y(t) = Be−γt/2 , (7.74)

onde B é uma constante. Inserindo este resultado na equação (7.72) obtemos uma equação
de primeira ordem não homogênea para Q(t).

d γ dQ γ γt/2
+ Q(t) = Be−γt/2 ⇐⇒ eγt/2 + e Q(t) = B. (7.75)
dt 2 dt 2
Note que o lado esquerdo da equação pode ser reescrito como
d γt/2
e Q(t) = B =⇒ eγt/2 Q(t) = A + Bt =⇒ Q(t) = (A + Bt)e−γt/2 , (7.76)
dt
que coincide com o resultado anterior. Este método é baseado no fato de que operadores
diferenciais com coeficientes constantes podem ser fatorados como polinômios.
Exercı́cios
(1) Resolva a equação diferencial

dy
= 2y + 1,
dx
sabendo que y(0) = 0.
dy
= y 2 sen (t),
dt
sabendo que y(0) = 1.

dv
= 1 − v2,
dt
sabendo que v(0) = 0.
Sugestão: use a propriedade

1 1 1 1
2
= + .
1−v 2 1−v 1+v
(4) A equação que descreve um circuito RC com um capacitor de capacitância C ligado

a um resistor com resistência R é dada por
dQ(t) Q(t)
R =− ,
dt C
onde Q(t) é a carga no capacitor e −dQ/dt é a corrente no circuito.
(a) Resolva a equação diferencial acima sabendo que no instante t = 0 a carga no

capacitor é Q0 .
(b) Mostre que toda a energia U0 = Q20 /2C armazenada no capacitor do problema
(6) em t = 0 é dissipada como calor no resistor.
Sugestão: a potência dissipada pelo resistor como calor é RI 2 = R(dQ/dt)2 .
Portanto, a energia total dissipada no resistor é
Z∞
Edis = RI 2 dt.
0
(5) Um circuito LR é constituı́do de um indutor com indutância L ligado em série a

um resistor com resistência R. A corrente I(t) que passa neste circuito satisfaz a
equação diferencial
dI(t)
R = −RI ,
dt
onde I(0) = I0 . Determine I(t).
(6) Uma massa m está presa a uma mola de constante elástica k. O sistema está imerso
num meio viscoso que oferece uma força de reação Fres = −ρv ao movimento da
massa. Escreva as soluções possı́veis da equação diferencial do sistema explicitando
em que condições elas ocorrem. Expresse seus resultados em função de m, k e ρ.
(7) Determine a solução estacionária do sistema massa mola descrito no problema (6)
quando
√ uma força F0 cos(ωt) é aplicada à massa. Para que valor de ω, quando
ρ << mk, a amplitude do movimento será máxima?
Apêndice A
Soluções dos exercı́cios
A.1 Limites
(1)
x2 − 5x + 6 (x − 2)(x − 3)
(a) lim = lim = lim (x − 3) = −1
x→2 x−2 x→2 x−2 x→2
(b) Substituição direta de x por 2

(c) Substituição direta de x por 0
y 2 − 25 (y + 5)(y − 5)
(d) lim = lim = lim (y + 5) = 10
y→5 y − 5 y→5 y−5 y→5
(2)
x−5 x−5 1 1
(a) lim = lim = lim =
x→5 x2 − 25 x→5 (x − 5)(x + 5) x→5 x + 5 10
2
x −x−2 (x − 2)(x + 1) x−2
(b) lim 2
= lim = lim ⇒ o limite não existe
x→1 x −1 x→1 (x + 1)(x − 1) x→1 x−1
x+5 x+5 1
(c) lim 2 = lim = lim ⇒ o limite não existe
x→5 x − 25 x→5 (x − 5)(x + 5) x→5 x − 5
5x3 + 8x2 x2 (5x + 8) 5x + 8 1

(e) lim 4 = lim = lim = −
x→0 3x − 16x2 x→0 x2 (3x2 − 16) x→0 3x2 − 16 2
1 1
(d) lim 2 = ∞ (Compare com lim que não existe. Neste caso, se x → 0
x→0 x x→0 x
por valores menores do que zero, 1/x → −∞, mas se x → 0 por valores maiores
do que zero, 1/x → ∞.)
83
APÊNDICE A. SOLUÇÕES DOS EXERCÍCIOS 84
A.2 Derivadas
dy1 (x) d 2
= (x + 5)8 = 8(x2 + 5)7 2x = 16x(x2 + 5)7
dx dx
dy2 (t) d
= cos(ωt + φ) = −ω sen (ωt + φ)
dt dt
dy3 (t) d
= [cos(ωt)]2 = −2ω cos(ωt) sen (ωt)
dt dt
dy4 (x) d
= [ sen (ax2 + bx)] = (2ax + b) cos(ax2 + bx)
dx dx
dy5 (t) d
= exp(−ωt) = −ω exp(−ωt)
dt dt
dy6 (x) d
= exp(ax2 + bx) = (2ax + b) exp(ax2 + bx)
dx dx
dy7 (x) d 2ax + b
= ln(ax2 + bx + c) = 2
dx dx ax + bx + c
dy8 (x) d 1 ax + b/2
= (ax2 + bx + c)1/2 = (2ax + b)(ax2 + bx + c)−1/2 = √
dx dx 2 ax2 + bx + c
dy8 (x) d 1 −ax
= (ax2 + c)−1/2 = − (2ax)(ax2 + c)−3/2 =
dx dx 2 (ax + c)3/2
2
A.3 Integrais
(1)
Z
x5
(2x + x4 + e−x ) dx = x2 + − e−x
5
Z
1
cos(ωt + φ) dt = sen (ωt + φ)
ω
Z
1 1
dx = ln(a + bx)
a + bx b
Z p 1
y 1 + y 2 dy = (1 + y 2 )3/2
3
Z
z
2
dz = ln(z 2 + 1)
z +1
Z
1 1 t
2 2
dt = arctan( ); fizemos a mudança de variável t = a tan θ.
t +a a a
Z
1 1 t−a 1 1 1 1
dt = ln ; usamos 2 = − .
t 2 − a2 2a t+a t − a2 2a t − a t + a
Z
x 1
2 2 3/2
dx = − √
(x + a ) x 2 + a2
Z
dx
ln(x)dx = x ln(x) − x; colocamos ln(x) = ln(x) e integramos por partes.
dx
(2)
O trabalho é dado por

Zx1 Zx1
kx2 |x1 kx21 kx20
W = F (x)dx = − kxdx = − =− +
2 | x0 2 2
x0 x0
(3)
A área entre as curvas é
Z1 Z1
2 x3 x2 |1 = 1 + 1 = 5
A= (y1 (x) − y2 (x))dx = (x + x)dx = + |0
3 2 3 2 6
0 0
(4)
Cálculo da velocidade
Z Z
ct2
v(t) = a(t)dt = ctdt = +k
2
A constante k é determinada com a condição inicial v(0) = v0
ct2
=⇒ k = v0 e v(t) = v0 + .
2
Alternativamente, podemos impor diretamente a condição inicial usando uma integral

definida:
Zt Zt
ct′ 2 |t ct2
′ ′
v(t) − v0 = a(t )dt = ct′ dt′ = = .
2 |0 2
0 0
Cálculo da posição
Z Z
ct2 ct3
x(t) = v(t)dt = v0 + dt = v0 t + + k̄.
2 6
Calculamos a constante k̄ impondo x(0) = x0
ct3
=⇒ k̄ = x0 e x(t) = x0 + v0 t + .
6
Como fizemos para a velocidade, podemos impor a condição inicial diretamente

Zt Z
′ ′ ct′ 2 ct′ 3 |t ct3
x(t) − x0 = v(t )dt = v0 + dt′ = v0 t′ + = v 0 t + .
2 6 |0 6
0
(5)
O trabalho aplicado sobre a partı́cula é

Z∞ Z∞
GM m GM m |∞ GM m
W =− F (r)dr = − 2
dr = | =
r r R T RT
RT RT
(6)
O impulso I de F (t) é
Zt2 Zt2 v(t

Z 2)
dv(t) v(t2 )
I= F (t)dt = m dt = mdv = mv || = mv(t2 ) − mv(t1 ).
dt v(t1 )
t1 t1 v(t1 )
A.4 Vetores
(1)
y
~u
2
−3 −2 −1 0 1 2 3
x
~v −1
−2
Figura A.1
(2)
1 1
(a) ~a = ~ı + ~
8 2
(b) by = −3
(c) ~a + ~b = 7~ı + 5~ − ~k
q p √
(d) |~b| = b2x + b2y + b2z = 32 + (−3)2 + (−1)2 = 19
(e) ~a · ~b = ax bx + ay by + az bz = 12 − 24 = −12
~a · ~b −12 −3
(f ) cos(θ) = = √ = √ ≈ −0, 3078 =⇒ θ ≈ 107, 9◦
|~a||~b| 4 95 95
(3)
d~y (t) d2 ~y (t)
(a) = −9t2~ı + 2~, = −18t~ı
dt dt2
d~r(t) ~ d2~r(t) ~
(b) = 0, =0
dt dt2
d~u(t) d2~u(t)
(c) = −Aω sen (ωt), 2
= −Aω 2 cos(ωt)
dt dt
(4)
Zt Zt
A 3
(a) ~v (t)dt = At2 dt ~ı = t ~ı
3
0 0
Zt Zt Zt
(b) ~v (t)dt = − ω sen (ωt)dt ~ı + ω cos(ωt)dt ~ = cos(ωt)~ı + sen (ωt) ~
0 0 0
(5)
O módulo da velocidade é constante.
d(~v · ~v ) ~ d~v
|~v | = const. =⇒ ~v · ~v = (const.)2 =⇒ = 0 =⇒ 2 · ~v = ~0 =⇒ ~a · ~v = 0.
dt dt
Se ~a · ~v = 0, ou ~a = 0 e o movimento é retilı́neo uniforme, ou ~a ⊥ ~v e o movimento é
circular uniforme.
A.5 Expansões em séries de potências

(1)

1 1+x x3 x5
f (x) = ln =x+ + + ... − 1 < x < 1.
2 1−x 3 5
(2)
x2 x4 x6
cosh(x) = 1 + + + + ...
2! 4! 6!
x3 x5 x7
senh (x) = x + + + + ...
3! 5! 7!
Comparando com as expansões de cos(x) e de sen (x) vemos que, ao contrário do que
ocorre nas expansões das funções trigonométricas, os sinais dos termos nas expansões das
funções hiperbólicas são todos positivos.
(3) e (4)
a
(1 ± δ)−2 = 1 ∓ 2δ + 3δ 2 ∓ 4δ 3 + 5δ 4 + . . . −1<δ ≡ < 1.
x
Para o dipolo vamos usar as expansões acima até ordem δ.

~ = q 1 1 q −2 −2

E − ~
ı = (1 − δ) − (1 + δ) ~ı
4πǫ0 (x − a)2 (x + a)2 4πǫ0 x2
q ~ ≈ q a ~ı.
≈ 4δ~
ı =⇒ E
4πǫ0 x2 πǫ0 x3
Para o quadrupolo vamos usar as expansões acima até ordem δ 2 .

~ q 1 2 1 q
E= 2
− 2+ 2
~ı = 2
(1 − δ)−2 − 2 + (1 + δ)−2 ~ı
4πǫ0 (x − a) x (x + a) 4πǫ0 x
2
q 2 ~ ≈ 3 q a ~ı.
≈ 6δ ~ı =⇒ E
4πǫ0 x2 2πǫ0 x4
(5)
Substituindo a solução q(t) = exp (pt) na equação do capacitor em série com um
resistor obtemos
dept ept ept 1 1
R =− =⇒ Rpept = − =⇒ Rp = − =⇒ p = −
dt C C C RC
t
=⇒ q(t) = exp(− )
RC
Como a equação é homogênea, a solução geral tem a forma
t
Q(t) = Aq(t) = A exp(− ) mas Q(0) = Q0 ⇒ A = Q0
RC
t
=⇒ Q(t) = Q0 exp(− ).
RC
A.6 Equações Diferenciais

(1)
Após separarmos as variáveis obtemos
dy 1 1 1
= dx =⇒ ln(2y + 1) = x + c′ =⇒ y(x) = Ce2x − ; y(0) = 0 ⇒ C = .
2y + 1 2 2 2
(2)
Após a separação de variáveis
dy 1 1
= sen (t)dt =⇒ = cos(t) + C =⇒ y = ; y(0) = 1 ⇒ C = 0.
y2 y cos(t) + C
(3)
Após a separação de variáveis, colocando a condição inicial v(0) = 0 no limite de

integração obtemos
Zv(t) Zt
dv 1 1 1
= dt =⇒ + dv ′ = dt′
1 − v2 2 1 + v′ 1 − v′
0 0
2t
1 + v(t) e −1 et − e−t
=⇒ ln = 2t =⇒ v(t) = 2t = t = tanh(t).
1 − v(t) e +1 e + e−t
(4)
(a) A equação do capacitor em série com o resistor pode ser integrada diretamente.
Colocando a condição inicial no limite de integração obtemos
Q(t)
Z Zt
dQ(t) Q(t) dQ′ 1 ′
R =− =⇒ =− dt
dt C Q′ RC
Q0 0

Q(t) t
=⇒ ln =− =⇒ Q(t) = Q0 e−t/RC
Q0 RC
(b) A energia dissipada é

Z∞
Edis = RI 2 dt.
0
O item (a) fornece
dQ(t) d Q0 −t/RC
I=− = − Q0 e−t/RC = e
dt dt RC
Substituindo este resultado na equação para Edis obtemos
Z∞
Q20 −2t/RC Q20 RC −2t/RC |∞ Q20
Edis = R e dt = − e |0 = .
R2 C 2 RC 2 2 2C
0
(5)
A equação do indutor em série com o resistor é
dI
L = −RI,
dt
separando as variáveis temos
dI R R R
= − dt =⇒ ln(I) = − t + c′ =⇒ I = C exp(− t) ; I(0) = I0 ⇒ C = I0 .
I L L L
(6)
A equação diferencial do sistema massa mola num meio viscoso é dada na equação
(7.40)
d2 x dx d2 x dx
m + ρ + kx = 0 =⇒ + γ + ω02 x = 0 (A.1)
dt2 dt dt2 dt
onde agora
ρ k
γ= e ω02 = . (A.2)
m m
A equação (A.1) tem a mesma forma da equação (7.40) do circuito RLC. Assim, as
soluções são as mesmas, mas com γ e ω02 dados pela equação (A.2).
√
(1) Amortecimento subcrı́tico para γ/2 < ω0 ⇔ ρ < 2 mk.
r
γ2
x(t) = Ae−γt/2 cos(ω1 t + φ) com ω1 = ω02 − .
4
√
(2) Amortecimento crı́tico para γ/2 > ω0 ⇔ ρ > 2 mk
r
−(γ/2−β)t −(γ/2+β)t γ2
x(t) = Ae + Be com β = − ω02 .
4
√
(3) Amortecimento supercrı́tico para γ/2 = ω0 ⇔ ρ = 2 mk
x(t) = (A + Bt)e−γt/2 .
(7)
A equação do oscilador amortecido e forçado é obtida adicionando-se à equação (7.40)
a força F = F0 cos ωt.
d2 x dx d2 x dx F0
m 2
= −kx − ρ + F 0 cos ωt =⇒ 2
+γ + ω02 x = cos ωt, (A.3)
dt dt dt dt m
onde
ρ k
γ= e ω02 = . (A.4)
m m
Comparando com a equação análoga do circuito RLC com uma fonte de corrente
alternada temos as correspondências x ↔ Q e F0 /m ↔ V0 /L. Portanto,

F0 γω
x(t) = p 2 cos(ωt + ϕ), ϕ = − arctan . (A.5)
m (ω0 − ω 2 )2 + γ 2 ω 2 ω02 − ω 2
A amplitude
F0
A= p . (A.6)
m (ω02 − ω 2 )2 + γ 2 ω 2
vai ser máxima quando o argumento da raiz for mı́nimo. Isto ocorre, no caso ρ <<
√
mk ⇔ γ << ω0 , para r
γ2
ωmáx = ω02 − . (A.7)
2
Como γ << ω0 , ωmáx ≈ ω0 e a amplitude do movimento pode se tornar muito grande.
Este fenômeno é chamado ressonância.
Referências Bibliográficas
[1] Hersch Moysés Nussenzveig, Curso de Fı́sica Básica - Vol. 1 – Mecânica, 4a edição,
2002, Editora Blucher, São Paulo.
[2] Hersch Moysés Nussenzveig, Curso de Fı́sica Básica - Vol. 3 – Eletromagnetismo,1a

edição, 5a reimpressão, 2006, Editora Blucher, São Paulo.
[3] Hamilton L. Guidorizzi, Um Curso de Cálculo - Vol. 1, Editora LTC, São Paulo,
2001.
[4] Paulo Boulos, Introdução ao Cálculo - Vol. 1 – Cálculo Diferencial, 1973 e Introdução
ao Cálculo - Vol. 2 – Cálculo Integral e Séries, 2a edição, 1983, Editora Blucher, São
Paulo.
[5] Michael Spivak, Calculus, Editora W. A. Benjamin, Inc., London, 1973.
[6] Murray R. Spiegel, Manual de Fórmulas e Tabelas Matemáticas, Editora McGraw-

Hill do Brasil, Ltda., 1973.
[7] I. S. Gradshteyn e I. M. Ryzhik, Table of Integrals, Series and Products, Sexta Edição,
Editora Academic Press, 2000.
[8] Wolfram Alpha em http://www.wolframalpha.com/input/?i=integral
91

Introducao Elementar As Tecnicas Do Calculo

Enviado por

Direitos autorais:

Formatos disponíveis

Introducao Elementar As Tecnicas Do Calculo

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Introducao Elementar As Tecnicas Do Calculo

Enviado por

Direitos autorais:

Formatos disponíveis

TERCEIRA EDIÇÃO

INTRODUÇÃO ELEMENTAR ÀS TÉCNICAS DO

Carlos E. I. Carneiro, Carmen P. C. Prado e Silvio R. A. Salinas

INSTITUTO DE FÍSICA DA UNIVERSIDADE DE SÃO PAULO

Carlos E. I. Carneiro, Carmen P. C. Prado e Silvio R. A. Salinas

INTRODUÇÃO ELEMENTAR ÀS TÉCNICAS DO

Este trabalho está licenciado com uma Licença

Capa: Carlos E. I. Carneiro

Carneiro, Carlos E. I..Prado, Carmen P. C..Salinas, Silvio R. A.

Introdução elementar às técnicas do calculo diferencial e integral.

1. CÁLCULO DIFERENCIAL E INTEGRAL (TEXTOS

Sobre os autores: Carlos E. I. Carneiro; Carmen P. C. Prado e Silvio R. A. Salinas são

6 Expansões em séries de potências. 56

7 Equações diferenciais simples 63

A Soluções dos exercı́cios 83

Prefácio da terceira edição

2.1 Limite de uma função

lim f (x) = L ≡ f (x0 ),

ou seja, o limite é igual ao valor da função em x0 .

(a) Com f (x) = 3x3 + 2x + 4, temos

lim f (x) = f (2) = 3 × 23 + 2 × 2 + 4 = 32,

(b) Com f (x) = 2 sen x + 3 cos 3x, temos

lim f (x) = 2 sen 0 + 3 cos 0 = 3,

Quando x → 1 o limite é claramente indefinido. Se fizermos x → 1 por valores maiores

2.2 Definição mais precisa de limite

A função f se aproxima do limite L próximo de x0 (limx→x0 f (x) = L) se

Dando o “limite de tolerância” ε = 0, 1 temos

2x21 = 18 − 0, 1 −→ x1 = 2, 991655... = 3 − 0, 008344...,

|f (x) − 18| < ε = 0, 1 quando 0 < |x − 3| < δ = 0, 008.

(1) Mostre que

(2) Calcule o valor dos seguintes limites:

3.1 Definição de derivada

∆f = f (x0 + ∆x) − f (x0 )

(a) Calcular a derivada de f (x) = 3x2 num certo ponto x0 .

∆f ≡ f (x0 + ∆x) − f (x0 ) = 3(x0 + ∆x)2 − 3x20 = 6x0 ∆x + 3∆x2

(b) Calcular a derivada de f (x) = Axn , onde A é uma constante e n é um número

Esses exemplos nos conduzem a uma regra de derivação importantı́ssima, que é

(c) Calcular a derivada de f (x) = cos x no ponto x0 .

Também é fácil obter uma regra para a função seno,

Tabela 3.1: Algumas derivadas fundamentais

O cosseno hiperbólico e o seno hiperbólico são dados por

3.2 Propriedades mais comuns das derivadas

(3) Dado o produto de funções, f (x) = f1 (x)f2 (x), temos

∆f = f1 (x0 + ∆x)f2 (x0 + ∆x) − f1 (x0 )f2 (x0 + ∆x) +

Em termos gerais, podemos escrever

(4) Dado o quociente de duas funções, f (x) = f1 (x)/f2 (x), temos

Faça um esforço para demonstrar essa propriedade.

∆f = f [y(x0 + ∆x)] − f [y(x0 )] .

Somando e subtraindo y(x0 ) também temos

∆f = f [y(x0 + ∆x) − y(x0 ) + y(x0 )] − f [y(x0 )] .

∆f f [y0 + ∆y] − f [y0 ] f [y0 + ∆y] − f [y0 ] ∆y

Na tabela 3.2, onde a, b e c são constantes, apresentamos um resumo destas proprie-

Tabela 3.2: Propriedades importantes das derivadas

(a) f (x) = Ax4 + Bx2 + Cx + D.

O cálculo da derivada é imediato,

(b) f (x) = (Ax4 + Bx2 + Cx + D)3 .

Basta fazer f (y) = y 3 , y = Ax4 + Bx2 + Cx + D e aplicar a “regra da cadeia”:

(c) f (t) = exp (at2 + b).

Basta fazer f (y) = exp (y) com y = at2 + b. Então