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 まで積分でπとなる。
- VB6
- VBScript
- Javascript
- PHP
- Java Applet Source
- Excel
- C