Física Matem Teoria e Aplicações PDF

PR REITORIA
DE GRADUAO
Programa de Apoio Produo de Material Didtico
Edson Sardella
FSICA-MATEMTICA:
TEORIA E APLICAES
So Paulo
2008
Pr-Reitoria de Graduao, Universidade Estadual Paulista, 2008.
1 reimpresso 2009
Sardella, Edson
S244f Fsica-Matemtica : teoria e aplicaes / Edson Sardella
So Paulo : Cultura Acadmica : Universidade Estadual Pau-
lista, Pr-Reitoria de Graduao, 2008.
212 p.
Inclui Bibliografia e ndice Remissivo
ISBN 978-85-98605-31-9
1. Fsica-Matemtica. I. Ttulo.
CDD 530.15
Ficha catalogrfica elaborada pela Coordenadoria Geral de Bibliotecas da Unesp

Universidade Estadual Paulista
Reitor
Marcos Macari
Vice-Reitor
Herman Jacobus Cornelis Voorwald
Chefe de Gabinete
Klber Toms Resende
Pr-Reitora de Graduao
Sheila Zambello de Pinho
Pr-Reitora de Ps-Graduao
Marilza Vieira Cunha Rudge
Pr-Reitor de Pesquisa
Jos Arana Varela
Pr-Reitoria de Extenso Universitria

Pr-Reitora
Maria Amlia Mximo de Arajo
Pr-Reitoria de Administrao
Pr-Reitor
Julio Cezar Durigan
Secretaria Geral
Secretria Geral
Maria Dalva Silva Pagotto
Cultura Acadmica Editora

Praa da S, 108 - Centro
CEP: 01001-900 - So Paulo-SP
Telefone: (11) 3242-7171
APOIO:
FUNDAO EDITORA DA UNESP

FUNDAO PARA O VESTIBULAR DA UNESP - VUNESP
CGB - COORDENADORIA GERAL DE BIBLIOTECAS
CONSELHO EDITORAL

Heraldo Luiz Giacheti
Joo Aristeu da Rosa

Miguel Ruiz

!
COMISSO EXECUTIVA
"# !

$$
%& # '
Klaus Schlnzen Junior
Leonor Maria Tanuri
APOIO TCNICO
$* $
+
%& ; <
PROJETO GRFICO
PROGRAMA DE APOIO
PRODUO DE MATERIAL DIDTICO
Considerando a importncia da produo de material didtico-

pedaggico dedicado ao ensino de graduao e de ps-graduao,
a Reitoria da UNESP, por meio da PrReitoria de Graduao
(PROGRAD) e em parceria com a Fundao Editora UNESP (FEU),
mantm o Programa de Apoio Produo de Material Didtico de
Docentes da UNESP, que contempla textos de apoio s aulas, material
audiovisual, homepages, softwares, material artstico e outras mdias,
sob o selo CULTURA ACADMICA da Editora da UNESP, dis-
ponibilizando aos alunos material didtico de qualidade com baixo
custo e editado sob demanda.
Assim, com satisfao que colocamos disposio da comuni-
dade acadmica esta obra, Fsica-Matemtica: teoria e aplicaes,
de autoria do Prof. Dr. Edson Sardella, da Faculdade de Cincias do
Campus de Bauru, esperando que ela traga contribuio no apenas
para estudantes da UNESP, mas para todos aqueles interessados no
assunto abordado.
Sumrio
Prefcio 12
Abreviaes 13
1 Frmulas e Funes Bsicas 14
1.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2 Funes Hiperblicas . . . . . . . . . . . . . . . . . . . . . . . 14
1.3 Relaes de Euler . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.4 Relaes entre Funes Trigonomtricas e Hiperblicas . . . . 15
1.5 Propriedades Bsicas da Funo ln x . . . . . . . . . . . . . . 15
1.6 Fatorial e Duplo Fatorial . . . . . . . . . . . . . . . . . . . . . 16
1.7 Regras de Derivao . . . . . . . . . . . . . . . . . . . . . . . 16
1.8 Regra de L'Hpital . . . . . . . . . . . . . . . . . . . . . . . . 16
1.9 Integrao por Partes . . . . . . . . . . . . . . . . . . . . . . 17
1.10 Notao de Vetores . . . . . . . . . . . . . . . . . . . . . . . . 17
2 Sries Innitas 18
2.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.2 Srie de Taylor . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.3 Srie de Maclaurin . . . . . . . . . . . . . . . . . . . . . . . . 20
2.4 Resto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.5 Srie de Taylor Revisitada . . . . . . . . . . . . . . . . . . . . 21
2.6 Critrio de Convergncia de d'Alembert . . . . . . . . . . . . 25
2.7 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3 Anlise Vetorial 30
3.1 Vetores e Escalares . . . . . . . . . . . . . . . . . . . . . . . . 30
3.2 lgebra dos Vetores . . . . . . . . . . . . . . . . . . . . . . . 31
3.3 Vetores Unitrios . . . . . . . . . . . . . . . . . . . . . . . . . 32
3.4 Sistema de Coordenadas Retangulares . . . . . . . . . . . . . 32
8 Sumrio
3.5 Produto Escalar . . . . . . . . . . . . . . . . . . . . . . . . . 34

3.6 Produto Vetorial . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.7 Lei dos Cossenos . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.8 Produtos Duplos . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.9 Projeo de um Vetor . . . . . . . . . . . . . . . . . . . . . . 37
3.10 Rotao de Coordenadas . . . . . . . . . . . . . . . . . . . . . 38
3.11 Gradiente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.12 Divergente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.13 Rotacional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.14 Exemplos de Campos Vetoriais . . . . . . . . . . . . . . . . . 45
3.15 Frmulas que Envolvem e . . . . . . . . . . . . . . . 46
3.16 Integral de Linha . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.17 Teorema do Divergente . . . . . . . . . . . . . . . . . . . . . . 51
3.18 Teorema de Stokes . . . . . . . . . . . . . . . . . . . . . . . . 53
3.19 Equaes de Maxwell . . . . . . . . . . . . . . . . . . . . . . . 53
3.19.1 Lei de Gauss . . . . . . . . . . . . . . . . . . . . . . . 53
3.19.2 Lei de Ampre . . . . . . . . . . . . . . . . . . . . . . 54
3.19.3 Lei de Faraday . . . . . . . . . . . . . . . . . . . . . . 55
3.19.4 Ausncia de Monopolo Magntico . . . . . . . . . . . 55
3.19.5 Lei de Ampre-Maxwell . . . . . . . . . . . . . . . . . 56
3.20 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4 Coordenadas Curvilneas 60
4.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
4.2 Elementos de Comprimento, rea, e Volume . . . . . . . . . . 62
4.3 Gradiente, Divergente, e Rotacional . . . . . . . . . . . . . . 63
4.4 Coordenadas Cilndricas . . . . . . . . . . . . . . . . . . . . . 66
4.5 Coordenadas Esfricas . . . . . . . . . . . . . . . . . . . . . . 67
4.6 Separao de Variveis . . . . . . . . . . . . . . . . . . . . . . 68
4.7 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
5 Equaes Diferenciais Ordinrias 75

5.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
5.2 Soluo Geral da EDH . . . . . . . . . . . . . . . . . . . . . . 76
5.3 Segunda Soluo . . . . . . . . . . . . . . . . . . . . . . . . . 77
5.4 Soluo Geral da EDNH . . . . . . . . . . . . . . . . . . . . . 78
5.5 Mtodo de Frobenius . . . . . . . . . . . . . . . . . . . . . . . 81
5.6 Equao Indicial . . . . . . . . . . . . . . . . . . . . . . . . . 87
5.7 Coecientes Constantes . . . . . . . . . . . . . . . . . . . . . 94
5.7.1 Razes Reais e Distintas . . . . . . . . . . . . . . . . . 94
5.7.2 Razes Reais e Iguais . . . . . . . . . . . . . . . . . . . 95
5.7.3 Razes Complexas . . . . . . . . . . . . . . . . . . . . 95
Sumrio 9
5.8 Solues Numricas . . . . . . . . . . . . . . . . . . . . . . . . 95

5.9 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
6 Variveis Complexas 102

6.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
6.2 Nmeros Complexos . . . . . . . . . . . . . . . . . . . . . . . 102
6.2.1 Operaes com Nmeros Complexos . . . . . . . . . . 102
6.2.2 Representao Grca de um Nmero Complexo . . . 103
6.3 Funes Complexas . . . . . . . . . . . . . . . . . . . . . . . . 105
6.4 Condies de Cauchy-Riemann . . . . . . . . . . . . . . . . . 105
6.5 Teorema Integral de Cauchy . . . . . . . . . . . . . . . . . . . 107
6.6 Frmula Integral de Cauchy . . . . . . . . . . . . . . . . . . . 111
6.7 Srie de Taylor . . . . . . . . . . . . . . . . . . . . . . . . . . 113
6.8 Pontos Singulares e Plos . . . . . . . . . . . . . . . . . . . . 115
6.9 Srie de Laurent . . . . . . . . . . . . . . . . . . . . . . . . . 115
6.10 Teorema dos Resduos . . . . . . . . . . . . . . . . . . . . . . 117
6.11 Aplicao: Integrais Trigonomtricas . . . . . . . . . . . . . . 120
6.12 Aplicao: Integrais Imprprias . . . . . . . . . . . . . . . . . 122
6.12.1 Integrais Imprprias de Funes Racionais . . . . . . . 123
6.12.2 Integrais Imprprias de Funes Racionais e Trigono-
mtricas . . . . . . . . . . . . . . . . . . . . . . . . . . 125
6.13 Aplicao: Integrais Imprprias com Plos . . . . . . . . . . . 126
6.14 Lema de Jordan . . . . . . . . . . . . . . . . . . . . . . . . . . 128
6.15 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
7 Sries de Fourier 133

7.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
7.2 Funes Peridicas . . . . . . . . . . . . . . . . . . . . . . . . 133
7.3 Funes Ortogonais . . . . . . . . . . . . . . . . . . . . . . . . 134
7.4 Srie de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . 136
7.5 Identidade de Parserval . . . . . . . . . . . . . . . . . . . . . 136
7.6 Funes Pares e mpares . . . . . . . . . . . . . . . . . . . . . 137
7.7 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
8 Problemas de Valores de Contorno 143

8.1 Equao de Difuso do Calor . . . . . . . . . . . . . . . . . . 143
8.1.1 Condies de Contorno Homogneas . . . . . . . . . . 143
8.1.2 Condies de Contorno No-Homogneas . . . . . . . 147
8.2 Equao da Onda . . . . . . . . . . . . . . . . . . . . . . . . . 148
8.3 Problemas de Dirichlet e de Neumann . . . . . . . . . . . . . 150
8.3.1 Equao de Laplace: Coordenadas Retangulares . . . 150
8.3.2 Equao de Laplace: Coordenadas Polares . . . . . . . 152
10 Sumrio
8.4 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155
9 Funo Delta de Dirac 157

9.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 157
9.2 Representaes da Funo Delta . . . . . . . . . . . . . . . . 158
9.3 Propriedades da Funo Delta . . . . . . . . . . . . . . . . . . 160
9.4 Funo de Heaviside . . . . . . . . . . . . . . . . . . . . . . . 162
9.5 Funo Delta em duas e trs Dimenses . . . . . . . . . . . . 163
9.6 Relao de Fechamento . . . . . . . . . . . . . . . . . . . . . 165
9.7 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 166
10 Transformadas de Fourier 168

10.1 Teorema de Fourier . . . . . . . . . . . . . . . . . . . . . . . . 168
10.2 Teorema de Convoluo . . . . . . . . . . . . . . . . . . . . . 171
10.3 Identidade de Parserval . . . . . . . . . . . . . . . . . . . . . 171
10.4 Transformada das Derivadas . . . . . . . . . . . . . . . . . . . 172
10.5 Transformada de Fourier Cosseno e Seno . . . . . . . . . . . . 172
10.6 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175
11 Teoria de Sturm-Liouville: Funes Ortogonais 177

11.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 177
11.2 Equao Diferencial Auto-adjunta . . . . . . . . . . . . . . . 178
11.3 Funes Ortogonais . . . . . . . . . . . . . . . . . . . . . . . . 180
11.4 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 181
12 Polinmios de Hermite 184

12.1 Introduo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 184
12.2 Oscilador Harmnico Simples Quntico . . . . . . . . . . . . . 184
12.3 Relaes de Recorrncia . . . . . . . . . . . . . . . . . . . . . 188
12.4 Funo Geradora . . . . . . . . . . . . . . . . . . . . . . . . . 189
12.5 Frmula de Rodrigues . . . . . . . . . . . . . . . . . . . . . . 191
12.6 Ortogonalidade . . . . . . . . . . . . . . . . . . . . . . . . . . 192
12.7 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 193
13 Problemas Adicionais 195

13.1 Sries Innitas . . . . . . . . . . . . . . . . . . . . . . . . . . 195
13.2 Anlise Vetorial . . . . . . . . . . . . . . . . . . . . . . . . . . 197
13.3 Coordenadas Curvilneas . . . . . . . . . . . . . . . . . . . . . 198
13.4 Equaes Diferenciais Ordinrias . . . . . . . . . . . . . . . . 199
13.5 Variveis Complexas . . . . . . . . . . . . . . . . . . . . . . . 201
13.6 Sries de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . 203
13.7 Equaes Diferenciais Parciais . . . . . . . . . . . . . . . . . . 204
Sumrio 11
13.8 Funo Delta de Dirac . . . . . . . . . . . . . . . . . . . . . . 205

13.9 Teoria de Sturm-Liouville . . . . . . . . . . . . . . . . . . . . 205
13.10Polinmios de Hermite . . . . . . . . . . . . . . . . . . . . . . 206
Bibliograa 207
ndice Remissivo 208
Prefcio
Este livro nasceu da experincia do autor em ministrar a disciplina Fsica-
Matemtica para o curso de Licenciatura em Fsica do Departamento de
Fsica, Campus Universitrio da UNESP de Bauru.
So vrios os livros textos sobre este assunto. Assim, perfeitamente
questionvel a confeco deste livro. Entretanto, os livros textos freqente-
mente contm uma enorme quantidade de informaes, o que muitas vezes
inviabiliza a abordagem de todos os sub-tpicos de um determinado tpico
especco numa disciplina de seis crditos num semestre. Foi diante desta
diculdade que o autor deparou-se com a necessidade de selecionar os prin-
cipais sub-tpicos de cada assunto sempre tendo como objetivo principal
proporcionar aos alunos de fsica os elementos essenciais para facilitar o
acompanhamento e compreenso mais efetivos das disciplinas prossiona-
lizantes do Curso de Licenciatura tais como Eletromagnetismo, Estrutura
da Matria, e Mecnica Quntica.
Em outras palavras, este livro foi concebido para ser um guia prtico
em sala de aula. Portanto, para complementar o contedo adquirido nas
aulas, recomenda-se fortemente a consulta de livros clssicos da disciplina,
tais como Fsica-Matemtica de Eugene Butkov, e Mathematical Methods
for Physicists de George Arfken.
Evitou-se tanto quanto possvel, demonstraes de teoremas. A nfase
dada nas aplicaes destes tanto para casos de interesse puramente aca-
dmico como para aqueles casos que aparecem na soluo de problemas de
fsica.
Como se trata de uma primeira tentativa de desenvolver um conjunto de
tpicos para a disciplina de Fsica-Matemtica, qualquer sugesto que vise
sua melhoria e aprimoramento ser bem vinda pelo autor.
Edson Sardella
Bauru, So Paulo, fevereiro de 2008.
Abreviaes
EDH Equao Diferencial Homognea
EDNH Equao Diferencial No-Homognea
Re(z) Parte real de um nmero complexo z
Im(z) Parte imaginria de um nmero complexo z
Arg(z) Argumento de um nmero complexo z
Res(f, zk ) Resduo de f (z) em z = zk
F {f (x)} Transformada de Fourier direta
F 1 {f (x)} Transformada de Fourier inversa
Captulo 1
Frmulas e Funes Bsicas
1.1 Introduo
Ao longo deste livro usaremos algumas funes elementares e importantes
frmulas matemticas. Embora a maioria delas de conhecimento do leitor,
com o objetivo de escrever um livro auto-contido, catalogamos as principais
funes elementares, suas propriedades fundamentais, bem como regras de
derivao e integrao, dentre outras.
1.2 Funes Hiperblicas

As funes hiperblicas so denidas por:
1 x
sinh x = (e ex ) ,
2
1 x
cosh x = (e + ex ) ,
2
sinh x
tanh x = ,
cosh x
coth x =
cosh x
sinh x
. (1.1)
Captulo 1: Frmulas e Funes Bsicas 15
1.3 Relaes de Euler

Existem importantes relaes entre as funes trigonomtricas e a exponen-
cial. Estas so conhecidas como relaes de Euler. Temos
1 i
cos = (e + ei ) ,
2
1 i
sin = (e ei ) , (1.2)
2i

onde i = 1 a unidade imaginria.
1.4 Relaes entre Funes Trigonomtricas e

Hiperblicas
sin(ix) = sinh(x) ,
cos(ix) = cosh(x) ,
tan(ix) = tanh(x) ,
cot(ix) = coth(x) . (1.3)
1.5 Propriedades Bsicas da Funo ln x

Embora de conhecimento do leitor, convm relembrar algumas propriedades
importantes da funo ln x. Estas propriedades so muito usadas na deter-
minao da segunda soluo de equaes diferenciais lineares de segunda
ordem. Temos
ln xa = a ln x ,
ln(xy) = ln x + ln y ,

x
ln = ln x ln y ,
y
eln x = x. (1.4)
Leonhard Euler (1707 - 1783) nasceu em Basilia, Sua, onde seu pai era ministro
religioso e possua alguns conhecimentos matemticos. Mesmo tendo perdido a viso
ainda jovem por causa de uma catarata, continuou suas intensas atividades de pesquisa
em matemtica.
16 Fsica-Matemtica: Teoria e Aplicaes
1.6 Fatorial e Duplo Fatorial

muito freqente o aparecimento de fatorial e duplo fatorial em problemas
de Fsica-Matemtica. Temos
n! = n (n 1) (n 2) . . . 2 1 ,
(2n)!! = (2n) (2n 2) . . . 4 2 ,
(2n + 1)!! = (2n + 1) (2n 1) . . . 3 1 .
A segunda (terceira) das relaes acima representa o produto somente dos

nmeros pares (mpares). Note que as denies acima seguem as seguintes
propriedades
(n + 1)! = (n + 1)n! ,
(2n)!! = 2n n! ,
(2n)!
(2n)!! = ,
(2n 1)!!
(2n + 1)!
(2n + 1)!! = . (1.5)
(2n)!!
1.7 Regras de Derivao

Derivada do produto e do quociente de duas funes:
d(f g) df dg
= g+f ,
dx dx dx
df dg
d f dx g f dx
= . (1.6)
dx g g2
Regra da cadeia:
df [g(x)] df dg
= . (1.7)
dx dg dx
1.8 Regra de L'Hpital

Considere duas funes f (x) e g(x) diferenciveis em a < x < b; seja x0 um
ponto neste intervalo tal que f (x0 ) = g(x0 ) = 0. Ento
Captulo 1: Frmulas e Funes Bsicas 17
f (x) f (x)
lim = lim , (1.8)
xx0 g(x) xx0 g (x)
onde a linha indica a primeira derivada. Se a indeterminao persistir,
necessrio a aplicao da regra mais de uma vez at que seja removida. A
regra de L'Hpital tambm se aplica a indeterminaes do tipo /.
1.9 Integrao por Partes

Seja u e v duas funes reais de uma varivel x. Ento

u dv = uv v du . (1.9)
1.10 Notao de Vetores

Para nalizar, importante chamar a ateno do leitor para a notao de
vetores que usaremos. Denotaremos um vetor por uma letra em negrito
como em ., e seu mdulo pela mesma letra, porm em itlico como em F .
Captulo 2
Sries Innitas
2.1 Introduo
Uma srie innita uma soma de innitos termos. Existem dois tipos de
srie, as sries numricas e as sries de funes. Neste Captulo estudaremos
apenas esta ltima. Primeiramente estudaremos a srie de Taylor. No
Captulo 6 veremos um novo tipo de srie conhecida como srie de Laurent, e
no Captulo 7 a srie de Fourier. O estudo das sries de funes importante
pois possui vrias aplicaes, tais como na soluo de equaes diferenciais
lineares de segunda ordem, sejam estas ordinrias ou parciais.
2.2 Srie de Taylor

Consideremos uma funo f (x). Suponhamos que f (x) e suas derivadas
de todas as ordens f (x), f (x), . . ., f (j) (x) so continuas no intervalo
b x c. Aqui usamos a notao
dj f (x)
f (j) (x) =
dxj
, f (0) (x) = f (x) . (2.1)
Consideremos a integral da j -sima derivada de f (x),

x
(j)
x
df (j1) (x)
f (x) dx = dx
a a dx
x

= f (j1) (x)
a
= f (j1) (x) f (j1) (a) . (2.2)
Captulo 2: Sries Innitas 19
Integrando novamente, encontramos
x x x
f (j) (x) dx dx = f (j1) (x) f (j1) (a) dx
a a a
(j2)
= f (x) f (j2) (a) (x a)f (j1) (a) ,
(2.3)
onde na passagem da primeira para a segunda linha usamos (2.2) recursi-

vamente.
Continuando, obtemos
x x x
f (j) (x) dx dx dx =
a a a
x
= f (j2) (x) f (j2) (a) (x a)f (j1) (a) dx
a
(x a)2 (j1)
= f (j3) (x) f (j3) (a) (x a)f (j2) (a) f (a) .
2
(2.4)
Continuando este processo, obtemos para a n-sima integral
x x x
f (j) (x) dx

dx dx =
a a a
nintegrais
(jn) (jn)
=f (x) f (a) (x a)f (jn+1) (a)
(x a)2 (jn+2) (x a)n1 (j1)
f (a) f (a) . (2.5)
2! (n 1)!
Tomando-se j = n, a equao acima toma a forma
x x x
f (n) (x) dx

dx dx =
a a a
nintegrais
(x a)2
= f (x) f (a) (x a)f (a) f (a)
2!
(x a)n1 (n1)
f (a) . (2.6)
(n 1)!
Resolvendo para f (x), nalmente encontramos

(x a)2
f (x) = f (a) + (x a)f (a) + f (a)
2!
(x a)n1 (n1)
+ + f (a) + Rn , (2.7)
(n 1)!
onde Rn conhecido como o resto da srie e denido por

x x x
Rn = f (n) (x) dx

dx dx . (2.8)
a a a
nintegrais
Quando a funo f (x) for tal que
lim Rn = 0 , (2.9)
n
ento a srie torna-se uma soma innita de termos. Temos
(x a)2
f (x) = f (a) + (x a)f (a) + f (a) +
2!

f (n) (a)
= (x a)n . (2.10)
n=0
n!
Esta srie conhecida como srie de Taylor.
2.3 Srie de Maclaurin

A expanso de uma funo em srie de Taylor em torno da origem conhe-
cida como srie de Maclaurin. Tomando-se a = 0 em (2.10), obtemos
x2
f (x) = f (0) + xf (0) + f (0) +
2!

f (n) (0) n
= x . (2.11)
n=0
n!
Publicado pelo matemtico ingls Brook Taylor (1685 - 1731) em 1715.

2.4 Resto
O resto (2.8) pode ser escrito numa forma mais prtica usando-se o teorema
do valor mdio do Clculo Integral. Se g(x) for uma funo contnua no
intervalo a x b, existe um ponto x = tal que
b
g(x) dx = (b a)g() . (2.12)
a
Aplicando este teorema para a integral mais interior de (2.8), temos
x x x
Rn = f (n) ()(x a) dx

dx dx
a a a
(n1)integrais
(x a) n
= f (n) () . (2.13)
n!
2.5 Srie de Taylor Revisitada

Nesta Seo mostraremos uma outra forma de obtermos a srie de Taylor.
Considere uma funo f (x) com derivadas contnuas de todas as ordens no
intervalo x h. Seja a seguinte integral
+h
I= f (x) dx = f ( + h) f () . (2.14)

Por meio da mudana de varivel x = + h t, podemos reescrever a

integral acima como
h
I= f ( + h t) dt . (2.15)
0
Integrando por partes, obtemos
h

tf ( + h t) dt
h
I = tf ( + h t)|0 +
0
h
= hf () + tf ( + h t) dt . (2.16)
0
Um nova integrao por partes nos leva seguinte expresso para a

integral

h2 h
t2
I = hf () + f () + f ( + h t) dt . (2.17)
2! 0 2!
Empregando este procedimento repetidas vezes e usando a equao
(2.14) obtemos a frmula geral de Taylor
h2 h3
f ( + h) = f () + hf () + f () + f () + . (2.18)
2! 3!
Podemos escrever a srie numa forma mais familiar usando = 0 e
h = x,
x2 x3
f (x) = f (0) + xf (0) + f (0) + f (0) + , (2.19)
2! 3!
que a srie de Maclaurin.
Exemplo 2.1 Expandir a funo exponencial f (x) = ex em uma srie
de Maclaurin.
Soluo: Primeiramente devemos calcular as derivadas da funo. Te-
mos, f (n) (x) = ex e conseqentemente f (n) (0) = 1. Substituindo as deriva-
das em (2.11), obtemos
x2 x3
ex = 1+x+ + +
2! 3!

xn
= . (2.20)
n=0
n!
No zemos nenhuma meno convergncia da srie. Para isto preci-

samos analisar o resto. Pela equao (2.13), temos que
xn
Rn = e
n!
n
x
ex , (2.21)
n!
pois x. Uma vez que x nito, podemos concluir que
xn
lim Rn = lim ex =0. (2.22)
n n n!
Portanto, a srie converge no intervalo < x < .
Exemplo 2.2 Determinar a expanso em srie de Maclaurin da funo
f (x) = ln(1 + x).
Soluo: Derivando, fcil mostrar que f (n) (x) = (1)n1 (n 1)!(1 +

n (n)
x) , e f (0) = (1)n1 (n 1)!. A substituio das derivadas em (2.11)
produz
x2 x3
ln(1 + x) = x +
2 3

xn
= (1)n1 . (2.23)
n=1
n
Usando (2.13), o resto pode ser escrito como
(n 1)! xn
Rn = (1)n1
(1 + )n n!
xn
, 0x, (2.24)
n
pois Rn |Rn |, e a funo 1/(1 + )n mxima quando = 0. Para
0 x 1, podemos garantir que
xn
lim Rn = lim =0. (2.25)
n n n
Com efeito, a srie converge no intervalo 0 x 1.

Exemplo 2.3 (Expanso Binomial) Expandir a funo f (x) = (1 +
x)m em uma srie de Maclaurin, onde m pertence ao conjunto dos nmeros
reais.
Soluo: Usando regras de derivao, podemos facilmente mostrar que
f (n) (x) = m(m 1)(m 2) (m n + 1)(1 + x)mn . Introduzindo as
derivadas em (2.11), encontramos
m(m 1) 2 m(m 1)(m 2) 3

f (x) = 1 + mx + x + x + . (2.26)
2! 3!
Neste caso, o resto dado por
m(m 1)(m 2) (m n + 1) xn
Rn =
(1 + )nm n!
xn
m(m 1)(m 2) (m n + 1) , 0 x , (2.27)
n!
pois 1/(1 + )nm mxima quando = 0. Supondo 0 x 1,
xn
lim Rn = lim m(m 1)(m 2) (m n + 1) =0. (2.28)
n n n!
Logo, a srie converge no intervalo 0 x 1.
Exemplo 2.4 A energia relativstica de uma partcula dada por
1/2
v2
E = mc2 1 2 . (2.29)
c
Obter a energia cintica da partcula no limite de baixa velocidade, isto ,
v/c 1.
Soluo: Comparando a ltima equao com (2.26), obtemos x =
v 2 /c2 e m = 1/2. Temos
2
1 v2 12 12 1 v2
E = mc2 1+ 2 + 2 +
2 c 2! c

2 1 v2 3 v4
= mc 1+ + + . (2.30)
2 c2 8 c4
O termo mc2 a energia de repouso, e a energia cintica dada por

E mc2 . Temos

1 3 v2
Ecinetica = mv 2 1 + + . (2.31)
2 4 c2
No limite de v/c 0, encontramos a expresso clssica Ecinetica = 12 mv 2 .
Exemplo 2.5 (Dipolo Eltrico) Um dipolo eltrico constitudo de
duas cargas pontuais de magnitude q como ilustrado na Fig. 2.1. Mostrar
que no limite de a/r 1, o potencial eltrico no ponto P dada por
p cos pr
(r, ) =
2
= 3 , (2.32)
r r
onde p = 2aq o momento de dipolo eltrico.
Soluo: Os vetores r1 e r2 esto relacionados com r atravs de r1 =
r ai e r2 = r + ai. Usando a lei dos cossenos podemos escrever
2 1/2
r1 = r + a2 2ar cos ,
2 1/2
r2 = r + a2 + 2ar cos . (2.33)
O potencial eltrico em P , no sistema CGS dado por

r2
r r1
q +q
x
a +a
Figura 2.1: Dipolo eltrico.
q q
(r, ) =
r1 r
2
q 1 1
= .(2.34)
r 1 +
a2 2ar cos 1/2
1+

a2 +2ar cos 1/2
r2 r2
Usando a expanso binomial, vem que
2
q 1 a 2ar cos
(r, ) = 1+ +
r 2 r2
2
1 a + 2ar cos
1+ +
2 r2

q 2a cos
= + . (2.35)
r r
Desprezando termos de ordem igual ou superior a (a/r)2 , obtemos a equao

(2.32).
2.6 Critrio de Convergncia de d'Alembert

Considere uma srie innita na sua forma geral
f (x) = a0 + a1 x + a2 x2 +

= an xn . (2.36)
n=0
O critrio de d'Alembert estabelece que, se o limite

an+1
lim = 1 , (2.37)
n an R
for nito, ento a srie converge no intervalo R < x < R. O valor de R
conhecido como raio de convergncia. O critrio no garante convergncia
nos extremos do intervalo, x = R.
Exemplo 2.6 Mostrar que a srie (2.20) converge no intervalo <
x < .
Soluo: Comparando (2.20) com (2.36), temos que an = 1/n!, e
an+1 /an = 1/(n + 1). Assim,

an+1 1
lim = lim =0, (2.38)
n an n n + 1
de onde se deduz que R = . Conseqentemente, a srie converge no

intervalo < x < .
Exemplo 2.7 Mostrar que a srie (2.23) converge no intervalo 1 <
x < 1.
Soluo: Comparando (2.23) com (2.36), obtemos an = 1/n, e
an+1 /an = n/(n + 1). Usando a regra de L'Hpital, temos que

an+1 n 1
lim = lim = lim =1. (2.39)
n an n n + 1 n 1
Da, obtemos R = 1. De maneira que a srie converge para todo |x| < 1.
Esta srie converge para x = 1, embora o critrio de d'Alembert no capaz
de detectar convergncia neste ponto.
Exemplo 2.8 Mostrar que a srie (2.26) converge no intervalo 1 <
x < 1.
Jean Le R. d'Alembert (1717 - 1783) abandonado quando pequeno nos degraus da
igreja de St. Jean Baptista de Rond, perto de Notre-Dame, em Paris, foi adotado por um
humilde casal. Mais tarde descobriu-se que seu pai era o general da artilharia Chevalier
Destouches e sua me a aristocrtica escritora Madame de Tenoim. Mas d'Alembert,
quando se tornou um matemtico famoso, preferiu ser reconhecido como lho de seus
pais adotivos.
Soluo: Procedendo de forma anloga aos casos anteriores, temos que

neste caso an = m(m 1)(m 2) (m n + 1)/n!. Usando as regras de
fatorial (1.5), tambm podemos escrever an = m!/n!(m n)!. Logo, temos
que

an+1 m! n!(m n)!

=
an (n + 1)!(m n 1)! m!

m! n!(m n)(m n 1)!
=
(n + 1)n!(m n 1)! m!

m n
= . (2.40)
n+1
Tomando o limite,

an+1 m n
lim
= lim =1. (2.41)
n an n n + 1
Ento a srie converge no intervalo 1 < x < 1.
2.7 Problemas
2.1 Mostrar que

x2n+1
sin x = (1)n ,
n=0
(2n + 1)!

x2n
cos x = (1)n .
n=0
(2n)!
Aplicando o critrio de d'Alembert, determinar os raios de convergn-

cia para ambos os casos.
2.2 Usando os resultados do Problema 2.1, e da srie da funo exponen-
cial, mostrar que
eix = cos x + i sin x ,

onde i a unidade imaginria.
2.3 Expandir (12tz +t2 )1/2 em potncias de t. Obter os coecientes de
t0 , t1 , e t2 . O resultado que voc acaba de encontrar so os chamados
polinmios de Legendre, P0 (z), P1 (z), P2 (z).
2.4 Obter os seguintes desenvolvimentos

