实验5复合辛普森公式 201226100108计自1201 一、实验目的 用复合辛普森公式计算积分dx ,使误差不超过-4 10 (注意所给积分特点, 做出相应的处理后再计算) 二、实验步骤 1.算法原理 复合辛普森原理: 将区间 lim此外,由于 中求积公系数均为正数,故知辛普森公式计算稳定。2.算法步骤 复合辛普森: 首先将区间 算法过程:这里将辛普森公式写为 Sn()函数,然后在 Solve()函数里依次计算 S1,S2,S4,S6.......当相 邻的精度小于eps 时退出循环,则S2n 保存结果。 三.程序代码 #include#include#define eps 1e-6 using namespace std; double f(double returnsqrt(1+cos(x)*cos(x)); }//被积函数 double Sn(double a,double b,double doublesum=0; for(int sum+=2*f(a+(2*k-1)*h);sum+=f(a+2*k*h); sum=(sum*2+f(a)-f(b))*h/3;return sum; }//辛普森公式 double Solve(double a,double doubleS1n=Sn(a,b,i); double S2n=Sn(a,b,2*i); while(fabs(S1n-S2n)>eps){ cout