[2-1] 定積分を計算するモジュールをインポートします。
In [1]:
import numpy as np
from scipy.integrate import quad
[2-2] 円弧を表す関数を定義します。
In [2]:
def circle(x):
y = np.sqrt(1-x*x)
return y
[2-3] 定積分を計算して、円の面積を求めます。
quadは、関数circleを0〜1の範囲で定積分した結果をresに保存します。（errは数値計算で発生した誤差の範囲を示します。）
In [3]:
res, err = quad(circle, 0, 1)
print res*4