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
-
前の記事
javascript lodashを使って2次元配列の行列を入れ替える 2021.09.09
-
次の記事
Ruby 文字を指定した文字に置換する 2021.09.09
コメントを書く