1. 算法原理简介
步1 将积分区间 [a,b] 分成 n 等分,分点为 xk=a+kh(k=0,1,⋯,n),其中 h=(b-a)/n。
步2 记区间 [xk,x(k+1)] 的中点为 x(k+1/2),在每个小区间上用辛普森公式,
则得到复化辛普森公式:

2. 应用实例
取 n=10,利用复化辛普森公式计算积分:

3. 程序代码
#include<iostream>
using namespace std;
double fun(double x)
{
return 4 / (1 + x *