qag qagi qk15

qk15() def qk15 ( f  , a  , b    ) 有限区間の積分 (固定分点) (15点ガウス・クロンロッド則) 目的 本ルーチンは   I = [a, b]におけるfの積分値(誤差推定付き) および   J = [a, b]におけるabs(f)の積分値 を求める. ここで, fはユーザー定義サブルーチンにより与えられる関数である. 積分は15点ガウス・クロンロッド則を用いて計算される. 戻り値 (result, abserr) result (float): 求められた I = [a, b]におけるfの積分値. abserr (float): 絶対誤差の推定値. 真の誤差に等しいかそれより大きい. 引数 [in] f 被積分関数f(x)を求めるユーザー定義サブルーチンで, 次のように定義すること. _CODE def f(x): return f(x)の関数値 _ENDCODE [in] a 積分区間の下限. [in] b 積分区間の上限. 出典 SLATEC (QUADPACK) 使用例 次の定積分を求める. ∫ 1/(1 + x^2) dx [0, 4] (= atan(4)) def f(x): return 1/(1 + x**2) def TestQk15(): a = 0 b = 4 s, abserr = qk15(f, a, b) print(s, abserr) 実行結果 >>> TestQk15() 1.3258176613637855 0.0014827239412162237