x2n+1
(a) arcsin x = x + n=1 (2n1)!! 2n+1 ,|x| 1 ,
(2n)!!
cos x x2 4x4 31x6
(b) e = e 1 2! + 4! 6! + , |x| < ,
x2n+1
(c) ln(x + x2 + 1) = x + n=1 (1)n (2n1)!!
(2n)!! 2n+1 , |x| 1 .
2.5 Expandir as funes cosh x e sinh x em srie de Taylor. Determine o

raio de convergncia para cada caso.
2.6 Mostrar que

x Bn n
= x ,
ex 1 n=0 n!
onde B0 = 1, B1 = 1/2, B2 = 1/6, B3 = 0, B4 = 1/30, etc. Os

coecientes da srie Bn so conhecidos como nmeros de Bernoulli.
2.7 A soma relativstica de duas velocidades u e v dada por
u+v
w= ,
1 + uv/c2
onde c a velocidade da luz no vcuo. Mostre que a expanso de w
em potncias de uv/c2 dada por

uv uv 2
w = (u + v) 1 2 + 2 .
c c
Note que no limite uv/c2 1 obtemos a expresso clssica de Galileu

w = u + v.
2.8 O deslocamento de uma partcula de massa de repouso m0 , resultando
de uma constante de fora m0 g ao longo do eixo y

2 1/2
c2 gt
y= 1+ 1 ,
g c
a qual inclui os efeitos relativsticos. Mostrar que a expanso de y em

potncias de gt/c dada por
2
1 2 1 gt
y = gt 1 + .
2 4 c
Note que no limite gt/c 1, y = gt2 /2, que a expresso clssica de

Galileu.
2.9 O perodo de um pndulo simples dado pela seguinte equao
/2 1/2
l M
T =4 1 sin2 sin2 d ,
g 0 2
onde M a amplitude mxima de oscilao. So dadas as seguintes
integrais
/2 /2
3
sin2 d = , sin4 d = .
0 4 0 16
Mostrar que

l 1 M 9 M
T = 2 1 + sin2 + sin4 + .
g 4 2 64 2
2.10 A teoria de radiao do corpo negro de Planck envolve a seguinte

integral:

x3
dx .
0 e 1
x
Mostrar que esta integral igual a 4 /15.

Sugesto: Primeiramente, usar a expanso binomial e mostrar que

(1 ex )1 = enx .
n=0
Segundo, usar os seguintes resultados

1 4
4
= , y 3 ey dy = 6 .
n=1
n 90 0
2.11 O uso da teoria relativstica de Dirac produz a seguinte frmula de

espectroscopia atmica
1/2
2 2
E = mc 1+ ,
(s + n |k|)2
onde s = (|k|2 2 )1/2 , k = 1, 2, 3, . . ., e 2 = Ze2 /c, onde Z
2 4
o nmero atmico. Expandir E em potncias de at ordem .
Captulo 3
Anlise Vetorial
3.1 Vetores e Escalares
Existem grandezas fsicas, tais como temperatura, presso, energia, den-
tre outras, que necessitam da especicao de apenas uma magnitude para
descreve-las. Estas quantidades so chamadas de grandezas escalares. Exis-
tem outras grandezas na fsica, tais como deslocamento, velocidade, campo
eltrico, dentre outras, que alm de uma magnitude, so caracterizadas por
uma direo e um sentido. Tais quantidades so denominadas por grandezas
vetoriais.
Representamos um vetor por um segmento de reta orientado AB dirigido
do ponto A para o ponto B , conforme ilustrado na Fig. 3.1. Para distinguir
um vetor de um escalar, o denotaremos por uma letra em negrito ou por
uma letra com uma seta em cima. Preferencialmente, usaremos a letra em
negrito. Assim, o vetor AB ser representado por A, e o seu mdulo (ou
magnitude) por uma letra em itlico. Assim, A representa o mdulo de A.
O mdulo de um vetor a medida do tamanho do segmento orientado.

Figura 3.1: Representao grca de um vetor.

Captulo 3: Anlise Vetorial 31
3.2 lgebra dos Vetores

Dois vetores so iguais se e somente se possuem o mesmo mdulo, a mesma
direo e o mesmo sentido. Na Fig. 3.1, A = B. Dizemos que dois vetores
so paralelos se tiverem a mesma direo.
O oposto de um vetor A um vetor que possui o mesmo mdulo, a
mesma direo, porm sentido oposto ao de A. Denotamos este por A
(ver Fig. 3.2).
Figura 3.2: Oposto de um vetor; multiplicao de um vetor por um escalar.
A multiplicao de um escalar por um vetor A resulta em um vetor

de mesma direo de A como ilustra a Fig. 3.2. O sentido de A ser o
mesmo de A se for positivo, e oposto de A se for negativo.
Dados A e B, determinamos a soma destes dois vetores reposicionando
os de modo que a origem de um deles coincida com a extremidade do outro
como ilustra a Fig. 3.3. O resultado um vetor C = A + B. Note que
esta regra de soma de vetores equivalente regra do paralelogramo (ver
Fig. 3.3).

Figura 3.3: Soma de vetores.
A subtrao de dois vetores, A B, equivalente a somar A com o

oposto de B (ver Fig. 3.4).
A regra do paralelogramo para adio de vetores pode ser generalizada
para mais de dois vetores. Esta regra conhecida como regra do polgono
(ver Fig. 3.5).
Temos as seguintes propriedades:

Figura 3.4: Subtrao de vetores.

Figura 3.5: Soma de vetores pela regra do polgono.
1. Comutativa da adio:
A+B=B+A, (3.1)
a qual pode ser vericada facilmente usando a regra do paralelogramo.
2. Associativa da adio:
A + (B + C) = (A + B) + C , (3.2)
a qual pode ser vericada facilmente usando a regra do polgono.
3.3 Vetores Unitrios

Denominamos por vetor unitrio aquele que possui mdulo exatamente igual
a um. Se A um vetor de mdulo A, ento a = A/A um vetor unitrio
de mesma direo e sentido de A.
3.4 Sistema de Coordenadas Retangulares

Considere trs vetores unitrios i, j, e k. Cada um deles dene o sentido
positivo dos eixos coordenados (x, y, z) (ver Fig. 3.6).
Os vetores unitrios so mutuamente perpendiculares. Assim, este sis-
tema de coordenadas ortogonal. Todo vetor A pode ser escrito como uma
combinao linear dos vetores unitrios (ver Fig. 3.7). Temos

Figura 3.6: Vetores unitrios i, j, k.
A = Ax i + Ay j + Az k , (3.3)
onde (Ax , Ay , Az ) so as componentes do vetor A. Assim, a soma de dois
vetores A e B ser dada pela soma de suas componentes. Este mtodo
algbrico de somar vetores muito mais prtico do que o mtodo grco.

Figura 3.7: Sistema de coordenadas retangulares.
O mdulo de um dado vetor u denido por

!
u = |u| = u2x + u2y + u2z . (3.4)
Um vetor muito especial que usaremos ao longo deste livro

" o vetor de
posio r = xi + yj + zk, cujo mdulo escreve-se como r = x2 + y 2 + z 2 ,
que a distncia de O ao ponto de coordenadas (x, y, z).
3.5 Produto Escalar

O produto escalar entre dois vetores u e v, denotado por u v, denido
por
u v = uv cos , (3.5)
onde o ngulo formado entre os dois vetores.
Com efeito, os vetores unitrios tm as seguintes propriedades:
ii=jj=kk=1,
ij=ik=jk=0. (3.6)
Note que o produto escalar resulta em um escalar e no em um vetor.

O produto escalar tem as seguintes propriedades:
1. Comutativa:
uv =vu. (3.7)
2. Distributiva:
u (v + w) = u v + u w . (3.8)
Usando-se (3.6) e a propriedade distributiva podemos reescrever o pro-

duto escalar como
u v = ux vx + uy vy + uz vz . (3.9)
2
Note tambm que (a) uu = u ; (b) uv = 0 se u e v so perpendiculares.
3.6 Produto Vetorial

Denota-se por uv o produto vetorial entre dois vetores u e v, cujo mdulo
dado por
|u v| = uv sin , 0, (3.10)
onde o menor ngulo formado entre os dois vetores u e v. Ento, se u

e v so paralelos, o produto vetorial entre eles identicamente nulo.
O vetor w = u v um vetor cuja direo perpendicular ao plano
formado por u e v. O sentido tal que u, v, e w formem um triedro direto.
O produto vetorial segue as seguintes propriedades:
1. Anti-Comutativa:
u v = v u . (3.11)
2. Distributiva:
u (v + w) = u v + u w . (3.12)
Usando-se a denio do produto vetorial, podemos mostrar que os ve-

tores unitrios tm as seguintes propriedades:
ii=jj=kk=0,
i j = k , i k = j , j k = i . (3.13)
Agora, usando estes resultados e a propriedade distributiva do produto

vetorial, podemos mostrar que
uv = (ux i + uy j + uz k) (vx i + vy j + vz k)
= (uy vz uz vy )i + (uz vx ux vz )j + (ux vy uy vx )k .
(3.14)
Esta equao tambm pode ser escrita mais convenientemente na forma

de um determinante

i j k

u v = ux uy uz . (3.15)

vx vy vz
O mdulo do produto vetorial pode ser facilmente obtido usando-se a

denio do produto escalar. Temos
(u v) (u v) = (uy vz uz vy )2 + (ux vz uz vx )2 + (ux vy uy vx )2

= u2y vz2 + u2z vy2 2uy uz vy vz
+u2z vx2 + u2x vz2 2ux uz vx vz

+u2x vy2 + u2y vx2 2ux uy vx vy
+u2x vx2 u2x vx2
+u2y vy2 u2y vy2
+u2z vz2 u2z vz2 , (3.16)
onde foram somados e subtrados seis termos sem alterar o resultado nal.
Os termos positivos podem ser fatorados na forma (u2x + u2y + u2z )(vx2 + vy2 +
vz2 ) = u2 v 2 e os negativos na forma (ux vx + uy vy + uz vz )2 = (u v)2 . Assim,
usando-se (3.5), obtemos
(u v) (u v) = u2 v 2 (u v)2
= u2 v 2 (uv cos )2
= u2 v 2 (1 cos2 )
= u2 v 2 sin2 . (3.17)
Extraindo a raiz quadrada, obtemos o resultado desejado (3.10).
3.7 Lei dos Cossenos

Seja w = u + v. Temos
w2 = (u + v) (u + v)
= uu+uv+vu+vv
= u2 + v 2 + 2u v
= u2 + v 2 + 2uv cos , (3.18)
onde, na ltima linha usamos a denio do produto escalar. Extraindo a

raiz quadrada, obtemos
"
|u + v| = u2 + v 2 + 2uv cos . (3.19)
3.8 Produtos Duplos

Os produtos escalar e vetorial podem se combinados de vrias maneiras,
tais como:
1. Produto misto:

ux uy uz

u (v w) = (w u) v = (u v) w = vx vy vz . (3.20)

wx wy wz
2. Duplo produto vetorial:
A (B C) = B(A C) C(A B) , (3.21)

conhecida como regra BAC CAB (l-se bac menos cab ).
O mdulo do produto misto igual ao volume de um paraleleppedo (ver

Fig. 3.8). A rea do paralelogramo formado pelos vetores B e C dada por
|B C|, ao passo que a altura h dada pela projeo do vetor A ao longo
do vetor unitrio n. Temos que o volume igual a
(A n)(|B C|) = A (|B C|n) = A (B C) . (3.22)

Se os vetores A, B e C no formam um triedro direto, ento A n < 0.
Para garantir que o resultado positivo adotamos o mdulo do produto
misto.

Figura 3.8:
3.9 Projeo de um Vetor

Consideremos um vetor u que forma um ngulo com um segundo vetor s
(ver Fig. 3.9). A projeo de u sobre s denida como sendo o mdulo de
u vezes o cosseno do ngulo entre os dois vetores:
us = u cos . (3.23)
Lembrando que u = ux i + uy j + uz k, a equao acima pode ser usada

recursivamente para obtermos
us = ux cos(i, s) + uy cos(j, s) + uz cos(k, s) , (3.24)

onde, por exemplo, (i, s) o ngulo formado entre (i e s).
z
s
1

u
y
O
Figura 3.9:
3.10 Rotao de Coordenadas

Considere uma rotao do sistema de coordenadas retangulares (x, y) de
um ngulo em torno do eixo z (ver Fig. 3.10). Denotemos por (x , y )
as coordenadas do sistema rotacionado. O nosso objetivo saber como se
relacionam as componentes de um vetor u nos dois sistemas de coordenadas.
Usando (3.24), as coordenadas do vetor u no sistema de coordenadas
(x , y ) sero dadas por
ux = ux cos(i , i) + uy cos(i , j)

= ux cos + uy cos(/2 )
= ux cos + uy sin , (3.25)
e

Figura 3.10:
uy = ux cos(j , i) + uy cos(j , j)

= ux cos(/2 + ) + uy cos
= ux sin + uy cos . (3.26)
Considere a generalizao para o caso tri-dimensional. Temos
ux = ux cos(i , i) + uy cos(i , j) + uz cos(i , k) ,

uy = ux cos(j , i) + uy cos(j , j) + uz cos(j , k) ,
uz = ux cos(k , i) + uy cos(k , j) + uz cos(k , k) . (3.27)
conveniente mudarmos a notao de (ux , uy , uz ) para (u1 , u2 , u3 ) e

denir amn = cos(m, n), onde (m, n) = 1, 2, 3. Assim, as equaes acima
podem ser reescritas como
u1 = a11 u1 + a12 u2 + a13 u3 ,

u2 = a21 u1 + a22 u2 + a23 u3 ,
u3 = a31 u1 + a32 u2 + a33 u3 . (3.28)
Os coecientes amn = cos(m, n) so conhecidos como cossenos diretores.

A equao anterior pode tambm ser escrita na forma matricial

u1 a11 a12 a13 u1
u2 = a21 a22 a23 u2 , (3.29)
u3 a31 a32 a33 u3
ou ainda na forma abreviada usando somatrio

3

um = amn un . (3.30)
n=1
Seja agora o caso particular u = (1, 0, 0) = i. Substituindo em (3.29),
obtemos para as novas componentes de u
u1 = a11 , u2 = a21 , u3 = a31 , (3.31)

de tal maneira que
i = a11 i + a21 j + a31 k . (3.32)

Analogamente, substituindo u por (0, 1, 0) ou (0, 0, 1) , obtemos
j = a12 i + a22 j + a32 k , (3.33)

e
k = a13 i + a23 j + a33 k . (3.34)

Suporemos que o novo sistema de coordenadas (x , y , z ) ortogonal,
isto , os vetores unitrios (i , j , k ) obedecem as mesmas propriedades que
as (3.6) para os vetores (i, j, k). Temos
i i = a211 + a221 + a231 = 1 ,

j j = a212 + a222 + a232 = 1 ,
k k = a213 + a223 + a233 = 1 ,
i j = a11 a12 + a21 a22 + a31 a32 = 0 ,
i k = a11 a13 + a21 a23 + a31 a33 = 0 ,
j k = a12 a13 + a22 a23 + a32 a33 = 0 . (3.35)
Portanto, a soma dos quadrados dos elementos dentro de uma coluna
igual a um, e a soma dos produtos dos elementos de diferentes colunas
igual a zero. Matrizes que obedecem esta propriedade so chamadas de
matrizes ortogonais .
As seis equaes (3.35) podem ser escritas em uma forma compacta como
segue
3

1, j=k
aij aik = jk = , (3.36)
k
0, j=
i=1
onde ij conhecido como delta de Kronecker.
3.11 Gradiente
Por convenincia, mudemos a notao de (x, y, z) para (x1 , x2 , x3 ). Seja
(x1 , x2 , x3 ) uma funo escalar da posio, por exemplo, temperatura,
densidade de energia, etc. O valor de , deve independer do sistema de
coordenadas,
(x1 , x2 , x3 ) = (x1 , x2 , x3 ) . (3.37)

Derivando parcialmente com relao a xi (i = 1, 2, 3), temos
(x1 , x2 , x3 ) (x1 , x2 , x3 )

=
xi xi
x1 x2 x3
= + +
x1 xi x2 xi x3 xi
3
xj
=
j=1
xj xi
3

= aij , (3.38)
j=1
xj
onde aij = xj /xi .

Na forma matricial tambm podemos escrever

x1 a11 a12 a13
x1

a21 a22 a23
. (3.39)
x2 = x2
a31 a32 a33

x3 x3
O vetor coluna que aparece em ambos os lados da equao acima

conhecido como vetor gradiente da funo ,

+j
= i +k , (3.40)
x y z
onde o operador nabla denido por

=i +j +k . (3.41)
x y z
Este operador diferencial s tem efeito quando aplicado sobre uma funo
escalar.
Uma interpretao geomtrica interessante do gradiente que se
(x, y, z) = C a equao de uma superfcie, ento normal a esta em
todos os seus pontos. Para ver isto, considere dois pontos P e Q de uma
superfcie arbitrria (ver Fig 3.11). Seja dr o deslocamento entre os dois
pontos. O deslocamento de P para Q causa uma mudana innitesimal em
= C dada por

d = dx + dy + dz
x y z

= i +j +k (idx + jdy + kdz)
x y z
= dr . (3.42)

Figura 3.11: Interpretao geomtrica do gradiente.
Por outro lado d = 0, logo perpendicular superfcie. O vetor

normal superfcie obtido dividindo-se o gradiente pelo seu mdulo, n =
/||.
Supor agora que P e Q esto localizados em superfcies diferentes, =
C1 e = C2 (ver Fig. 3.12). Considerando os pontos innitesimalmente
prximos, um deslocamento de um ponto ao outro causa uma mudana
dada por
d = C2 C1
=
= dr . (3.43)
Q
= C 2 > C1
dr
P
= C1
Figura 3.12: Outra interpretao geomtrica do gradiente.
Esta mudana mxima quando o vetor paralelo a dr. Portanto,

o gradiente indica o direo de maior crescimento da funo escalar .
Exemplo " 3.1 Uma superfcie esfrica de raio r dada pela equao
(x, y, z) = x2 + y 2 + z 2 = r. Determinar o vetor normal a qualquer
ponto da esfera.
Soluo: Temos as seguintes derivadas
x x
= " = ,
x x2 + y2 + z2 r
y y
= " = ,
y x2 + y2 + z2 r
y z
= " = . (3.44)
z x2 + y2 + z2 r
Substituindo na equao do gradiente (3.40), encontramos

xi + yj + zk r
= = . (3.45)
r r
O vetor acima unitrio e portanto no h necessidade de dividi-lo pelo
seu mdulo.
3.12 Divergente
Existe algumas formas de combinar o operador nabla com funes escalares
e vetoriais. Na Seo anterior vimos que este operador aplicado a uma
funo escalar resulta no gradiente. Podemos tambm combinar com
uma funo vetorial atravs de um produto escalar. O divergente de uma
funo vetorial denido por
ux uy uz
u= + + . (3.46)
x y z
Note que o divergente de um vetor um escalar.
Quando u tem valor positivo em um determinado ponto, dizemos
que o campo vetorial tem uma vazo lquida para fora da vizinhana deste
ponto. Este o caso do campo eltrico produzido por uma carga pontual.
Dependendo do sinal da carga, as linhas de campo ou convergem para, ou
divergem do ponto onde ela se encontra. O inverso ocorre com o campo
magntico, pois neste caso as linhas de campo se fecham em si mesmas,
isto , o campo tem circuitao (ou vorticidade). Neste ltimo caso, no h
vazo: todas as linhas de campo que adentram uma determina superfcie
fechada, iro emergir em algum outro ponto da mesma superfcie.
3.13 Rotacional
Tambm podemos combinar o operador nabla com um vetor atravs de um
produto vetorial. O rotacional denido por

i j k

u = x
y

z . (3.47)
ux uy uz
Ao contrrio do divergente, se o campo tem vorticidade, ento o rotaci-

onal deste campo no nulo.
3.14 Exemplos de Campos Vetoriais
0.8
0.6
0.4
0.2
0
y
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.5 0 0.5 1
x
Figura 3.13: Campo vetorial A = xi + yj.
Um exemplo tpico de campo vetorial com divergente no nulo e rotacio-

nal identicamente nulo A = xi + yj, cuja representao grca encontra-se
ilustrada na Fig. 3.13.
Um segundo exemplo onde ocorre exatamente o oposto do exemplo an-
terior A = yi + xj, cujo grco mostrado na Fig. 3.14.
Pode ocorrer de o campo vetorial no possuir simultaneamente vor-
ticidade e divergncia. Como um caso tpico, temos o campo vetorial
A = yi + xj (ver Fig 3.15).
0.8
0.6
0.4
0.2
0
y
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.5 0 0.5 1
x
Figura 3.14: Campo vetorial A = yi + xj.
3.15 Frmulas que Envolvem e

Supondo que as derivadas parciais de f , g, u, e v existem, ento so vlidas
as seguintes identidades
(f + g) = f + g , (3.48)
(f g) = f g + gf , (3.49)
(f u) = f u + u f , (3.50)
(f u) = f ( u) + (f ) u , (3.51)
(f ) = 0 , (3.52)
( u) = 0 , (3.53)
0.8
0.6
0.4
0.2
0
y
-0.2
-0.4
-0.6
-0.8
-1
-1 -0.5 0 0.5 1
x
Figura 3.15: Campo vetorial A = yi + xj.
(u v) = (v )u + (u )v + v ( u) + u ( v) , (3.54)
(u v) = v ( u) u ( v) , (3.55)
(u v) = u( v) v( u) + (v )u (u )v , (3.56)
2f 2f 2f
(f ) = 2 f = + + . (3.57)
x2 y 2 z 2
Esta ltima conhecida como o Laplaciano de f . Tambm temos
( u) = ( u) 2 u , (3.58)
2 (f g) = f 2 g + g2 f + 2(f ) (g) . (3.59)

Exemplo 3.2 As foras centrais podem ser escritas como F = f (r)r.
Calcular (f (r)r).
Soluo: Usando a identidade (3.51), temos que
(f (r)r) = f r + f r . (3.60)
Usando a denio de rotacional (3.47) fcil mostrar que r = 0.
Com efeito (f (r)r) = f r. Uma vez que f depende apenas do
mdulo de r, temos
f f f
f = i +j +k
x y z
df r df r df r
= i +j +k . (3.61)
dr x dr y dr z
Similarmente ao Exemplo 3.1, obtemos f = (r/r)df /dr. Portanto,

r df
(f (r)r) = r=0. (3.62)
r dr
Exemplos de fora central so o da fora gravitacional e da fora de
Coulomb, em que f (r) = C/r3 , onde C = Gm1 m2 para o primeiro caso e
C = q1 q2 /40 para o segundo caso.
Exemplo 3.3 Calcular r.
Soluo: Pela denio (3.46), temos que
x y z
r= + + =1+1+1=3. (3.63)
x y z

Exemplo 3.4 Calcular 2 g(r), onde g uma funo que depende
somente do mdulo de r.
Soluo: Do Exemplo 3.2, vem que

r dg
(g(r)) = . (3.64)
r dr
Agora, usando a identidade (3.50) e o resultado do Exemplo 3.3, ob-
temos

1 dg 1 dg
(g(r)) = r+r
r dr r dr

3 dg r d 1 dg
= +r
r dr r dr r dr
2
2 dg d g
= + 2 . (3.65)
r dr dr
Se g(r) = rn , ento
(g(r)) = n(n + 1)rn2 . (3.66)

Esta equao ser nula somente se n = 0 ou n = 1. Ou seja, o potencial
de Coulomb g(r) = 1/r satisfaz a equao de Laplace 2 g = 0.
3.16 Integral de Linha

Consideremos (x, y, z) uma funo escalar e u(x, y, z) uma funo vetorial.
Basicamente, existem trs formas de integral de linha. So elas,

dr ,
C
u dr ,
C
u dr ,
C
(3.67)
onde dr um elemento de comprimento de um contorno (aberto ou fechado)

C de integrao.
O mtodo de resoluo de uma integral de linha consiste em reduzi-la a
integrais ordinrias as quais podemos resolver pelos mtodos usuais.
Vejamos alguns exemplos.
Exemplo 3.5 Integrar a funo (x, y) = xy desde a origem at o
ponto (1, 1), seguindo os caminhos C1 e C2 indicados na Fig. 3.16.
Soluo: Atravs do caminho C1 , temos que
1 1
dr = i dx + j dy
C 0,y=0 0,x=1
1 1
= i 0 dx + j y dy
0 0
1
= j . (3.68)
2
Atravs de C2 temos
1 1
dr = i dx + j dy
C 0,y=x 0,y=x
(1, 1)
C2
x
C1 (1, 0)
Figura 3.16: O caminho C1 : y = 0 entre (0, 0) e (1, 0), e x = 1 entre (1, 0)

e (0, 1); o caminho C2 : y = x entre os pontos (0, 0) e (1, 1).
1 1
= i x2 dx + j x2 dx
0 0
1 1
= i +j . (3.69)
3 3

Exemplo 3.6 Integrar a funo F(x, y) = (x y)i + (x + y)j atravs
da parbola que une os pontos (0, 0) e (1, 1) como ilustrado na Fig. 3.17.
(1, 1)
Figura 3.17: O caminho C dado por y = x2 .
Soluo: Ao longo do caminho C , y = x2 . Temos

1 1
F dr = (x y) dx + (x + y) dy
C 0,y=x2 0,y=x2
1 1
= (x x2 ) dx + 2x(x + x2 ) dx
0 0
1
x2 x3 x4
= + +
2 3 2 0
11
= . (3.70)
6

3.17 Teorema do Divergente

Seja V um volume limitado por uma superfcie fechada S . Suponhamos que
as componentes de um vetor u tenham derivadas parciais contnuas em S e
no seu interior. Ento
)
u n dS = u dV , (3.71)
S V
onde n um vetor normal superfcie em qualquer ponto de S . Do lado

esquerdo da equao acima temos uma integral de superfcie e do lado direito
uma integral de volume.
Para demonstrar este teorema, consideremos primeiro um elemento de
volume V . Consideremos duas faces deste volume como ilustrado na
Fig. 3.18. O uxo de u atravs das faces 1 e 2 so dados por
(valor de ux na face 1)yz ,

(valor de ux na face 2)yz . (3.72)
O uxo lquido ento dado por
[(valor de ux na face 2) (valor de ux na face 1)]yz = ux yz .

(3.73)
Repetindo-se o mesmo raciocnio para as outras faces obtemos o uxo
total atravs de todas as faces
Figura 3.18: A face 1 corresponde ao plano x = 0 e a face 2 paralela a este

plano.
ux uy uz
ux yz + uy xz + uz xy = ( + + )xyz .
x y z
(3.74)
Na forma diferencial esta equao pode ser reescrita como

ux uy uz
+ + dxdydz = udV . (3.75)
x y z
Pela denio de uxo,

u nS = uV , (3.76)
onde a soma refere-se a todas as faces do volume V .
Imaginemos agora que o volume V possa ser subdividido em pequenas
unidades innitesimais de volume Vi . Usando a ltima equao acima,
temos

u dV = uVi
V i

= u nSi
i
)
= u n dS . (3.77)
S
3.18 Teorema de Stokes

Seja S uma superfcie bilateral limitada por uma curva fechada C . O teo-
rema de Stokes estabelece que
)
u n dS = u dr . (3.78)
S C
Deixaremos de lado a demonstrao deste teorema.
3.19 Equaes de Maxwell

3.19.1 Lei de Gauss
Esta lei arma que o campo eltrico E produzido por uma certa distribuio
de carga pode ser determinado pela equao
)
q
E n dS = , (3.79)
S 0
onde q a carga total no interior de S . Pelo teorema do divergente (3.71),

temos
)
E n dS = E dV . (3.80)
S V
Por outro lado, a carga total pode ser expressa em termos da densidade de
carga atravs de

q= dV . (3.81)
V
Substituindo as duas ltimas equaes em (3.79), obtemos

E dV = 0 . (3.82)
V 0
Uma vez que V completamente arbitrrio, s podemos garantir que a
integral acima seja nula se e somente se
Esta lei foi primeiramente proposta por Carl Friedrich Gauss (1777 - 1855), nascido
na cidade de Brunswick, Alemanha. Trabalhou em diversos campos da matemtica e da
fsica dentre eles a teoria dos nmeros, geometria diferencial, magnetismo, astronomia e
tica.

E= . (3.83)
0
Esta a lei de Gauss na forma diferencial.
3.19.2 Lei de Ampre

O campo magntico B no interior e nas vizinhanas de uma distribuio de
corrente pode ser calculado pela lei de Ampre, a qual estabelece que
)
B dr = 0 I , (3.84)
C
onde I a corrente total no interior de C .

Do teorema de Stokes (3.78), temos que
)
B dr = B n dS . (3.85)
C S
Agora, a corrente total no interior de C em termos da densidade de

corrente J dada por

I= J n dS . (3.86)
S
Combinando as duas equaes acima com a lei de Ampre (3.84), temos

que

( B 0 J) n dS = 0 . (3.87)
S
Sendo S uma superfcie inteiramente arbitrria, encontramos
B = 0 J . (3.88)
Esta equao a lei de Ampre na forma diferencial. Veremos mais adiante

que esta lei est incompleta. necessrio adicionar um termo a mais do
lado direito da equao acima de tal forma a torna-la consistente com a lei
de conservao da carga.
Esta lei foi primeiramente enunciada por Andr Marie Ampre (1775 - 1836), nascido
em Polemieux-Le-Mont-d'Or, uma herdade prxima a Lyon, Frana. Devido grande
contribuio de Ampre eletrodinmica, Maxwell o apelidou de o Newton da Eletrici-
dade.
3.19.3 Lei de Faraday

O campo eltrico pode ser gerado por um campo magntico que varia no
tempo. A lei de Faraday! estabelece que
)
d
E dr = B n dS . (3.89)
C dt S
Usando o teorema de Stokes no lado esquerdo da lei de Faraday, encontramos

B
E+ n dS = 0 . (3.90)
S t
Com efeito
B
E= , (3.91)
t
que a lei de Faraday escrita na forma diferencial.
3.19.4 Ausncia de Monopolo Magntico

Assumindo que as linhas de campo magntico tem uma vorticidade, o uxo
de campo magntico atravs de qualquer superfcie fechada anula-se,
)
B n dS = 0 . (3.92)
S
Pelo teorema do divergente (3.71), temos que

)
B n dS = B dV = 0 . (3.93)
S V
Conseqentemente
B = 0. (3.94)
Esta equao signica que na natureza no h monopolo magntico.

! Michael Faraday (1791 - 1867) nasceu em Newington Butts, Surrey, em Londres,
Inglaterra. Considerado por muitos como o maior fsico experimental da histria da
cincia. De origem pobre (o seu pai trabalhava como ferreiro), logo cedo foi trabalhar em
uma livraria como transportador e encadernador. Nas horas vagas dedicava-se leitura.
Por isso, considerado por muitos como um autodidata.
3.19.5 Lei de Ampre-Maxwell

J mencionamos anteriormente que a lei de Ampre no descreve satisfato-
riamente a conservao da carga. Vejamos por que. Da equao (3.88), e
usando a identidade (3.53), temos que
( B) = 0 J = 0 . (3.95)
Este resultado est em desacordo com a equao de continuidade que ex-
pressa a conservao da carga,

J+ =0. (3.96)
t
Agora, usando a lei de Gauss (3.83), temos que

E E
J + 0 = J + 0 =0. (3.97)
t t
Vemos ento que esta equao e (3.95) so incongruentes. Este apa-
rente paradoxo pode ser resolvido adicionando-se um termo a mais na lei
de Ampre de maneira que a equao da continuidade seja satisfeita,
B = 0 J + 0 Jd . (3.98)
A densidade de corrente de deslocamento Jd pode ser determinada obser-
vando-se que
( B) = 0 (J + Jd ) = 0 . (3.99)
Comparando esta equao com (3.97), obtemos
E
Jd = 0
. (3.100)
t
"
Assim, as quatro equaes fundamentais de Maxwell do eletromagne-
tismo podem ser escritas como

E ,=
0
B = 0,
B
E = ,
t
E
B = 0 J + 0 0 . (3.101)
t
" James Clerk Maxwell (1831 - 1879) nasceu em Edimburgo, Esccia. Entre os anos
de 1860 e 1865 trabalhou no Kings College, de Londres, onde elaborou a teoria do ele-
tromagnetismo sintetizada nas quatro clebres equaes que levam o seu nome.
3.20 Problemas
3.1 Ache o volume do paraleleppedo de arestas A = 3i j, B = j + 2k,
C = i + 5j + 4k.
3.2 Determine a de modo que 2i 3j + 5k e 3i + aj 2k sejam perpendi-

culares.
3.3 Se A = i + j, B = 2i 3j + k, C = 4j 3k, calcule o produto vetorial

triplo A (B C) (a) usando a denio de produto vetorial, (b)
usando a relao BAC-CAB.
3.4 So dados os vetores
bc
a = ,
a (b c)
ca
b = ,
a (b c)
ab
c = .
a (b c)
Supor a (b c) = 0. Mostrar que
x y = xy ,
onde (x, y = a, b, c). Este resultado importante no estudo de redes

