牛顿-科特斯公式∫abf(x)dx≈¿(b−a)∑i=0nCi(n)f(xi)¿科特斯(Cotes)系数Ci(n),特点:Cotes系数仅取决于n和i,可通过查表得到
与被积函数f(x)及积分区间[a,b]均无关
n=1:C0(1)=12,C1(1)=12为梯形求积公式∫abf(x)dx≈b−a2[f(a)+f(b)]梯形求积公式的几何意义:用梯形面积近似代替曲边梯形的面积梯形公式的余项为−(b−a)312f''(η)代数精度=1-0
5n=2:C0(2)=16,C1(2)=23,C2(2)=16Simpson求积公式(为抛物线求积公式)∫abf(x)dx≈b−a6[f(a)+4f(a+b2)+f(b)]辛普森公式的余项为−b−a180(b−a2)4f(4)(η)代数精度=3n=4:科特斯(Cotes)求积公式(五点公式)∫abf(x)dx≈b−a90[7f(x0)+32f(x1)+12f(x2)+32f(x3)+7f(x4)]xi=a+i⋅h,h=(b−a)/4柯特斯公式的余项为−2(b−a)495(b−a4)6f(6)(η)柯特斯公式具有5次代数精度科特斯系数具有以下特点:(1)∑i=0nCi(n)=1(2)Ci(n)=Cn−i(n)(3)当n8时,出现负数,稳定性得不到保证
而且当n较大时,由于Runge现象,收敛性也无法保证
一般不采用高阶的牛顿-科特斯求积公式
当n7时,牛顿-科特斯公式是稳定的
当n为偶数时,牛顿-科特斯公式至少有n+1阶代数精度
牛顿-柯特斯公式的舍入误差只是函数值误差的(b−a)倍复化求积公式特点当积分区间[a,b]的长度较大,而节点个数n+1固定时直接使用牛顿-柯特斯公式余项将会较大而如果增加节点个数,n即+1增加时当n>8时,公式的舍入误差又很难得到控制此时,使用复化方法,即将积