Pythonで数値積分 〜フーリエ級数展開を例に、てかscipy.integrate.quad()かわいいよscipy.integrate.quad() - naoya_t@hatenablogはいはいまたグラフ描きたいだけのエントリですよ。てか数値積分が簡単に出来ちゃうSciPy(・∀・)ｲｲﾈ!! が from scipy.integrate import quad y, abserr = quad(f, a, b)で計算できちゃうのです。（y が積分値、abserrは推定誤差の絶対値） 基礎編 → を期待。 # -*- coding: utf-8 -*- from scipy.integrate import quad y, abserr = quad(lambda x:x*x*x, 0, 2) print "∫_0^2 x^3 dx = %g ± %g" % (y, abserr)→ ∫_0^2 x^3 dx = 4 ± 4.44089e-14 正解。推定誤差も小さい。 相手はSciPy様なので人間に正解とか言われてもおそらくあれですね。 → を期待。 # -*- cod
