Sympyで定積分を計算する

前回の記事の続き

環境

• Python 3.5
• Sympy 1.4

コード

from sympy import integrate, Symbol

x = Symbol('x')
expr = 'x**2'
# 不定積分
print(integrate(expr, x))
# 定積分 (var, 下限, 上限)
print(integrate(expr, (x, 2, 4)))

以下実行結果

x**3/3
56/3

載ってた0から10の不定積分も求めてみる

from sympy import integrate, Symbol
from sympy.plotting import plot

x = Symbol('x')
expr = 'x'
print(integrate(expr, (x, 0, 10)))
plot(expr, (x, 0, 10))
=> 50

なるほど、三角形の面積として求めても10*10/2=50になるので体積が等しくなる。