台形公式を用いた数値積分を行います.

$ f(x) = x^2 - 2x + 2 $ の積分値を台形公式によって求めます.

• 積分区間は $ [a, b] $ とします.
• 実際に計算する式は以下の通りです. ここで $ h $ は $ h = frac{b-a}{n} $ です.

\[

   \int_{a}^{b} f(x) dx =
     \frac{h}{2} \left\{f(a) + f(b)\right\} + h\sum_{i=1}^{n-1} f(a+ih)

\]

このサンプルプログラムは www.gfd-dennou.org/library/dcmodel にて入手できます.