Buscar

Integração_Numérica_Material

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes
Você viu 3, do total de 33 páginas

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes
Você viu 6, do total de 33 páginas

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes
Você viu 9, do total de 33 páginas

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você também pode ser Premium ajudando estudantes

Prévia do material em texto

Professores: 
Erica Boizan Batista 
Glauber Marcio Silveira Pereira 
Introdução à Integração Numérica 
Unidade 1. 
Apresentação da ideia geométrica 
TÓPICOS 
O que é a integral de uma função? 
Soma de Riemann 
 
Regra do Trapézio e Regra de Simpson 
 
Quadratura de Gauss 
Erro 
O que é uma integral? 
 
Uma integral pode ser entendida 
geometricamente como a área ou 
região compreendida entre o 
gráfico de uma função e o eixo x em 
um plano cartesiano. 
 
Neste curso trabalharemos apenas 
com integrais definidas, o que 
significa que estaremos focados em 
uma determinada seção da área, 
delimitada por retas x=a e x=b. 
Neste caso denotamos a integral 
por: Imagem elaborada pelos autores no software Geogebra 
Imagine agora um função f qualquer, 
definida em um intervalo [a,b]. Na prática 
geralmente utilizamos a primitiva da 
função f para calcular sua integral, mas 
podemos ter alguns casos em que isso não 
é possível. 
 
Por exemplo, essa função f pode ter 
primitiva desconhecida, ou seja, não 
sabemos como calcular diretamente sua 
integral. Outra possibilidade é que a 
primitiva seja muito difícil de ser calculada. 
Também é possível, dependendo do caso, 
que a função estudada seja conhecida 
apenas em alguns pontos. Então como 
fazer para encontrar a área que queremos? 
 
Uma alternativa é tentar encontrar uma 
aproximação para a área procurada usando 
formas geométricas cujas áreas nós já 
conhecemos. 
Imagem retirada do site Pixabay: https://pixabay.com/pt/ 
illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando 
-3978394/; autor(a): Peggy_Marco 
MOTIVAÇÃO 
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/illustrations/ponto-de-interroga%C3%A7%C3%A3o-perguntando-3978394/
https://pixabay.com/pt/users/Peggy_Marco-1553824/
Perceba que esse é um raciocínio 
semelhante ao que usamos quando 
queremos calcular a área de um polígono 
irregular. 
 
Neste caso nós subdividimos nosso 
polígono em áreas conhecidas, e depois de 
calcular as áreas das formas menores basta 
somar tudo, certo? 
 
Apesar do nosso problema ser um pouco 
mais complexo a estratégia de resolução 
não é muito diferente. Imagem elaborada pelos autores no software Geogebra 
Somas de Riemann 
à esquerda 
Imagens elaboradas pelos autores no software Geogebra. 
A ideia da Soma de Riemann é fazer a 
aproximação da área sob o gráfico de uma 
função usando retângulos. Em particular, a 
Soma de Riemann à esquerda define a altura 
desses retângulos (geralmente de bases de 
igual largura) usando o valor e f no ponto da 
extremidade esquerda da sua base. 
Somas de Riemann 
à direita 
Imagens elaboradas pelos autores no software Geogebra. 
A Soma de Riemann à direita define a altura 
desses retângulos (geralmente de bases de 
igual largura) usando o valor e f no ponto da 
extremidade direita da sua base. 
Somas de Riemann 
no ponto médio 
Imagens elaboradas pelos autores no software Geogebra. 
A Soma de Riemann no ponto médio define a 
altura de desses retângulos (geralmente de 
bases de igual largura) usando o valor e f no 
ponto médio da sua base. 
Regra do Trapézio 
Imagem elaborada pelos autores no software Geogebra. 
A ideia desta regra é pegar os dois 
pontos extremos da parte do gráfico 
que desejamos calcular a área e traçar 
um seguimento de reta, construindo 
assim um trapézio. 
Regra do Trapézio 
Composta 
Imagens elaboradas pelos autores no software Geogebra. 
De forma semelhante à Soma de Riemann 
podemos subdividir o intervalo [a,b] no eixo x 
em intervalos (geralmente de mesmo tamanho) 
e usar seus extremos para construir trapézios, 
que servirão para aproximar a área que 
queremos. 
Regra de Simpson 
Imagem elaborada pelos autores no software Geogebra. 
A ideia desta regra é, ao invés de retas, 
usar polinômios que passem por pontos 
do gráfico, e calcular a integral desses 
polinômios para conseguir uma 
aproximação da área que queremos. 
Regra de Simpson 
Composta 
Imagens elaboradas pelos autores no software Geogebra. 
De forma semelhante à Regra do Trapézio 
Composta podemos subdividir o intervalo 
[a,b] no eixo x em intervalos (geralmente de 
mesmo tamanho) e usar vários polinômios, 
que servirão para aproximar a área que 
queremos. 
Calma! A palavra erro assusta um pouco, 
mas o erro aqui faz parte do processo de 
aproximação, ele é a diferença entre o 
valor obtido (aproximado) e o valor exato 
da integral. 
 
