◆ qk51()

有限区間の積分 (固定分点) (51点ガウス・クロンロッド則)

目的: 本ルーチンは
I = [a, b]におけるfの積分値(誤差推定付き) および
J = [a, b]におけるabs(f)の積分値
を求める. ここで, fはユーザー定義サブルーチンにより与えられる関数である.
積分は51点ガウス・クロンロッド則を用いて計算される.

引数

[in]	f	被積分関数f(x)を求めるユーザー定義サブルーチンで, 次のように定義すること. double f(double x) { return f(x)の計算値 }
[in]	a	積分区間の下限.
[in]	b	積分区間の上限.
[out]	result	求められた I = [a, b]におけるfの積分値.
[out]	abserr	絶対誤差の推定値. 真の誤差に等しいかそれより大きい.
[out]	resabs	求められた J = [a, b]におけるabs(f)の積分値.
[out]	resasc	求められた [a, b]におけるabs(f - I/(b - a))の積分値.