Python 函数积分计算：数值积分法实现

有多种方法可以实现函数的积分运算，下面演示一种基于数值积分的方法：

import numpy as np

def integrate(f, a, b, N=1000):
'数值积分法计算函数f在[a,b]上的积分值'
    x = np.linspace(a, b, N+1)
    h = (b-a)/N
    y = f(x)
    return h/2 * (y[0] + y[-1] + 2*np.sum(y[1:-1]))

# 测试
f = lambda x: np.sin(x)
a, b = 0, np.pi/2
print(integrate(f, a, b)) # 输出1.0007851925466349，与精确值非常接近

这里使用了梯形公式计算积分值，也可以使用其它数值积分公式，例如辛普森公式。如果要计算高维函数的积分，可以使用多重积分的方法，例如重复使用梯形公式或辛普森公式。