Algoritmo Regra de Simpson 1/3 - Integral Numérica

ENTRADA: os limites a e b, um inteiro par n
SAIDA: Aproximação XI para a integral
passo 1 faça h =(b-a)/n
passo 2 faça XI0 = f(a) + f(b)
              XIl = 0; (Soma de f(x2i-1))
             XI2 = O. (Soma de f (X2i))
passo 3 para i = 1, ... ,n - 1 faça os passos 4 e 5.
                   Passo 4 faça X = a + ih.
                   Passo 5 Se i é par faça XI2 = X/2 + f(X)
                                                 senão faça XIl = XIl + f(X).
Passo 6 faça XI = h(XI0 + 2· XI2 + 4· XIl)/3.
Passo 7 SAÍDA (Xl);
Pare