使用Matplotlib库的Simpson方法

Question

帮助，请使用matplotlib库绘制图表。 +告诉我如何使用Runge方法在此处进行检查。这是原始代码。

import numpy as np

def simps(f,a,b,N=50):

    if N % 2 == 1:
        raise ValueError("N must be an even integer.")
    dx = (b - a) / N
    x = np.linspace(a, b, N + 1)
    y = f(x)
    S = dx / 3 * np.sum(y[0:-1:2] + 4 * y[1::2] + y[2::2])
    return S

print(simps(lambda x : np.sqrt(16-x**2), -2,2,10))