Simpson 1/3
resultado
função a ser integrada
limite inferior de integração
limite superior de integração
número de pontos [n+1], ou número do intervalos[n]
Trata-se do caso mais simples de Fórmula de Newton-Cotes fechada
res = cpc_simpson_13( fx, 0, 0.8, 4 );