Integral by Simpson's Formula

 関数 f(x) を x=a から x=b まで積分するには
b から a までを偶数個に分割しその幅を h とする。h=(b-a)/2n
i=0 から2n-2 まで1つおきに和 s を計算する。
    x=a+ih   s=s+f(x)+4f(x+h)+f(x+2h)
この s を h/3 倍すると求める積分値となる。
左図は f(x)=4/(1+x2) の 0 から 1 まで積分でπとなる。
 
  1. VB6
  2. VBScript
  3. Javascript
  4. PHP
  5. Java Applet  Source
  6. Excel
  7. C