python辛普森积分

Python中辛普森积分的实现可以通过定义一个函数来完成。下面是一个示例代码：

import numpy as np

def simpson(f, a, b, n=100):
    """
    f: 要积分的函数
    a: 积分下界
    b: 积分上界
    n: 划分区间个数，默认为100
    """
    if n % 2 == 1:
        raise ValueError("n must be an even number")
    
    h = (b - a) / n
    x = np.linspace(a, b, n + 1)
    y = f(x)
    
    return h/3 * (y + y[-1 + 4*np.sum(y