cristalinas em Fsica do Estado Slido.
3.5 Prove que a rea de um paralelogramo de lados ) e * dado por

|A B|.
*
3.6 Se A = (3x 6yz)i + (2y + 3xz)j + (1 4xyz)k, calcular C A dr de
(0,0,0) a (1,1,1) ao longo dos seguintes percursos C : (a) x = t, y = t2 ,
z = t3 ; (b) segmentos de reta de (0,0,0) a (0,0,1), deste a (0,1,1), e
deste a (1,1,1).
*
3.7 Se F = (x2 y 2 )i + 2xyj, calcular C A dr ao longo da curva C no
plano xy , dada por y = x2 x do ponto (1,0) a (2,2).
3.8 Usando os teoremas do divergente e de Stokes, caso conveniente, cal-

cule as seguintes integrais.
+
(a) u n dS , onde u = x3 i + y 3 j + z 3 k e S a esfera de raio R com
centro na origem.
+
(b) u n dS , onde u = x5 i + y 5 j + z 5 k e S a esfera do tem (a).
+
(c) C u dr, onde u = 3yi + 3xj + k, e C o crculo x2 + y 2 = 1
situado no plano z = 2.
3.9 Mostrar que

1
r n dA = V ,
3 S
onde V o volume limitado pela superfcie fechada S .

3.10 Se B = A, mostrar que

B n dS = 0 ,
S
para uma superfcie qualquer fechada S .

3.11 Sendo t = iy + jx, mostrar que
)
1
t dr = A ,
2 C
onde A a rea limitada pelo caminho fechado C .

3.12 As foras conservativas podem sempre ser escritas como menos o gra-
diente de uma funo escalar
F = .
Mostrar que o trabalho realizado por esta fora independe da trajet-

ria da partcula.
Analogamente, o campo eltrico dito conservativo se E = .
Mostrar que

2 = ,
0
onde a densidade de carga. Esta conhecida como equao de
Poisson. Se = 0, a equao conhecida como equao de Laplace.
3.13 A energia magntica (por unidade de volume) armazenada nos campos

eletromagnticos dada por
1 1
U= 0 E E + BB.
2 20
Usando as equaes de Maxwell, mostrar que
U
=JE+S,
t
onde S o vetor de Poyting e dado por
1
S= EB.
0
3.14 Supor que os campos eletromagnticos dependem apenas da coorde-

nada z e do tempo. Usando as equaes de Maxwell, mostrar que
By Ex
= 0 0 ,
z t
Ex By
= ,
z t
Bx Ey
= 0 0 ,
z t
Ey Bx
= .
z t
Usando estas equaes, mostrar que as componentes dos campos sa-

tisfazem a equao da onda
2 f (z, t) 1 2 f (z, t)
= ,
z 2 c2 t2

onde c = 1/ 0 0 a velocidade da luz no vcuo.
Captulo 4
Coordenadas Curvilneas
4.1 Introduo
Com muita freqncia, encontramos problemas em fsica que no possuem
uma simetria retangular. Portanto, conveniente mudarmos do sistema
de coordenadas retangulares para o sistema de coordenadas correspondente
simetria do problema em questo. Neste Captulo aprenderemos como
encontrar uma correspondncia entre os pontos do sistema de coordenadas
retangulares (x, y, z) e os pontos do sistema de coordenadas curvilneas que
denotaremos por (q1 , q2 , q3 ).
Na forma genrica, as equaes transformativas so escritas como
x = x(q1 , q2 , q3 ) ,
y = y(q1 , q2 , q3 ) ,
z = z(q1 , q2 , q3 ) . (4.1)
Iremos impor que o novo sistema de coordenadas seja ortogonal. Para
isto, necessrio encontrarmos os vetores unitrios que denam os senti-
dos positivos do eixos q1 , q2 , e q3 . Vejamos como encontrar estes vetores.
Usando regras de derivao parcial, podemos escrever
x x x
dx = dq1 + dq2 + dq3 ,
q1 q2 q3
y y y
dy = dq1 + dq2 + dq3 ,
q1 q2 q3
z z z
dz =
q1
dq1 +
q2
dq2 +
q3
dq3 . (4.2)
Captulo 4: Coordenadas Curvilneas 61
Estas trs equaes tambm podem ser escritas abreviadamente como
3
xi
dxi = dqj . (4.3)
j=1
qj
Primeiramente, consideremos q2 e q3 constantes. Ento,
dr = (dxi + dyj + dzk)q2 ,q3

x y z
= i+ j+ k dq1
q1 q1 q1

r
= dq1 . (4.4)
q1 q2 ,q3

O vetor r
q1 tangencial curva descrita por r com q2 e q3 cons-
q2 ,q3
1 o vetor unitrio na direo do eixo q1 , podemos
tantes. Denotando por a
escrever

1 r 1 x y z
1 =
a = i+ j+ k , (4.5)
h1 q1 q2 ,q3 h1 q1 q1 q1
onde

r 2 2 2
x y z
h1 = = + + , (4.6)
q1 q2 ,q3 q1 q1 q1
conhecida como mtrica ou fator de escala.
Repetindo-se o mesmo procedimento para os eixos q2 e q3 encontramos
para os outros vetores unitrios e mtricas,

1 r 1 x y z
2 =
a = i+ j+ k , (4.7)
h2 q2 q1 ,q3 h2 q2 q2 q2

1 r 1 x y z
a
3 = = i+ j+ k , (4.8)
h3 q3 q1 ,q2 h3 q3 q3 q3

r 2 2 2
x y z
h2 = = + + , (4.9)
q2 q1 ,q3 q2 q2 q2

r 2 2 2
x y z
h3 = = + + . (4.10)
q3 q1 ,q2 q3 q3 q3
Tambm podemos escrever na forma compacta,
3
1 xj
a
i = x
j , (4.11)
hi j=1 qi
3 2
xj
h2i = , (4.12)
j=1
qi
onde (
x1 , x 3 ) = (i, j, k).
3 , x
Como veremos adiante, as mtricas assumem um papel fundamental nas
novas representaes no sistema de coordenadas curvilneas.
Usando o fato que x i x
j = ij , e a equao (4.11), temos que
, 3
- , 3
-
1 xk 1 xl
i a
a j = x
k x
l ,
hi qi hj qj
k=1 l=1
3 3
1 xk xl
= k x
x l ,
hi hj qi qj
k=1 l=1
3 3
1 xk xl
= kl ,
hi hj qi qj
k=1 l=1
3
1 xk xk
= . (4.13)
hi hj qi qj
k=1
i
Se assumirmos que os novos vetores unitrios so ortogonais, ou seja, a aj =
ij , ento
3
1 xk xk
= ij . (4.14)
hi h j qi qj
k=1
Assim, Aij = (xi /qj )/hj uma matriz ortogonal.
4.2 Elementos de Comprimento, rea, e Vo-

lume
O comprimento de um elemento de arco ds dado por
ds2 = dr dr
= (dx)2 + (dy)2 + (dz)2
3

= dx2i .
i=1
(4.15)
Usando a equao (4.3), temos que
3 3 3
xi xi
ds2 = dqj dqk
i=1 j=1
q j q k
k=1
3 3
, 3 -
xi xi
= dqj dqk
j=1 i=1
qj qk
k=1
3
3

= hj hk jk dqj dqk
j=1 k=1
3

= h2j dqj2 , (4.16)
j=1
onde, na passagem da segunda para a terceira linha usamos a relao de

ortogonalidade (4.14).
O ltimo resultado sugere que no espao de coordenadas curvilneas o
elemento de volume um paraleleppedo de lados h1 dq1 , h2 dq2 , e h3 dq3 .
Conseqentemente, o elemento de volume em coordenadas curvilneas pode
ser escrito como
dV = h1 h2 h3 dq1 dq2 dq3 . (4.17)
Por m, o elemento de rea de qualquer uma das faces do paraleleppedo

ser dada por dij = hi hj dqi dqj , com i = j .
4.3 Gradiente, Divergente, e Rotacional

Seja (x, y, z) uma funo escalar com derivadas parciais contnuas. Usando
a regra da cadeia, temos que
x y z
= + +
q1 x q1 y q1 z q
1
x y z
= i+ j+ k i+ j+ k
x y z q1 q1 q1

r
= ()
q1 q2 ,q3
= () (h1 a
1 ) . (4.18)
De forma anloga, podemos mostrar que

= () (h2 a
2 ) , (4.19)
q2

= () (h3 a
3 ) . (4.20)
q3
As trs ltimas equaes tm solues nicas e so dadas por
1 1 1
=a
1 + a
2 + a
3 . (4.21)
h1 q1 h2 q2 h3 q3
A m de obter a representao do divergente em coordenadas curvil-
neas usaremos o mesmo procedimento que na Seo 3.17. Imaginemos um
elemento de volume em coordenadas curvilneas como ilustrado na Fig. 4.1.
O uxo lquido de um campo vetorial u atravs das faces 1, 2, e 3 dado
por
(u1 h2 h3 )
(u1 h2 h3 )q2 q3 = q1 q2 q3 ,
q1
(u2 h1 h3 )
(u2 h1 h3 )q1 q3 = q1 q2 q3 ,
q2
(u3 h1 h2 )
(u3 h1 h2 )q1 q2 = q1 q2 q3 . (4.22)
q3
Na forma diferencial o uxo lquido total pode ser escrito como

(u1 h2 h3 ) (u2 h1 h3 ) (u3 h1 h2 )
+ + dq1 dq2 dq3 = udV . (4.23)
q1 q2 q3
Por outro lado, usando a expresso do elemento de volume em coordenadas
curvilneas (4.17), obtemos

Figura 4.1: Elemento de volume em coordenadas curvilneas.

1 (u1 h2 h3 ) (u2 h1 h3 ) (u3 h1 h2 )
u= + + . (4.24)
h 1 h2 h3 q1 q2 q3
A determinao do rotacional em coordenadas curvilneas um pouco

mais complexa e a deixaremos de lado. Temos que

h1 a 3
1 1 h2 a
2 h3 a
u= . (4.25)
h1 h2 h3 q1 q2 q3
h1 u1 h2 u2 h3 u3
Da denio do Laplaciano (3.57), e das representaes do gradiente e
do divergente em coordenadas curvilneas (4.21) e (4.24), respectivamente,
temos que

1 h2 h3 h1 h3
2 = () = +
h1 h2 h3 q1 h1 q1 q2 h2 q2

h1 h2
+ . (4.26)
q3 h3 q3
Nas prximas duas Sees iremos aplicar o formulrio acima para dois
casos especcos de grande utilidade em fsica.
4.4 Coordenadas Cilndricas

Equaes transformativas:
x = cos , y = sin , z=z, (4.27)

onde 0, 0 2 , z (ver Fig. 4.2).
* 0
0
z
r
Figura 4.2: Coordenadas cilndricas.
Vetor de posio:
r = cos i + sin j + zk
= 0 + zk . (4.28)
Mtricas:
h1 = 1 , h 2 = , h3 = 1 . (4.29)
Elemento de volume:
dV = dddz . (4.30)
Gradiente:
1
= 0 + 0 +k . (4.31)
z
Divergente:
1 1 u uz
u= (u ) + + . (4.32)
z
Rotacional:

0 k
1 0
u =

z
.
(4.33)

u u uz
Laplaciano:

2 1 1 2 2
= + 2 +
2 z 2
2
1 1 2
2
= + + + . (4.34)
2 2 2 z 2
4.5 Coordenadas Esfricas

Equaes transformativas:
x = r sin cos , y = r sin sin , z = r cos , (4.35)

onde r 0, 0 , 0 2 (ver Fig. 4.3).
Vetor de posio:
r = r sin cos i + r sin sin j + r cos k

= rr0 . (4.36)
Mtricas:
h1 = 1 , h2 = r , h3 = r sin . (4.37)
Elemento de volume:
dV = r2 sin drdd . (4.38)

ro
(x, y, z) o
r o
(x, y, 0)
Figura 4.3: Coordenadas esfricas.
Gradiente:
1 1
= r0 + 0 + 0 . (4.39)
r r r sin
Divergente:

1 2 u
u= 2 sin (r ur ) + r (sin u ) + r . (4.40)
r sin r
Rotacional:

r0 r 0 r sin 0
1
u= 2 . (4.41)
r sin r
ur ru r sin u
Laplaciano:

2 1 2 1 1 2
= 2 r + 2 sin + 2 2 . (4.42)
r r r r sin r sin 2
4.6 Separao de Variveis

No que segue abaixo, listamos algumas das principais equaes diferenciais
parciais da fsica.
Equao de Laplace:
2 = 0 . (4.43)
Equao de Poisson:

2 = . (4.44)
0
Equao de difuso independente do tempo:
2 k 2 = 0 . (4.45)
Equao de difuso dependente do tempo:
1
2 = . (4.46)
a2 t
Equao de Schrdinger:
2 2
+ V = E . (4.47)
2m
A incgnita nas equaes acima . Ao resolvermos as equaes acima,
tentamos uma forma de transferir a soluo das equaes diferenciais par-
ciais para a soluo de equaes diferenciais ordinrias mais simples. Este
objetivo pode ser alcanado, usando-se o mtodo de separao de variveis.
Para ilustrar como este mtodo funciona, consideremos a equao de difu-
so independente do tempo com o sinal mais. A idia separar a equao
em trs equaes diferenciais ordinrias e independentes cada qual em uma
nica varivel. Escrevemos
(x, y, z) = X(x)Y (y)Z(z) . (4.48)

Substituindo esta equao em (4.45), obtemos
d2 X d2 Y d2 Z
YZ + XZ + XY + k 2 XY Z = 0 , (4.49)
dx2 dy 2 dz 2
ou ainda
1 d2 X 2 1 d2 Y 1 d2 Z
= k . (4.50)
X dx2 Y dy 2 Z dz 2
Uma vez que as coordenadas (x, y, z) so independentes, ambos os lados

da equao acima devem ser iguais a uma constante. Temos
1 d2 X
= l2 (4.51)
X dx2
2
1 d Y 1 d2 Z
k 2 = l2 . (4.52)
Y dy 2 Z dz 2
Rearranjando a ltima das duas equaes, tambm podemos escrever
1 d2 Y 2 2 1 d2 Z
= k + l . (4.53)
Y dy 2 Z dz 2
Igualmente ao caso anterior, ambos os membros da equao acima so
iguais se ambos forem constantes. Temos
1 d2 Y 2
2 1 d Z
= m , = n2 , (4.54)
Y dy 2 Z dz 2
onde k 2 = l2 + m2 + n2 .
Para cada trinca de valores (l, m, n) teremos uma soluo lmn . Como
a equao (4.45) linear (ver Captulo 5) em , uma combinao linear de
solues tambm ser uma soluo da equao diferencial,

(x, y, z) = almn lmn (x, y, z) . (4.55)
lmn
As constantes (l, m, n), e almn so determinadas de acordo com as con-

dies de contorno impostas sobre a funo (x, y, z). Para o momento,
o que importante notar, que transformamos o problema de resolver
uma equao diferencial parcial na soluo de trs equaes diferenciais
ordinrias.
Este ltimo exemplo aplica-se a problemas com simetria retangular. Su-
por agora que a simetria seja esfrica. Neste caso, usamos o Laplaciano em
coordenadas esfricas (4.42). A equao de difuso independente do tempo
pode ser escrita como

1 2 1 1 2
r + 2 sin + + k 2 = 0 . (4.56)
r2 r r r sin r2 sin2 2
As condies de contorno so especicadas pelos valores da funo e/ou suas deri-
vadas nas fronteiras de um determinado domnio (para maiores detalhes, ver Captulo
8).
Usando o mtodo de separao de variveis, podemos escrever
(r, , ) = R(r)()() . (4.57)

Substituindo esta equao em (4.56), encontramos

1 d2 2 2 2 1 d 2 dR 1 d d
= r sin k + r + sin .
d2 r2 R dr dr r2 sin d d
(4.58)
Analogamente ao caso de coordenadas retangulares, assumimos (r, , )
como coordenadas independentes. Temos que
1 d2
+ m2 = 0 , (4.59)
d2

1 d 2 dR QR
r + k2 R 2 = 0 , (4.60)
r2 dr dr r
2
1 d d m
sin + Q = 0 , (4.61)
sin d d sin2
onde m e Q so duas constantes de separao. Mais uma vez, vemos que
a soluo de uma equao diferencial parcial se reduz soluo de trs
equaes diferenciais ordinrias.
Consideremos o caso particular k = 0, o qual corresponde equao de
Laplace (4.43). Obtemos

1 d 2 dR QR
r 2 =0. (4.62)
r2 dr dr r
Supor que a soluo seja do tipo R(r) = Arn , onde A uma constante.
Substituindo na equao diferencial para R, vem que
[n(n + 1) Q]rn2 = 0 . (4.63)

Assumindo r = 0, encontramos Q = n(n + 1).
Assim, a equao para , com m = 0, toma a forma

1 d d
sin + n(n + 1) = 0 . (4.64)
sin d d
A mudana de varivel x = cos , e o uso da regra da cadeia nos leva a
equao
d2 d
(1 x2 ) 2x + n(n + 1) = 0 . (4.65)
dx2 dx
Esta equao conhecida como equao de Legendre. No prximo Captulo

veremos como resolve-la.
4.7 Problemas
4.1 A transformao de coordenadas retangulares para coordenadas cilin-
dro-parablicas denida pelas equaes x = 12 (u2 v 2 ), y = uv ,
z = z . (a) Calcule as mtricas h1 , h2 e h3 . (b) Prove que a matriz
Aij = h1j x
qj ortogonal. (c) Escrever o gradiente, o divergente, o
i
Laplaciano e o rotacional em coordenadas cilindro-parablicas.

4.2 Considere as equaes transformativas de coordenadas abaixo:
x = a cosh u cos v ,
y = a sinh u sin v ,
z = z.
(a) Calcule as mtricas h1 , h2 e h3 . (b) Determinar a matriz Aij =

1 xi
hj qj e mostrar que ortogonal. (c) Escrever o gradiente, o diver-
gente, o Laplaciano e o rotacional no novo sistema de coordenadas.
4.3 As equaes transformativas de coordenadas retangulares xyz para
esfricas so dadas por x = r sin cos , y = r sin sin , z = r cos .
Monte a matriz Bij = h1i qij . Calcule a matriz inversa de B .
x
4.4 Usando as equaes (4.5-4.10) e o resultado do Problema 4.3, mos-

tre que os vetores unitrios nos dois sistemas de coordenadas esto
relacionados pelas equaes:
i = r0 sin cos + 0 cos cos 0 sin ,

j = r0 sin sin + 0 cos sin + 0 cos ,
k = r0 cos 0 sin .
4.5 Mostrar que (/x, /y, /z) em coordenadas esfricas so dados

por
1 sin
= sin cos + cos cos ,
x r r r sin
1 cos
= sin sin + cos sin + ,
y r r r sin
1
= cos sin .
z r r
Sugesto : Calcular a inversa da matriz xj

qi e multiplicar o resultado
pela matriz coluna (/r, /, /).
4.6 O operador momento angular denido por L = i(r ), onde i
a unidade imaginria e = h/2 , sendo h a constante de Planck. (a)
Usando os resultados do Problema 4.5, mostrar que Lz dado por

Lz = i x y = i .
y x
(b) Mostrar que

L+ = Lx + iLy = ei + i cot ,

L = Lx iLy = ei i cot .

Em Mecnica Quntica, L+ e L so conhecidos como operadores de
levantamento e abaixamento, respectivamente.
4.7 Considere a equao de Laplace em coordenadas polares-cilndricas:

1 1 2 2
+ + =0.
2 2 z 2
Usando o mtodo de separao de variveis, mostre que a soluo
desta equao se reduz soluo de trs equaes diferenciais lineares
ordinrias dadas por
d2 Z
l2 Z = 0 ,
dz 2
d2
+ m2 = 0 ,
d2

1 d dR 2 m2
+ l 2 R = 0.
d d
Fazendo a mudana de varivel x = l, mostrar que a equao radial

pode ser escrita como

1 d dR m2
x + 1 2 R=0. (4.66)
x dx dx x
Esta equao conhecida como equao diferencial de Bessel. No

Captulo 5 iremos resolver esta equao.
Captulo 5
Equaes Diferenciais
Ordinrias
5.1 Introduo
Uma equao diferencial ordinria e linear de segunda ordem tem a seguinte
forma geral
d2 y
A(x)
dx2
+ B(x)
dy
dx
+ C(x)y = D(x) , (5.1)
onde A(x), B(x), C(x) e D(x) so funes conhecidas; y(x) a funo
incgnita.
Esta equao tambm pode ser escrita na forma
L y(x) = D(x) , (5.2)
onde
d2
L = A(x)
dx 2
+ B(x)
d
dx
+ C(x) (5.3)
um operador diferencial. Este operador tem as seguintes propriedades
L [y1 (x) + y2 (x)] = L y1 (x) + L y2 (x) ,

L [ay(x)] = aL y(x) , (5.4)
onde a uma constante. Operadores que possuem estas propriedades so
chamados de operadores lineares.
A equao (5.1) tambm pode ser escrita na forma
d2 y dy
2
+ P (x) + Q(x)y = R(x) , (5.5)
dx dx
onde
B(x)
P (x) = , (5.6)
A(x)
C(x)
Q(x) = , (5.7)
A(x)
D(x)
R(x) = . (5.8)
A(x)
Se R(x) = 0, a equao diferencial chamada de equao diferencial

homognea (EDH). Caso contrrio, denominamos esta equao de equao
diferencial no-homognea (EDNH). Primeiro vamos estudar a EDH.
5.2 Soluo Geral da EDH

Considere a EDH
d2 y dy
+ P (x) + Q(x)y = 0 . (5.9)
dx2 dx
Se y1 (x) soluo de (5.9), ento y = C1 y1 (x), tambm uma soluo,
onde C1 uma constante arbitrria. Se y = y1 (x), e y = y2 (x) so solues
de (5.9), ento y = C1 y1 (x) + C2 y2 (x), tambm uma soluo de (5.9).
Estas armaes so verdadeiras por causa que (5.9) linear em y e suas
derivadas. Em outras palavras, (5.9) segue as propriedades (5.4).
Um par de solues y = y1 (x) e y = y2 (x) so denominadas linearmente
independentes se a igualdade
C1 y1 + C2 y2 = 0 , (5.10)
onde os C 's so constantes, for verdadeira se e somente se C1 = C2 = 0.
Se as solues forem linearmente dependentes, y1 = ky2 , onde k uma
constante que no depende de x. Logo, podemos escrever y1 = ky2 , e
y1 y2
= , (5.11)
y1 y2
ou ainda
Captulo 5: Equaes Diferenciais Ordinrias 77
y1 y2 y1 y2 = 0 . (5.12)

A funo
W (x) = y1 (x)y2 (x) y1 (x)y2 (x) , (5.13)

conhecida como Wronskiano. Uma condio necessria e suciente para
que um par de solues sejam linearmente independentes que

y y2
W (x) = 1 = 0 . (5.14)
y1 y2
Concluindo, se y = y1 (x) e y = y2 (x) um par de solues linearmente
independentes, ento a soluo geral da EDH dada por
y = C1 y1 (x) + C2 y2 (x) , (5.15)

onde C1 e C2 so duas constantes arbitrrias.
5.3 Segunda Soluo

Agora aprenderemos como construir uma soluo y2 (x) dado que conhece-
mos uma primeira soluo y1 (x). Derivando (5.13), obtemos
W = y1 y2 + y1 y2 y1 y2 y1 y2

= y1 [P (x)y2 Q(x)y2 ] y2 [P (x)y1 Q(x)y1 ]
= P (x)(y1 y2 y1 y2 ) . (5.16)
A expresso entre parnteses exatamente o Wronskiano. Temos

dW
= P (x)W . (5.17)
dx
Esta equao tambm pode ser escrita na forma
dW
= P (x) dx . (5.18)
W
Integrando ambos os lados da equao acima, encontramos

ln W = P (x) dx + C , (5.19)
onde C uma constante de integrao. O Wronskiano pode ser facilmente

obtido
*
W (x) = W0 e P (x) dx
, (5.20)
onde W0 = eC .
Derivando o quociente y2 /y1 , vem que

d y2 y1 y2 y1 y2
=
dx y1 y12
W
= . (5.21)
y12
Integrando esta ltima equao, encontramos

*
y2 (x) e P (x) dx
= dx , (5.22)
y1 (x) [y1 (x)]2
onde, sem perda de generalidade, tomamos W0 = 1. Alm disso, todas as
constantes de integrao foram ignoradas. Levar em conta estas constantes,
apenas muda os valores das constantes C1 e C2 da soluo geral da EDH.
Portanto, nada de novo acrescentado soluo geral.
5.4 Soluo Geral da EDNH

Se u(x) uma soluo particular de (5.5), ento a soluo geral da EDNH
dada por
y = C1 y1 (x) + C2 y2 (x) + u(x) . (5.23)

Veremos agora como determinar u(x). Conhecendo-se um par de solu-
es da EDH, possvel determinarmos uma soluo particular da EDNH.
Para isto, usamos o mtodo de variao das constantes. A idia escrever a
soluo particular da EDNH na mesma forma que a soluo geral da EDH,
substituindo-se as constantes arbitrrias C1 e C2 por duas funes v1 (x) e
v2 (x). Temos
u(x) = v1 (x)y1 (x) + v2 (x)y2 (x) . (5.24)

O prximo passo a determinao de v1 e v2 . Derivando u(x), temos
u = v1 y1 + v1 y1 + v2 y2 + v2 y2 . (5.25)

Agora, note que as funes v1 e v2 no so nicas, pois existe uma funo
arbitrria g tal que

y1
v1 (x)y1 (x) + v2 (x)y2 (x) = (v1 + g)y1 + v2 g y2 . (5.26)
y2
Portanto, podemos escolher v1 e v2 convenientemente tal que
v1 y1 + v2 y2 = 0 . (5.27)

Substituindo esta equao em (5.25), obtemos
u = v1 y1 + v2 y2 . (5.28)

Derivando-se novamente, vem que
u = v1 y1 + v2 y2 + v1 y1 + v2 y2 . (5.29)

Substituindo (5.24), (5.28), e (5.29) em (5.5), vem que
v1 y1 + v2 y2 + v1 y1 + v2 y2

+P (x)(v1 y1 + v2 y2 ) + Q(x)(v1 y1 + v2 y2 ) = R , (5.30)
ou ainda
v1 y1 + v2 y2 + v1 [y1 + P (x)y1 + Q(x)y1 ]

+v2 [y2 + P (x)y2 + Q(x)y2 ] = R . (5.31)
Os termos entre colchetes se anulam pois y1 e y2 so solues da EDH.

Assim, juntamente com (5.27) obtemos uma segunda equao para as deri-
vadas de v1 e v2 . Temos
v1 y1 + v2 y2 = R . (5.32)

As equaes (5.27) e (5.32) podem ser facilmente resolvidas para as
derivadas. Obtemos
Ry2
v1 = , (5.33)
W
Ry1
v2 = . (5.34)
W
Integrando, nalmente obtemos

y2 (x)R(x)
v1 (x) = dx ,
W (x)

y1 (x)R(x)
v2 (x) = dx . (5.35)
W (x)
As constantes de integrao podem ser tomadas todas nulas, sem perda

de generalidade, conforme justicado anteriormente.
Exemplo 5.1 Considere a seguinte EDNH
d2 y(x)
2 y(x) = ex , (5.36)
dx2
onde uma constante no nula. Determinar a segunda soluo.
Soluo: Por substituio direta na EDH, podemos vericar que
y1 (x) = ex soluo. Comparando esta equao com (5.5), encontra-
mos P (x) = 0 e Q(x) = 2 . Conseqentemente, usando (5.22), temos
que

y2 (x) 1
= e2x dx = e2x . (5.37)
ex 2
ou ainda
y2 (x) = ex , (5.38)
onde a constante multiplicativa foi ignorada. Como j mencionamos anteri-
ormente, esta constante no acrescenta nada de novo soluo. Usando-se
(5.13), obtemos W (x) = 2. De posse do par de solues da EDH e do
Wronskiano, podemos determinar a soluo particular da EDNH. De (5.35),
temos que

ex ex x
v1 (x) = dx = , (5.39)
2 2
x x
e e 1
v2 (x) = dx = 2 e2x . (5.40)
2 2
Substituindo estes resultados em (5.24), vem que

1 1
u(x) = x ex . (5.41)
2 2

5.5 Mtodo de Frobenius

Consideremos a EDH na forma geral
d2 y dy
A(x) + B(x) + C(x)y = 0 , (5.42)
dx2 dx
onde A(x), B(x), e C(x) so polinmios em x. Denominamos x = a um
ponto ordinrio da EDH se A(a) = 0, e ser um ponto singular em caso
contrrio.
Se x = 0 um ponto ordinrio, ento a EDH pode ser resolvida em srie
nas vizinhanas de x = 0 como

y(x) = cj xj
j=0
= C1 {srie} + C2 {srie} , (5.43)
onde C1 e C2 so duas constantes arbitrrias. As duas sries so linearmente

independentes e ambas convergem nas vizinhanas de x = 0.
Um ponto x = a dito singular regular se a EDH pode ser escrita na
forma
d2 y R1 (x) dy R2 (x)
2
+ + y=0, (5.44)
dx (x a) dx (x a)2
onde R1 (x) e R2 (x) admitem expanso em srie de Taylor nas vizinhanas
de x = a. Em caso contrrio, x = a chamado de um ponto singular
irregular.
Uma forma alternativa de vericar se um ponto singular x = a regular
ou irregular testar os limites
lim (x a)P (x) ,

xa
lim (x a)2 Q(x) . (5.45)
xa
Se estes limites existirem, dizemos que x = a um ponto singular regular

e, em caso contrario, irregular.
Se x = a um ponto singular irregular, a EDH no possui uma soluo
em forma de srie nas vizinhanas deste ponto. Se x = a um ponto singular
regular e, R1 (x) e R2 (x) puderem ser desenvolvidas em srie de Taylor nas
vizinhanas deste ponto, ento a EDH admite uma soluo em forma de
srie

y(x) = cj xj+s , (5.46)
j=0
onde c0 = 0, s e todos os demais coecientes da srie so constantes a

determinar de tal modo a satisfazer a EDH.
Este procedimento para resolver as equaes diferenciais lineares de se-
gunda ordem conhecido como mtodo de Frobenius.
Exemplo 5.2 Resolver a EDH
(1 + x2 )y + xy y = 0 , (5.47)
usando o mtodo de Frobenius.
Soluo: Temos A(x) = 1 + x2 . Esta funo no possui nenhuma raiz
real, e portanto no tem ponto singular. Uma vez que A(0) = 0, ento
x = 0 um ponto ordinrio da EDH. Com efeito, a soluo da EDH admite
uma soluo em forma de srie como em (5.43). Temos

y(x) = cj xj ,
j=0

y (x) = cj jxj1 ,
j=0

y (x) = cj j(j 1)xj2 . (5.48)
j=0
Substituindo na EDH, obtemos

(1 + x2 ) cj j(j 1)xj2 + x cj jxj1 cj xj = 0 , (5.49)
j=0 j=0 j=0
que ainda pode ser escrita como

cj j(j 1)xj2 + cj j(j 1)xj + cj jxj cj xj = 0 . (5.50)
j=0 j=0 j=0 j=0
A primeira soma tambm pode ser escrita como

Ferdinand Georg Frobenius (1849 - 1917), nasceu em Berlim, Alemanha.

cj j(j 1)xj2 = cj j(j 1)xj2
j=0 j=2

= cj+2 (j + 2)(j + 1)xj . (5.51)
j=0
Substituindo em (5.50), encontramos

{cj+2 (j + 2)(j + 1) + cj [j(j 1) + j 1]} xj = 0 . (5.52)
j=0
Igualando a zero todas os coecientes das potncias de x, encontramos

j1
cj+2 = cj , j 0 . (5.53)
j+2
Esta equao conhecida como relao de recorrncia. Note que para j = 1
obtemos c3 = 0. Logo, para todo j mpar, temos cj = 0, exceto c1 . Se j for
par, temos que
1 1 1
c2 =c0 , c4 = c2 = c0 , . . . . (5.54)
2 4 8
Para j par (j = 2k ), temos que
2k 3
c2k = c2k2
2k
(2k 3)(2k 5)
= c2k4
2k(2k 2)
..
.
(2k 3)!!
= (1)k+1 c0 , k 2 . (5.55)
2k k!
A soluo geral da EDH dada por

1 2 k+1 (2k 3)!! 2k
y(x) = c0 1 + x + (1) x + c1 x . (5.56)
2 2k k!
k=2

c
Note que limj j+2 j1
cj = limj j+2 = 1. Portanto, a srie converge
para todo |x| < 1.