Ao contrário da forma clássica de calcular 
integrais (usando a primitiva), que nos 
permite obter soluções exatas para o 
cálculo da integral, os métodos numéricos 
produzem, em geral, apenas soluções 
aproximadas. 
 
Por isso, antes da utilização de qualquer 
método numérico é necessário decidir qual 
a precisão dos cálculos que você precisa, já 
que métodos diferentes apresentam erros 
diferentes. 
Imagem retirada do site Pixabay: https://pixabay.com/ 
pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/; 
autor(a): OpenClipart-Vectors 
Erro 
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/vectors/alerta-parar-aviso-%C3%ADcone-erro-146730/
https://pixabay.com/pt/users/OpenClipart-Vectors-30363/
https://pixabay.com/pt/users/OpenClipart-Vectors-30363/
https://pixabay.com/pt/users/OpenClipart-Vectors-30363/
Quadratura de Gauss 
Imagem elaborada pelos autores no software Geogebra. 
A quadratura de Gauss é um método 
que permite escolhermos o tamanho 
das subdivisões que faremos do 
intervalo [a,b] no eixo x, de forma a 
obtermos erros muito pequenos. 
Ela pode ser aplicada a qualquer um dos 
métodos apresentados anteriormente, 
e para o caso em que f é um polinômio 
ela garante erro igual a zero. 
Introdução à Integração Numérica 
Resumo da apresentação 
 
Apresentamos aqui a ideia geométrica 
dos métodos numéricos mais usados 
para aproximação de uma integral 
definida a partir de somas finitas. 
 
No próximo materialveremos a 
formalização desses métodos. 
CRÉDITOS 
Autores 
Colaboração 
Glauber Marcio Silveria Pereira 
Graduado em Matemática 
Mestre em Bioestatística 
Doutor em Estatística 
 
Erica Boizan Batista 
Graduada em Matemática 
Mestre em Matemática 
Doutora em Matemática 
 
REFERÊNCIAS 
BARROSO, L. C., Cálculo Numérico. S. Paulo: HARBRA, 1987. 
 
BURDEN, R. L.; FAIRES, J. D. Análise Numérica. Cengage Learning, Tradução da 8. Ed. Americana, 
2008. 
 
FRANCO, N. B. Cálculo Numérico. Pearson Prentice Hall, 2007. 
 
HOHEBWARTER, Markus. GeoGebra clássico 6: Software de geometria dinâmica. Disponível em: 
<https://www.geogebra.org/download>. Acesso em: 23 de agosto de 2020. 
 
REAMAT - Cálculo Numérico. REAMAT, 2020. Cálculo Numérico - Versão GNU Octave. Disponível em: 
<https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html>. Acesso em: 09 de junho de 
2020. 
 
RUGGIERO, M. A. G., Cálculo Numérico: Aspectos Teóricos e Computacionais, São Paulo, McGraw-
Hill, 1998. 
 
https://www.geogebra.org/download
https://www.geogebra.org/download
https://www.geogebra.org/download
https://www.geogebra.org/download
https://www.geogebra.org/download
https://www.geogebra.org/download
https://www.geogebra.org/download
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
 
