よくpythonを使って積分計算をすることがあるので、その時のためにテンプレートを用意しておきます。
目次
モジュールのインポート
Sympyを使用します。
import sympy as sp 不定積分
以下の不定積分の計算結果をターミナルに出力するプログラムを書いてみます。
$$\int tdt$$
import sympy as sp
sp.init_printing(use_unicode=True)
t = sp.symbols("t")
f = t
int = sp.integrate(f, t)
print(int)
>>>t**2/2定積分
閉区間における定積分の値を取得するプログラムを書いてみます。
$$\int_{0}^{10}e^t dt$$
import sympy as sp
sp.init_printing(use_unicode=True)
t = sp.symbols("t")
f = sp.exp(t)
int= sp.integrate(f, (t, 0, 10)).evalf()
print(int)
>>>22025.4657948067広義積分
以下のガウス積分の値を求めてみます。
$$\int_{-\infty}^{\infty}e^{-t^2}dt$$
import sympy as sp
sp.init_printing(use_unicode=True)
t = sp.symbols("t")
infty_p= str('oo')
infty_n = str('-oo')
f = sp.exp(-t**2)
int= sp.integrate(f, (t, infty_n, infty_p)).evalf()
print(int)
>>>1.77245385090552