2x2 y xy + (x2 + 1)y = 0 , (5.57)
Soluo: Neste caso A(x) = 2x2 , e x = 0 uma raiz desta funo. Logo
x = 0 um ponto singular regular da EDH, pois esta pode ser escrita na
forma
1 (x2 + 1)
y
y + y=0. (5.58)
2x 2x2
Assim, a soluo da EDH admite uma soluo em forma de srie como
em (5.46). Temos

y(x) = cj xj+s ,
j=0

y (x) = cj (j + s)xj+s1 ,
j=0

y (x) = cj (j + s)(j + s 1)xj+s2 . (5.59)
j=0
Substituindo na EDH, vem que

2x2 cj (j + s)(j + s 1)xj+s2 x cj (j + s)xj+s1
j=0 j=0

+(x2 + 1) cj xj+s = 0 , (5.60)
j=0

2cj (j + s)(j + s 1)xj+s + (1)cj (j + s)xj+s
j=0 j=0

+ cj xj+s+2 + cj xj+s = 0 . (5.61)
j=0 j=0

A terceira soma tambm pode ser escrita como j=2 cj2 xj+s . Con-
seqentemente,
c0 [2s(s 1) s + 1] xs + c1 [2(s + 1)s (s + 1) + 1] xs+1

+ {cj [2(j + s)(j + s 1) (j + s) + 1] + cj2 } xj+s = 0 . (5.62)
j=2
Todos os coecientes da srie devem se anular. Temos
c0 [2s(s 1) s + 1] = 0,
c1 [2(s + 1)s (s + 1) + 1] = 0. (5.63)
A primeira destas equaes conhecida como equao indicial. Supor

c0 = 0. Ento
2s(s 1) s + 1 = 2s2 3s + 1 = 0 , (5.64)
cujas razes so dadas por s = 1/2 e s = 1. Substituindo estes valores

2
na segunda das equaes (5.63), vem que s = 1/2 c1 2(1/2) +2 1/2 =
[1] c1 = 0 c1 = 0. Analogamente, temos que s = 1 c1 2(1) + 1 =
[3] c1 = 0 c1 = 0. Se em algum destes dois casos, o termo entre colchetes
fosse exatamente zero, ento c1 seria no nula.
Para todo j 2, da equao (5.62), obtemos
cj [2(j + s)(j + s 1) (j + s) + 1] + cj2 = 0 , (5.65)
ou ainda
1
cj = cj2 . (5.66)
(j + s)(2j + 2s 3) + 1
Para s = 1/2, esta equao toma a forma
1
cj = cj2 (5.67)
j(2j 1)
A partir desta relao de recorrncia geramos todos os coecientes da
srie,
Do ingls, indicial. Alguns autores traduzem este termo por indicadora. J em outros
livros este termo no recebeu traduo permanecendo como equao indicial. Neste livro
usamos esta ltima opo.
1 1
c2 = c0 = c0 ,
23 6
1 1
c4 = c2 = c0 ,
47 168
..
. (5.68)
Uma vez que c1 = 0, temos cj = 0 para j mpar. Substituindo os coecientes

na srie, com s = 1/2, encontramos

y1 (x) = cj xj+1/2
j=0

x2 x4
= c0 x 1 + . (5.69)
6 168
Igualmente para s = 1, obtemos
1
cj = cj2 , (5.70)
j(2j + 1)
1 1
c2 = c0 = c0 ,
25 10
1 1
c4 = c2 = c0 ,
49 360
..
. (5.71)
Substituindo os coecientes na srie, com s = 1, obtemos uma segunda

soluo

y2 (x) = cj xj+1
j=0

x2 x4
= c0 x 1 + . (5.72)
10 360
As duas solues y1 (x) e y2 (x) so linearmente independentes. A soluo
geral da EDH uma combinao linear destas duas funes.
Para ambos
cj
os casos, pode-se vericar facilmente que limj cj2 = 0 . Portanto, as
sries convergem para todo < x < .
5.6 Equao Indicial

As solues da equao indicial podem ser classicadas em trs casos:
1. as razes so distintas, no diferindo por um inteiro;
2. as razes so iguais;
3. as razes so distintas, porm diferindo por um inteiro.
O primeiro caso est ilustrado no Exemplo 5.3. Quando as razes s1 e

s2 da equao indicial so iguais, as solues correspondentes so idnticas.
Neste caso as duas solues linearmente independentes so obtidas por
y1 (x) = y(x, s1 ) ,

y(x, s)
y2 (x) = . (5.73)
s s=s1
Quando as razes s1 < s2 da equao indicial diferirem por um inteiro, a

maior das duas, s2 , d uma soluo, enquanto que a menor, s1 , poder dar
ou no. No caso de no dar uma soluo, ento fazemos c0 = b0 (s s1 ) e
denotemos a nova funo por y(x, s). A solues linearmente independentes
sero dadas por
y1 (x) = y(x, s1 ) ,

y(x, s)
y2 (x) = . (5.74)
s s=s1
Vejamos alguns exemplos ilustrativos de como o mtodo funciona.

xy + y y = 0 , (5.75)
Soluo: A equao diferencial tambm pode ser escrita na forma
1 x
y + y 2y = 0 . (5.76)
x x
Portanto, x = 0 um ponto singular regular da EDH. Substituindo as sries
(5.59) na EDH, obtemos

x cj (j + s)(j + s 1)xj+s2 + cj (j + s)xj+s1
j=0 j=0

cj xj+s = 0 , (5.77)
j=0

cj [(j + s)(j + s 1) + (j + s)] xj+s1 + (1)cj xj+s = 0 . (5.78)
j=0 j=0
Fazendo a mudana de ndice j (j + 1) na primeira soma, temos que

c0 [s(s 1) + s] xs1 + cj+1 [(j + 1 + s)(j + s) + (j + 1 + s)] xj+s
j=0

+ (1)cj xj+s = 0 ,
j=0
(5.79)
ou ainda

. /
c0 s2 xs1 + cj+1 (j + s + 1)2 cj xj+s = 0 . (5.80)
j=0
Supondo c0 = 0 obtemos a equao indicial s2 = 0, cuja raiz dupla

dada por s = 0. Da equao (5.80) obtemos a relao de recorrncia
1
cj+1 = cj . (5.81)
(s + j + 1)2
Atravs desta relao podemos gerar todos os coecientes da srie. Com
s = 0, temos que
c1 = c0 ,
1 1
c2 = c0 = c0 ,
22 (2!)2
1 1 1
c3 = 2
c1 = 2 3 c0 = c0 ,
3 2 3 (3!)2
..
.
1
cn = c0 . (5.82)
(n!)2
Uma soluo ento dada por

1 2 1 3
y1 (x) = y(x, s)|s=0 = c0 1 + x + x + x + . (5.83)
(2!)2 (3!)2
Da primeira das equaes (5.59), temos para y2 (x)

y(x, s)
y2 (x) = ,
s
s=0

j dcj
= ln x cj (0)x + xj , (5.84)
j=0 j=0
ds s=0
onde cj (0) = cj (s = 0), cujos valores so dados por (5.82). De (5.81), para
s arbitrrio, encontramos os seguintes coecientes
1
c1 = c0 ,
(s + 1)2
1
c2 = c0 ,
(s + 1)2 (s + 2)2
1
c3 = c0 ,
(s + 1)2 (s + 2)2 (s + 3)2
..
.
1
cn = c0 . (5.85)
(s + 1)2 (s + 2)2 (s + 3)2 (s + n)2
Embora envolvendo uma extensa lgebra, fcil mostrar que

dc1
= 2c0 ,
ds s=0

dc2 2 1
= 1+ c0 ,
ds s=0 (2!)2 2

dc3 2 1 1
= 1+ + c0 . (5.86)
ds s=0 (3!)2 2 3
Substituindo (5.82) e (5.86) em (5.84), nalmente obtemos

1 2 1 3
y2 (x) = c0 ln x 1 + x + x + x +
(2!)2 (3!)2

x2 1 x3 1 1
2c0 x + 1+ + 1+ + + .
(2!)2 2 (3!)2 2 3
(5.87)

x2 y + 2xy + (x2 2)y = 0 , (5.88)

Soluo: Neste caso A(x) = x2 , e x = 0 uma raiz desta funo. Logo
x = 0 um ponto singular regular da EDH, pois esta pode ser escrita na
forma
2 (x2 2)
y + y + y=0. (5.89)
x x2
Assim, a soluo da EDH admite uma soluo em forma de srie como
em (5.46). Substituindo (5.59) na EDH, vem que

x2 cj (j + s)(j + s 1)xj+s2 + 2x cj (j + s)xj+s1
j=0 j=0

2
+(x 2) cj xj+s = 0 , (5.90)
j=0

cj (j + s)(j + s 1)x j+s
+ 2cj (j + s)xj+s
j=0 j=0

+ cj xj+s+2 + (2)cj xj+s = 0 . (5.91)
j=0 j=0

A terceira soma tambm pode ser escrita como j=2 cj2 xj+s . Con-
seqentemente,
c0 [s(s 1) + 2s 2] xs + c1 [(s + 1)s + 2(s + 1) 2] xs+1

+ {cj [(j + s)(j + s 1) + 2(j + 2) 2] + cj2 } xj+2 = 0 . (5.92)
j=2
Todos os coecientes da srie devem se anular. Temos
c0 [s(s 1) + 2s 2] = 0 ,
c1 [(s + 1)s2(s + 1) 2] = 0 . (5.93)
Supor c0 = 0. Ento
s(s + 1) + 2s 2 = (s + 2)(s 1) = 0 , (5.94)

cujas razes so dadas por s = 1 e s = 2. Neste caso as razes da equao
indicial diferem por um inteiro. O valor de s = 2 pode dar ou no uma
soluo. Substituindo os valores de s na segunda das equaes (5.93), vem
que s = 1 c1 [2 + 4 2] = 2c1 = 0 c1 = 0. Analogamente temos
s = 2 c1 [2 2 2] = 2c1 = 0 c1 = 0. Se em algum destes dois
casos, o termo entre colchetes fosse exatamente zero, ento c1 seria no nula.
Para todo j 2, da equao (5.92), obtemos
cj [(j + s)(j + s 1) + 2(j + 2) 2] + cj2 = 0 , (5.95)

ou ainda
1
cj = cj2 . (5.96)
(j + s 1)(j + s + 2)
A partir desta relao de recorrncia, para s = 1 geramos todos os
coecientes da srie
1 1
c2 = c0 = c 0 ,
25 10
1 1
c4 = c2 = c0 ,
47 280
..
. (5.97)
Uma vez que c1 = 0, temos cj = 0 para j mpar. Substituindo os coecientes

na srie, com s = 1, encontramos

y1 (x) = cj xj+1
j=0

x2 x4
= c0 x 1 +
10 280

1 2k
= c0 x 1 + 3 k
(1) k x . (5.98)
2 k!(2k + 3)!!
k=1
Igualmente, para s = 2, obtemos
1 1
c2 = c0 = c0 ,
(1) 2 2
1 1
c4 = c2 = c0 ,
14 8
1 1
c6 = c4 = c0 ,
36 144
..
. (5.99)
Note que neste caso pudemos gerar todos os coecientes sem diculda-
des. Assim, no precisamos regularizar a soluo como veremos no prximo
exemplo. Substituindo os coecientes na srie, com s = 2, obtemos uma
segunda soluo

y2 (x) = cj xj2
j=0

1 x2 x4 x6
= c0 2 1 + +
x 2 8 144

1 x2 k+1 1 2k
= c0 2 1 + + (1) x .
x 2 (2k 3)!!(2k)!!
k=2
(5.100)
As duas solues y1 (x) e y2 (x) so linearmente independentes. A solu-

o geral da EDH uma combinao linear destas duas funes. Pode-se
cj
vericar facilmente que limj cj2 = 0. Portanto, a srie de y1 (x) con-
verge para todo < x < . O mesmo ocorre no caso da srie de y2 (x)
exceto para o caso x = 0.

xy 3y + xy = 0 , (5.101)

Soluo: fcil ver que x = 0 um ponto singular regular. Substi-
tuindo as sries de y , y e y na EDH, encontramos
c0 [s(s 4)] xs1 + c1 [(s 3)(s + 1)] xs

+ j=0 [(j + s 4)(j + s)cj + cj2 ] xj+s1 = 0 . (5.102)
As razes da equao indicial so s = 0 e s = 4. Ambas as razes anulam

a constante c1 e portanto acarretar a nulidade de todos os coecientes
mpares.
A relao de recorrncia dada por
1
cj = cj2 . (5.103)
(j + s 4)(j + s)
Note que esta relao, para s = 4, d valores nitos para todos os coeci-
entes. No entanto, quando s = 0, c4 . Seguindo a regra, substitumos
c0 por b0 (s 0) = b0 s. Obtemos para y

s 1
y(x, s) = b0 x s
s
x2 + x4
(s 2)(s + 2) (s 2)(s + 2)(s + 4)
1
x6
(s 2)(s + 2)2 (s + 4)(s + 6)

1 8
+ x . (5.104)
(s 2)(s + 2)2 (s + 4)2 (s + 6)(s + 8)
O restante da soluo anlogo ao Exemplo 5.5, e portanto apenas

apresentamos as duas solues linearmente independentes

1
y1 (x) = y(x, 0) = c0 x4
224

1 6 1 8
+ x x .
2 22 4 6 2 22 42 6 8
(5.105)

y(x, 0) 1 1
y2 (x) = = y1 (x) ln x + 1 + 2 x2 + 5 x4
s 2 2 2!
s=0
1 1 1
6 1+ + x6
2 3!1! 2 3

1 1 1 1 1 8
8 1+ + + + x . (5.106)
2 4!2! 2 3 4 2
5.7 Coecientes Constantes

A equao
y + by + cy = 0 , (5.107)
onde b e c so duas constantes, a forma mais geral de uma equao di-

ferencial linear homognea de coecientes constantes. Uma vez que, pelo
menos uma das constantes no nula, vemos claramente que y , y e y
so linearmente dependentes. A nica funo cuja primeira derivada que
no linearmente independente de si prpria a exponencial. Desta forma,
escrevemos a soluo de (5.107) como
y(x) = epx . (5.108)
Substituindo em (5.107), vem que
(p2 + bp + c)epx = 0 p2 + bp + c = 0 , (5.109)
o qual conhecido como polinmio caracterstico , cujas razes podem ser:

(a) reais e distintas; (b) reais e iguais; (c) complexas, sendo que uma ser o
conjugado da outra. Sejam (p1 , p2 ) duas razes do polinmio caracterstico.
Discutimos abaixo cada caso especco.
5.7.1 Razes Reais e Distintas

Se p1 = p2 , ento a soluo geral da EDH ser
y(x) = C1 ep1 x + C2 ep2 x . (5.110)

5.7.2 Razes Reais e Iguais

Se b2 4c = 0, o polinmio caracterstico possui apenas uma soluo p1 =
p2 = p = b/2. Ento, uma soluo da EDH dada por y1 (x) = ebx/2 .
Substituindo em (5.22), obtemos y2 (x) = xebx/2 . Logo a soluo geral da
EDH
y(x) = (C1 + C2 x)ebx/2 . (5.111)
5.7.3 Razes Complexas

Sejam p = i duas razes complexas do polinmio caracterstico. Usando
a frmula de Euler, obtemos a soluo geral da EDH
y(x) = C1 e(+i)x + C2 e(i)x

= ex [B1 cos(x) + B2 sin(x)] , (5.112)
onde B1 = C1 + C2 e B2 = C1 iC2 .
Exemplo 5.7 Resolver a equao do oscilador harmnico simples
y (x) + 2 y(x) = 0 . (5.113)

Soluo: O presente caso corresponde ao de duas razes complexas
p = i . Usando = 0 e = acima, temos que
y(x) = B1 cos(x) + B2 sin(x) . (5.114)

5.8 Solues Numricas

Em geral, os problemas de fsica e engenharia no possuem solues exatas.
Se estivermos interessados em solues analticas, necessrio usar mto-
dos perturbativos . Em geral estes mtodos funcionam da seguinte forma.
Primeiramente identicamos o termo da equao diferencial o qual torna a
soluo no-exata; este o termo perturbativo. Em seguida, consideramos
a equao diferencial na ausncia deste termo, e determinamos as solues
exatas. A seguir, expandimos a soluo da equao diferencial original como
uma combinao linear das solues exatas no-perturbadas. Finalmente,
determinamos os coecientes da expanso. Em geral, somente os poucos
primeiros coecientes da expanso podem ser determinados, o que torna a
soluo muito restrita.
Devido grande complexidade dos clculos envolvidos nos mtodos de

perturbao, e devido as suas limitaes, procuramos por solues num-
ricas, as quais so exatas. Com o advento dos computadores digitais com
processadores de alta velocidade, nas ltimas dcadas cresceu enormemente
o interesse por mtodos de soluo numrica de equaes diferenciais. Os
fabricantes de software tambm tem desenvolvido rotinas que implementam
estes mtodos. O estudo dos mtodos numricos esto alm dos objetivos
deste livro. Nesta Seo, veremos como resolver uma equao diferencial or-
dinria usando o programa MATLAB.! Para isto, consideremos a seguinte
EDH:
xy + y + xy = 0 . (5.115)
Note que x = 0 um ponto singular regular. Do ponto de vista numrico,
isto exige cuidado especial para se evitar overow. Usemos as seguintes
denies: y1 = y , y2 = y , dy1 = y = y2 , e dy2 = y . Ento temos,
y
2
dy2 = + y1 . (5.116)
x
O cdigo a seguir implementa o clculo da funo e suas derivadas. Salve
este cdigo em alguma rea de seu computador com o nome EDO.m."
function dy = EDO (x,y)

y1 = y(1); % funo y(x)
y2 = y(2); % primeira derivada de y(x)
dy1 = y2;
%%%
dy2 = -(y2./x+y1); % calcula a segunda
% derivada de y(x)
%%%
dy = [dy1;dy2]; % constri o vetor dy
% cujas componentes so as derivadas
Nesta mesma rea, crie um script contendo os comandos listados abaixo

e o salve com o nome Bessel.m. A seguir, execute o script. Alternativamente,
podemos executar os comandos um de cada vez no prompt do MATLAB,
sem necessidade de criar um script.
Os programas Mathematica e Maple tambm podem ser boas alternativas para esta
!
tarefa.
" Os arquivos do MATLAB devem ter a extenso m. Os valores da funo e sua
primeira derivada sero armazenados nas componentes O e O do vetor O, respecti-
vamente. O texto a frente do smbolo % interpretado pelo MATLAB como comentrio.
Portanto ignorado na execuo do cdigo.
xmin = 0.000001; xmax = 10;

y10 = 1; % valor da funo em x=0
y20 = 0; % valor da derivada em x=0
[x,y] = ode45(@EDO, [xmin xmax], [y10; y20]);
% o comando ode45 do MATLAB chama a funo EDO
% e a resolve entre xmin e xmax
n = size(y,1); % o comando size determina
% o tamanho da matriz y
figure; % abre um frame
plot(x,y(1:n,1),'k'); % constri o grfico
% de y(x) versus x
xlabel('x','FontSize',14);
ylabel('y(x)','FontSize',14);
figure; % abre um novo frame
plot(x,y(1:n,2),'k'); % constri o grfico
% da derivada de y(x) versus x
xlabel('x','FontSize',14);
ylabel('y^{\prime}(x)','FontSize',14);
Os resultados da soluo numrica da EDH (5.115) so exibidos na
Fig. 5.1.
1 0.4
0.3
0.2
0.1
0.5
0
y (x)
y(x)
-0.1
-0.2
0
-0.3
-0.4
-0.5
-0.5 -0.6
0 1 2 3 4 5 6 7 8 9 10 0 1 2 3 4 5 6 7 8 9 10
x x
Figura 5.1:
5.9 Problemas
5.1 Considere a equao diferencial linear de segunda ordem
3 3
y y + 2 y = 2x 1 .
x x
Mostrar que y1 (x) = x soluo da EDH associada. Faa isto substi-

tuindo y1 e suas respectivas derivadas na EDH. Determine a segunda
soluo da EDH y2 (x). A partir desta soluo determinar a soluo
particular da EDNH.
5.2 Repita o exerccio anterior para as seguintes equaes diferenciais:
(a) xy (2x + 1)y + (x + 1)y = (x2 + x 1)e2x ,

(b) (x + 3)y (2x + 7)y + 2y = (x + 3)2 ex ,
(c) (x sin x + cos x)y x(cos x)y + (cos x)y = x ,
sendo que para o caso (a) y1 (x) = ex ; para o caso (b) y1 (x) = 2x + 7;
e para o caso (c) y1 (x) = x.
5.3 Resolver as seguintes equaes diferenciais usando o mtodo de Fro-
benius:
(a) y x2 y y = 0 ,
(b) y + (x 1)y + y = 0 .
5.4 Tomando-se c0 = 1/3, mostre que y1 (x) do Exemplo 5.5 torna-se

sin x cos x
j1 (x) = ,
x2 x
que conhecida com funo esfrica de Bessel# de primeira ordem.
Analogamente, se c0 = 1, mostre que y2 (x) torna-se
cos x sin x
n1 (x) = ,
x2 x
que conhecida como funo esfrica de Neumann de primeira ordem.
5.5 Resolver as seguintes equaes diferenciais usando o mtodo de Fro-
benius:
(a) 3xy + 2y + x2 y = 0 ,
(b) 2(x2 + x3 )y (x 3x2 )y + y = 0 ,
(c) xy + y + x2 y = 0 .
Em cada um dos casos acima determinar os pontos singulares da EDH

(se houver) e classica-los em regular ou irregular.
# Friedrich Wilhelm Bessel (1784 - 1846). Matemtico, astrnomo e geodesista alemo.
5.6 Considere a equao diferencial linear homognea de segunda ordem

de Legendre$
(1 x2 )y 2xy + n(n + 1)y = 0 .
Resolver esta equao usando o mtodo de Frobenius. Siga as seguin-

tes etapas:
Mostre que x = 0 um ponto ordinrio da EDH. Assim podemos
supor que a soluo da EDH tem a forma

y(x) = cj xj .
j=0
Mostrar que a relao de recorrncia dada por
(j n)(j + n + 1)
cj+2 = cj .
(j + 1)(j + 2)
Mostrar que a soluo ser dada por y(x) = c0 y1 (x) + c1 y2 (x) =
F1 (x) + F2 (x), onde
n(n + 1) 2 n(n 2)(n + 1)(n + 3) 4

y1 (x) = 1 x + x
2! 4!
(n 1)(n + 2) 3 (n 1)(n 3)(n + 2)(n + 4) 5
y2 (x) = x x + x
3! 5!
Ambas as sries convergem para todo |x| < 1, qualquer que seja n.
Supor n inteiro e positivo. Ento as solues sero polinomiais.
Encontrar c0 e c1 tal que F1 (1) = F2 (1) = 1. Nestas condies,
mostrar que as solues sero dadas por
P0 (x) = 1,
P1 (x) = x,
1
P2 (x) = (3x2 1) ,
2
1
P3 (x) = (5x3 3x) ,
2
..
.
$ Adrien-Marie Legendre (1752 - 1833). Matemtico francs com importantes contri-
buies para a estatstica, teoria dos nmeros, lgebra abstrata e anlise matemtica.
Estas funes so conhecidas como polinmios de Legendre (ver Pro-

blema 2.3).
Note que os polinmios de Legendre podem ser representados na forma
geral
l
[2]
1 (2l 2j)!
Pl (x) = l (1)j xl2j ,
2 j=0 j!(l 2j)!(l j)!
onde [ 2l ] indica o maior inteiro menor ou igual a 2l .
Vericar que esta frmula verdadeira para l = 0, 1, 2, 3.
Outra representao til dos polinmios de Legendre a chamada

frmula de Rodrigues
1 dl 2
Pl (x) = (x 1)l .
2l l! dxl
Vericar que esta frmula verdadeira para l = 0, 1, 2.
5.7 Considere a equao diferencial linear homognea de segunda ordem

de Bessel
x2 y + xy + (x2 n2 )y = 0 .
Resolver esta equao usando o mtodo de Frobenius. Siga as seguin-

tes etapas:
Mostrar que x = 0 um ponto singular regular da EDH. Desta

forma, a soluo da EDH pode ser escrita como

y(x) = cj xj+s .
j=0
Mostrar que a equao indicial tem duas razes dadas por s = n.

Considere o caso s = n. Mostrar que neste caso a relao de
recorrncia dada por
1
cj = cj2 .
j(j + 2n)
Mostrar que a soluo correspondente a s = n dada por

1 1
y1 (x) = xn
1 x2 + 4 x4
22 (n + 1) 2 2!(n + 1)(n + 2)

1 6
6 x + .
2 3!(n + 1)(n + 2)(n + 3)
Multiplicando e dividindo a equao acima por 2n n!, mostrar que

y1 (x) = 2n n!Jn (x), onde

1 x 2j+n
Jn (x) = (1)j ,
j=0
j!(j + n)! 2
a funo de Bessel de ordem n. Note que y2 (x) pode ser obtida a

partir de y1 (x) substituindo-se n por n. Se n for um inteiro posi-
tivo y2 (x) tem signicado? E se n for um inteiro negativo y1 (x) tem
signicado?
5.8 Resolver a equao diferencial linear homognea
y 3y + 2y = 0 .
Captulo 6
Variveis Complexas
6.1 Introduo
O surgimento dos nmeros complexos tem origem na diculdade em deter-
minar razes reais de determinadas equaes polinomiais. A mais simples
delas x2 + 1 = 0. Matemticos do sculo XVI , tais como Cardano e
Bombelli, deram os primeiros passos para encontrar uma soluo para esse
problema. Eles perceberam que era necessrio introduzirum nmero, o qual
chamamos de unidade imaginria e denido por i = 1.
6.2 Nmeros Complexos

Sejam a e b dois nmeros reais. Denimos um nmero complexo z como
sendo z = a + ib. Dado um nmero complexo z = a + ib, denotamos a
por parte real de z e b por parte imaginria de z: a = Re(z), b = Im(z).
Denotamos um nmero complexo de imaginrio puro se sua parte real for
nula, isto , z = ib.
Chamamos de conjugado de um nmero complexo z = a + ib, e o deno-
tamos por z, o nmero cuja parte real igual a Re(z) e a parte imaginria
igual ao oposto de Im(z). Assim, z = a ib.
6.2.1 Operaes com Nmeros Complexos
Adio e Subtrao
A soma (subtrao) de dois nmeros complexos z1 = a + ib e z2 = c + id,
obtida somando-se (subtraindo-se) as partes reais e imaginrias. Assim,
Captulo 6: Variveis Complexas 103
z1 z2 = (a c) + i(b d) . (6.1)
Multiplicao
Sejam dois nmeros complexos z1 = a + ib e z2 = c + id. Usando-se a
propriedade distributiva, e que i2 = 1, podemos escrever
z1 z2 = (ac bd) + i(ad + bc) . (6.2)

Ento, segue desta propriedade que z z = a2 + b2 .
Diviso
A diviso entre nmeros complexos pode ser obtida usando-se as proprie-
dades acima. Temos
z1 a + ib
=
z2 c + id
a + ib c id
=
c + id c id
ac + bd bc ad
= +i 2 . (6.3)
c2 + d2 c + d2
Denimos o mdulo
de um nmero complexo z = a + ib, e o denotamos
por |z|, como |z| = a2 + b2 . Note que z z = |z|2
6.2.2 Representao Grca de um Nmero Complexo

A partir de agora usaremos a terminologia varivel complexa
para enfatizar
que as partes real e imaginria do nmero complexo so variveis ao invs
de constantes. Assim, usaremos a notao z = x + iy .
Uma varivel complexa z = x+iy tambm pode ser representada por um
par ordenado (x, y). Desta forma, podemos representar z no plano xy como
ilustrado na Fig. 6.1, tambm conhecido como plano complexo
ou diagrama
de Argand.
Usando as relaes trigonomtricas em um tringulo retngulo, podemos
escrever
x = r cos , y = r sin , (6.4)

"
onde r = x2 + y 2 = |x + iy| e = Arg(z) = arctan(y/x) o argumento
de z . Decorre que
y (x, y)
r

x
Figura 6.1: Representao grca de um nmero complexo.
z = x + iy = r(cos + i sin ) . (6.5)

Esta conhecida como forma polar (ou trigonomtrica ) de um nmero
complexo, onde (r, ) so as coordenadas polares. Usando o resultado do
Problema 2.2, podemos escrever
z = rei , (6.6)
conhecida como frmula de Euler.
Seja z1 = x1 + iy1 = r1 (cos 1 + i sin 1 ) e z2 = x2 + iy2 = r2 (cos 2 +
i sin 2 ). Usando-se a frmula de Euler, podemos facilmente mostrar as
seguintes frmulas
z1 z2 = r1 r2 [cos(1 + 2 ) + i sin(1 + 2 )] ,
z1 r1
= [cos(1 2 ) + i sin(1 2 )] ,
z2 r2
n
zn = [r(cos + i sin )] = rn [cos(n) + i sin(n)] ,
1/n
z 1/n = [r(cos + i sin )]

1/n + 2k + 2k
= r cos + i sin ,
n n
k = 0, 1, 2, . . . , n 1 . (6.7)
A ltima destas equaes conhecida como frmula de Moivre.

Deixa-se como exerccio para o leitor a demonstrao dessas frmulas.
Note que o valor de z 1/n no unvoco, isto , pode assumir mltiplos
valores para cada valor de z .
Fica claro pela representao grca de um nmero complexo que este
pode ser identicado com um vetor no plano xy . Portanto, vale a desigual-
dade triangular
|z1 | |z2 | |z1 + z2 | |z1 | + |z2 | . (6.8)

A demonstrao simples,
|z1 + z2 |2 = (z1 + z2 )(z1 + z2 )

= |z1 |2 + |z2 |2 + 2Re(z1 z2 )
= |z1 |2 + |z2 |2 + 2|z1 ||z2 | cos(1 2 )
|z1 |2 + |z2 |2 + 2|z1 ||z2 |
= (|z|1 + |z|2 )2 . (6.9)
6.3 Funes Complexas

Se a cada valor de uma varivel complexa z corresponder um nico ou mais
valores de uma segunda varivel complexa w, ento dizemos que w uma
funo de uma varivel complexa z , isto w = f (z). Esta funo dita
uniforme se a cada valor de z corresponder um nico valor de w = f (z),
e multiforme ou plurvoca em caso contrrio. Por exemplo, f (z) = z 2 =
(x+iy) = x2 y 2 +2ixy . Em geral podemos escrever f (z) = u(x, y)+iv(x, y),
onde u e v so as partes real e imaginria de f (z).
6.4 Condies de Cauchy-Riemann

Analogamente ao caso de uma varivel real, deni-se a derivada de uma
funo complexa f (z) como segue
df (z) f (z + z) f (z) f (z)

= lim = lim . (6.10)
dz z0 z z0 z
Se o limite existe e independente do caminho usado para se aproximar
do ponto z , ento dizemos que f (z) analtica em z .
Considere os dois caminhos simples da Fig. 6.2. Aproximando-se do
ponto z0 pelo caminho horizontal, y = 0 e z = x. Temos
f (z) u + iv
lim = lim
z0 z x0 x
u v
= lim + i lim
x0 x x0 x
u v
= +i . (6.11)
x x
x 0 z0
y 0
x 0
y 0
Figura 6.2: Dois caminhos independentes aproximando-se do ponto z0 .
Analogamente, aproximando-se de z pelo caminho vertical, obtemos
f (z) u + iv
lim = lim
z0 z y0 iy
v u
= lim i lim
y0 y y0 y
v u
= i . (6.12)
y y
A condio para que f (z) seja analtica em z a de que os dois limites

sejam idnticos. Temos
u v v u
= , = . (6.13)
x y x y
Estas so conhecidas como equaes de Cauchy-Riemann. Elas so con-
dies necessrias para que f (z) seja analtica em um determinado domnio
D do plano complexo.
Exemplo 6.1 Mostrar que f (z) = z 2 analtica em todo plano
complexo xy .
Soluo: As partes real e imaginria so dadas por u(x, y) = x2 y2 e
v(x, y) = 2xy . Temos para as derivadas parciais,
u v
= 2x =
x y
v u
= 2y = . (6.14)
x y
Com efeito, f (z) = z 2 analtica em todo plano complexo, pois no

existem valores de x e/ou de y que restringem a validade das condies de
Cauchy-Riemann.
Exemplo 6.2 Mostrar que f (z) = z no analtica em qualquer
subdomnio do plano complexo xy .
Soluo: Neste caso, u(x, y) = x e v(x, y) = y. Temos
u v
= 1 = = 1
x y
v u
=0= , (6.15)
x y
para todo (x, y). Logo f (z) = z no analtica em qualquer ponto do