Profa. Érica Boizan Batista 
Prof. Glauber Márcio 
Silveira Pereira 
 
 
Unidade 2. 
Integração Numérica: 
formalização dos 
conceitos 
 
 
Integração Numérica1 
1. Introdução. 
Uma função 𝑓: [𝑎, 𝑏] → ℝ pode ter primitiva desconhecida, ou seja, que não 
sabemos diretamente sua integral, ∫ 𝑓(𝑥)𝑑𝑥
𝑏
𝑎
 . Outra possibilidade é que a primitiva seja 
muito difícil de ser calculada. Também é possível, dependendo do caso, que a função 
estudada seja conhecida apenas em alguns pontos. No estudo de Cálculo Numérico existem 
métodos que fazem aproximações da primitiva através de somas finitas de áreas. A ideia 
destes métodos é subdividir o domínio da função [𝑎, 𝑏] e calcular de forma aproximada 
∫ 𝑓(𝑥)𝑑𝑥 ≅
𝑏
𝑎
∑ 𝑎𝑖𝑓(𝑥𝑖)
𝑛
𝑖=0
, com 𝑥𝑖 o 𝑖-ésimo ponto, 𝑎𝑖 o 𝑖-ésimo peso. 
Neste curso apresentamos os métodos mais conhecidos, são eles Soma de Riemann; 
Newton Cotes simples e composta, com regra de Simpson e regra dos trapézios, e o último 
método a Quadratura de Gauss. 
 
2. Soma de Riemann 
O primeiro método apresentado é conhecido pelos alunos que já cursaram ou 
estudaram Cálculo diferencial e Integral. O cálculo da integral definida limitada por um 
intervalo [𝑎, 𝑏] é a área limitada pelas retas 𝑦 = 0, 𝑥 = 𝑎, 𝑥 = 𝑏 e uma função 𝑓: [𝑎, 𝑏] →
ℝ , como vemos no gráfico abaixo: 
 
Gráfico 1: Área limitada pelas retas 𝑦 = 0, 𝑥 = 𝑎, 𝑥 = 𝑏 e a função 𝑓. 
 
Fonte: Autoria própria. 
 
1 Érica B. Batista doutora em Matemática e Glauber M. S. Pereira doutor em Estatística. 
 
 
 
Desta forma a Soma de Riemann trata-se de uma soma finita de 𝑛 áreas de 
retângulos na tentativa de termos uma aproximação da área já descrita. Primeiro 
subdividimos o domínio da função 𝑓 , o intervalo [𝑎, 𝑏], em 𝑛 ∈ ℕ subintervalos não 
necessariamente equidistantes, ou seja, que não tem necessariamente a mesma distância. 
Definimos os pontos de um conjunto finito contido no domínio que se chama uma partição 
de [𝑎, 𝑏], sendo os pontos definidos da seguinte forma 𝑎 = 𝑥0 ≤ 𝑥1 ≤ ⋯ ≤ 𝑥𝑛 = 𝑏 , Seja 
ℎ𝑖 = 𝑥𝑖 − 𝑥𝑖−1, os pontos são 𝑥𝑖 = 𝑎 + 𝑖ℎ𝑖 para cada 𝑖 = 0,1, … , 𝑛. Note que ℎ𝑖 será a base 
do 𝑖-ésimo retângulo. Agora só falta definirmos as alturas dos retângulos. As alturas podem 
ser definidas de diversas formas. Podemos tomar o máximo ou mínimo 𝑓(𝑥) no 𝑖-ésimo 
intervalo por exemplo. Aqui apresentamos três possibilidades, tomamos 𝑓(𝑥𝑖) e definimos 
a soma de Riemann a esquerda (SRE), tomamos 𝑓(𝑥𝑖) e definimos a soma de Riemann a 
esquerda (SRD), tomamos 𝑓(
𝑥𝑖−𝑥𝑖−1
2
) e definimos a soma de Riemann ponto médio (SRPM). 
Mostraremos graficamente os três métodos. 
 
