• The code example computes \( \pi \) using the trapezoidal rules.
• The trapezoidal rule

$$ I=\int_a^bf(x) dx\approx h\left(f(a)/2 + f(a+h) +f(a+2h)+\dots +f(b-h)+ f(b)/2\right). $$

Click on this link for the full program.