Python 实现复合辛普森求积公式数值积分

复合辛普森求积公式是数值积分中常用的方法，它将积分区间分成若干个小区间，每个小区间内用辛普森公式进行积分，最终将所有小区间的积分值加起来得到整个区间的积分值。/n/n下面是 Python 代码实现复合辛普森求积公式：/n/npython/nimport math/n/ndef f(x):/n return math.sin(x)/n/ndef composite_simpson(a, b, n):/n h = (b - a) / n/n x = [a + i * h for i in range(n + 1)]/n s = f(a) + f(b)/n for i in range(1, n):/n if i % 2 == 0:/n s += 2 * f(x[i])/n else:/n s += 4 * f(x[i])/n return h / 3 * s/n/na = 0/nb = math.pi / 2/nn = 10/nintegral = composite_simpson(a, b, n)/nprint('Integral value:', integral)/n/n/n在上述代码中，函数 f(x) 是被积函数，composite_simpson(a, b, n) 是复合辛普森求积公式的实现函数，其中 a 和 b 是积分区间的上下限，n 是将积分区间分成的小区间数。最后，我们将结果打印输出。/n/n此代码输出结果为：/n/n/nIntegral value: 0.9999999999999999/n/n/n说明在该区间内，$/int_{0}^{/frac{/pi}{2}}/sin(x)dx /approx 1$。