A Regra de Simpson calcula uma aproximação para
$\displaystyle \int_a^b f(x)dx)$ pela fórmula
$$\displaystyle \int_a^b f(x)dx \approx
\frac{h}{3}\left[f(x_0)+4\sum_{i=1}^{\frac{n}{2}} f(x_{2i-1})+
2\sum_{i=1}^{\frac{n-2}{2}} f(x_{2i})+ f(x_n)\right]. $$

Continue lendo o arquivo pdf com o procedimento em Python.