数値積分
#
被積分関数をクラス
GSL_Function
で定義
f = GSL_Function.new{ |x|
exp(-x*x)
←
exp(-
x
2
)
}
#
オブジェクト
f
に積分用メソッド
integration_qng
を起動
#
引数は積分の下端、上端、絶対精度、相対精度
#
戻り値は
[
積分値、誤差、関数値の評価回数、積分関数の返したステータス
]
#
を要素とする配列
ans = f.integration_qng(-10.0, 10.0, 0.0, 1.0e-7)
p ans
出力:
[1.772453851, 1.967820304e-14, 87, 0]
√
p