NEW ARTICLE
MariaDB クエリの結果をcsvに出力する
  1. HOME
  2. python
  3. python scipyで定積分を計算する

python scipyで定積分を計算する

python scipyで定積分を計算する

pythonで、ライブラリscipyのintegrate.quadを使用して、定積分を計算するサンプルコードを記述してます。pythonのバージョンは3.8.5を使用してます。

目次

環境

  • OS windows10 pro 64bit
  • python 3.8.5

scipyインストール

scipysをインストールされていない方は、pipでインストールしておきます。

pip install scipy

# Successfully installed SciPy-1.5.4

integrate.quad使い方

integrate.quadを使用すると、定積分を計算することが可能です。

from scipy import integrate

ix, err = integrate.quad(関数, 積分区間開始, 積分区間終了)

以下は、xの2乗+1の積分区間[0,5]を計算するサンプルコードとなります。

from scipy import integrate
 
# 関数
y = lambda x: x**2+1
 
# integrate.quad(関数,積分区間[0, 5])
ix, err = integrate.quad(y, 0, 5)
 
print('積分値:', ix)
# 46.66666666666667

print('推定誤差：', err)
# 5.181040781584064e-13