【數據分析day05】積分求圓周率

積分

• 就是求曲線下的面積
  

導包： from scipy import integrate

X2 + Y2 = 1,半徑是1
pi×r**2,只要求得面積—>pi

f = lambda x: (1 - x **2)  ** 0.5

首先畫一個圓

x = np.linspace(-1, 1, 1000)
y = (1 - x **2)  ** 0.5
plt.figure(figsize=(8,8))
plt.plot(x, y, x, -y)

使用scipy.integrate進行積分，調用quad()方法
注意：integrate.quad(）返回的是一個元組，裏面有兩個元素，
前面的是半圓的面積（也就是半個pi,因爲r=1）後面的是誤差

from scipy import integrate
#half_pi是半圓的面積（也就是半個pi,因爲r=1），deviation 是誤差
half_pi, deviation = integrate.quad(f, -1, 1)
half_pi * 2

對比系統中的pi

np.pi