plano xy . Este um caso tpico de uma funo que contnua mas no
necessariamente analtica.
Exemplo 6.3 Mostrar que se f (z) = u(x, y) + iv(x, y) uma funo
analtica, ento u e v satisfazem a equao de Laplace.
Soluo: Derivando ambos os membros da primeira das equaes de
Cauchy-Riemann (6.13) com relao a x e a segunda com relao a y , ob-
temos
2u 2v 2v 2u
= , = 2 . (6.16)
x2 xy yx y
Ento, somando-se
2u 2u
+ 2 =0. (6.17)
x2 y
De maneira anloga podemos demonstrar que 2 v = 0.
6.5 Teorema Integral de Cauchy

Antes de enunciar e provar o teorema integral de Cauchy, precisamos denir
a integral no plano complexo. Seja f (z) uma funo uniforme e contnua
num domnio D. Seja uma curva C unindo dois pontos z0 = x0 + iy0 e
z0 = x0 + iy0 no interior de D. Denimos a integral complexa como sendo
Augustin Louis Cauchy (1789 - 1857) nasceu em Paris, Frana. Laplace e Langrange,
freqentadores da casa dos pais de Cauchy, logo cedo perceberam o interesse do jovem
pela matemtica.
(x0 ,y0 )
f (z) dz = (u + iv)(dx + idy)
C (x0 ,y0 )
(x0 ,y0 ) (x0 ,y0 )
= (udx vdy) + i (vdx + udy) . (6.18)
(x0 ,y0 ) (x0 ,y0 )
Conseqentemente, o clculo da integral complexa se reduz ao clculo

de integrais reais de linha.
Teorema 6.1 (Teorema Integral de Cauchy, domnio simplesmente conexo) Se-

ja f (z) uma funo analtica num domnio D limitado por um contorno
fechado C , inclusive. Ento

f (z) dz = 0. (6.19)
C
Demonstrao: A demonstrao do teorema pode ser feita via teorema

de Stokes visto no Captulo 3. Conforme visto acima

f (z) dz = (udx vdy) + i (vdx + udy) . (6.20)
C C C
Mostraremos separadamente que as duas integrais acima anulam-se iden-

ticamente. Seja o vetor bidimensional V = Vx i+Vy j. Pelo teorema de Stokes
(3.78), temos que
)
V dr = V n dS
C D

Vy Vx
(Vx dx + Vy dy) = dS . (6.21)
C D x y
Demonstremos o teorema em duas etapas. Primeiro empregamos a igual-

dade acima tomando-se Vx = u e Vy = v . Obtemos

v u
(udx vdy) = + =0, (6.22)
C D x y
onde na ltima passagem usamos a segunda das equaes de Cauchy-Rie-
mann (6.13). Segundo, usamos o teorema de Stokes com Vx = v e Vy = u,

u u
(vdx + udy) = =0, (6.23)
C D x y
onde agora usamos a segunda das equaes de Cauchy-Riemann (6.13). Isto

naliza a demonstrao do teorema.
Seja f (z) uma funo analtica num domnio D. Este chamado de sim-
plesmente conexo se e somente se f (z) analtica em qualquer subdomnio
de D. Do contrrio, D chamado de multiplamente conexo (ver Fig. 6.3).
O teorema integral de Cauchy s valido para domnios simplesmente co-
nexos. Veremos agora a verso deste teorema para domnios multiplamente
conexos.
y y
C
D D
x x
(a) (b)
Figura 6.3: (a) Domnio simplesmente conexo; (b) domnio multiplamente

conexo. As partes hachuradas da gura representam as regies onde a fun-
o analtica e em branco onde no analtica.
Teorema 6.2 (Teorema Integral de Cauchy, domnio multiplamente conexo)

Seja f (z) uma funo analtica limitada por dois contornos fechados C1
e C2 , com C2 C1 . Ento
) )
f (z) dz = f (z) dz . (6.24)
C1 C2
Demonstrao: Consideremos os caminhos C1 e C2 como ilustrado na

Fig. 6.4. Pelo teorema integral de Cauchy para domnios simplesmente
conexos, podemos escrever
)
f (z) dz = 0 . (6.25)
AP QA B T SRBA
Por outro lado, a integral pode ser subdividida em vrios caminhos,

Q
A B
C1 C2
A B
P
S R
Figura 6.4:
)
f (z) dz = f (z) dz + f (z) dz
AP QA B T SRBA AP QA A B

+ f (z) dz + f (z) dz .
B T SRB BA
(6.26)
Suponhamos que os caminhos A B e AB so paralelos e de mesmo

comprimento. Isto no uma limitao na validade do teorema, pois se a
distncia entre os dois segmentos se aproxima de zero, nossa aproximao
se torna exata. Usando propriedade de integrais de linha, temos que

f (z) dz = f (z) dz . (6.27)
A B BA
Ento,
)
f (z) dz = f (z) dz + f (z) dz
AP QA B T SRBA AP QA B T SRB
) )
= f (z) dz + f (z) dz
C2 C1
= 0, (6.28)
de onde se obtm
) )
f (z) dz = f (z) dz . (6.29)
C1 C2
6.6 Frmula Integral de Cauchy

Nesta Seo aplicaremos o teorema integral de Cauchy para a deduo da
importante frmula integral de Cauchy. Considere f (z) uma funo anal-
tica num domnio D limitado por um contorno fechado C . Temos
)
n! f (z)
f (n) (z0 ) = dz , (6.30)
2i C (z z0 )n+1
lembrando que

(n) dn f (z)
f (z0 ) = , (6.31)
dz n z=z0
com z0 no interior de D.
Note que f (z) analtica em D, mas f (z)/(z z0 )n+1 no .
Provaremos essa frmula para o caso n = 0. Temos
)
1 f (z)
f (z0 ) = dz . (6.32)
2i C (z z0 )
Seguindo os caminhos da Fig. 6.5, e usando o teorema integral de Cauchy,
podemos escrever
) )
f (z) f (z)
dz = dz . (6.33)
C (z z0 ) C (z z0 )
A equao do crculo centrado em z0 dada por z z0 = ri . Diferen-

ciando, com r xo, dz = irei d. Portanto,
) )
f (z) f (z0 + rei )
dz = irei d
C (z z0 ) C rei
)
= i f (z0 + rei ) d . (6.34)
C
Tomando-se o limite de r 0, provamos o teorema para n = 0,

C
z0
Figura 6.5:
) )
f (z)
dz = if (z0 ) d
C (z z0 ) C
= 2if (z0 ) . (6.35)
Demonstremos o teorema para o caso n = 1, partindo da expresso da

derivada de f (z),
f (z0 + z) f (z0 )
f (z0 ) = lim . (6.36)
z0z
Usando recursivamente o resultado do teorema para o caso n = 0, temos
que
)
1 f (w)
f (z0 + z) = dw
2i C (w z0 z)
)
1 f (w)
f (z0 ) = dw . (6.37)
2i C (w z0 )
Subtraindo uma equao da outra, obtemos
)
f (z0 + z0 ) f (z0 ) 1 1 1 1
= f (w) dw
z 2i z C w z0 z w z0
)
1 f (w)
= dw . (6.38)
2i C (w z0 z)(w z0 )
Levando ao limite de z 0,
)
1 f (w)
f (z0 ) = dw , (6.39)
2i C (w z0 )2
que equivalente a (6.30) para n = 1. A demonstrao pode ser feita para
n arbitrrio usando o mtodo de induo nita: admitimos que o teorema
vale para n e ento provamos que vale para n + 1. Isto deixado para o
leitor.
Exemplo 6.4 Mostrar que
)
1 1 1, n=1
dz = , (6.40)
2i C (z z0 )n 0 , n = 1
onde C um contorno fechado incluindo o ponto z0 .
Soluo: Reescrevemos a frmula integral de Cauchy (6.30) com n tro-
cado por n 1,
)
(n 1)! f (z)
f (n1) (z0 ) = dz . (6.41)
C (z z0 )
2i n
Tomando o caso particular f (z) = 1, a integral acima ser no nula se

e somente se n = 1, do contrrio todas as derivadas se anulam, pois f (z)
uma funo constante.
6.7 Srie de Taylor

Considere uma funo analtica e uniforme em um domnio D limitado por
um contorno fechado C . Seja um contorno fechado no interior de D
contendo os pontos z e z0 , e z sobre (ver Fig. 6.6). Ento, para qualquer
z satisfazendo |z z0 | < |z z0 |, temos a frmula integral de Cauchy
)
1 f (z )
f (z) = dz . (6.42)
2i C (z z)
Agora, reescrevemos o integrando como
1 1 1
= = , (6.43)
z z z z0 (z z0 )
(z z0 ) 1 zz0
z z0
e usamos a expanso binomial para escrever

z

z0 z
Figura 6.6:
n
1 z z0
= , (6.44)
1 zz0
z z0 n=0
z z0
a qual vlida para |z z0 | < |z z0 |. Juntando os termos, temos

1 (z z0 )n
= . (6.45)
z z n=0 (z z0 )n+1

Inserindo dentro da frmula integral de Cauchy,

)

1 f (z )
f (z) = (z z0 )n dz . (6.46)
2i C n=0 (z z0 )n+1
Agora, invertemos a ordem de integrao e soma,
)
1 f (z )
f (z) = z )n+1
(z z0 )n dz
n=0
2i C (z 0

f (n) (z0 )
= (z z0 )n , (6.47)
n=0
n!
onde na passagem da primeira para a segunda linha usamos a frmula in-

tegral de Cauchy. Essa a srie de Taylor para funes complexas.
6.8 Pontos Singulares e Plos

Seja f (z) uma funo analtica num domnio D, exceto num ponto z = z0 .
Denominamos z0 de ponto singular isolado de f (z).
Considere funes que tenham a forma
(z)
f (z) = , (6.48)
(z z0 )n
onde (z) analtica em qualquer ponto do domnio D, incluindo o ponto
z0 . Ento f (z) tem um ponto singular isolado em D, z = z0 . Denominamos
este ponto de plo de ordem n.
Se n = 1, denominamos a singularidade de plo simples, e se n = 2 de
plo duplo.
6.9 Srie de Laurent

Conforme vimos anteriormente, uma funo analtica nas vizinhanas de
um ponto z0 pode ser expandida em uma srie de Taylor em torno deste
ponto. A srie de Laurent uma generalizao da srie de Taylor para os
casos em que z0 um ponto singular isolado.
A srie de Laurent em torno do ponto z0 uma soma de duas sries
de potncias, uma consistindo de potncias positivas de z z0 e outra de
potncias negativas,

f (z) = aj (z z0 )j + aj (z z0 )j , (6.49)
j=0 j=1
a qual pode ser escrita abreviadamente como

f (z) = aj (z z0 )j . (6.50)
j=
Considere uma funo f (z) com um ponto singular isolado z0 , isto ,

pode ser escrita como em (6.48). Uma vez que (z) analtica em D,
podemos escrever

(z) = aj (z z0 )j , (6.51)
j=0
onde aj = (j) (z0 )/j!. Combinando as equaes, obtemos
a0 a1 an1
f (z) = + + + + an + an+1 (z z0 ) + .
(z z0 )n (z z0 )n1 z z0
(6.52)
Renomeando os coecientes, {a0 an , a1 an+1 , . . .,an1 a1 ,
an a0 , an+1 a1 , . . .}, encontramos
an an+1 a1
f (z) = + + + + a0 + a1 (z z0 ) +
(z z0 )n (z z0 )n1 z z0

= aj (z z0 )j . (6.53)
j=n
Tomando-se o limite n , obtemos a srie de Laurent. Assim, uma

funo que possui um ponto singular isolado, admite uma expanso em srie
de Laurent em torno deste ponto. Se a srie de Laurent em torno de z0 tem
innitos termos, chamamos z0 de singularidade essencial.
Exemplo 6.5 Expandir a funo f (z) = 1/(z z0 ) em uma srie de
Taylor.
Soluo: Para |z| < |z0 |, temos
1 1 zj
= = , (6.54)
z z0 z0 1 z z0j+1
z0
e para |z0 | < |z|,
1 1 zj
0
= = . (6.55)
z z0 z 1 z0
z
z j+1

Exemplo 6.6 Expandir a funo f (z) = 1/z(z 1) em uma srie de
Laurent em torno de z = 0.
Soluo:
1 1 1
= = 1 + z + z2 + z3
z(z 1) z(1 z) z
1
= 1 z z2 . (6.56)
z

Captulo 6: Variveis Complexas %
Exemplo 6.7 Expandir a funo f (z) = z/(z + 1)(z + 2) em uma

srie de Laurent em torno de z = 1 e z = 2.
Soluo: Para expandir em torno de z = 1 faamos a mudana de
varivel w = z + 1,
z w1 w1
= = 1 w + w2 w3 +
(z + 1)(z + 2) w(w + 1) w
1
= + 2 2w + 2w2 2w3 +
w
1
= + 2 2(z + 1) + 2(z + 1)2 2(z + 1)3 + .
z+1
(6.57)
Procedendo de forma anloga, a expanso em torno do ponto z = 2
dada por
z 2
= + 1 + (z + 2) + (z + 2)2 + (z + 2)3 + . (6.58)
(z + 1)(z + 2) z+2

6.10 Teorema dos Resduos

As Sees anteriores foram uma preparao para chegarmos parte mais
importante deste Captulo que o teorema dos resduos. Este teorema
muito importante para o clculo de integrais denidas que normalmente
no podem ser resolvidas usando os mtodos convencionais.
Primeiramente comeamos com a denio de resduos. Consideremos
o comportamento de uma funo analtica ao redor de um ponto singular
isolado z0 . Assim, a funo admite uma expanso de Laurent,
an an+1 a1
f (z) = + + + +a0 +a1 (z z0 )+ . (6.59)
(z z0 )n (z z0 )n1 z z0
Multiplicando os dois membros da equao por (z z0 )n , temos
(z z0 )n f (z) = an + an+1 (z z0 ) + + a1 (z z0 )n1

+a0 (z z0 )n + a1 (z z0 )n+1 + . (6.60)
Agora, derivando (n 1)-vezes,
dn1 n!
[(z z0 )n f (z)] = (n 1)!a1 + a0 (z z0 )
dz n1 1!
(n + 1)!
+ a0 (z z0 )2 + , (6.61)
2!
e levando ao limite z z0 , obtemos
1 dn1
a1 = lim [(z z0 )n f (z)] . (6.62)
zz0 (n 1)! dz n1
Este coeciente da expanso de Laurent conhecido como resduo de
f (z) no plo de ordem n, z = z0 .
Teorema 6.3 (Teorema dos Resduos) Considere uma funo analtica

f (z) em um domnio D limitado por um contorno fechado C , exceto em
um determinado nmero de pontos singulares isolados zk , k = 1, 2, . . . , N .
Ento
)
N
f (z) dz = 2i Res(f, zk ) , (6.63)
C k=1
onde Res(f, zk ) o resduo de f (z) em z = zk , no sentido da equao (6.62).
Demonstrao: Faamos a demonstrao do teorema para o caso de

um nico ponto singular isolado. Integrando a expanso de Laurent (6.53)
no contorno fechado C , vem que
) ) )
1 1
f (z) dz = an dz + a n+1 dz +
C (z z0 ) C (z z0 )
n n1
C
) ) )
1
+a1 dz + a0 dz + a1 (z z0 ) dz + .
C z z0 C C
(6.64)
+
Usando a equao (6.40), todas as integrais se anulam, exceto C (z
z0 )1 dz . Obtemos
)
f (z) dz = 2ia1 . (6.65)
C
Se o domnio D conter mais de um ponto singular isolado, ento

)
f (z) dz = 2i(a1 + b1 + c1 + )
C

N
= 2i Res(f, zk ) . (6.66)
k=1
Exemplo 6.8 Calcular

)
z2
dz , (6.67)
C (z 2)(z 2 + 1)
onde C dado pelo crculo |z| = 3.
Soluo: As razes do denominador do integrando so dadas por z = 2
e z = i. Fatorando, encontramos

1 1 1 1 1
= = . (6.68)
z2 + 1 (z + i)(z i) 2i z i z + i
Logo, todos os plos so simples. Da denio de resduo (6.62),
z2 4
a1 = lim (z 2) = ,
z2 (z 2)(z 2 + 1) 5
z2 1 2i
b1 = lim (z i) = ,
zi (z 2)(z 2 + 1) 10
z2 1 + 2i
c1 = lim (z + i) = . (6.69)
zi (z + 1)(z 2 + 1) 10
Portanto, pelo teorema dos resduos (6.63), temos que
)
z2 4 1 2i 1 + 2i
dz = 2i + + = 2i . (6.70)
C (z 2)(z 2 + 1) 5 10 10

)
1
dz , (6.71)
C z(z + 2)3
Soluo: O plo z = 0 simples e z = 2 de ordem trs. Tomamos
somente o resduo de z = 0, pois somente este plo est no interior do
domnio em considerao. Temos
1 1
a1 = lim z = . (6.72)
z0 z(z + 2)3 8
Ento,
)
1 1 i
3
dz = 2i = . (6.73)
C z(z + 2) 8 4

)
1
dz , (6.74)
C (z 1)(z + 3)2
Soluo: O plo z = 1 simples e z = 3 duplo. Os dois esto no
interior do domnio. Temos
1 1
a1 = lim (z 1) 2
= , (6.75)
z1 (z 1)(z + 3) 16
e

d 2 1
b1 = lim (z + 3)
z3 dz (z 1)(z + 3)2
1
= lim
z3 (z 1)2
1
= . (6.76)
16
Para a integral, temos que
)
1 1 1
2
dz = 2i =0. (6.77)
C (z 1)(z + 3) 16 16

6.11 Aplicao: Integrais Trigonomtricas

A primeira aplicao do teorema dos resduos que consideraremos o clculo
de integrais trigonomtricas do tipo
2
R(cos , sin ) d , (6.78)
0
onde R uma funo racional de seus argumentos e nita no intervalo

0 2 . O nosso objetivo transformar a integral em uma integral
complexa tal que possamos aplicar o teorema dos resduos. O crculo de
raio unitrio dado por z = ei . Lembrando que

ei + ei 1 1
cos = = z+ ,
2 2 z

ei ei 1 1
sin = = z . (6.79)
2i 2i z
Tambm temos, dz = izd, de onde se obtm d = iz .

dz
Juntando todos os
termos,
)
2
1 z + z1 z z1
R(cos , sin ) d = R , dz , (6.80)
0 C iz 2 2i
onde C o crculo de raio unitrio centrado na origem, e portanto a integral
pode ser calculada como segue
2
N
R(cos , sin ) d = 2 Res(f, zk ) , (6.81)
0 k=1
onde f (z) a funo racional

1 z + z1 z z1
f (z) = R , . (6.82)
z 2 2i
Exemplo 6.11 Calcular a integral
2
1
d , (6.83)
0 1 2p cos + p2
onde p real e p = 1.
Soluo: A restrio em p garante que o denominador do integrando
sempre no nulo, logo podemos ir adiante. A funo racional dada por
(6.82)
1 1 1
f (z) = 1 = 2 , (6.84)
z 1 2p z+ z
+ p2 pz (1 + p2 )z + p
2
e portanto tem dois plos z = 1/p e z = p. Fatorando, temos que

1
pz 2 (1 + p2 )z + p = p(z p)(z ) . (6.85)
p
Se 0 p < 1, somente z = p est no interior do domnio em questo; se
p > 1, ento somente z = 1/p est dentro do domnio. No primeiro caso,
temos que

1 1
Res(f, p) = lim (z p)
zp p (z p)(z p1 )
1
= , (6.86)
1 p2
e no segundo,

1 1 1 1
Res(f, ) = lim (z )
p z1/p p p (z p)(z p1 )
1
= . (6.87)
1 p2
Em qualquer dos casos podemos escrever o resultado de forma compacta
como 1/|1 p2 |, e para a integral temos que
2
1 2
2
d = . (6.88)
0 1 2p cos + p |1 p2 |

6.12 Aplicao: Integrais Imprprias

Seja f (x) uma funo contnua em 0 x . A integral imprpria de
f (x) denida por
R
f (x) dx = lim f (x) dx . (6.89)
0 R 0
No uso do teorema dos resduos, em geral mais conveniente considerar
a integral imprpria no intervalo simtrico com relao a origem,
R
f (x) dx = lim f (x) dx , (6.90)
R R
desde que f (x) tambm seja contnua em x 0.

6.12.1 Integrais Imprprias de Funes Racionais

Ilustraremos o clculo de integrais imprprias de funes racionais atravs
de um exemplo.

1
2 + a2
dx , (6.91)
0 x
onde a real e a = 0.
Soluo: Como uma tentativa consideremos a integral complexa equi-
valente substituindo a varivel real x pela varivel complexa z ,
)
1
2 + a2
dz , (6.92)
C z
onde C o semicrculo da Fig. 6.7. Os plos do integrando so z = ia;
fatorando, z 2 + a2 = (z ia)(z + ia), logo os plos so simples. Somente
z = ia est dentro do semicrculo, e pelo teorema dos resduos, obtemos
)
1 1
dz = 2i lim (z ia) = . (6.93)
C z2 +a 2 zia (z ia)(z + ia) a
+ia R CR
R +R x
ia
Figura 6.7:
Por outro lado, a integral complexa pode ser separada em duas outras
integrais,
) R
1 1 1
dz = dx + dz , (6.94)
C z 2 + a2 R x2 + a2 CR z 2 + a2
de onde se obtm
R
1
dx + IR = , (6.95)
R x2 +a 2 a
onde

1
IR = dz . (6.96)
CR + a2 z2
Agora, provaremos que IR 0, quando R . Usando a propriedade de
integrais complexas, podemos escrever

f (z) dz |f (z)||dz|

M |dz| = M L

(6.97)
onde M = max|f (z)| o mximo valor de |f (z)| e L o comprimento do
caminho . Pela desigualdade triangular (6.8), obtemos
|z 2 + a2 | |z 2 | a2 = R2 a2 , (6.98)
e portanto
1 1
2 . (6.99)
|z 2 + a2 | R a2
Combinando os resultados, conclumos que
R
|IR | . (6.100)
R2 a2
Tomando o limite de raio innito, vem que
R
lim |IR | = lim =0. (6.101)
R R R 2 a2
Finalmente obtemos

1
dx = , (6.102)
x2 +a2 a
ou ainda

1
dx = . (6.103)
0 x2 +a 2 2a

6.12.2 Integrais Imprprias de Funes Racionais e Tri-

gonomtricas
Igualmente ao caso anterior, ilustraremos o clculo de integrais imprprias
de funes racionais e trigonomtricas atravs de um exemplo.
Consideraremos integrais do tipo

cos(ax)
f (x) dx . (6.104)
sin(ax)
A estratgia de clculo destas integrais substituir a varivel real x pela
complexa z em f (x), mas no nas funes trigonomtricas, pois o seno e o
cosseno de uma varivel complexa no so funes limitadas. Do contrrio
a integral ir divergir. Alternativamente, podemos considerar a integral
complexa
)
f (z)eiaz dz , (6.105)
C
onde C dado pela Fig. 6.7.

cos(bx)
2 2
dx . (6.106)
x + a
Soluo: Neste caso a integral complexa

)
eibz
dz . (6.107)
C z 2 + a2
Analogamente ao Exemplo 6.12, temos que
)
eibz eibz
dz = 2i lim (z ia) = eab . (6.108)
C z2+a 2 zia (z ia)(z + ia) a
Agora,
) R
eibz eibx
dz = dx + IR , (6.109)
C z 2 + a2 R x2 + a2
onde

eibz
IR = dz . (6.110)
CR z 2 + a2
Na Seo 6.14 provaremos que IR 0, quando R . Logo,

eibx
dx = eab . (6.111)
x2+ a2 a
O uso da frmula de Euler eibx = cos(bx) + i sin(bx), resulta em

cos(bx) sin(bx)
2 2
dx + i 2 2
dx = eab . (6.112)
x + a x + a a
A segunda integral do membro esquerdo anula-se, pois o integrando mpar,
de onde se obtm

cos(bx)
2 2
dx = eab . (6.113)
x + a a
Alternativamente, podemos escrever

cos(bx) ab
dx = e . (6.114)
0 x2 + a2 2a

6.13 Aplicao: Integrais Imprprias com P-

los
Existem casos onde a funo possui plos localizados ao longo da reta real.
Este tipo de cenrio exige um tratamento especial que exemplicamos a
seguir.

sin
dx . (6.115)
0 x
Soluo: Note que x = 0 um ponto singular isolado. No entanto, a
singularidade removvel, pois no limite x 0 o integrando nito. Ento
o integrando contnuo para todo x. Seguimos a mesma idia do Exemplo
6.13 partindo da integral
)
eiz
dz . (6.116)
C z
No entanto, agora o caminho a ser seguido o da Fig. 6.8, pois z = 0
est no caminho de integrao, e como no podemos integrar passando pela
singularidade, fazemos um desvio em z = 0 por um semicrculo de raio r.
R CR
Cr
r
R r +r +R x
Figura 6.8:
O integrando no tem plos no interior do domnio, e pelo teorema

integral de Cauchy para domnios simplesmente conexos, temos que
r R
eiz eix eiz eiz eix
dz = dx + dz + dz = 0 . dx +
C z R x Cr r z CR z x
(6.117)
A integral real ser obtida nos limites R e r 0. Veremos na
Seo 6.14 que a integral no semicrculo maior se anula. Consideremos a
integral no semicrculo menor; Cr descrito por z = rei , com 0 , e
dz = irei d. Temos
0 i
eiz ere
dz = irei d , (6.118)
Cr z rei
e levando ao limite r 0,
0 i
eiz ere
lim dz = lim irei d
r0 Cr z r0 rei
0 i
ere
= lim irei d
r0 rei
0
= i d

= i . (6.119)
Reunindo os resultados,

eix
dx i = 0 , (6.120)
x
& Fsica-Matemtica: Teoria e Aplicaes
e usando a frmula de Euler, obtemos

cos x sin x
x
dx + i
x
dx = i . (6.121)

A primeira integral se anula pois o integrando uma funo mpar; o inte-

grando da segunda integral uma funo par, ento

sin x
x
dx = .
2
(6.122)
0

6.14 Lema de Jordan

Seja a integral

IR = f (z)eiaz dz , (6.123)
CR
onde CR o semicrculo de raio R da Fig. 6.7. Suponhamos que |f (z)|

seja limitado por um mximo, isto , |f (z)| Rc , onde c uma constante
k
positiva e k 1. Usando (6.97), temos que

|IR |
= f (z)e dz
iaz
CR
|f (z)||eiaz | |dz|
CR

c iaz
k
|e | |dz| . (6.124)
CR R
Substituindo z = Rei = R cos + iR sin e dz = iRei d,

c iR cos R sin
|IR | k
|e ||e | |iRei d|
0 R

c R sin
= k
e Rd
0 R

cR
=
Rk 0
eR sin d . (6.125)
Para encontrar um limite superior para a integral acima, primeiramente
notemos que sin = sin( ) para 0 ,
/2
eR sin d = 2 eR sin d . (6.126)
0 0
Agora, notamos que sin 2/ para todo 0 /2 (ver Fig. 6.9).
Ento

cR
|IR | eR sin d
Rk 0

2cR /2 2R/
e d
Rk 0
c
= k
1 eR . (6.127)
R
1
sin

2
Figura 6.9:
Tomando o limite de R , e supondo k 1, obtemos

c
lim |IR | = lim k
1 eR = 0 . (6.128)
R R R
As funes dos Exemplos 6.13 e 6.14 satisfazem as condies do lema
de Jordan pois, |1/(z 2 +a2 )| c/R2 , onde c = 1/(1a2 /R2 ) para o primeiro
caso, e |1/z| = 1/R para o segundo caso.
6.15 Problemas
6.1 Usando a primeira das relaes de Euler (1.2), mostrar que

1 p cos x
pn cos(nx) = .
n=0
1 2p cos x + p2
6.2 Determinar as funes analticas w(z) = u(x, y) + iv(x, y) para os

casos (a) u(x, y) = x3 3xy 2 , e (b) v(x, y) = ey sin x.
6.3 Considere a funo
1
f (z) = z + .
z
Determine o domnio em que f (z) analtica.
6.4 Mostre que a funo f (z) = x2 + iy no analtica para z arbitrrio.
* 2+i
Por conseguinte, espera-se que a integral 0 f (z) dz dependa do
caminho. Ilustre isso, calculando a integral sobre trs caminhos suaves
escolhidos por voc e ligando os pontos z1 = 0 e z2 = 2 + i.
6.5 Se f (z) = u(x, y) + iv(x, y) analtica no domnio D, provar que as
famlias de curvas de um parmetro u(x, y) = C1 e v(x, y) = C2 so
famlias ortogonais.
Sugesto: usando as condies de Cauchy-Riemann, mostrar que o
produto da inclinao das famlias de curvas -1.
6.6 Usando a frmula integral de Cauchy, mostrar que
)
2i , n = 1
(z z0 ) dz =
n
C 0 , n = 1 ,
onde C um contorno fechado circundando o ponto z = z0 no sentido
anti-horrio.
6.7 Usando a frmula integral de Cauchy, mostrar que
)
1
z mn1 dz = mn .
2i C
onde C um contorno fechado circundando a origem no sentido anti-
horrio.
6.8 Usando o teorema de Cauchy, a frmula integral de Cauchy, ou suas
conseqncias, calcule as seguintes integrais, todas tomadas sobre o
crculo |z| = 2:
+ cos z + sin z + ez
(a) z dz , (b) z dz , (c) z1 dz ,
+ 2z 2 +3z1
+ 2z
+ sin z
(d) z1+i dz , (e) z 2 9 dz , (f ) z2 dz ,
+ ez
+ sin z
(g) (z1)2 dz , (h) zn dz .
6.9 Desenvolva as seguintes funes em srie de Taylor ou de Laurent

vlidas nas vizinhanas dos pontos indicados. Determine as regies
em que os desenvolvimentos so vlidos.
1
(a) f (z) = cos z (z = 0) , (b) f (z) = (z 2 +1) (z = 0) ,
z 1
(c) f (z) = z 2 1 (z = 2) , (d) f (z) = z 2 +1 (z = 1) ,
ez
(e) f (z) = z(z 2 +1) (z = 0) , (f ) f (z) = cosh z (z = i) .
6.10 Determinar os resduos nos respectivos plos das seguintes funes
2z + 3 z3 ezt z
(a) , (b) , (c) , (d) 2 .
z2 4 z 3 + 5z 2 (z 2)3 (z + 1)2
+ z2
6.11 Calcular C (z+1)(z+3) dz , onde C uma curva fechada, envolvendo
todos os plos.
6.12 Se C uma curva fechada envolvendo z = i, mostrar que

)
1 zezt 1
dz = t sin t .
2i C (z 2 + 1)2 2
6.13 Usando o teorema dos resduos, calcule as seguintes integrais sobre os

contornos indicados:
)
1
(a) dz ,
C z2 + 4
onde C o crculo |z 2i| = 1;

)
cosh z
(b) dz ,
C z(z 2 + 1)
onde C o crculo |z| = 2;

)
z1
(c) dz ,
C z2 + iz + 2
onde C um caminho fechado envolvendo todos os plos.

6.14 Usando o teorema dos resduos, calcule as seguintes integrais:

* 2
(a) d
= 2 (|| < 1) ,
0 1+ cos 1 2
* 2 cos 3
(b) 0 54 cos d = 12 ,
* d
(c) = ,
0 1+sin2 2
* 2 2a
(d) 0
d
(a+b sin )2 = (a2 b2 )3/2
(a > |b|) ,
* 2 135(n1)
(e) 0
cosn d = 246n 2 .
6.15 Escolhendo apropriadamente o caminho de integrao, use o teorema

dos resduos para calcular as seguintes integrais:
* + dx
(a) x4 +a4
= a3 2
,
* x2
(b) 0 x6 +1 dx = 6 ,
* + dx 3
(c) (x2 +1)3
= 8 ,
* + cos kx kb
(d) (xa)2 +b2
dx = be cos ka (k, b > 0) ,
* x sin x
(e) 0 x2 +1 dx = 2e ,
* +
cos x eb ea
(f ) (x2 +a2 )(x2 +b2 ) dx =
a2 b2 b a (a = b) .
Captulo 7
Sries de Fourier
7.1 Introduo
J vimos no Captulo 2 que uma funo contnua com derivadas contnuas
de todas as ordens pode ter uma representao em forma de srie de Taylor
em torno de um determinado ponto. Neste Captulo, veremos um outro
tipo de srie conhecida como srie de Fourier. Se uma funo for peridica,
esta pode ter uma representao em forma de uma srie de Fourier. A srie
de Fourier tem vrias aplicaes. Para ilustrar sua utilidade iremos aplica-
la para o clculo de somas numricas. Mas a srie de Fourier tem outras
utilidades tais como na soluo de problemas de valores de contorno (ver
Captulo 8).
7.2 Funes Peridicas

Uma funo f (x) dita peridica de perodo T se, para todo x, satiszer
f (x) = f (x + T ) , (7.1)
onde T uma constante.
As Figs. 7.1 e 7.2 so exemplos de funes peridicas. No caso da funo
de onda quadrada temos que
Jean Baptiste Joseph Fourier (1768 - 1830) nasceu em Auxerre, Frana. rfo aos
8 anos, Fourier foi colocado no Colgio Militar, dirigido pelos beneditinos. Ao estudar a
conduo de calor foi levado a criar um novo tipo de desenvolvimento em srie, diferente
da srie de Taylor, empregando funes peridicas em vez de potncias, e que recebeu
seu nome.

