实验四复合辛普森公式.docVIP

实验四 复合辛普森公式 一、实验目的 1、理解复合辛普森公式的方法原理； 2、利用复合辛普森公式计算积分。 二、实验步骤 1.算法原理 复合辛普森原理： 将区间划分为n等分，在每个子区间上采用辛普森公式，若记则得 记 称为复合辛普森求积公式，其余项为 于是当时，与复合梯形公式相似有 易知误差阶为，收敛性是显然的，实际上，只要则可得到收敛性，即 此外，由于中求积公系数均为正数，故知辛普森公式计算稳定。 2.算法步骤 复合辛普森： 首先将区间划分为n等分，在每个子区间上采用辛普森公式，若记则得 3.程序流程图 三.程序代码 #include #include #include #define eps 0.000001 //计算精度 float func double x //% func-- 求积函数（函数文件） while x 0 return 1; double y; y sin x /x ; return y ; double bsimpson double a,double b //% a, b -- 求积区间 int i,n; double h,p,e,s; double t1,t2,s1,s2,x; n 1; h b-a; t1 h* func a +func b /2.0; s1 t1; /*用t1代替s1*/ e eps+1.0; whi