2.1. Soma de Riemann a esquerda (SRE): 
Temos a seguinte aproximação ∫ 𝑓(𝑥)𝑑𝑥 ≅ 
𝑏
𝑎
∑ ℎ𝑖𝑓(𝑥𝑖)
𝑛−1
𝑖=0
= ℎ1𝑓(𝑥1) + ℎ2𝑓(𝑥2) +
⋯ + ℎ𝑛−1𝑓(𝑥𝑛+1). Vemos graficamente o método com 𝑓: [0,4] → ℝ, 𝑓(𝑥) = sin 𝑥 + 1 para 
𝑛 = 2, 𝑛 = 4 e 𝑛 = 5 : 
 
Gráfico 2: SRE para 𝑛 = 2, 𝑛 = 4 e 𝑛 = 5. 
 
Fonte: Autoria própria. 
 
2.2. Soma de Riemann a direita (SRD): 
Temos a seguinte aproximação ∫ 𝑓(𝑥)𝑑𝑥 ≅
𝑏
𝑎
∑ ℎ𝑖𝑓(𝑥𝑖+1)
𝑛−1
𝑖=0
. Vemos graficamente 
o método com 𝑓: [0,4] → ℝ, 𝑓(𝑥) = sin 𝑥 + 1 para 𝑛 = 2, 𝑛 = 4 e 𝑛 = 6 : 
 
Gráfico 3: SRD para 𝑛 = 2, 𝑛 = 4 e 𝑛 = 6. 
 
Fonte: Autoria própria. 
 
 
 
2.3. Soma de Riemann ponto médio (SRPM): 
Temos a seguinte aproximação ∫ 𝑓(𝑥)𝑑𝑥 ≅
𝑏
𝑎
∑ ℎ𝑖𝑓(𝜉𝑖)
𝑛−1
𝑖=0
, sendo 𝜉𝑖 = 
𝑥𝑖−𝑥𝑖−1
2
, para 
cada 𝑖 = 0,1, … , 𝑛. graficamente o método com 𝑓: [0,4] → ℝ, 𝑓(𝑥) = sin 𝑥 + 1 para 
𝑛 = 2, 𝑛 = 4 e 𝑛 = 6 : 
 
Gráfico 4: SRPM para 𝑛 = 2, 𝑛 = 4 e 𝑛 = 6. 
 
Fonte: Autoria própria. 
 
3. Método Newton Cotes. 
Método com o intervalo [𝑎, 𝑏] , ou seja, o domínio da função 𝑓 é dividido em 
𝑛 subintervalos equidistantes. Assim temos 𝑥0 = 𝑎; 𝑥𝑛 = 𝑏 como estremos do intervalo e 
ℎ =
𝑏−𝑎
𝑛
 como passo e os pontos 𝑥𝑖 = 𝑎 + 𝑖ℎ para cada 𝑖 = 0,1, … , 𝑛. Nos baseamos no 
polinômio interpolador de Lagrange 𝑃𝑛 = ∑ 𝐿𝑖𝑓(𝑥𝑖) =
𝑛
𝑖=0
𝐿0𝑓(𝑥0) + ⋯ + 𝐿𝑛𝑓(𝑥𝑛), sendo 
𝐿𝑖 = ∏
𝑥−𝑥𝑗
𝑥𝑖−𝑥𝑗
𝑛
𝑗=0:𝑗≠𝑖 . Temos a aproximação do cálculo ∫ 𝑓(𝑥)𝑑𝑥 = 
𝑏
𝑎
∑ 𝑎𝑖𝑓(𝑥𝑖) +
𝑛
𝑖=0
𝐸(𝑓), 
em que 𝑎𝑖 = ∫ 𝐿𝑖(𝑥)𝑑𝑥,
𝑏
𝑎
 erro 𝐸(𝑓) = 
1
(𝑛+1)!
 ∫ ∏ (𝑥 − 𝑥𝑖)
𝑛
𝑖=0 𝑓
(𝑛+1)(𝜉(𝑥))𝑑𝑥 
𝑏
𝑎
 , para 
𝜉(𝑥) ∈ (𝑎, 𝑏). As duas principais regras são a dos trapézios e de Simpson, em que é feita 
uma interpolação (aproximação da função pelo polinômio de Lagrange). 
 