1 , 1 < x < 0 ,
f (x) = (7.2)
1, 0<x<1,
onde T = 2. No caso da funo dente-de-serra temos
f (x) = |x| , 2 x 2 , (7.3)

onde T = 4.
f (x)
+1
3 2 1 +1 +2 +3 x
1
Figura 7.1: Funo de onda quadrada.
f (x)
6 4 2 +2 +4 +6 x
Figura 7.2: Funo dente-de-serra.
7.3 Funes Ortogonais

Considere um conjunto de funes {m (x) , m = 0, 1, 2, . . .}. Este con-
junto dito ortogonal se obedecer a propriedade
b
m (x)n (x) dx = mn , (7.4)
a
Captulo 7: Sries de Fourier 135
onde as funes so consideradas normalizadas. O conjunto de funes que

segue esta propriedade chamado de linearmente independente. Assim,
qualquer funo f (x) no intervalo a x b pode ser escrito como uma
combinao linear das 's. Temos
+

f (x) = cm m (x) , (7.5)
m=
Esta conhecida como srie de Fourier generalizada. Os coecientes

podem ser obtidos multiplicando-se ambos os lados da equao (7.5) por
n (x) e integrar o resultado no intervalo a x b. Obtemos
b
cm = f (x)m (x) dx . (7.6)
a
0 Exemplo 7.1
Mostre que1 o conjunto de funes
1
m (x) = L cos L x , m = 0, 1, 2, . . . , com L x L, ortogo-
m
nal.
Soluo: Denotemos por Imn a integral
L m n
1
Imn = cos x cos x dx . (7.7)
L L
L L
Usando a identidade trigonomtrica
1
cos A cos B = [cos(A B) + cos(A + B)] , (7.8)
2
com A e B escolhidos convenientemente, encontramos
L
1 (m n) (m + n)
Imn = cos x + cos x dx
2L L L L

1 L (m n)
= sin x
2L (m n) L
L
L (m + n)
+ sin x . (7.9)
(m + n) L L
Se m = n, ento a integral identicamente nula.

Considere agora o caso m = n,
L
1 2m
Imm = 1 + cos x dx
2L L L
L
1 L (m + n)
= x+ sin x
2L 2m L L
= 1. (7.10)
Ento Imn = mn .
Podemos
0 mostrar de forma anloga1 que o conjunto de

funes m (x) = 1L sin m
L x , m = 1, 2, . . . , com L x L, tam-
bm ortogonal.
7.4 Srie de Fourier

Seja uma funo f (x) denida no intervalo L x L e fora deste por
f (x + 2L) = f (x), isto , o perodo da funo T = 2L. A srie de Fourier
de f (x) denida por
m m
a0
f (x) = + am cos x + bm sin x . (7.11)
2 m=1
L L
Os coecientes da srie podem ser obtidos

multiplicando-se
m ambos os
lados da equao (7.11) por L1 cos mL x ou 1
L sin L x , conforme conve-
niente, e integrando o resultado de L a L. Em seguida, usamos a ortogo-
nalidade destas funes conforme visto no Exemplo 7.1. Obtemos
L m
1
am = f (x) cos x dx , m = 0, 1, 2, 3, . . . ,
L L L
L m
1
bm = f (x) sin x dx , m = 1, 2, 3, . . . . (7.12)
L L L
A srie (7.11) um caso particular da srie de Fourier generalizada (7.5).

Embora importante, deixaremos de lado os aspectos de convergncia da srie
de Fourier.
7.5 Identidade de Parserval

Supor que a srie de Fourier correspondente a uma funo f (x) converge
para esta funo no intervalo L x L. Desta forma, f (x) pode ser
representada pela srie (7.11). Multiplicando ambos os lados desta equao
por L1 f (x), e integrando o resultado, temos que
L L
1 2 a0 1
f (x) dx = f (x) dx
L L 2 L L

a0

L m
1
+ am x dx
f (x) cos

L L L
m=1

am

1
L m
+bm f (x) sin x dx
L L

L

bm

a20
= + (a2m + b2m ) . (7.13)
2 m=1
A equao

1 L
a20
f 2 (x) dx = + (a2m + b2m ) , (7.14)
L L 2 m=1
conhecida como identidade de Parserval.
7.6 Funes Pares e mpares

A srie de Fourier das funes que possuem esta propriedade pode ser simpli-
cada consideravelmente. Antes de prosseguirmos, vejamos uma importante
propriedade das funes pares e mpares quando integradas num intervalo
simtrico com relao origem.
Se f (x) uma funo par, isto , f (x) = f (x), ento
L 0 L
f (x) dx = f (x) dx + f (x) dx
L L 0
0 L
= f (x) dx + f (x) dx
L 0
L L
= f (x) dx + f (x) dx
0 0
L
= 2 f (x) dx . (7.15)
0
Por outro lado, se f (x) uma funo mpar, isto , f (x) = f (x),
ento
L 0 L
f (x) dx = f (x) dx + f (x) dx
L L 0
0 L
= f (x) dx + f (x) dx
L 0
L L
= f (x) dx + f (x) dx
0 0
= 0. (7.16)
m
Se f (x) uma funo par, ento f (x) cos m
L x par e f (x) sin L x
mpar. Logo, de (7.12) e (7.11), conclumos que
L m
1
am = f (x) cos x dx , m = 0, 1, 2, 3, . . . , (7.17)
L L L
bm = 0, (7.18)
a0

m
f (x) = + am cos x . (7.19)
2 m=1
L

Se f(x) uma funo mpar, ento f (x) cos m L x mpar e
f (x) sin m
L x par. Logo, de (7.12) e (7.11), conclumos que
am = 0, (7.20)
L m
1
bm = f (x) sin x dx , m = 1, 2, 3, . . . , (7.21)
L L L
m
f (x) = bm sin x . (7.22)
m=1
L
Exemplo 7.2 Determinar a srie de Fourier de f (x) = x2 , para

x .
Soluo: Neste caso a funo par. Ento, somente os am 's so no
nulos. Temos

2
am = x2 cos(mx) dx
0

1 1 2 2x
= x2 2 sin(mx) + 2 cos(mx)
m m m 0
m
(1)
= 4 , (7.23)
m2
para todo m = 0. Se m = 0,

2
2 2
a0 = x2 dx = . (7.24)
0 3
Inserindo os resultados na equao (7.19) encontramos

2 (1)m
x2 = +4 cos(mx) , (7.25)
3 m=1
m2
para todo x .
fcil mostrar que, tomando-se x = na equao acima, obtemos o
valor da soma numrica

1 2
= . (7.26)
m=1
m2 6

Exemplo 7.3 Achar a srie de Fourier da funo dente-de-serra da
Fig. 7.2.
Soluo: A funo par, com efeito bm = 0, e
2 mx
am = x cos dx
0 2
mx mx 2
2 4
= x sin + 2 2 cos
m 2 m 2 0
4
= [(1) 1] ,
m
(7.27)
m2 2
para todo m = 0. Se m = 0, ento
2
a0 = x dx = 2 . (7.28)
0
Substituindo os coecientes na srie de Fourier (7.19), temos que

4 [(1)m 1] mx
|x| = 1+ cos
2 m=1 m2 2

8 1 (2m + 1)x
= 1 2 cos , (7.29)
m=0 (2m + 1)2 2
para todo 2 x 2.
Exemplo 7.4 Usando a identidade de Parserval e o Exemplo 7.3,
mostrar que

1 4
= . (7.30)
m=1
m4 90
Note que este resultado foi usado na soluo do Problema 2.10.
Soluo: Pela identidade de Parserval (7.14),
2 2
1 22 16 [(1)m 1]
x2 dx = + ,
2 2 2 m=1
m4 4

8 64 1
= 2+ ,
3 4 m=0 (2m + 1)4

1 4
= . (7.31)
m=0
(2m + 1)4 96
Seja agora a soma

1
S= . (7.32)
m=1
m4
Temos

1 1
S = +
m=0
(2m + 1)4 m=1 (2m)4

4 1 1
= + 4
96 2 m=1 m4
4 S
= + , (7.33)
96 16
de onde se obtm o valor desejado S = 4 /90.
7.7 Problemas
7.1 Seja a funo

0 , t 0 ,
f (t) =
sin t , 0 t .
Esta a sada de uma meia-onda reticadora simples. tambm uma

aproximao do efeito trmico solar que produz mars na atmosfera.
Mostrar que o desenvolvimento em srie de Fourier de f (t) dado por
1 1 2 cos t
f (t) = + sin t .
2 n=2,4,6,...
n2 1
7.2 A funo dente-de-serra (ou onda triangular) dada por

x , x 0 ,
f (x) =
x, 0x.
Mostrar que a srie de Fourier de f (x) dada por
4 cos nx
f (x) = .
2 n=1,3,5,...
n2
7.3 Usando a srie de Fourier da onda triangular do Problema 7.2, mostrar

que

1 2
= .
n=1
(2n 1)2 8
7.4 Usando
f (x) = x4 , x ,
mostrar que

1 4
4
= ,
n=1
n 90
e que

(1)n+1 7 4
= .
n=1
n4 720
7.5 Mostre que a srie de Fourier de f (x) = | sin x| no intervalo x

, dada por

2 4 cos(2x) cos(4x) cos(6x)
| sin x| = + 2 + 2 + .
22 1 4 1 6 1
Usando este resultado, mostrar que o nmero pode ser calculado

usando a srie

1 1 1
=1+2 + .
2 22 1 42 1 62 1
7.6 Desenvolvendo f (x) = cosh(ax) em srie de Fourier, mostre que

sinh(a) 2a sinh(a) (1)n
cosh(ax) = + cos(nx) ,
a n=1
n 2 + a2
para todo x . Usando este resultado, mostrar que

1 cosh(a) 1
= 2 ,
n=1
n2 + a2 2a sinh(a) 2a
para todo a = 0.
7.7 Mostre que
nx
2 (1)m b2 m na mx
sinh = sinh sin ,
b m=1 b2 m2 + a2 n2 b a
para todo x .
Captulo 8
Problemas de Valores de
Contorno
8.1 Equao de Difuso do Calor
8.1.1 Condies de Contorno Homogneas
Neste Captulo no discutiremos a teoria geral das equaes diferenciais
parciais. Ao contrrio, aprenderemos como resolve-las atravs de alguns
exemplos relacionados com alguns fenmenos fsicos importantes. Comece-
mos pela equao de difuso do calor a qual dada por
2u
2
x2
=
u
t
, (8.1)
onde u a temperatura ao longo de uma barra metlica de comprimento 2L.
Suporemos que a face lateral da barra esteja isolada, de tal modo que no
haja uxo de calor ao longo da direo radial, mas sim somente ao longo
da direo paralela ao eixo da barra. A constante est relacionada com a
condutividade trmica, a densidade especca, e o calor especco da barra
(ver Fig. 8.1).
Suponhamos que a distribuio inicial de temperatura seja conhecida,
u(x, 0) = f (x), L x L . (8.2)
Esta equao tambm conhecida como condio inicial.
Suponhamos tambm que as extremidades da barra sejam mantidas s
temperaturas xas T1 e T2 . Estas so conhecidas como condies de con-
torno. Sem perda de generalidade, podemos considerar as condies homo-
x=0
x
2L
x = L x = +L
Figura 8.1: Barra metlica mantida s temperaturas T1 e T2 xas nas ex-

tremidades esquerda e direita, respectivamente.
gneas T1 = T2 = 0. Como veremos mais adiante, a soluo da equao

(8.1) com condies de contorno no-homogneas mapeada em cima da
soluo para o caso homogneo. Assim, temos que
u(L, t) = 0, t > 0 . (8.3)
Este problema pode ser considerado como um problema de valor de

contorno no plano xt (ver Fig. 8.2).
u(L, t) = 0 u(L, t) = 0
x
L u(x, 0) = f (x) +L
Figura 8.2:
Para resolver a equao de difuso do calor usamos o mtodo de sepa-

rao de variveis. Temos
u(x, t) = X(x)T (t) . (8.4)

Substituindo em (8.1), obtemos
Captulo 8: Problemas de Valores de Contorno 145

X 1 T
= 2 , (8.5)
X T
onde a linha denota a derivada com relao a x e o ponto a derivada com
relao a t. Uma vez que (x, t) so duas variveis independentes, os dois
lados da equao acima devem ser uma constante. Temos
X
= 2 ,
X

1 T
= 2 . (8.6)
2 T
Podemos escrever as equaes acima como

X + 2 X = 0 , 2 2
T + T = 0 . (8.7)
A constante de separao foi escolhida convenientemente como 2 , de
tal modo que a equao para X parea-se com a equao de um oscilador
harmnico simples. Consideremos o primeiro caso = 0. Temos, X = 0
cuja soluo dada por X(x) = k1 x + k2 . As constantes de integrao k1
e k2 podem ser obtidas usando-se as condies de contorno (8.3), ou seja,
X(L) = 0. Temos
k1 L + k2 = 0 ,
k1 L + k2 = 0 , (8.8)
cuja soluo a trivial k1 = k2 = 0.

Considere agora o caso = 0. A soluo geral para a parte espacial
dada por
X(x) = k1 cos(x) + k2 sin(x) . (8.9)

Impondo as condies de contorno, vem que
k1 cos(L) + k2 sin(L) = 0 ,
k1 cos(L) k2 sin(L) = 0 . (8.10)
Em princpio deveramos ter escolhido uma constante de separao arbitrria po-
sitiva, negativa ou nula. No difcil mostrar que no existe soluo para X(x) que
satisfaa as condies de contorno para o caso < 0. Logo, consideramos somente os
casos = 2 0.
O sistema de equaes acima possui soluo no-trivial se e somente se

cos(L) sin(L)

cos(L) sin(L) = 2 cos(L) sin(L) = 0 , (8.11)

de onde se obtm L = n ou L = n + 12 , com n = 1, 2, 3, . . .. Valores
negativos de n no geram solues novas, pois estas diferem apenas por
um sinal das solues com n positivo. Consideremos o caso = n/L.
Substituindo na primeira das equaes (8.10), com sin(L) = 0 e cos(L) =
(1)n , obtemos k1 = 0. Logo,
nx
X(x) = k2 sin . (8.12)
L
Igualmente, usando a expresso de na segunda das equaes (8.7),
obtemos
n 2
2
T + T =0, (8.13)
L
cuja soluo dada por
2
(n/L)2 t
T (t) = Ae . (8.14)
Combinando as solues (8.12) e (8.14), com k2 = A = 1, sem perda de
generalidade, encontramos a soluo fundamental
2 2
nx
un (x, t) = e (n/L) t sin , n = 1, 2, 3, . . . . (8.15)
L
Determinamos a soluo geral usando o princpio de superposio, ou
seja, expressando a soluo geral como um combinao linear das solues
fundamentais,

u(x, t) = bn un (x, t)
n=1

2
nx
(n/L)2 t
= bn e sin . (8.16)
n=1
L
As constantes bn 's so determinadas usando-se as condies de contorno.

Temos

nx
u(x, 0) = f (x) = bn sin . (8.17)
n=1
L
Usando a teoria das sries de Fourier do Captulo 7, obtemos

L n
1
bn = f (x) sin x dx . (8.18)
L L L
A srie para u(x, t) converge devido ao fator exponencial negativo. Em
geral podemos armar que
lim u(x, t) = 0 . (8.19)

t
8.1.2 Condies de Contorno No-Homogneas

Consideremos agora as condies de contorno sicamente corretas, pois em
caso contrrio, no poderia haver troca de calor entre as extremidades da
barra metlica. Temos
u(L, t) = T1 , u(L, t) = T2 , t > 0 . (8.20)

Suponhamos que no estado estacionrio a distribuio de temperatura
seja dada por
v(x) = lim u(x, t) . (8.21)

t
Ento, pela equao de difuso do calor (8.1), a soluo estacionria

satisfaz
2 v (x) = 0 , L x L , (8.22)
cuja soluo geral dada por v(x) = k1 x + k2 . Das condies de contorno
(8.20), vemos que v(x) satisfaz as seguintes condies
v(L) = T1 , v(L) = T2 , (8.23)

de onde obtemos k1 = (T2 T1 )/2L e k2 = (T1 + T2 )/2. Logo,
T2 T1 T 1 + T2
v(x) = x+ . (8.24)
2L 2
A soluo geral pode ser obtida superpondo a soluo estacionria (8.24)
com uma distribuio transiente w(x, t). Temos
u(x, t) = v(x) + w(x, t) , (8.25)

onde w(x, t) uma funo a ser determinada. Substituindo esta equao na
equao de difuso do calor (8.1), e usando (8.22), encontramos
2w w
2
2
= . (8.26)
x t
As condies de contorno para a soluo transiente so homogneas,
pois, de (8.20) e (8.23) fcil vericar que
w(L, t) = u(L, t) v(L) = 0 . (8.27)

Ento, a soluo da equao (8.26) dada por (8.16), exceto pelo fato
de que bn agora determinado pela condio inicial
w(x, 0) = u(x, 0) v(x) = f (x) v(x) , (8.28)

onde f (x) uma funo dada. Portanto, a soluo geral da equao de
difuso do calor satisfazendo as condies de contorno ser dada por
T2 T1 T 1 + T2
nx
u(x, t) = x+ + bn sin , (8.29)
2L 2 n=1
L
onde agora
n
1 L
T 2 T1 T 1 + T2
bn = f (x) x sin x dx . (8.30)
L L 2L 2 L
8.2 Equao da Onda

Como um segundo exemplo, consideremos a equao da onda
2u 2u
= a2 , (8.31)
x2 t2
sujeitas s condies de contorno conforme ilustradas na Fig. 8.3.
Note que agora temos duas condies inicias, pois a equao diferencial
parcial de segunda ordem no tempo.
Igualmente ao caso anterior, usamos o mtodo de separao de variveis
escrevendo
u(x, t) = X(x)T (t) . (8.32)

Substituindo esta equao em (8.31) obtemos duas equaes diferenciais
ordinrias independentes,

X + 2 X = 0 , 2 2
T +a T = 0 . (8.33)
u(L, t) = 0 u(L, t) = 0
x
L u(x, 0) = f (x) +L

u(x,t)
t = g(x)
t=0
Figura 8.3:
A soluo para a primeira das equaes (8.33) a mesma que para o

caso da equao de difuso do calor. Temos
nx
X(x) = k2 sin . (8.34)
L
As funes cos(at) e sin(at) so duas solues linearmente indepen-
dentes da segunda das equaes (8.33). Portanto, temos duas solues fun-
damentais,
nx
nat
un (x, t) = sin cos ,
L L
nx
nat
vn (x, t) = sin sin . (8.35)
L L
Usando o princpio de superposio obtemos a soluo geral de (8.31),

nx
nat

nat

un (x, t) = sin bn cos + an sin . (8.36)
n=1
L L L
Para completar a soluo, falta ainda determinar as constantes bn e cn .

Para isto usamos as condies iniciais. Temos

nx
u(x, 0) = f (x) = bn sin ,
n=1
L
na nx
u(x, t)
= g(x) = cn sin . (8.37)
t t=0 n=1
L L
Por meio da teoria de Fourier, podemos determinar as constantes. Temos
L n
1
bn = f (x) sin x dx ,
L L L
na L n
1
cn = g(x) sin x dx . (8.38)
L L L L
8.3 Problemas de Dirichlet e de Neumann

Suponhamos um determinado problema fsico, podendo ou no haver depen-
dncia temporal, descrito por alguma equao diferencial parcial de segunda
ordem nas coordenadas espaciais. O problema de Dirichlet aquele em que
as condies sobre a funo incgnita so conhecidas em todos os contornos
da geometria em considerao. Por outro lado, o problema de Neumann!
aquele em que os valores das derivadas normais da funo incgnita so
conhecidas nos contornos.
Como um exemplo, discutiremos a soluo da equao de Laplace em
coordenares retangulares e polares.
8.3.1 Equao de Laplace: Coordenadas Retangulares

Seja a equao de Laplace no plano xy ,
2u 2u
+ 2 =0, (8.39)
x2 y
sujeita s condies de contorno de Dirichlet num retngulo conforme ilus-
trado na Fig. 8.4.
Johann Peter Gustav Lejeune Dirichlet (1805 - 1859), nascido em Dren, Alemanha.
Notabilizou-se por importantes contribuies teoria dos nmeros e pelo melhoramento
na demonstrao no teorema de Fourier.
! John Von Neumann (1903 - 1957), matemtico hngaro de origem judaica. Participou
da construo do ENIAC (Electrical Numerical Integrator and Calculator ), o primeiro
computador digital eletrnico da era moderna.
+b u(x, b) = 0
u(a, y) = 0
+a
x
a
u(a, y) = f (y)
u(x, b) = 0 b
Figura 8.4:
Empregando o mtodo de separao de variveis, escrevemos
u(x, y) = X(x)Y (y) . (8.40)

Substituindo na equao de Laplace (8.39), obtemos duas equaes di-
ferenciais ordinrias e independentes,
X 2 X = 0 , Y + 2 Y = 0 . (8.41)
As condies de contorno da Fig. 8.4 implicam em Y (b) = 0. Assim,
similarmente aos casos anteriores, encontramos a seguinte soluo para a
segunda das equaes (8.41)
ny
Y (y) = k2 sin . (8.42)
b
As funes ex so duas solues linearmente independentes da primei-
ra das equaes (8.41). Portanto
X(x) = j1 ex + j2 ex . (8.43)
Atravs da condio de contorno X(a) = 0, facilmente obtida da
Fig. 8.4, obtemos j1 = j2 e2a . Substituindo j1 de volta ltima equao
acima, encontramos
X(x) = j2 sinh[(x + a)] , (8.44)

onde j2 = j2 ea .
Sem perda de generalidade, podemos tomar k2 = j2 = 1. Conseqente-

mente, obtemos a seguinte soluo fundamental
ny
n(x + a)
un (x, y) = sinh sin . (8.45)
b b
Recorrendo ao princpio de superposio, obtemos nalmente a soluo
geral para a equao de Laplace,

ny
n(x + a)
u(x, y) = cn sinh sin . (8.46)
n=1
b b
Para determinar as constantes cn 's, usamos a condio de contorno

u(a, y) = f (y). Encontramos
L ny
2na 1
sinh cn = f (y) sin dy . (8.47)
b L L L
8.3.2 Equao de Laplace: Coordenadas Polares

Considere agora a equao de Laplace em coordenadas polares no plano,
2 u 1 u 1 2u
+ + 2 2 =0. (8.48)
r 2 r r r
Suponhamos o mesmo problema de Dirichlet, porm com as condio de
contorno sobre um crculo conforme ilustrado na Fig. 8.5.
Alm da condio de contorno da gura, tambm iremos impor que
u(0, ) nita, e que u(r, ) seja peridica de perodo 2 .
Como usual, empregamos o mtodo de separao de variveis
u(r, ) = R(r)() . (8.49)
Introduzindo esta equao em (8.48) obtemos duas equaes diferenciais

ordinrias e independentes,
+ = 0 , r2 R + rR R = 0 , (8.50)
onde uma constante de separao.

Resolveremos a equao de Laplace em trs etapas para o domnio r<
a. Primeiramente supor = 2 < 0. A soluo geral da primeira das
equaes (8.50) ser dada por
() = j1 e + j2 e . (8.51)
u(a, ) = f ()
Figura 8.5:
A funo () assim escrita no satisfaz a condio de periodicidade, a

menos que j1 = j2 = 0. Logo, se for negativa, obtemos a soluo trivial,
a qual no sicamente signicativa.
Em segundo, supor que = 0. Com efeito, = 0, cuja soluo
geral = j1 + j2 . Impondo a condio de periodicidade ( + 2) =
(), conclumos que j2 = 0; de tal maneira que () uma constante.
Com = 0, de (8.50) obtemos r2 R + rR = 0. Esta equao pode ser
facilmente integrada e obtermos R = k1 + k2 ln r. Supondo R nita em
r = 0, devemos ter k2 = 0. Conseqentemente, u(r, ) uma constante, e a
soluo fundamental ser
u0 (r, ) = 1 . (8.52)
Finalmente, consideremos o caso = 2 > 0. A soluo geral das

equaes (8.50) ser dada por
= j1 cos() + j2 sin() , R = k1 r + k2 r . (8.53)
A condio de periodicidade para a soluo da parte angular nos fornece

cos() {j1 [cos(2) 1] + j2 sin(2)}

+ sin() {j2 [cos(2) 1] j1 sin(2)} = 0 . (8.54)
Usando o fato que cos() e sin() so duas funes linearmente in-

dependentes, seus coecientes na equao acima devem se anular identica-
mente. Temos
j1 [cos(2) 1] + j2 sin(2) = 0 ,
j2 [cos(2) 1] j1 sin(2) = 0 . (8.55)
A soluo no-trivial do sistema de duas equaes e duas incgnitas

existe se o determinante da matriz dos coecientes anular-se. Temos

[cos(2) 1] sin(2)
=0, (8.56)
sin(2) [cos(2) 1]
cuja soluo dada por cos(2) = 1. Ou seja, = n = 0 um nmero
inteiro no nulo. Sem perda de generalidade podemos considerar somente
positivo. Da equao (8.53) podemos observar que a soluo da parte radial
tem uma singularidade em r = 0. No entanto, supondo R nita neste ponto,
devemos ter k1 = 0. Assim, obtemos as seguintes solues fundamentais
un (r, ) = rn cos(n) ,
vn (r, ) = rn sin(n) . (8.57)
Finalmente, a soluo geral da equao de Laplace pode ser obtida

tomando-se uma combinao linear das solues fundamentais. De (8.52) e
(8.57), obtemos

c0 n
u(r, ) = + r [cn cos(n) + kn sin(n)] . (8.58)
2 n=1
Por ltimo, devemos determinar os coecientes da srie. Usando a con-

dio de contorno u(a, ) = f () e a teoria de Fourier, encontramos

n 1
a cn = f () cos(n) d , n = 0, 1, 2, 3, . . . ,

1
an kn = f () sin(n) d , n = 1, 2, 3, . . . . (8.59)

A soluo para o domnio r > a pode ser determinada de forma anloga.

Encontramos

c0 n
u(r, ) = + r [cn cos(n) + kn sin(n)] , (8.60)
2 n=1
onde

1
an cn = f () cos(n) d , n = 0, 1, 2, 3, . . . ,

1
an kn = f () sin(n) d , n = 1, 2, 3, . . . . (8.61)

u(x,t) u(x,t)
x =0 x =0
x=L x=L
x
L u(x, 0) = f (x) +L
Figura 8.6:
8.4 Problemas
8.1 Considere uma barra com as extremidades isoladas. Neste caso no
h passagem de calor pelas extremidades. As condies de contorno
esto ilustradas na Fig. 8.6. Resolver a equao de difuso do calor
(8.1) empregando as condies de contorno especicadas no digrama
xt.
8.2 Determine o deslocamento u(x, t) de uma corda elstica de compri-

mento 2L, xa em ambos os extremos, que posta em movimento a
partir de sua posio de equilbrio reta com velocidade inicial g de-
nida por
+b u(x, b) = g(x)
u(a, y) = 0
+a x
a
u(a, y) = 0
u(x, b) = 0 b
Figura 8.7:
A
L (x + L) , L x L/2 ,
L x , L/2 x L/2 ,
A
g(x) =
A
L (L x) , L/2 x L ,
onde A uma constante. Suponha

agora que a condio inicial fosse
substituda por g(x) = A2 sin xL , L x L. Voc esperaria
diferena signicativa entre os dois resultados para o deslocamento?
Em caso armativo, em quais pontos do domnio em questo?
8.3 Considere o problema de Dirichlet para um retngulo. Resolva a equa-

o de Laplace usando as condies de contorno ilustradas na Fig. 8.7.
8.4 Seja o problema de Dirichlet em um crculo de raio a. Usando os

resultados deduzidos na Seo 8.3.2, mostre que a soluo da equao
de Laplace tambm pode ser escrita na forma

1
r 2 a2
u(r, ) = f () d ,
2 r2 2ar cos( ) + a2
onde r < a.
Sugesto : Use o resultado do Problema 6.1.
Captulo 9
Funo Delta de Dirac
9.1 Introduo
Antes de iniciarmos o estudo da funo delta de Dirac, recordemos o conceito
de funo. Seja X um conjunto de elementos de nmeros reais simbolizados
pela varivel x. Seja tambm um segundo conjunto de nmeros reais Y .
Uma funo f (x) uma aplicao ou transformao de X em Y . A cada
valor de x em X , associamos um valor f (x) em Y ; f (x) pode associar
um valor de X a um ou mais valores de Y . O conjunto X chamado de
domnio da funo e Y de sua imagem. Este conceito pressupe que ambos
os conjuntos sejam formados por nmeros nitos.
A funo delta de Dirac denida por

x = 0 ,
(x) =
0,
, x=0.
(9.1)
Note que esta funo no se enquadra dentro do conceito de funo
estabelecido anteriormente, pois singular no ponto x = 0. Por esta razo,
a funo delta de Dirac freqentemente chamada de funo imprpria.
Conforme veremos mais adiante, essa funo no deve ser entendida no
sentido usual de funo, mas sim como uma distribuio. Isto signica
dizer que a funo delta deve ser tal que

(x) dx = 1 . (9.2)

A funo delta foi primeiramente proposta por Paul Dirac no incio do

Este conceito pode facilmente ser extendido para o conjunto de nmeros complexos.
sculo passado. Ele inventou uma engenhosa forma de lidar com innitos,
o que inicialmente causou muito desconforto aos matemticos. Porm, mais
tarde, Schwartz! desenvolveu uma nova teoria matemtica, chamada de
teoria das distribuies, onde se provou que a proposio de Dirac estava
correta e a dava consistncia.
Existem vrias representaes da funo delta que obedecem a proprie-
dade fundamental (9.2). Na realidade, a funo delta de Dirac uma funo
extremamente pontiaguda em torno de um determinado ponto da reta (ver
Fig. 9.1). Pulsos so exemplos tpicos de funes delta, como por exemplo
uma fora aplicada a um corpo em um tempo innitesimalmente pequeno.
y
(x)
Figura 9.1: Esboo da funo delta de Dirac.
9.2 Representaes da Funo Delta

Seja a funo denida como segue
Paul Adrien Maurice Dirac (1902 - 1984) nasceu em Bristol, Inglaterra. Filho de
pai suo com me britnica. Seu pai, de natureza extremamente autoritria, lhe en-
sinou francs e o encaminhou para as cincias exatas. Paul Dirac foi responsvel por
importantes contribuies Mecnica Quntica, como a descoberta da antimatria.
! Laurent-Mose Schwartz (1915 - 2002), importante matemtico francs do sculo XX,
por muito tempo teve de esconder seu verdadeiro nome para evitar a perseguio nazista
aos judeus durante a segunda guerra mundial.
Captulo 9: Funo Delta de Dirac 159

0 , |x| > 1/n ,
n (x) = (9.3)
2 , |x| < 1/n ,
n
para n = 1, 2, 3, . . .. Na forma imprpria podemos escrever
(x) = lim n (x) . (9.4)

n
Agora, note que esta representao da funo delta satisfaz a propriedade
(9.2), pois
1/n
n
lim n (x) dx = lim dx
n n 1/n 2
n 1/n
= lim x
n 2 1/n
= 1. (9.5)
Existem inmeras outras representaes da funo delta. As mais nor-
malmente usadas so
1 n
(x) = lim , (9.6)
1 + n2 x2
n
1 sin(nx)
(x) = lim , (9.7)
n x
n 2 2
(x) = lim en x . (9.8)
n
A seqncia {n (x), n = 1, 2, 3, . . .} conhecida como seqncia delta. A
seqncia (9.3) ilustrada na Fig. 9.2. Note que quanto maior o valor de
n, maior a intensidade do pico da funo n (x) em torno do ponto x = 0
e menor a sua largura. Usando o teorema dos resduos (ver Seo 6.10),
podemos mostrar que

1 n
dx = 1 , (9.9)
1 + n2 x2

1 sin(nx)
dx = 1 . (9.10)
x
Uma forma equivalente representao (9.7) e de grande utilidade
+
1
(x) = eikx dk , (9.11)
2
(ver Problema 9.4).
x
1 21 51 1 1
5 2 1
Figura 9.2: Seqncia delta da equao (9.3).
9.3 Propriedades da Funo Delta

Seja f (x) uma funo contnua nas vizinhanas de um ponto da reta x = a.
A propriedade fundamental da funo delta

(x a)f (x) dx = f (a) . (9.12)

Assim, a funo delta atua como um ltro, ou seja, ela seleciona apenas o
valor f (a) de f (x). Na Fig. 9.3 ilustramos geometricamente esta importante
propriedade.
Esta propriedade pode ser facilmente demonstrada usando-se o teorema
da mdia do clculo integral. Temos

