Python可视化 函数积分图

def func(x):
    return -(x-2)*(x-8)+40

x = np.linspace(0, 10)
y = func(x)

fig, ax = plt.subplots()
plt.plot(x, y, 'r', linewidth=2)
a = 2
b = 9
ax.set_xticks([a, b])
ax.set_yticks([])
ax.set_xticklabels(['$a$', '$b$'])  #给2 9 打上a b的标签
plt.figtext(0.9, 0.1, '$x$')
plt.figtext(0.1, 0.9, '$y$')
plt.show()

ix = np.linspace(a, b) #取曲线上的点
iy = func(ix)

ixy = zip(ix, iy)  #拼接成坐标对
verts = [(a, 0)] + list(ixy)+[(b, 0)]
poly = Polygon(verts, facecolor='0.9', edgecolor='0.5')
ax.add_patch(poly)
plt.figtext(0.9, 0.1, '$x$')
plt.figtext(0.1, 0.9, '$y$')

x_math = (a+b)*0.5
y_math = 35

plt.text(x_math, y_math, r'$\int_a^b(-(x-2)*((x-8)+40)$', horizontalalignment='center') #自动对齐
plt.ylim(ymin=25)
plt.show()