3.1 Regra do trapézio: 
A ideia desta regra é pegar os dois pontos extremos da parte do gráfico que 
desejamos calcular a área e traçar um seguimento de reta, construindo assim um trapézio. 
Vemos graficamente o método com 𝑓: [0,4] → ℝ, 𝑓(𝑥) = sin 𝑥 + 1: 
 
 
 
 
 
 
 
 
 
Gráfico 5: Regra do Trapézio. 
 
Fonte: Autoria própria. 
 
Matematicamente temos aproximação da integral com 𝑥0 = 𝑎, 𝑥1 = 𝑏, ℎ = 𝑏 − 𝑎,
𝑃1(𝑥) = 𝐿0𝑓(𝑥0) + 𝐿1𝑓(𝑥1) =
𝑥−𝑥1
𝑥0−𝑥1
𝑓(𝑥0) +
𝑥−𝑥0
𝑥1−𝑥0
𝑓(𝑥1) , temos ∫ 𝑓(𝑥)𝑑𝑥 = 
𝑏
𝑎
ℎ
2
(𝑓(𝑎) +
𝑓(𝑏)) −
𝑏−𝑎
12
ℎ2𝑓′′(𝜉(𝑥)). 
 
3.2 Regra de Simpson: 
 
A ideia desta regra é pegar os dois pontos extremos e o ponto médio da parte do 
gráfico que desejamos calcular a área e traçar uma curva polinomial, cuja área será mais fácil 
de calcular e nos dará uma aproximação para a área que queremos. Vemos graficamente o 
método com 𝑓: [0,4] → ℝ, 𝑓(𝑥) = sin 𝑥 + 1: 
 
Gráfico 6: Regra de Simpson. 
 
Fonte: Autoria própria. 
 
Matematicamente temos aproximação da integral com 𝑥0 = 𝑎, 𝑥2 = 𝑏, ℎ =
𝑏−𝑎
2
, 𝑥1 = 𝑥0 + ℎ, 𝑃2(𝑥) =
(𝑥−𝑥1)(𝑥−𝑥2)
(𝑥0−𝑥1)(𝑥0−𝑥2)
𝑓(𝑥0) +
(𝑥−𝑥0)(𝑥−𝑥2)
(𝑥1−𝑥0)(𝑥1−𝑥2)
𝑓(𝑥1) +
(𝑥−𝑥0)(𝑥−𝑥1)
(𝑥2−𝑥0)(𝑥2−𝑥1)
𝑓(𝑥2), temos ∫ 𝑓(𝑥)𝑑𝑥 = 
𝑏
𝑎
ℎ
2
(𝑓(𝑎) + 𝑓(𝑏)) −
𝑏−𝑎
12
ℎ2𝑓′′(𝜉(𝑥)). 
 
 
 
3.3 Integração Composta: 
Para obter melhorprecisão de aproximação subdividimos o intervalo [𝑎, 𝑏] em 𝑛 
subintervalos, com 𝑛 ∈ ℕ, e aplicamos a cada subintervalo as regras já apresentadas. 
 
3.3.1 Regra do trapézio Composta: 
Temos aproximação da integral com 𝑥0 = 𝑎, 𝑥𝑛 = 𝑏, ℎ =
𝑏−𝑎
2
, 𝑥𝑗 = 𝑥0 + 𝑗ℎ, para 
𝑗 = 0,1, … , 𝑛, temos ∫ 𝑓(𝑥)𝑑𝑥 = 
𝑏
𝑎
ℎ
2
(𝑓(𝑎) + 2 ∑ 𝑓(𝑥𝑗) +
𝑛−1
𝑗=1
 𝑓(𝑏)) −
𝑏−𝑎
12
ℎ2𝑓′′(𝜉(𝑥)) . 
Vemos graficamente o método com 𝑓: [0,4] → ℝ, 𝑓(𝑥) = sin 𝑥 + 1 para 𝑛 = 2, 𝑛 = 4 e 𝑛 = 6 : 
 
Gráfico 7: Regra do trapézio composta para 𝑛 = 2, 𝑛 = 4 e 𝑛 = 6. 
 
 
Fonte: Autoria própria. 
 
