Python数值积分：使用梯形公式计算任意函数积分

我们可以使用数值积分方法来实现任意函数的积分运算。以下是一种基于梯形公式的数值积分方法的Python代码：

def integrate(f, a, b, n=1000):
    '梯形公式的数值积分方法'
    :param f: 被积函数
    :param a: 积分下限
    :param b: 积分上限
    :param n: 积分区间数量
    :return: 积分结果
    dx = (b - a) / n
    x = [a + i * dx for i in range(n + 1)]
    y = [f(x_i) for x_i in x]
    return (dx / 2) * (y[0] + y[-1] + 2 * sum(y[1:-1]))


# 示例：计算 f(x) = x^2 在 [0, 1] 区间的积分
def f(x):
    return x ** 2

result = integrate(f, 0, 1)
print(result)

输出结果为：

0.33333350000000004

注意：数值积分方法的精度受到区间数量 n 的影响，一般来说 n 越大，积分结果越接近真实值，但计算时间也会增加。在实际应用中需要根据实际情况选择合适的 n 值。