scipy quad 定积分

• 画单位圆

import numpy as np
from matplotlib import pyplot as plt
from scipy import integrate 

# 画圆
f = lambda x : (1 - x ** 2)**0.5
x = np.linspace(-1, 1, 100)
y = f(x)
plt.figure(figsize=(5, 5))
plt.plot(x,y)
plt.plot(x,-y)

• scipy.integrate.quad 进行定积分 求pi 即圆周率

# 使用scipy.integrate.quad进行定积分 
integrate.quad(f, -1, 1) # 积分函数 + 积分区间