3.3.2 Regra de Simpson Composta: 
Temos aproximação da integral com 𝑥0 = 𝑎, 𝑥𝑛 = 𝑏, ℎ =
𝑏−𝑎
2
, 𝑥𝑗 = 𝑥0 + 𝑗ℎ, para 
𝑗 = 0,1, … , 𝑛, temos ∫ 𝑓(𝑥)𝑑𝑥 = 
𝑏
𝑎
ℎ
3
(𝑓(𝑎) + 2 ∑ 𝑓(𝑥2𝑗)
𝑛/2−1
𝑗=1
+ 4 ∑ 𝑓(𝑥2𝑗−1)
𝑛/2
𝑗=1
+
 𝑓(𝑏)) −
𝑏−𝑎
180
ℎ4𝑓(4)(𝜉(𝑥)). Vemos graficamente o método com 𝑓: [1,4] → ℝ, 𝑓(𝑥) = sin 𝑥 + 1 
para 𝑛 = 2, 𝑛 = 4 e 𝑛 = 6 : 
 
Gráfico 8: Regra de Simpson composta para 𝑛 = 2, 𝑛 = 4 e 𝑛 = 6. 
 
 
Fonte: Autoria própria. 
 
 
 
4 Quadratura de Gauss. 
Método que otimiza a área em função da oscilação da função 𝑓. Segue graficamente 
a regra dos trapézios com duas funções como exemplo comparando o método Newton 
Cotes com a quadratura de Gauss. 
Na aproximação da quadratura de Gauss da forma composta 
∫ 𝑓(𝑥)𝑑𝑥 ≅
𝑏
𝑎
∑ 𝑐𝑖𝑓(𝑥𝑖)
𝑛
𝑖=1
 o método tem 2𝑛 parâmetros, que são os nós 𝑥1, … , 𝑥𝑛 e os 
coeficientes 𝑐1, … , 𝑐𝑛 . 
Definimos o polinômio de Legendre {𝑃0, … , 𝑃𝑛} com as propriedades para 𝑃𝑗(𝑥): 
I. Para cada 𝑗 ∈ ℕ, é polinômio mônico, ou seja, o coeficiente dominante é 1; 
II. ∫ 𝑃(𝑥)𝑃𝑗(𝑥)𝑑𝑥 = 0
1
−1
 , ou seja, ortogonal para polinômio com grau menor que 𝑗; 
III. Seja o produto escalar 〈𝑓, 𝑔〉 = ∫ 𝑓(𝑥)𝑔(𝑥)𝑑𝑥
1
−1
 e definimos seguindo a 
sequencia 𝑃0(𝑥) = 1, 𝑃1(𝑥) = (𝑥 − 𝛼0)𝑃0(𝑥) = 𝑥, 𝑃𝑗+1(𝑥) = (𝑥 − 𝛼𝑗)𝑃𝑗(𝑥) −
𝛽𝑗𝑃𝑗−1(𝑥) , para cada 𝑗 = 2,3, … , em que 𝛼𝑗 =
〈𝑥𝑃𝑗(𝑥),𝑃𝑗(𝑥)〉
〈𝑃𝑗(𝑥),𝑃𝑗(𝑥)〉
=
∫ 𝑥(𝑃𝑗(𝑥))
2𝑑𝑥
1
−1
∫ (𝑃𝑗(𝑥))
2𝑑𝑥
1
−1
 e 
