Fórmula de Simpson
Em análise numérica, a fórmula de Simpson (em nome de Thomas Simpson, um matemático inglês) também conhecida como regra de Simpson é uma forma de se obter uma aproximação da integral definida
A regra de Simpson baseia-se em aproximar a integral definida pela área sob arcos de parábola que interpolam a função.
Derivação da fórmula de Simpson
[editar | editar código-fonte]Vamos considerar o intervalo [a,b] com b>a e seja m o ponto central do intervalo (veja a figura ao lado). Se os pontos (a,f(a)), (m,f(m)) e (b,f(b)) não forem colineares, ou seja, se não existir uma reta que passe pelos três pontos, existirá uma única parábola com eixo vertical que passa por todos esses três pontos. De fato, a equação de qualquer parábola com eixo vertical tem a forma , onde é um polinômio quadrático.
A fórmula de Simpson faz uma aproximação de pelo polinômio de grau no máximo 2 que admite o mesmo valor de em a, b, e no ponto central . Pode-se utilizar interpolação por polinômios de Lagrange para encontrar uma expressão para essa função polinomial.
Segue, através de um cálculo simples, que:
O erro na aproximação da integral por meio da fórmula de Simpson é dado pela seguinte expressão:
com e um número entre e . Esta expressão do termo de erro significa que o método de Simpson é preciso (isto é, o termo de erro é zero) para qualquer polinomial inferior ou igual a 3 graus. Além disso, este método é a ordem 4 para qualquer função de quatro vezes continuamente diferenciável sobre [a, b].
Fórmula de Simpson Composta
[editar | editar código-fonte]Vemos que a fórmula de Simpson fornece uma boa aproximação se o intervalo de integração for pequeno, o que não acontece na maior parte dos casos. A solução óbvia é dividir o intervalo de integração em intervalos menores, aplicar a fórmula de Simpson para cada um destes e somar os resultados. Deste modo obtemos a fórmula de Simpson composta:
onde é o número de partes em que o intervalo foi dividido com par, igual ao comprimento de cada sub-intervalo e para , em particular, e
A estrutura da expressão entre parênteses: e ocorrem com coeficiente ; os com par ocorrem com coeficiente ; e os com ímpar ocorrem com coeficiente . Alternativamente, pode-se reescrever a expressão da seguinte forma:
A regra de Simpson composta admite um erro de truncamento que pode ser estimado conforme a seguinte expressão:
Onde é o comprimento do "passo", dado por
Devido à sua simplicidade de implementação e sua alta precisão, este método é utilizado pela maioria das calculadoras para cálculos aproximados de integrais função explícitas.
Exemplo
[editar | editar código-fonte]Utilizando a regra de Simpson repetida com para calcular um valor aproximado para a integral
Solução: Temos .
- ⇒
- ⇒
- ⇒
- .
Pela Regra de Simpson obtemos o resultado:
- ≈
Implementando a regra de Simpson 1/3 composta utilizando octave:
function [x,fx,v] = Simpson(a,b,n,func)
## 2013 Lesliê cardoso da silva (c)
h = (b-a)/n;
x = a:((b-a)/n):b;
func = inline(func);
for i = 1:length(x)
fx(i) = func(x(i));
end
v=fx(1)+fx(length(x));
for i= 2:length(x)-1
if mod(i,2) == 0
v+= 4*fx(i);
else
v+= 2*fx(i);
end
end
v = h/3*v;
endfunction
O exemplo acima poderia ser executado da seguinte forma no octave: [x,fx,v] = Simpson(0,2,4,'1/(1+x^4)')
Veja também
[editar | editar código-fonte]- ↑ «Confira este exemplo e faça outros com O Monitor». omonitor.io. Consultado em 1 de abril de 2016
Fontes
[editar | editar código-fonte]- Ruggiero, Lopes. Calculo Numérico. Aspectos teóricos e computação. único 2 ed. [S.l.]: Pearson
- Sperandio, Mendes e Monken (2003). Calculo Numérico. características matemáticas e computacionais dos métodos numéricos. [S.l.]: Pearson
Ligações externas
[editar | editar código-fonte]- Formules de Newton-Cotes - math-linux.com (em inglês)
- «mathworld.wolfram.com» (em inglês)