(x a)f (x) dx = lim n (x a)f (x) dx
n
(x)
f (a)
f (x)
x
a
Figura 9.3:
a+1/n
n
= lim f (x) dx
a1/n 2
n
n 2 1 1
= lim f () , a a +
n 2 n n n
= lim f ()
n
= f (a) . (9.13)
Outras propriedades importantes da funo delta so:
(x) = (x) , (9.14)

1
(ax)f (x) dx = f (0) , (9.15)
|a|

1
(x2 a2 )f (x) dx = [f (a) + f (a)] , a > 0 . (9.16)
2a
As duas primeiras propriedades podem ser provadas sem diculdade.
Provemos a terceira,

2 2
(x a )f (x) dx = ((x a)(x + a))f (x) dx

a+
= ((x a)(x + a))f (x) dx
a
a
+ ((x a)(x + a))f (x) dx
a+
a+
= (2a(x + a))f (x) dx
a
a
+ (2a(x a))f (x) dx . (9.17)
a+
Agora, usando as propriedades (9.15) e (9.12) em seqncia, obtemos o

resultado desejado.
Essa propriedade pode ser generalizada para situaes em que o argu-
mento da funo delta uma funo. Seja y(x) uma funo que possua
razes simples nos pontos {xi , i = 1, 2, . . . , N }. Suponhamos que as deriva-
das de y(x) existam nos pontos xi . Ento

N
f (x )
(y(x))f (x) dx = i . (9.18)
dy
i=1
dx x=xi
9.4 Funo de Heaviside

A funo de Heaviside," tambm conhecida como funo escada, denida
por

1 , x > x ,
(x x ) = (9.19)
0 , x < x .
O esboo do grco desta funo encontra-se na Fig. 9.4.

Existe uma importante relao entre as funes delta e de Heaviside.
Considere a integral
" Oliver Heaviside (1850 - 1925), engenheiro eletricista, matemtico e fsico britnico,
nasceu em Camden Town, cidade nos arredores de Londres. Dentre muitas contribuies
importantes para a cincia, cou famoso por simplicar e tornar til as equaes de
Maxwell.
(x x )
x
x
Figura 9.4: Funo de Heaviside.

d d
f (x ) (x x ) dx

= f (x )(x x ) dx
dx dx

x
d
= f (x ) (x x ) dx
dx

=1

f (x ) (x x ) dx
x

=0
x
d
= f (x ) dx
dx
= f (x)

= f (x )(x x) dx , (9.20)

de onde se obtm
d
(x x) = (x x ) . (9.21)
dx
9.5 Funo Delta em duas e trs Dimenses

A funo delta (9.1) pode ser facilmente generalizada para dimenses supe-
riores. Em trs dimenses, temos que

0, r = 0 ,
(r) = (9.22)
, r=0,
onde r um vetor de coordenadas (x, y, z). Similarmente ao caso unidimen-
sional (9.2), temos que

(r) dx dy dz = 1 , (9.23)

de onde se obtm, em coordenadas retangulares,
(r) = (x)(y)(z) . (9.24)

No caso bidimensional, o grco da funo delta se parece com a super-
fcie de um sino, conforme ilustrado na Fig. 9.5.
40
30
(x,y)
20
10
0
0.2
0.1 0.2
0 0.1
-0.1 0
-0.1
-0.2 -0.2
y x
Figura 9.5: Visualizao da funo delta em duas dimenses; no grco

acima usamos a representao (x, y) = limn n2 exp(n2 (x2 + y 2 ))/ .
A funo delta pode facilmente ser escrita em outros sistemas de co-

ordenadas por meio da teoria de coordenadas curvilneas do Captulo 4.
Por exemplo, em coordenadas cilndricas, usando a equao (4.30) para o
elemento de volume, temos que
2
(r r ) dx dy dz = (r r ) d d dz .
0 0
(9.25)
A integral acima s poder ser identicamente 1, se e somente se,

1
(r r ) = ( )( )(z z ) . (9.26)

9.6 Relao de Fechamento

A funo delta tem uma estreita relao com um conjunto de funes orto-
normais. Enunciaremos esta relao em forma de um teorema.
Teorema 9.1 Seja {n (x)}, sendo n um inteiro, um conjunto de funes
ortonormais. Isto ,
b
m (x)n (x) dx = mn . (9.27)
a
Ento, uma representao da funo delta em termos destas funes ser

(x x ) = n (x)n (x ) . (9.28)
n
Demonstrao: Uma vez que {n (x)} um conjunto de funes line-

armente independentes, uma funo arbitrria poder ser escrita como uma
combinao linear destas. Temos

f (x) = cn n (x) . (9.29)
n
Multiplicando os dois membros da equao por m (x) e integrando o
resultado, vem que
b b
f (x)m (x) dx = cn m (x)n (x) dx
a n a

= cn mn
n
= cm . (9.30)
Substituindo cn de volta em (9.29), obtemos
, -
b
f (x) = f (x )n (x ) dx n (x)
n a
, -
b

= f (x ) n (x)n (x ) dx . (9.31)
a n
Usando a propriedade fundamental (9.12), vemos que a igualdade acima

s pode ser satisfeita, se e somente se, a funo delta for escrita como em
(9.28). Abreviadamente, dizemos que {n } forma um conjunto completo de
funes.
9.7 Problemas
9.1 Mostre que uma distribuio de cargas pontuais pode ser escrita em
termos de uma funo delta:

N
(r) = qi (r ri ) ,
i=1
onde ri a posio da i-sima carga qi .

Sugesto : Mostre que a integral da densidade de carga resulta na
carga total.
9.2 Considere as seguintes seqncias delta:
1 n 1 sin(nx) n 2 2
n (x) = 2 2
, n (x) = , n (x) = en x .
1+n x x
Desenhe os grcos destas funes para n = 1, 2, 3. Verique que,
quanto maior for o valor de n, mais estas funes se localizam em
torno de x = 0.
9.3 Considere a representao (9.8) da funo delta. Mostre que os pontos,
onde os o valores da funo correspondem
metade do mximo em
x = 0, so dados por x = ln 2/n. A distncia entre estes dois
pontos a largura do pico. Note que quando n , a largura
tende a zero e o mximo ao innito, porm a integral da funo em
toda reta vale um. Faa um esboo desta funo demonstrando estas
caractersticas.
9.4 Mostre que a seqncia delta
1 sin(nx)
n (x) = ,
x
tambm pode ser escrita como
n
1
n (x) = eikx dk .
2 n
Sugesto : Calcule a integral explicitamente.

9.5 Usando integrao por partes, mostre que

f (x) (x) dx = f (0) .

9.6 Tomado (0) = 1/2, mostre gracamente que
1
(x) = lim ,
k 1 + e2kx
uma boa aproximao para a funo de Heaviside.
9.7 Mostre que a representao da funo delta em coordenadas esfricas

dada por
1
(r r ) = (r r )(cos cos )( ) .
r2

9.8 Mostre que {n (x) = einx/L / 2L , n = 0, 1, 2, . . .} forma um
conjunto completo de funes ortonormais no intervalo L x L.
E, conseqentemente,

1 in(xx )/L
(x x ) = e .
2L n=
Captulo 10
Transformadas de Fourier
10.1 Teorema de Fourier
O teorema de Fourier pode ser enunciado da seguinte forma.
Teorema 10.1 Seja f (x) uma funo tal que satisfaa as seguintes condi-
es:
1. A srie de Fourier para f (x) converge,

*
2. A integral
|f (x)| dx < ,
ento temos as relaes

1
f (x) = F (q)eiqx dq ,
2

F (q) = f (x)eiqx dx . (10.1)

Demonstrao: Faremos a demonstrao do teorema de forma heurs-

tica. Primeiramente lembramos que a srie de Fourier dada por
m m
f (x) =
a0
2
+ am cos
L
x + bm sin
L
x . (10.2)
m=1
Usando as relaes de Euler (1.2), temos que

Captulo 10: Transformadas de Fourier 169

, imx - , imx -
a0 imx
e L + e L e L e L
imx
f (x) = + am + bm
2 m=1
2 2i

1
a0 1 imx 1 imx
= + am + bm e L + am bm e L .
2 m=1
2 i i
(10.3)
Substituindo m por m na segunda soma, podemos reescrever a equao

acima como
1
a0 1 mx 1 imx
f (x) = + (am ibm ) e L + (am + ibm ) e L .
2 m=1
2 m=
2
(10.4)
A equao acima pode ser escrita abreviadamente na forma

1
f (x) = Fm eimx/L , (10.5)
2 m=
onde

am ibm , m1,
Fm = a0 , m=0, (10.6)

am + ibm , m 1 .
Se f (x) uma funo conhecida, podemos determinar os coecientes Fm
em termos de uma integral envolvendo esta funo. Com este objetivo, pri-
meiramente multiplicamos ambos os lados da equao (10.5) por einx/L /L
e integramos o resultado entre L e L. Temos
L
L
1 1 1
f (x)einx/L dx = Fm ei(nm)x/L dx ,
L L 2 m=
L L
L
1
ei(nm)x/L dx = mn . (10.7)
L L
Logo, os coecientes so dados por

L
1
Fm = f (x)eimx/L dx . (10.8)
L L
Agora, substituindo os coecientes de volta na srie de Fourier dada por

(10.5), encontramos

, -
1 1 L
imx /L
f (x) = f (x )e dx eimx/L
2 m=
L L

1 L
1
im(x x)
= f (x ) e dx . (10.9)
2 L L m=
O ltimo passo da demonstrao no rigorosa, o que a torna heurstica.

Tomemos o limite L , de tal forma que q = L
0, e q = m
L torna-se
uma varivel contnua. Assim, temos que

1
f (x) = f (x ) eiq(x x) dq dx
2

1
= F (q)eiqx dx , (10.10)
2
onde

F (q) = f (x)eiqx dx . (10.11)

mais freqente nos livros textos encontrarmos o par de transformadas

de Fourier na forma

1
f (x) = F (q)eiqx dq ,
2

1
F (q) = f (x)eiqx dx , (10.12)
2
com as quais lidaremos daqui em diante.
Chamamos F (q) de transformada de Fourier direta de f (x), e de f (x) a
transformada de Fourier inversa de F (q). Simbolicamente temos
F (q) = F {f (x)} , f (x) = F 1 {F (q)} . (10.13)

As integrais (10.1) existem se satiszerem a segunda condio do teorema
enunciado acima. Para ver isto, basta notar que

f (x)e iqx
dx |f (x)e iqx
| dx = |f (x)| dx < . (10.14)

10.2 Teorema de Convoluo

Teorema 10.2 Se F (q) e G(q) so as transformadas de Fourier de f (x) e
g(x) respectivamente, ento

F (q)G(q)eiqx dx = f (y)g(x y) dy

= f (x y)g(y) dy

= 2f g , (10.15)
onde

1
f g = f (x y)g(y) dy , (10.16)
2
usualmente conhecido como produto de convoluo. Simbolicamente po-

demos escrever F 1 {F (q)G(q)} = f g , isto , F (q)G(q) = F {f g}.
Demonstrao: Este teorema pode ser demonstrado em poucas linhas

1 iq(xy)
f (y)g(x y) dy = f (y) G(q)e dq dy
2

1
= G(q)eiqx iqy
f (y)e dy dq
2

= F (q)G(q)eiqx dx . (10.17)

10.3 Identidade de Parserval

Se F (q) e G(q) so as transformadas de Fourier de f (x) e g(x), respectiva-
mente, ento

dq =
F (q)G(q) f (x)
g (x) dx , (10.18)

onde as barras denotam o complexo conjugado. Usando a representao

( 9.11) da funo delta, podemos facilmente demonstrar a identidade de
Parserval,

1
dq
F (q)G(q) = f (x)eiqx dx
2

1
g(y)eiqy dy dq
2

1
= f (x) g (y) eiq(xy) dq dx dy
2

= f (x) g (y)(x y) dx dy

= f (x)g (x) dx . (10.19)

10.4 Transformada das Derivadas

A transformada das derivadas muito importante no uso de soluo de
equaes diferenciais. Se F (q) a transformada de Fourier de f (x), ento
F {f (n) (x)} = (iq)n F (q) . (10.20)

Provemos para o caso n = 1. Por denio, temos que

1 df (x) iqx
F {f (x)} = e dx
2 dx

1
iqx 1
= f (x)e f (x)iqeiqx dx
2 2
= iqF {f (x)} , (10.21)
onde usamos o fato que f (x) 0 quando x . A demonstrao para
o caso n = 1 pode ser feita usando-se o mtodo de induo.
10.5 Transformada de Fourier Cosseno e Seno

As transformadas de Fourier cosseno e seno, respectivamente, so denidas
por
7
2
Fc {f (x)} = f (x) cos(qx) dx , (10.22)
0
7
2
Fs {f (x)} = f (x) sin(qx) dx . (10.23)
0
Note que se f (x) for uma funo par, ento F {f (x)} = Fc {f (x)}, e se f (x)
for mpar, ento F {f (x)} = Fs {f (x)}.
Exemplo 10.1 Determinar a transformada de Fourier da funo
f (x) = ea|x| , para todo a > 0.
Soluo: Lembrando que |x| = x se x > 0, e que |x| = x se x < 0,
temos que
0
1 (a+iq)x (a+iq)x
F (q) = e dx + e dx
2 0

1 1 1
= +
2 a + iq a + iq
7
2 a
= . (10.24)
a2 + q 2

Exemplo " 10.2 Determinar a transformada de Fourier inversa da
funo F (q) = 2/ sin(aq)/q para todo a > 0.
Soluo:

1 sin(aq) iqx
f (x) = e dq
q

1 sin(aq)
= cos(qx) dq
q

1 sin[(a + x)q] sin[(a x)q]
= dq + dq .
2 q q
(10.25)
Usando o teorema de Cauchy, mostramos anteriormente que a primeira
das integrais acima dada por para todo x > a e a segunda por para
todo x < a. Logo f (x) = 1 para todo |x| < a, e zero para todos os outros
casos.
Exemplo 10.3 Determinar a transformada de Fourier inversa da
funo H(q) = (2/)ab/(q 2 + a2 )(q 2 + b2 ).
Soluo: H duas maneiras de determinar a transformada inversa. Uma
seria usando o teorema dos resduos. A outra, talvez mais simples, seria
usando o teorema de convoluo. Seguiremos esse ltimo caminho. Note
que H(q) = F (q)G(q) um produto de duas transformadas de Fourier,
onde " "
F (q) = 2/a/(q 2 + a2 ) e G(q) = 2/b/(q 2 + b2 ). Assim, usando o
resultado do Exemplo 10.1 e (10.15), temos que

1
h(x) = ea|y| eb|xy| dy
2
0
1 ay b|xy| ay b|xy|
= e e dy + e e dy
2 0

1 ay b|x+y| ay b|xy|
= e e dy + e e dy (10.26)
,
2 0 0
onde na passagem da segunda para a terceira linha zemos a mudana de

varivel y y e em seguida invertemos os limites de integrao. Note que
f (x) invariante sob a troca de x por x. Com efeito, basta considerarmos
o caso x > 0. Temos
x
1 bx (a+b)y bx
h(x) = e e dy + e e(ba)y dy
2 0 0

bx (a+b)y
+e e dy
x
2 1 bx
= 2 2
ae beax . (10.27)
(a b )

Exemplo 10.4 Determinar a soluo da equao diferencial h (x) +
2 h = f (x) usando as condies de contorno h() = 0.

Soluo: Usando a transformada de Fourier da derivada, e transfor-
mando Fourier ambos os lados da equao diferencial, encontramos
q 2 H(q) + 2 H(q) = F (q) , (10.28)

de onde obtemos H(q) = F (q)/(2 q 2 ). Para obter a transformada de
Fourier inversa, usamos o teorema de convoluo,
7
1
h(x) = f (y) sin(|x y|) dy , (10.29)
2
onde usamos o resultado

7
1 1 1
F { 2 }= sin(x) , (10.30)
q2 2
a qual pode ser obtida facilmente usando o teorema dos resduos.
10.6 Problemas
10.1 Calcular a transformada de Fourier cosseno das seguintes funes
1
(a) f (x) = x2 +a2 ,
(b) f (x) = eax ,

1 , |x| < a ,
(c) f (x) =
0 , |x| > a ,
com a > 0 em todos os casos.
10.2 Calcular as transformadas de Fourier seno das seguintes funes
x
(a) f (x) = x2 +a2 ,
(b) f (x) = xeax ,
sin(ax)
(c) f (x) = x ,
com a > 0 em todos os casos.
10.3 Calcule a transformada de Fourier inversa das seguintes funes

!
2 2aiq
(a) F (q) = q 2 +a2 ,

(b) F (q) = 2 (q + a) ,
1
(c) F (q) = |q| ,
com a > 0 para o caso (a).
10.4 A equao de Fermi unidimensional para a difuso de nutrons em

algum meio (tal como o grate) dada por
2 q(x, ) q(x, )
2
= .
x
Na equao acima, q o nmero de nutrons que decaem abaixo de

alguma dada energia por segundo por unidade de volume. O tempo
de Fermi , a medida da perda de energia.
Se q(x, 0) = S(x), correspondendo a uma fonte plana de nutrons em

x = 0, emitindo S nutrons por unidade de rea por segundo, mostrar
que
2
ex /4
q(x, ) = S .
4
Sugesto: Substituir q pela sua transformada de Fourier inversa

1
q(x, ) = Q(k, )eikx dk .
2
Obtenha ento uma expresso para Q substituindo a equao acima

na equao diferencial.
Captulo 11
Teoria de Sturm-Liouville:
Funes Ortogonais
11.1 Introduo
A equao de Sturm-Liouville uma equao diferencial de segunda ordem
ordinria e linear que tem a forma

d
dx
p(x)
dy
dx
+ q(x)y = w(x)y , (11.1)
onde p(x), q(x) e w(x) so funes conhecidas e contnuas em um determi-

nado intervalo a x b. A constante determinada pelas condies de
Dirichlet ou de Neumann nos extremos do intervalo. Os valores de para
os quais as solues da equao diferencial so no-triviais, so conhecidos
como auto-valores e as correspondentes y(x) como auto-funes. Nem sem-
pre as equaes a auto-valores e auto-funes aparecem como em (11.1).
necessrio algumas manipulaes algbricas para coloca-las na forma de
Sturm-Liouville. Dizemos ento que a equao diferencial pode ser escrita
na forma auto-adjunta. Conforme veremos mais adiante, se a equao di-
ferencial auto-adjunta, ento as solues formam um conjunto de funes
ortogonais.
Jacques Charles Franois Sturm (1803 - 1855), matemtico francs, nascido em Ge-
neva; Joseph Liouville (1809 - 1882), matemtico francs, nascido em Saint-Omer.
11.2 Equao Diferencial Auto-adjunta

Sob determinadas condies, nesta Seo veremos como transformar uma
EDH em uma equao diferencial na forma auto-adjunta. Seja a EDH
escrita na forma geral
d2 y dy
+ B(x)
A(x) + C(x)y = 0 . (11.2)
dx2 dx
Conforme j vimos no Captulo 5, a EDH tambm pode ser reescrita
como
L y(x) = 0 , (11.3)
onde L o operador diferencial linear denido por (5.3). O operador ad-
junto denido por
d2 [A(x)y(x)] d[B(x)y(x)]
Ly(x) = + C(x)y(x) . (11.4)
dx2 dx
Sejam u e v duas solues da EDH. O operador diferencial L dito
auto-adjunto se satiszer a condio
b b
uL v dx = v Lu dx . (11.5)
a a
O operador L em geral no auto-adjunto. No entanto, impondo deter-

minadas condies sobre os coecientes das derivadas possvel transforma-
lo em um operador auto-adjunto. Usando integrao por partes, temos que
b b
b
d2 v
b
dv
uL v dx = uA 2 dx + uB dx + uCv dx
a a dx a dx a
b b
dv d(Au) dv
= uA dx
dx a a dx dx
b b
b d(Bu)
+ uBv|a v dx + vCu dx
a dx a
b b b
dv d(Au) d2 (Au)
= uA v + v dx
dx a dx a a dx2
b b
b d(Bu)
+ uBv|a v dx + vCu dx
a dx a
b
d2 (Au) d(Bu)
= v dx + Cu dx
a dx2 dx
Captulo 11: Teoria de Sturm-Liouville: Funes Ortogonais 179
+ u(A B)v|a + A(uv vu )|a .

b b
(11.6)
Se as solues satiszerem as condies de contorno de Dirichlet ou de Neu-

mann homogneas e
A (x) = B(x) , (11.7)

ento L auto-adjunto. Se os coecientes das primeira e segunda deri-
vadas no satiszerem esta condio, mesmo assim podemos realizar uma
transformao na EDH que a torna em uma auto-adjunta. Primeiramente
multiplicamos a EDH pelo fator
1 * B(x)
dx
e A(x) , A(x) = 0 , x , (11.8)
A(x)
onde a integral indenida. Temos
* B(x)
dx dy 2
B(x) * B(x)
dx dy C(x) * B(x)
dx
e A(x) + e A(x) + e A(x) y=0. (11.9)
dx2 A(x) dx A(x)
Agora, note que
* *
d B(x)
dx dy
B(x)
dx dy 2
B(x) * B(x)
dx dy
e A(x) =e A(x) + e A(x) . (11.10)
dx dx dx2 A(x) dx
Assim, a equao (11.9) pode ser reescrita na forma

d dy C(x)
p(x) + p(x)y = 0 , (11.11)
dx dx A(x)
onde
* B(x)
dx
p(x) = e A(x) . (11.12)
Escrita na forma (11.11), fcil ver que os novos coecientes das primeira
e segunda derivadas obedecem a equao (11.7), e assim auto-adjunta.
Escolhendo q(x), w(x) e apropriadamente, podemos reescrever (11.11)
na forma de Sturm-Liouville (11.1), de tal modo que
C(x)
p(x) = w(x) q(x) . (11.13)
A(x)
A equao de Sturm-Liouville usualmente escrita como
Por razes bvias, em problemas prticos, desconsideramos a constante de integrao.
L y(x) + w(x)y(x) = 0 , (11.14)
onde

d d
L = p(x) q(x) . (11.15)
dx dx
11.3 Funes Ortogonais

Suponhamos que y(x) satisfaa as condies de contorno de Dirichlet ou
de Neumann nos extremos do intervalo fechado a x b. As condies
de contorno podem dar origem a innitos valores de , conforme vimos no
Captulo 8. Consideremos duas auto-funes (ym (x), yn (x)) da equao de
Sturm-Liouville com correspondentes auto-valores (m , n ), m = n. As
duas funes satisfazem simultaneamente a equao (11.1),

d dym
p(x) + q(x)ym = n w(x)ym ,
dx dx

d dyn
p(x) + q(x)yn = m w(x)yn . (11.16)
dx dx
Multiplicando a primeira das duas equaes acima por yn , a segunda
por ym , subtraindo os resultados e integrando no intervalo a x b,
encontramos
b
d dyn d dym
ym p(x) yn p(x) dx
a dx dx dx dx
a
+(m n ) w(x)ym (x)yn (x) dx = 0 . (11.17)
b
Usando integrao por partes e as condies de contorno de Dirichlet

ou de Neumann homogneas, podemos facilmente mostrar que a primeira
integral se anula identicamente. Assim, temos que
a
(m n ) w(x)ym (x)yn (x) dx = 0 . (11.18)
b
Uma vez que m = n , obtemos

a
w(x)ym (x)yn (x) dx = 0 . (11.19)
b
Conseqentemente, {ym (x)}, com m inteiro, um conjunto de funes

ortogonais com relao funo peso w(x).
Exemplo 11.1 Determinar a forma auto-adjunta da equao diferen-
cial de Hermite
y 2xy + 2y = 0 , (11.20)
onde uma constante.
Soluo: No caso em questo, A(x) = 1, B(x) = 2x e C(x) = 2;
substituindo em (11.12), vem que
* 2
(2x) dx
p(x) = e = ex . (11.21)
Por meio de (11.13), temos que
2 x2 2
w(x) q(x) = e = 2ex , (11.22)
1
2
de onde obtemos q(x) = 0, = 2 e w(x) = ex . Portanto, substituindo
os resultados em (11.1), nalmente obtemos

d x2 dy 2
e + 2ex y = 0 . (11.23)
dx dx

11.4 Problemas
11.1 Considere as seguites equaes diferenciais ordinrias, lineares, e de
segunda ordem.
Equao de Legendre:
(1 x2 )y 2xy + l(l + 1)y = 0 .
Equao de Legendre associada:

m2
(1 x2 )y 2xy + l(l + 1) y=0.
1 x2
Equao de Chebyshev:
(1 x2 )y xy + n2 y = 0 .
Equao de Bessel:
x2 y + xy (x2 n2 )y = 0 .
Equao de Laguerre:
xy + (1 x)y + y = 0 .
Equao de Laguerre associada:
xy + (k + 1 x)y ( k)y = 0 .
Equao de Hermite:
y 2xy + 2y = 0 .
Equao do oscilador harmnico simples (OHS):
y + 2 y = 0 .
Usando a teoria de Sturm-Liouville, construa a Tab. 11.1.

Tabela 11.1:
Equao p(x) q(x) w(x)
de Legendre 1 x2 0 l(l + 1) 1
2
de Legendre associada 1 x2 1x
m
2 l(l + 1) 1

de Chebyshev 1 x2 0 n2 1
1x2
2
de Bessel x nx a2 x
de Laguerre xex 0 ex
de Laguerre associada xk+1 ex 0 k xk ex

2 2
de Hermite ex 0 2 ex
do OHS 1 0 2 1
Captulo 12
Polinmios de Hermite
12.1 Introduo
Neste Captulo estudaremos os polinmios de Hermite dentro de um con-
texto fsico onde eles emergem naturalmente. Para isto escolhemos a soluo
da equao de Schrdinger para o oscilador harmnico simples.
12.2 Oscilador Harmnico Simples Quntico

Considere a equao de Schrdinger unidimensional para uma partcula de
massa m presa a um potencial V (),
2 d2 ()

2m d 2
+ V ()() = E() , (12.1)
onde E a energia total da partcula. Estudaremos o caso simples, porm
muito importante, do oscilar harmnico, onde V () = 12 k 2 , sendo k a
constante de mola. Temos
2 d2 () 1 2

2m d 2
+ k () = E() .
2
(12.2)
conveniente fazer uma mudana de varivel de tal forma que as gran-
dezas fsicas tornem-se adimensionais. Seja, = x, onde = (2 /mk)1/4 ;
note que tem dimenso de comprimento. Desta forma, usando a regra da
cadeia, encontramos
d2 () 1 d2 (x)
d()
d
=
1 d(x)
dx
,
d 2
=
2 dx2
. (12.3)
Captulo 12: Polinmios de Hermite 185
Substituindo estes resultados em (12.2), obtemos a equao a auto-valores

e auto-funes,
+ x2 = , (12.4)
"
onde = 2E/ , sendo que = k/m a freqncia natural de oscilao
do oscilador harmnico simples.
Suporemos que a funo de onda de quadrado somvel , isto , (x)
0, quando x . Antes de resolver a equao (12.4), busquemos as
solues assintticas que satisfaam as condies de contorno. Na regio
assinttica, os termos dominates da equao so
+ x2 = 0 , (12.5)
2
cuja soluo aproximada dada por (x) = ex /2 . Isto pode ser vericado
por substituio direta da soluo aproximada no membro esquerdo da EDH
acima. Temos
2
+ x2 = ex /2
0 , quando x . (12.6)
Assim, para satisfazer as condies de contorno, escrevemos a soluo
de (12.4) como
2
(x) = ex /2
y(x) , (12.7)
onde y(x) uma funo a ser determinada. Temos
2
(x) = [y (x) xy(x)] ex /2 ,
2
(x) = y (x) 2xy (x) y(x) + x2 y(x) ex /2 . (12.8)
Introduzindo estes resultados na equao (12.4), somos conduzidos a
y (x) 2xy (x) + ( 1)y(x) = 0 . (12.9)

Esta equao conhecida como equao diferencial de Hermite (ver
Tab. 11.1). Note que a EDH regular x. Usaremos o mtodo de Frobenius
para resolve-la,
Charles Hermite (1822 - 1901), nasceu em Dieuze, Frana. Embora tenha tido uma
decincia fsica de nascena que o obrigou a usar bengala por toda a vida, isto no
desenvolveu nenhum tipo de complexo em sua personalidade. Por causa de sua decincia,
em 1842 foi dispensado da Escola Politcnica. Ironicamente, em 1846 voltou para a
mesma Escola como examinador de admisso.

y(x) = an xn , y (x) = nan xn1 , y (x) = n(n 1)an xn2 .
n=0 n=0 n=0
(12.10)
Substituindo na EDH (12.9), vem que

[(n + 1)(n + 2)an+2 + ( 2n 1)] xn = 0 , (12.11)
n=0
de onde se obtm a relao de recorrncia
2n + 1
an+2 = an , n 0 . (12.12)
(n + 1)(n + 2)
Os coecientes pares esto relacionados com a0 e os mpares com a1 .
Isto produz duas solues linearmente independentes. Agora, note que
an+2 2n + 1
lim = lim
n an n(n + 1)(n + 2)
1
= lim
l l + 1
= 0, (12.13)
onde usamos a regra de L'Hpital e l = n + 1/2. Logo a srie converge para
todo < x < . No entanto, note que 1/(n + 1) a razo entre os
coecientes adjacentes da srie da funo
2 x2 x4 x6
ex = 1+ + + +
1! 2! 3!

x2n
= , (12.14)
n=0
n!
pois
an+2 n!
lim = lim
n an (n + 1)!
n
1
= lim
n n + 1
= 0. (12.15)
2
Isto signica que o comportamento assinttico de y(x) do tipo ex para
x . Este comportamento de y(x) viola as condies de contorno, pois
2
(x) ex /2 , a qual tende ao innito quando x . Todavia, as
solues podem ser regularizadas truncando a srie para um determinado
valor de n = m. Temos
2m + 1 = 0 , (12.16)
m
de onde obtemos a quantizao da energia E = 2 + 1 .
Assim, as duas solues linearmente independentes sero dadas pelos
seguintes polinmios
y1 (x) = a0 + a2 x2 + a4 x4 + + a2m x2m ,

y2 (x) = a1 x + a3 x3 + a5 x5 + + a2m+1 x2m+1 , (12.17)
onde, em qualquer dos casos

2(n m)
an+2 = an , 0 n < m . (12.18)
(n + 1)(n + 2)
A seguir, geramos todas as solues em forma de polinmios para cada
valor de m. Se m for par, usamos somente y1 (x) e se m for mpar, usa-
mos somente y2 (x), pois a condio de truncamento no pode ser satisfeita
simultaneamente para um determinado valor de m. Temos
1. m = 0: a0 = 0, a2 = a4 = = 0, y1 (x) = a0 = a0 H0 (x), onde

H0 (x) = 1;
2. m = 1: a1 = 0, a3 = a5 = = 0, y2 (x) = a1 x = 2 H1 (x),
a1
onde
H1 (x) = 2x;
3. m = 2: a0 = 0, a2 = 2a0 , a4 = a6 = = 0, y1 (x) = a0 2a0 x2 =
a20 H2 (x), onde H2 (x) = 4x2 2;
4. m = 3: a1 = 0, a3 = 23 a1 , a5 = a7 = = 0, y2 (x) = a1 x 23 a1 x3 =
a121 H3 (x), onde H3 (x) = 8x3 12x;
e assim por diante. As funes {Hn (x), n 0} so conhecidas como polin-

mios de Hermite. Eles podem ser representados de forma genrica como
n
[2]
(1)l
Hn (x) = n! (2x)n2l , (12.19)
l!(n 2l)!
l=0
onde [ n2 ] o maior inteiro menor ou igual a n2 .

Finalmente, a soluo fundamental para a equao (12.4) pode ser escrita
como
2
n (x) = ex Hn (x) . (12.20)
A soluo geral de (12.4) obtida pelo princpio de superposio. Mais
adiante
* veremos como normalizar a funo de onda de tal modo que

m (x) n (x) dx = mn . Antes, iremos estudar uma srie de propri-
edades dos polinmios de Hermite.
12.3 Relaes de Recorrncia

Derivando-se a equao diferencial para os polinmios de Hermite
d2 Hn (x) dHn (x)

2x + 2nHn (x) = 0 , (12.21)
dx2 dx
obtemos

d2 dHn (x) d dHn (x) dHn (x)
2x + 2(n 1) =0. (12.22)
dx2 dx dx dx dx
Ento, Hn (x) soluo da equao de Hermite com n substitudo por

(n 1). Conseqentemente, Hn (x) e Hn1 (x) no so linearmente indepen-
dentes,
Hn (x) = CHn1 (x) , (12.23)

onde C uma constante a ser determinada. Usando (12.19), temos que
n
[2]
(1)l 2(n 2l)
n! (2x)n2l1 = C(n 1)!
l!(n 2l)!
l=0
[ n1
2 ]
(1)l
(2x)n2l1 . (12.24)
l!(n 2l 1)!
l=0
Igualando os coecientes das potncias de x, determinamos a constante