𝛽𝑗 =
〈𝑃𝑗(𝑥),𝑃𝑗(𝑥)〉
〈𝑃𝑗−1(𝑥),𝑃𝑗−1(𝑥)〉
=
∫ (𝑃𝑗(𝑥))
2𝑑𝑥
1
−1
∫ (𝑃𝑗−1(𝑥))
2𝑑𝑥
1
−1
. 
Exemplo: 𝑃2(𝑥) = (𝑥 − 𝛼2)𝑃1(𝑥) − 𝛽2𝑃0(𝑥) = (𝑥 − 0)𝑥 −
1
3
= 𝑥2 −
1
3
. 
O polinômio de Legendre é usado no principal resultado da quadratura de Gauss. 
Teorema: Sejam 𝑥1, … , 𝑥𝑛 raízes de 𝑃𝑛(𝑥), 𝑛-ésimo polinômio de Legendre e definimos os 
coeficientes 𝑐𝑖 = ∫ ∏
𝑥−𝑥𝑗
𝑥𝑖−𝑥𝑗
𝑛
𝑗=1:𝑗≠𝑖 𝑑𝑥
1
−1
, para cada 1, … , 𝑛. Se 𝑃(𝑥) é um polinômio qualquer 
de grau menor que 2𝑛, então ∫ 𝑃(𝑥)𝑑𝑥 =
1
−1
∑ 𝑐𝑖𝑃(𝑥𝑖)
𝑛
𝑖=1
. 
Observação: Para um domínio qualquer [𝑎, 𝑏] podemos aplicar os resultados apresentados 
nesta seção usando a seguinte mudança de variável: 𝑥 =
1
2
((𝑏 − 𝑎)𝑡 + 𝑎 + 𝑏) ⟺ 𝑡 =
2𝑥−𝑎−𝑏
𝑏−𝑎
, logo ∫ 𝑓(𝑥)𝑑𝑥 = 
𝑏
𝑎
∫ 𝑓 (
1
2
((𝑏 − 𝑎)𝑡 + 𝑎 + 𝑏))
𝑏−𝑎
2
𝑑𝑡 
1
−1
. 
 
5 Referências. 
BARROSO, L. C., Cálculo Numérico. S. Paulo: HARBRA, 1987. 
BURDEN, R. L. ; FAIRES, J. D.. Análise Numérica, Cengage Learning, Tradução da 8. Ed. 
Americana, 2008. 
FRANCO, N. B. Cálculo Numérico. Pearson Prentice Hall, 2007. 
 
 
http://poca.ufscar.br/ 
HOHEBWARTER, Markus. GeoGebra clássico 6: Software de geometria dinâmica. Disponível 
em: <https://www.geogebra.org/download>. Acesso em: 23 de agosto de 2020. 
RUGGIERO, M. A. G., Cálculo Numérico: Aspectos Teóricos e Computacionais, São Paulo, 
McGraw-Hill, 1998. 
REAMAT - Cálculo Numérico. REAMAT, 2020. Cálculo Numérico - Versão GNU Octave. 
Disponível em: <https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html>. Acesso 
em: 09 de junho de 2020. 
 
 
 
http://poca.ufscar.br/
https://www.geogebra.org/download
https://www.ufrgs.br/reamat/CalculoNumerico/livro-oct/in.html
 
 
http://poca.ufscar.br/ 
 
http://poca.ufscar.br/
28/10/2022 10:26 Questionário: Revisão da tentativa
https://cursos.poca.ufscar.br/mod/quiz/review.php?attempt=621445&cmid=3096 1/7
Página inicial Meus cursos Integra_Num Atividade Questionário
Iniciado em quinta-feira, 22 set 2022, 07:25
Estado Finalizada
Concluída em quinta-feira, 22 set 2022, 07:27
Tempo
empregado
2 minutos 36 segundos
Avaliar 10,00 de um máximo de 10,00(100%)
Questão 1
Correto
Atingiu 1,00 de 1,00
Questão 2
Correto
Atingiu 1,00 de 1,00
Considere uma função ƒ definida em um intervalo [a,b]. O método que usa a soma de retângulos para aproximar a área
definida pela integral é:
Escolha uma opção:
Regra de Romberg 
Regra de Simpson 
Soma de Riemann 
Regra dos Trapézios 
∫abf(x)dxf(x)dx∫ ba
Sua resposta está correta.
Considere uma função ƒ definida em um intervalo [a,b]. O método que usa polinômios para aproximar a área definida pela
integral é: 
Escolha uma opção:
Soma de Riemann à direita 
Soma de Riemann à esquerda 
Regra de Romberg 
Regra de Simpson
∫abf(x)dxf(x)dx∫ ba
Sua resposta está correta.
https://cursos.poca.ufscar.br/
https://cursos.poca.ufscar.br/course/view.php?id=73&section=0
https://cursos.poca.ufscar.br/course/view.php?id=73&section=3
https://cursos.poca.ufscar.br/mod/quiz/view.php?id=3096
28/10/2022 10:26 Questionário: Revisão da tentativa
https://cursos.poca.ufscar.br/mod/quiz/review.php?attempt=621445&cmid=3096 2/7
Questão 3
Correto
Atingiu 1,00 de 1,00
Questão 4
Correto
Atingiu 1,00 de 1,00
Considere uma função ƒ definida em um intervalo [a,b]. Qual é o método que subdivide o intervalo [a,b] de forma equidistante
e usa o polinômio de Lagrange? 
Escolha uma opção:
Newton Cotes 
Soma de Riemann 
Quadratura de Gauss 
Método de Taylor 
Sua resposta está correta.
Considere uma função ƒ definida em um intervalo [a,b]. Utilizando o método Newton Cotes, qual é a regra que utiliza apenas o
polinômio de Lagrange P (x). 
Escolha uma opção:
Regra de Simpson 
Regra de Romberg 
Soma de Riemann no ponto médio 
Regra dos Trapézios 
1
Sua resposta está correta.
28/10/2022 10:26 Questionário: Revisão da tentativa
https://cursos.poca.ufscar.br/mod/quiz/review.php?attempt=621445&cmid=3096 3/7
Questão 5
Correto
Atingiu 1,00 de 1,00
Observe a figura abaixo e indique qual regra está sendo utilizada para encontrar uma aproximação de . 
 
Escolha uma opção:
Soma de Riemann no ponto médio 
Soma de Riemann à direita 
Soma de Riemann à esquerda 
Regra dos Trapézios 
∫abf(x)dxf(x)dx∫ ba
Sua resposta está correta.
28/10/2022 10:26 Questionário: Revisão da tentativa
https://cursos.poca.ufscar.br/mod/quiz/review.php?attempt=621445&cmid=3096 4/7
Questão 6
Correto
Atingiu 1,00 de 1,00
Questão 7
Correto
Atingiu 1,00 de 1,00
Observe a figura abaixo e indique qual método está sendo utilizado para encontrar uma aproximação de . 
 
Escolha uma opção:
Quadratura adaptativa com regra do Trapézio 
Quadratura de Gauss com regra do Trapézio 
Newton Cotes com regra do Trapézio 
Método de Taylor com regra do Trapézio 
∫abf(x)dxf(x)dx∫ ba
Sua resposta está correta.
De acordo com o material apresentado, o que significa erro em termos de Integração Numérica? 
Escolha uma opção:
Usar o método errado 
O pesquisador errar o cálculo 
É a diferença entre o valor obtido (aproximado) e o valor exato da integral 
O software errar o cálculo 
Sua resposta está correta.
28/10/2022 10:26 Questionário: Revisão da tentativa
https://cursos.poca.ufscar.br/mod/quiz/review.php?attempt=621445&cmid=3096 5/7
Questão 8
Correto
Atingiu 1,00 de 1,00
Questão 9
Correto
Atingiu 1,00 de 1,00
Qual é o método que para polinômios garante erro igual a zero? 
Escolha uma opção:
Newton Cotes 
Quadratura de Gauss 
Soma de Riemann à direita 
Regra do Trapézio 
Sua resposta está correta.
Algum dos métodos apresentados neste curso não usa soma finita de áreas para aproximação do cálculo de integrais? Qual? 
Escolha uma opção:
Quadratura de Gauss 
Regra de Simpson 
Soma de Riemann no ponto médioTodos os métodos apresentados usam soma finita de áreas 
Sua resposta está correta.
28/10/2022 10:26 Questionário: Revisão da tentativa
https://cursos.poca.ufscar.br/mod/quiz/review.php?attempt=621445&cmid=3096 6/7
Questão 10
Correto
Atingiu 1,00 de 1,00
Assinale qual das alternativas apresenta o método de Simpson Composto.
Escolha uma opção:
Sua resposta está correta.
28/10/2022 10:26 Questionário: Revisão da tentativa
https://cursos.poca.ufscar.br/mod/quiz/review.php?attempt=621445&cmid=3096 7/7

Outros materiais