C . Por exemplo, para l = 0, temos que
1 1
n! 2n = C(n 1)! C = 2n . (12.25)
n! (n 1)!
e
Hn (x) = 2nHn1 (x) . (12.26)
Esta relao de recorrncia pode ser usada para obter outras relaes im-
portantes. Vejamos um exemplo. Derivando (12.26), encontramos Hn (x) =

2nHn1 (x). Levando este resultado e (12.26) na equao de Hermite
(12.21), temos que

2nHn1 (x) 2x[2nHn1 (x)] + 2nHn (x) = 0 , (12.27)
de onde se deduz que

d
Hn (x) = 2xHn1 (x) Hn1 (x) = 2x Hn1 (x) , (12.28)
dx
como tambm
Hn+1 (x) = 2xHn (x) 2nHn1 (x) . (12.29)

A relao de recorrncia (12.28) pode ser usada recursivamente para
obtermos
2
d
Hn (x) = 2x Hn2 (x)
dx
3
d
= 2x Hn3 (x)
dx
..
. n
d
= 2x H0 (x) . (12.30)
dx
Assim, os polinmios de Hermite podem ser determinados por aplicao
sucessiva do operador diferencial (2x d/dx) ao polinmio de mais baixa
ordem H0 (x).
12.4 Funo Geradora

Considere uma seqncia de polinmios {pn (x), n = 0, 1, 2, . . .}. A idia
fundamental da teoria da funo geradora determinar uma funo g(x, t)
cujos coecientes da expanso desta em potncia de t sejam dados pelas
funes pn (x). Formalmente, temos que

pn (x) n
g(x, t) = t . (12.31)
n=0
n!
Iremos aplicar esta relao para os polinmios de Hermite e determinar

g(x, t) que os geram,

Hn (x) n
g(x, t) = t . (12.32)
n=0
n!
Derivando os dois membros com relao a t, vem que

g(x, t) Hn (x) n1
= nt
t n=1
n!

Hn+1 (x)
= (n + 1)tn
n=0
(n + 1)!

Hn+1 (x) n
= t . (12.33)
n=0
n!
Por outro lado, usando (12.29), temos que

g(x, t) Hn (x) n Hn1 (x) n
= 2x t 2 n t
t n=0
n! n=1
n!

Hn (x) n+1
= 2xg(x, t) 2 (n + 1) t
n=0
(n + 1)!

Hn (x) n
= 2xg(x, t) 2t t
n=0
n!
= (2x 2t)g(x, t) . (12.34)
Agora, a funo geradora pode ser facilmente obtida por integrao do

resultado acima. Temos

dg
= (2x 2t) dt = 2xt t2 + C(x) , (12.35)
g
donde
2
g(x, t) = h(x)e2xtt . (12.36)
A constante de integrao h(x) = eC(x) pode ser determinada usando-se

a expanso

h(x) 1 + (2xt t2 ) + = H0 (x) + H1 (x)t + ,
h(x) + h(x)(2x)t + = H0 (x) + H1 (x)t + , (12.37)
donde, igualando os coecientes das potncias de t, encontramos h(x) =

H0 (x) = 1, ou ainda, h(x)(2x) = H1 (x) = 2x h(x) = 1. Finalmente
obtemos para a funo geradora
2
g(x, t) = e2xtt . (12.38)
12.5 Frmula de Rodrigues

Por meio da funo geradora podemos deduzir uma srie de propriedades dos
polinmios de Hermite. Discutiremos uma delas que a obteno da frmula
de Rodrigues. Para isto, primeiramente expandimos a funo geradora em
uma srie de Taylor,

1 n g(x, t) Hn (x) n
g(x, t) = n
t n
= t , (12.39)
n=0
n! t t=0 n=0
n!
donde obtemos que

n g(x, t)
Hn (x) = . (12.40)
tn t=0
Antes de calcular a derivada acima, notemos que a funo geradora

2 2
tambm pode ser escrita como g(x, t) = ex (tx) . Assim,
n (tx)2 n
n 2
n
e = (1) n
e(tx) . (12.41)
t x
Portanto,

n x2 n (tx)2
Hn (x) = (1) e e
xn t=0
2 d
n 2
Hn (x) = (1)n ex ex . (12.42)
dxn
Esta equao conhecida como frmula de Rodrigues para os polinmios
de Hermite.
12.6 Ortogonalidade
Conforme vimos no Captulo 11, a teoria de Sturm-Liouville nos garante que
os polinmios de Hermite formam um conjunto completo de funes, pois
a equao diferencial de Hermite pode ser escrita na forma auto-adjunta.
Temos

2
ex Hm (x)Hn (x) dx = 0 , m = n . (12.43)

Consideremos agora o caso m = n,

2 2 dn x2
In = ex [Hn (x)] dx = Hn (x)(1)n e dx , (12.44)
dxn
onde usamos a frmula de Rodrigues (12.42). Usando integrao por par-

tes, a frmula de Rodrigues novamente e a relao de recorrncia (12.26),
encontramos

dn1 x2 dHn (x) dn1
In = n
(1) Hn (x) n1 e (1)n n1 dx
dx dx dx

2 2 2
= ex Hn (x)Hn1 (x) + 2n ex [Hn1 (x)] dx .

(12.45)
Ora, mas o primeiro termo do segundo membro se anula identicamente,

pois a funo exponencial tende mais rapidamente a zero do que os polin-
mios de Hermite ao innito. Assim, In = 2nIn1 . Usando esta frmula
recursivamente, temos que
In = [2n][2(n 1)]In2
= [2n][2(n 1)][2(n 2)]In3
..
.
= [2n][2(n 1)][2(n 2)] [2.1] I0 = 2n n!I0 ,

nprodutos
(12.46)
* 2
onde I0 =
ex dx = . Finalmente obtemos In = 2n n! . Denin-
do
1 2
n (x) = ex Hn (x) , (12.47)
2n n!
a relao de ortogonalidade pode ser escrita como

m (x)n (x) dx = mn . (12.48)

A funo (12.47) de fato a soluo fundamental da equao (12.4) com

= 2n + 1.
12.7 Problemas
12.1 Mostre que
(2m)!
H2m (0) = (1)m ,
m!
H2m+1 (0) = 0.
2
12.2 Mostre que et cos(2xt) a funo geradora dos polinmios de Hermite
2
de ordem par e que et sin(2xt) a funo geradora dos polinmios de
Hermite de ordem mpar.
12.3 Mostre que os polinmios de Hermite tm a seguinte representao
integral
)
2 1 2
Hn (x) = ex n! tn1 e(tx) dt ,
2i C
onde C um caminho fechado em torno da origem.

Sugesto : Multiplicar ambos os lados da identidade

2xtt2 Hn (x) n
e = t ,
n=0
n!
por tm1 . Pela frmula integral de Cauchy, podemos escrever

)
1 1 dm n
tnm1 dt = t = mn .
2i C n! dtm
Usando este resultado proceda com a demonstrao.

12.4 Usando a formula integral de Cauchy, e usando o resultado do Pro-

blema 12.3, obtenha a frmula de Rodrigues
2 dn x2
Hn (x) = (1)n ex e .
dxn
12.5 Usando a frmula de Rodrigues para n = 0, 1, 2, 3, obtenha os polin-

mios de Hermite determinados na Seo 12.2.
12.6 Mostre que Hn (x) = (1)n Hn (x), ou seja, H2m (x) uma funo
par e H2m+1 (x) uma funo mpar.
12.7 Usando a relao de recorrncia (12.29) demonstre a frmula

2 1 1
x Hn (x) = n(n 1)Hn2 (x) + n + Hn (x) + Hn+2 (x) .
2 4
Captulo 13
Problemas Adicionais
13.1 Sries Innitas
1. A integral elptica completa de segunda espcie dada por
/2
E(m) = (1 m sin2 )1/2 d .
0
Usar a expanso binomial para mostrar que

8 2 2
1 m 13 m2
E(m) = 1
2 2 1 24 3
2 9
135 m3
.
246 5
Sugesto: Usar o resultado

/2
(2n 1)!!
sin2n d = .
0 (2n)!! 2
2. A funo beta incompleta dada pela integral

x
Bx (p, q) = tp1 (1 t)q1 dt .
0
Usando a expanso binomial, mostrar que


1 1q (1 q)(2 q) 2
Bx (p, q) = xp + x+ x +
p p+1 2!(p + 2)

(1 q)(2 q) (n q) n
+ x + .
n!(p + n)
A expanso acima s vale para os casos 0 x < 1, p > 0, e q > 0 (se

x = 1).
3. Usando a expanso binomial, mostrar que

tanh x = 1 + 2 (1)n e2nx ,
n=1
para todo x > 0.

4. Dado que
1
dx
= ,
0 1 + x2 4
mostrar que o numero pode ser representado pela srie
1 1 1 1 1
= 1 + + + (1)n + .
4 3 5 7 9 2n + 1
5. Usar a expanso binomial para mostrar que
x
= x + 2x2 + 3x3 +
(1 x)2

= nxn .
n=1
6. Considere um corpo em queda livre sujeito resistncia do ar. A

velocidade do ponto material em funo do tempo dada por
mg
v(t) = 1 ebt/m ,
b
onde b o atrito, m a massa do corpo em queda livre e g a acelerao
da gravidade local. Mostrar que
Captulo 13: Problemas Adicionais 197
1 bg 2
v(t) = gt + t + .
2m
Note que no limite bt/m 1, obtemos o movimento uniformemente

variado, isto , v(t) = gt. Por outro lado, quando bt/m 1, o
movimento uniforme, pois a velocidade constante e vale v(t) =
mg/b.
13.2 Anlise Vetorial

+
1. Calcular A n dS , onde A = xz 2 i + (x2 y z 3 )j + (2xy + y 2 z)k e S
uma superfcie fechada limitada por uma esfera de raio a.
+
2. Calcular C A dr, onde A = yi + xj + k e C um quadrado de lado
1 com centro na origem.
3. Provar que F = (2xz 3 + 6y)i + (6x 2yz)j + (3x2 z 2 y 2 )k um campo

de foras conservativo, isto , mostrar que F = 0. Determinar a
funo escalar (x, y, z) tal que F = .
4. A induo magntica B relaciona-se com o potencial vetor A por

B = A. Pelo teorema de Stokes
)
B n dS = A dr .
Mostrar que cada lado desta equao invariante sob a transformao

de calibre A A + , onde uma funo escalar.
5. Mostrar que se B constante, a equao
1
A= Br,
2
satisfaz B = A.
6. Seja uma funo denida num volume V limitado por uma superfcie
fechada S . Supor que satisfaz a equao de Laplace. Mostrar que
a integral sobre qualquer superfcie fechada dentro de V da derivada
normal n ser zero.
13.3 Coordenadas Curvilneas

1. Considere a transformativas de coordenadas abaixo:
x = cos ,
y = sin ,
1 2
z= ( 2 ) .
2
Determinar a matriz de transformao Aij = h1i xqj e mostrar que

i
ortogonal. Calcule as mtricas h1 , h2 e h3 . Escrever o gradiente, o

divergente, o Laplaciano e o rotacional no novo sistema de coordena-
das.
2. Repita o exerccio anterior para as equaes transformativas de coor-

denadas abaixo:
u2 + v 2
x = ,
2
u2 v 2
y = ,
2
z = z.
3. Se e so coordenadas polares, e A, B , n constantes quaisquer,

prove que U = n (A cos n + B sin n) satisfaz a equao de Laplace:
2 U = 0.
4. Se
2 cos + 3 sin3 cos

V = ,
r2
mostrar que
6 sin cos (4 5 sin2 )

2 V = .
r4
13.4 Equaes Diferenciais Ordinrias

1. Seja a EDNH
d2 y dy
x2 2x + (x2 + 2)y = x3 ex .
dx2 dx
Se y1 (x) = x cos x soluo da EDH associada, mostrar que a segunda
soluo dada por y2 (x) = x sin x, e que uma soluo particular da
EDNH dada por u(x) = 12 xex .
2. Seja a EDNH
d2 y dy 1 x2
(1 + x2 ) 2x + 2y = .
dx2 dx x
Se y1 (x) = x2 1 soluo da EDH associada, mostrar que a segunda
soluo dada por y2 (x) = x, e que uma soluo particular da EDNH
dada por u(x) = x ln x.
3. Seja a EDNH
d2 y dy
x2 x(2x + 3) + (x2 + 3x + 3)y = (6 x2 )ex .
dx2 dx
Se y1 (x) = x3 ex soluo da EDH associada, mostrar que a segunda
soluo dada por y2 (x) = xex , e que uma soluo particular da
EDNH dada por u(x) = ex (x2 + 2).
4. Resolver as seguintes equaes diferenciais usando o mtodo de Fro-

benius:
(a) 2x2 y xy + (x2 + 1)y = 0 ,

(b) x2 y + xy + x2 y = 0 .
5. A equao hipergeomtrica de Gauss tem a forma
x(1 x)y (x) + [c (a + b + 1)x]y (x) aby(x) = 0 .
Mostre que x = 0 e x = 1 so dois pontos singulares regulares da

equao diferencial. Usando o mtodo de Frobenius, mostrar que uma
soluo desta equao dada por
ab x a(a + 1)b(b + 1) x2
y(x) = 1 + + + ,
c 1! c(c + 1) 2!
com c = 0, 1, 2, 3, . . . .
Sugesto: Use apenas a raiz s = 0 da equao indicial.
6. A equao hipergeomtrica conuente de Gauss tem a forma
xy (x) + (c x)y (x) ay(x) = 0 .
Mostre que x = 0 um ponto singular regular da equao diferen-

cial. Usando o mtodo de Frobenius, mostrar que uma soluo desta
equao dada por
y(x) = M (a, c; x)
ax a(a + 1) x2
= 1+ + +
c 1! c(c + 1) 2!
a(a + 1) (a + n) xn
+ + .
c(c + 1) (c + n) n!
Mostre que M (1, 1; x) = ex .

Sugesto : Use apenas a raiz s = 0 da equao indicial.
7. A equao de Bessel modicada de ordem zero tem a forma
xy (x) + y (x) xy(x) = 0 .
Mostre que x = 0 um ponto singular regular da equao diferen-

cial. Usando o mtodo de Frobenius, mostrar que uma soluo desta
equao dada por
x2 x4 x6 x2n
y(x) = 1 + + + + + .
[2!!]2 [4!!]2 [6!!]2 [(2n)!!]2
Sugesto : Use apenas a raiz s = 0 da equao indicial.

8. Considere a equao diferencial homognea abaixo
(1 x2 )y (x) xy (x) + n2 y(x) = 0 .
Mostre que x = 1 so pontos singulares regulares da equao dife-

rencial. Usando o mtodo de Frobenius, mostrar que duas solues
linearmente independentes desta equao so dadas por

(12 n2 ) 3 (32 n2 ) 5
y1 (x) = c1 x + x + x + ,
3! 5!
2 2 2 2
n 2 (n 2 )n 4
y2 (x) = c0 1 x + x
2! 4!

(n2 42 )(n2 22 )n2 6
x + .
6!
Ajuste as constantes c0 e c1 tal que y1 (1) = y2 (1) = 1 e mostre que as

solues so polinomiais se n for um inteiro positivo:
T0 (x) = 1 ,
T1 (x) = x ,
T2 (x) = 2x2 1 ,
T3 (x) = 4x3 3x ,
e assim por diante. Estas funes so conhecidas como polinmios de

Chebyshev.
Sugesto : considere somente a raiz s = 0 da equao indicial.
13.5 Variveis Complexas

1. Escolhendo apropriadamente o caminho de integrao, use o teorema
dos resduos para calcular as seguintes integrais:
* + x sin(ax)
(a) (x2 +b2 ) dx = eab , (a > 0, b > 0) ,
* x sin x (ea eb )
(b) 0 (x2 +a2 )(x2 +b2 ) dx = 2(b2 a2 ) , (a = b) ,
* cos(ax) a
(c) (x)2 + 2
dx = e cos(a) , (a > 0 , > 0 , > 0) ,
*
x sin(ax) ab/ 2 ab
(d) dx =
0 x4 +b4 2b2 e sin 2
, (a > 0 , b > 0) .
2. Considere a funo
eaz
f (z) = ,
1 + ez
no domnio R x R, 0 y 2 , ou seja um retngulo de
lados 2R e 2 (ver Fig. 13.1). Mostre que z = i um plo de f (z).
Expandir f (z) em uma srie de Laurent em torno deste ponto, isto ,
mostrar que

eaz eaz (z i) (z i)2
= 1 + .
1 + ez (z i) 2 12
Sugesto : Expandir 1 + ez em srie de Taylor em torno do ponto
z = i . Em seguida, escolhendo w apropriadamente, use 1/(1 + w) =
1 w + w2 , para calcular 1/(1 + ez ).
y
R + 2i R + 2i
R r R x
Figura 13.1:
3. Usando o resultado da exerccio anterior, e o teorema dos resduos,

mostrar que
+
eax
dx = , 0 < a < 1.
1 + ex sin(a)
Sugesto : Primeiramente mostrar que
) , -
+R R
eaz eax eax
dz = lim dx ei2a dx
C 1 + ez R R 1 + ex R 1 + ex
+ ax
e
= (1 ei2a ) dx ,
1 + ex
e que ao longo dos ramos verticais x = R, a integral anula-se no
limite R .
13.6 Sries de Fourier

1. Considere a funo f (x) = x para 2 < x < 2. Mostre que a srie de
Fourier de f (x) nesse intervalo dada por

4 (1)n+1 nx
x= sin .
n=1 n 2
Usando a identidade de Parserval, mostrar que

1 2
= .
n=1
n2 6
2. Considere a funo

1, 1 < x < 0 ,
f (x) =
1 , 0<x<1.
Mostre que a srie de Fourier de f (x) nesse intervalo dada por

4 sin((2n + 1)x)
f (x) = .
n=0 2n + 1
Usando a identidade de Parserval, mostrar que

1 2
= .
n=0
(2n + 1)2 8
13.7 Equaes Diferenciais Parciais

1. Considere a seguinte equao diferencial parcial linear de segunda or-
dem

2 2h 2h
+ 2 +h=0,
x2 y
sujeita s condies de contorno de Dirichlet da Fig. 13.2. Na gura,

H uma constante. Mostre que a soluo desta equao dada por

(1)n n x cosh(n y)
h(x, y) = 2H cos( ) ,
n=0
n a cosh(n b)
"
onde n = (n + 1/2) , e n = 1/2 + n2 /a2 .
+b h(x, b) = H
h(a, y) = 0
+a
x
a
h(a, y) = 0
h(x, b) = H b
Figura 13.2:
Sugesto : Use o mtodo de separao de variveis,

h(x, y) = X(x)Y (y) e mostre que
X 1 Y
= 2 =.
X Y
Considere trs casos distintos: (a) = 2 ( positivo); (b) = 0; (c)

= 2 ( negativo). Mostre que para os casos (a) e (b), usando as
condies de contorno, obtemos a soluo trivial para X(x). Para o
caso ?, usando as condies de contorno, mostre que a soluo para

X(x) no-trivial se e somente se = n/a ou = (n + 1/2)/a.
Use apenas a segunda soluo para . Mostre ento que
n x
Xn (x) = cos( ),
a
onde foi ignorado a constante de integrao. Usando as condies de
contorno, mostre que a soluo para Y (y) dada por
cosh(n y)
Yn (y) = H .
cosh(n b)
Na construo desta soluo importante o uso da identidade

sinh(2c) = 2 sinh c cosh c.
Use o princpio de superposio para obter a soluo formal

n x cosh(n y)
h(x, y) = H cn cos( ) .
n=1
a cosh(n b)
Use o fato de que { 1a cos(n x/a) , n = 0, 1, 2, 3, . . .}, forma um con-

junto de funes linearmente independentes para obter
cn = 2(1)n /n .
13.8 Funo Delta de Dirac

1. Determine a srie de Fourier da funo delta de Dirac no intervalo
L x L.
2. No Problema 13.7.1 usamos o fato que
{ 1a cos(n x/a) , n = 0, 1, 2, 3, . . .} forma um conjunto completo de
funes ortogonais. Mostre que isto verdadeiro para todo a < x <
a. Qual deve ento ser a representao da funo delta em termos
destas funes?
13.9 Teoria de Sturm-Liouville

1. Considere a equao de Gegenbauer,
(1 x2 )y 2(1 + )xy + n(n + 2 + 1)y = 0 ,

onde n um inteiro e real.

Esta equao pode ser escrita na forma auto-adjunta de Sturm-Liou-
ville? Em caso armativo, obtenha a forma auto-adjunta da equao.
13.10 Polinmios de Hermite

1. Usando a funo geradora dos polinmios de Hermite, mostre que

2
ex Hm (x)Hn (x) dx = 2n n! mn .

* 2
Sugesto : Usar o resultado
ey dy = .
Referncias Bibliogrcas
[1] GRADSHTEYN, I. S.; RYZHIK, I. M. Tables of integrals, series, and
products . San Diego: Academic Press, 1994.
[2] SPIEGEL, M. R. Manual de frmulas e tabelas matemticas . So Paulo:
McGraw-Hill, 1968.
[3] ARFKEN, G. Mathematical methods for physicists . New York: Acade-
mic Press, 1970.
[4] SPIEGEL, M. R. Clculo avanado . So Paulo: McGraw-Hill, 1972.
[5] BUTKOV, E. Fsica matemtica . Rio de Janairo: Guanabara-Koogan,
1988.
[6] CHURCHILL, R. V. Variveis complexas e aplicaes . So Paulo:
McGraw-Hill, 1975.
[7] VILA, G. Variveis complexas e aplicaes . Rio de Janeiro: LTC,
2000.
ndice Remissivo
Condies Elemento de rea em coordenadas cur-

de Cauchy-Riemann, 105 vilneas, 63
de contorno, 70, 143 Equao
no-homogneas, 147 da continuidade, 56
problemas de Dirichlet e de da onda, 59, 148
Neumann, 150 de Bessel, 74, 182
Critrio de Convergncia de d'Alem- de Chebyshev, 181
bert, 25 de difuso, 69
de difuso do calor, 143
Delta de Kronecker, 40 de Hermite, 181, 182
Densidade de corrente, 54 de Laguerre, 182
de deslocamento, 56 de Laguerre associada, 182
Dipolo eltrico, 24 de Laplace, 69
Distribuio transiente, 147 em coordenadas polares-cilndricas,
Divergente 73
de uma funo vetorial, 44 de Legendre, 181
em coordenadas cilndricas, 67 de Legendre associada, 181
em coordenadas esfricas, 68 de Poisson, 69
em coordenas curvilneas, 65 de Schrdinger, 69
Domnio unidimensional, 184
multiplamente conexo, 109 de Sturm-Liouville, 177
simplemente conexo, 109 diferencial homognea, 76
Duplo Fatorial
de nmero par, 16 diferencial no-homognea, 76
de nmero mpar, 16 soluo geral de, 78
Duplo produto vetorial, 37 hipergeomtrica conuente de Gauss,
200
Elemento de comprimento em coor- hipergeomtrica de Gauss, 199
denadas curvilneas, 62 indicial, 85, 87
Elemento de volume raizes da, 87
em coordenadas cilndricas, 66 Equaes
em coordenadas curvilneas, 63 de Cauchy-Riemann, 106
em coordenadas esfricas, 67 de Maxwell, 53
ndice Remissivo 209
transformativas, 60 de Moivre, 104

Estado estacionrio, 147 de Rodrigues
Expanso para os polinmios de Hermite,
binomial, 23 191
da funo exponencial, 22 para os polinmios de Legen-
da funo logaritmo, 22 dre, 100
integral de Cauchy, 111
Fatorial, 16
Fora central, 47 Gradiente
de Coulomb, 48 de uma funo escalar, 41
gravitacional, 48 em coordenadas cilndricas, 67
Funo em coordenadas esfricas, 68
analtica, 105 em coordenas curvilneas, 64
beta incompleta, 195 interpretao geomtrica, 41
de Bessel, 101
de Heaviside, 162 Identidade de Parserval, 137
de onda quadrada, 133 Integral
de quadrado somvel, 185 de linha, 49
delta de Dirac, 157 elptica completa de segunda es-
em coordenadas cilndricas, 164 pcie, 195
em coordenadas esfricas, 167 por partes, 17
em trs dimenses, 163
largura de, 159 Laplaciano
propriedade fundamental, 160 em coordenadas cilndricas, 67
representaes de, 158 em coordenadas esfricas, 68
dente-de-serra, 134 em coordenas curvilneas, 65
efrica de Bessel, 98 Laplaciano de uma funo escalar, 47
efrica de Neumann, 98 Lei
imprpria, 157 de Ampre, 54
par, 137 de Ampre-Maxwell, 56
peso, 181 de Faraday, 55
mpar, 138 de Gauss, 53
Funes dos cossenos, 36
complexas Lema de Jordan, 128
plurvoca, 105
uniforme, 105 Mtodo
hiperblicas, 14 de Frobenius, 81
ortogonais, 134, 181 de separao de variveis, 69
peridicas, 133 Mtricas
trigonomtricas, 15 em coordenadas cilndricas, 66
Frmula em coordenadas esfricas, 67
de Euler, 15, 104 Maple, 96
210 ndice Remissivo
Mathematica, 96 Ponto
MATLAB, 96 ordinrio, 81
Matriz singular, 81, 115
ortogonal, 62 irregular, 81
Matriz ortogonal, 40 regular, 81
Monopolo magntico, 55 Princpio de superposio, 146
Mtricas, 61 Produto
escalar, 34
Nmeros propriedade comutativa, 34
complexos, 102 propriedade distributiva, 34
adio e subtrao de, 102 misto, 37
conjugado de, 102 vetorial, 34
desigualdade triangular, 104 mdulo de, 35
diagrama de Argand, 103 propriedade anti-comutativa, 35
diviso de, 103 propriedade distributiva, 35
forma polar de, 104 Projeo de um vetor sobre um se-
mdulo de, 103 gundo vetor, 37
multiplicao de, 103
parte imaginria de, 102 Regra
parte real de, 102 de derivao, 16
de Bernoulli, 28 de L'Hpital, 16
do paralelogramo, 31
Operador Relao de recorrncia, 83
adjunto, 178 Resto da srie de Taylor, 21
auto-adjunto, 178 Rotacional
de abaixamento, 73 de uma funo vetorial, 44
de levantamento, 73 em coordenadas cilndricas, 67
diferencial L , 75, 178 em coordenadas esfricas, 68
momento angular, 73 em coordenas curvilneas, 65
nabla , 41 Rotao de coordenadas
frmulas que envolvem o, 46 em trs dimenses, 39
no plano, 38
Plo
de ordem n, 115 Srie
duplo, 115 de Fourier, 136
simples, 115 generalizada, 135
Perodo de um pndulo simples, 29 de Laurent, 116
Polinmio caracterstico, 94 de Maclaurin, 20
Polinmios de Taylor, 18, 81, 115
de Chebyshev, 201 de Taylor revisitada, 21
de Hermite, 187 Seqncia delta, 159
de Legendre, 27, 100 Singularidade essencial, 116
ndice Remissivo 211
Sistemas de coordenadas
cilndricas, 66
curvilneas, 60
esfricas, 67
ortogonais, 32, 40
retangulares, 32, 60
Teorema
de Stokes, 53, 108
do divergente, 51
dos resduos, 118
integral de Cauchy
domnio multiplamente conexo,
109
domnio simplesmente conexo,
108
Varivel complexa, 103, 105

Vetor
componentes de um, 33
de posio, 34
em coordenadas cilndricas, 66
em coordenadas esfricas, 67
de Poyting, 59
mdulo de um, 30, 33
multiplicaao de um escalar por
um, 31
oposto de um, 31
Vetores
lgebra dos, 31
e escalares, 30
notao de, 17
propriedade associativa, 32
propriedade comutativa, 32
soma de, 31
subtrao de, 31
unitrios, 32
Volume de um paraleleppedo, 37
em coordenadas curvilneas, 63
Vorticidade, 44, 45, 55
Wronskiano, 77
Sobre o Autor
Edson Sardella professor do Departamento de Fsica da Faculdade de Ci-

ncias do Campus da UNESP de Bauru desde 1995. Tambm j foi professor
da Faculdade de Engenharia de Ilha Solteira do Campus da UNESP de Ilha
Solteira de 1985 a 1992. Foi graduado em Licenciatura em Fsica pela Uni-
versidade Federal de So Carlos em 1981. Em 1985 adquiriu o ttulo de
Mestre em Fsica pelo Instituto de Fsica Gleb Wataghin da Universidade
Estadual de Campinas. Ao nal de 1993 doutorou-se em Fsica pela Uni-
versidade de Manchester, Inglaterra. Sua rea de atuao em pesquisa
Fsica da Matria Condensada, onde possu vrios trabalhos em Supercon-
dutividade.
PR REITORIA
DE GRADUAO

Física Matem Teoria e Aplicações PDF

Enviado por

Direitos autorais:

Formatos disponíveis

Física Matem Teoria e Aplicações PDF

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

Física Matem Teoria e Aplicações PDF

Enviado por

Direitos autorais:

Formatos disponíveis

PR REITORIA

Inclui Bibliografia e ndice Remissivo

Ficha catalogrfica elaborada pela Coordenadoria Geral de Bibliotecas da Unesp

Pr-Reitoria de Extenso Universitria

Cultura Acadmica Editora

FUNDAO EDITORA DA UNESP

" #  ! 

Considerando a importncia da produo de material didtico-

3.5 Produto Escalar . . . . . . . . . . . . . . . . . . . . . . . . . 34

5 Equaes Diferenciais Ordinrias 75

5.8 Solues Numricas . . . . . . . . . . . . . . . . . . . . . . . . 95

6 Variveis Complexas 102

7 Sries de Fourier 133

8 Problemas de Valores de Contorno 143

8.4 Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155

9 Funo Delta de Dirac 157

10 Transformadas de Fourier 168

11 Teoria de Sturm-Liouville: Funes Ortogonais 177

12 Polinmios de Hermite 184

13 Problemas Adicionais 195

13.8 Funo Delta de Dirac . . . . . . . . . . . . . . . . . . . . . . 205

1.2 Funes Hiperblicas

1.3 Relaes de Euler

1.4 Relaes entre Funes Trigonomtricas e

1.5 Propriedades Bsicas da Funo ln x

1.6 Fatorial e Duplo Fatorial

A segunda (terceira) das relaes acima representa o produto somente dos

1.7 Regras de Derivao

1.8 Regra de L'Hpital

1.9 Integrao por Partes

1.10 Notao de Vetores

2.2 Srie de Taylor

Integrando novamente, encontramos

onde na passagem da primeira para a segunda linha usamos (2.2) recursi-

Continuando este processo, obtemos para a n-sima integral

Resolvendo para f (x), nalmente encontramos

onde Rn conhecido como o resto da srie e denido por

Quando a funo f (x) for tal que

ento a srie torna-se uma soma innita de termos. Temos

Esta srie conhecida como srie de Taylor.

2.3 Srie de Maclaurin

 Publicado pelo matemtico ingls Brook Taylor (1685 - 1731) em 1715.

Aplicando este teorema para a integral mais interior de (2.8), temos

2.5 Srie de Taylor Revisitada

Por meio da mudana de varivel x = + h t, podemos reescrever a

Integrando por partes, obtemos

Um nova integrao por partes nos leva seguinte expresso para a

No zemos nenhuma meno convergncia da srie. Para isto preci-

Soluo: Derivando, fcil mostrar que f (n) (x) = (1)n1 (n 1)!(1 +

Usando (2.13), o resto pode ser escrito como

Com efeito, a srie converge no intervalo 0 x 1.

m(m 1) 2 m(m 1)(m 2) 3

O termo mc2 a energia de repouso, e a energia cintica dada por

O potencial eltrico em P , no sistema CGS dado por

Figura 2.1: Dipolo eltrico.

Usando a expanso binomial, vem que

Desprezando termos de ordem igual ou superior a (a/r)2 , obtemos a equao

2.6 Critrio de Convergncia de d'Alembert

O critrio de d'Alembert estabelece que, se o limite

"# !

Resolvendo para f (x), nalmente encontramos

onde Rn conhecido como o resto da srie e denido por

ento a srie torna-se uma soma innita de termos. Temos

Esta srie conhecida como srie de Taylor.

Publicado pelo matemtico ingls Brook Taylor (1685 - 1731) em 1715.

No zemos nenhuma meno convergncia da srie. Para isto preci-

Note que no limite uv/c2 1 obtemos a expresso clssica de Galileu

Note que no limite gt/c 1, y = gt2 /2, que a expresso clssica de

Figura 3.1: Representao grca de um vetor.

O mdulo de um dado vetor u denido por

Usando-se a denio do produto vetorial, podemos mostrar que os ve-

onde, na ltima linha usamos a denio do produto escalar. Extraindo a

Os coecientes amn = cos(m, n) so conhecidos como cossenos diretores.